Table of Contents
Fetching ...

Modelling multiwavelength afterglows of the VHE-GRB population

Monica Barnard, Ankur Ghosh, Jagdish C. Joshi, Soebur Razzaque

TL;DR

The paper presents a uniform one-zone forward-shock modelling approach for five VHE-detected GRBs, using NAIMA with MCMC to fit multiwavelength afterglow SEDs and assess ISM vs wind environments. It finds SSC-dominated VHE emission with negligible external Compton, generally favors ISM for several bursts and wind for the brightest, and reveals a population-level anti-correlation between $E_{ m k,iso}$ and $\epsilon_B$, while noting radio flux overpredictions in some epochs. The study constrains emission sites and microphysical parameters, highlighting that VHE GRBs probe low $\epsilon_B$ environments and dense media, and suggests that more complex, multi-zone or structured-jet models may be required for some events. The results provide benchmarks for current and upcoming facilities (CTA, next-gen γ-ray observatories) to refine our understanding of particle acceleration and radiation in extreme relativistic jets.

Abstract

The recent detection of very high energy (VHE, $\gtrsim$ 100 GeV) $γ$-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 $γ$ rays. Previous studies interpreted afterglow of VHE GRBs using a range of frameworks, including single- and multi-zone jet configurations, synchrotron radiation from forward and reverse shocks, synchrotron self-Compton (SSC) processes, as well as hadronic emission processes. We have modeled 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σ$ region of the $E_{\rm k,iso}$ - $ε_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 $γ$-ray observatories.

Modelling multiwavelength afterglows of the VHE-GRB population

TL;DR

The paper presents a uniform one-zone forward-shock modelling approach for five VHE-detected GRBs, using NAIMA with MCMC to fit multiwavelength afterglow SEDs and assess ISM vs wind environments. It finds SSC-dominated VHE emission with negligible external Compton, generally favors ISM for several bursts and wind for the brightest, and reveals a population-level anti-correlation between and , while noting radio flux overpredictions in some epochs. The study constrains emission sites and microphysical parameters, highlighting that VHE GRBs probe low environments and dense media, and suggests that more complex, multi-zone or structured-jet models may be required for some events. The results provide benchmarks for current and upcoming facilities (CTA, next-gen γ-ray observatories) to refine our understanding of particle acceleration and radiation in extreme relativistic jets.

Abstract

The recent detection of very high energy (VHE, 100 GeV) -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 rays. Previous studies interpreted afterglow of VHE GRBs using a range of frameworks, including single- and multi-zone jet configurations, synchrotron radiation from forward and reverse shocks, synchrotron self-Compton (SSC) processes, as well as hadronic emission processes. We have modeled 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 region of the - 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 -ray observatories.

Paper Structure

This paper contains 22 sections, 18 equations, 26 figures, 2 tables.

Figures (26)

  • Figure 1: The SEDs for GRB 180720B at several epochs, as indicated in the legend, for the ISM (left) and wind (right) density scenarios. For each epoch the synchrotron (dashed), SSC (dashed dotted), EC (dotted) and total (solid) emission are modelled and fitted to the multiwavelength data, i.e., optical (crosses), Swift-XRT (squares), Swift-BAT and Fermi-LAT (shaded regions), and H.E.S.S. (circles). For the 10 hr epoch the X-ray data is an extrapolation to the time window of the H.E.S.S. observation. The data were collected from 2018GCN.22977....1SAbdalla2019Fraija_2019Ronchi2020 and https://www.swift.ac.uk/.
  • Figure 2: The SEDs for GRB 190114C at several epochs, as indicated in the legend, for both the ISM (left) and wind (right) density scenarios. For each epoch the synchrotron (dashed), SSC (dashed dotted), EC (dotted) and total (solid) emission are modelled and fitted to the multiwavelength data, i.e., optical (crosses), Swift-XRT (diamonds and shaded regions), Swift-BAT (squares), Fermi-LAT (stars), and MAGIC (circles). The data were collected from 2019Natur.575..459M and https://www.swift.ac.uk/.
  • Figure 3: The SEDs for GRB 190829A at several epochs, as indicated in the legend, for the ISM (left) and wind (right) density scenarios. For each epoch the synchrotron (dashed), SSC (dashed dotted), EC (dotted) and total (solid) emission are modelled and fitted to the multiwavelength data, i.e., radio (triangles), optical (crosses), Swift-XRT (squares and shaded region), Fermi-LAT upper limit, and H.E.S.S. (circles). The data were collected from 2019GCN.25589....1D2019GCN.25676....1L2020MNRAS.496.3326RHESS:2021dbzHu2021 and https://www.swift.ac.uk/.
  • Figure 4: The SEDs for GRB 201216C at several epochs, as indicated in the legend, for the ISM density (left) and the stellar wind (right) scenarios. For each epoch the synchrotron (dashed), SSC (dashed dotted), EC (dotted) and total (solid) emission are modelled and fitted to the multiwavelength data, i.e., radio (triangles), optical (crosses), Swift-XRT (diamonds), and MAGIC (circles). The data were collected from 2020GCN.29066....1I2022MNRAS.513.1895R and Abe2024.
  • Figure 5: The SEDs for GRB 221009A at several epochs, as indicated in the legend, for the ISM density (left) and the stellar wind (right) scenarios. For each epoch the synchrotron (dashed), SSC (dashed dotted), and total (solid) emission are modelled and fitted to the multiwavelength data, i.e., radio (triangles), optical (crosses), Swift-XRT (diamonds), GBM (squares), AGILE (stars), and LHAASO (circles). The data were collected from 2023NatAs...7..986B2023Sci...380.1390Loconnor2023SciA....9I1405O2023ApJ...956L..23T and Banerjee2024. The EC was omitted for figure clarity.
  • ...and 21 more figures