Abstract
The recent detection of very high energy (VHE, greater than or similar to 100 GeV) gamma-ray emission from gamma-ray bursts (GRBs) has provided new insights into afterglow physics. Understanding the temporal and spectral evolution of VHE GRBs requires detailed modelling of multiwavelength observations spanning radio to VHE gamma-rays. Previous studies interpreted afterglow of VHE GRBs using a range of frameworks, including single-and multizone jet configurations, synchrotron radiation from forward and reverse shocks, synchrotron self-Compton (SSC) processes, as well as hadronic emission processes. We have modelled five long-duration VHE GRBs-GRB 180720B, GRB 190114C, GRB 190829A, GRB 201216C, and GRB 221009A; using the NAIMA code and modifications to it. The results from our analysis indicate that SSC is the dominant VHE emission mechanism, with negligible contribution from external Compton. Most VHE GRBs are well described by the forward shock model in a spherical jet configuration, where constant density interstellar medium is preferred over wind medium. Additionally, we find that VHE GRBs tend to occur in environments with lower magnetic fields and higher ambient medium densities. Interestingly, VHE GRBs lie at the edge of the 3(sigma) region of the E-k,E- iso-is an element of(B) correlation observed in other energetic GRBs. Our model slightly over predicts the radio fluxes, indicating that a more complicated modelling might be required in some cases. These findings provide crucial constraints on VHE GRB emission sites and mechanisms and serve as a benchmark for future observations and theoretical studies in the era of CTA and next-generation gamma-ray observatories.