An Afterglow Study of the "New Year's Burst" GRB 220101A

TL;DR

This study analyzes the broadband afterglow of GRB 220101A, an extremely energetic burst at $z=4.618$ with $E_{ m 3iso}\sim3.6\times10^{54}$ erg, by combining X-ray, optical/IR, sub-mm, and radio data from $10^4\lesssim\Delta T\lesssim10^7$ s. It employs afterglowpy to model both Top-hat and Gaussian jets in an ISM-like circumburst medium, with two electron-participation fractions $\xi=0.1$ and $1.0$, and complements this with analytical fits; a Markov-chain Monte Carlo framework is used to constrain $E_{ m K,iso}$, $p$, $\theta_j$, $\theta_{ m obs}$, $n_0$, $\epsilon_e$, $\epsilon_B$, and extinction $A_B$. The results show a distinct jet break near $t_j\sim8.6$–9 days, and a remarkably steep post-break slope $\alpha_{ m post}\approx -2.99\pm0.10$, which is significantly steeper than the spectral-inferred $p\sim2.0$–$2.6$ and cannot be easily explained by simple lateral expansion; both jet geometries require extremely low circumburst densities, $n_0\lesssim10^{-4}$ cm$^{-3}$ (with radio data driving the low densities), while the energy budget remains modest, $E_{ m K,iso}\sim10^{51}$–$10^{52}$ erg. Excluding the radio data raises $n_0$ by orders of magnitude but then overpredicts the radio flux, a tension also seen in other LAT-detected GRBs. The authors discuss a wind-evacuation (cavity) scenario as a natural way to obtain such low densities and highlight potential limitations of the standard model, including the need for time-dependent microphysics or more complex hydrodynamics to fully explain the radio behavior. Overall, GRB 220101A provides a stringent test of jet physics, the circumburst environment, and the completeness of afterglow modelling in the era of high-energy LAT GRBs.

Abstract

We present a detailed broadband afterglow study of GRB 220101A ($10^4\lesssimΔT\lesssim10^7$ s) combining multi-wavelength data from soft X-rays until 6 GHz. The afterglow light curves in both X-ray and optical show distinct steepening around $\sim9$ days, followed by a sharp post-break decay index of $\sim2.99\pm0.10$. We fit the light curves using the afterglow modelling package \texttt{afterglowpy} for both Top-hat and Gaussian jets for different values of the electronic participation fraction $ξ$ from 0.01 to 1.0 and find that, although the radio behavior is well described by the $ξ=1.0$ case, the required circumburst medium (CBM) densities are very low, $<10^{-4}$ cm$^{-3}$. However, the resulting energy requirements are modest, $\sim10^{52}$ erg, with an electron energy distribution (EED) index $p\sim2.05$. Similar results are also obtained from an analytic model fit to the light curve, except the predicted $p$ is higher, $\sim2.40$. The observed post-break decay index of $2.99$ is at least 5$σ$ away from $p$, which is one of the steepest observed so far. We also find that when ignoring the radio observations, the CBM density is raised by a few orders of magnitude to $\sim0.01$ cm$^{-3}$ for $ξ=1.0$, still far from the expected ISM density ($>1$ cm$^{-3}$) of GRB environments, which are highly star forming regions. Similarly low ISM densities have been seen in modeling of other LAT GRBs as well, especially ones with reverse-shock features (e.g., GRBs 130427A, 160509A and 160625B), thereby hinting at either an issue with the standard model or possible evacuated cavities where GRBs explode.

Figures (8)

• Figure 1: HST F125W image of the GRB 220101A. The two blue lines mark the position of the transient. Left : at $\sim38$ days; Middle: at $\sim236$ days; Right: background subtracted image at $\sim38$ days.
• Figure 2: Figure showing the I band and the X-ray points, which were sampled the most among all wavelengths, and the results of the broken power law fit using Equation \ref{['eqn:bpl']}.
• Figure 3: Radio to submm light curves plotted separately (scaled accordingly) for better clarity. While no clear smooth/broken power law behaviour is evident like optical and X-rays, a general change in trend after 10 days is evident for the $6-10$ GHz light curves.
• Figure 4: Afterglowpy top-hat jet fits to the light curves for $\xi=0.1$ and $\xi=1.0$. Both the optical/IR and the radio data are fit better with the higher $\xi$. The grey points show the unbinned X-ray data.
• Figure 5: Afterglowpy spectral evolution for $\xi=0.1$ and $\xi=1.0$ with the top-hat jet. The evolution of the break frequencies is considerably different between the two models, and the $\xi=1.0$ model fits the X-ray evolution better. The kink in the figures is due to the Heaviside-like (step function in frequency) extinction correction applied to the optical bands.