Gamma-ray burst afterglow scaling relations for the full blast wave evolution

We demonstrate that gamma-ray burst afterglow spectra and light curves can be calculated for arbitrary explosion and radiation parameters by scaling the peak flux and the critical frequencies connecting different spectral regimes. Only one baseline calculation needs to be done for each jet opening angle and observer angle. These calculations are done numerically using high-resolution relativistic hydrodynamical afterglow blast wave simulations which include the two-dimensional dynamical features of expanding and decelerating afterglow blast waves. Any light curve can then be generated by applying scaling relations to the baseline calculations. As a result, it is now possible to fully fit for the shape of the jet break, e.g. at early time X-ray and optical frequencies. In addition, late-time radio calorimetry can be improved since the general shape of the transition into the Sedov-Taylor regime is now known for arbitrary explosion parameters so the exact moment when the Sedov-Taylor asymptote is reached in the light curve is no longer relevant. When calculating the baselines, we find that the synchrotron critical frequency and the cooling break frequency are strongly affected by the jet break. The synchrotron break temporal slope quickly drops to the steep late time Sedov-Taylor slope, while the cooling break first steepens then rises to meet the level of its shallow late time asymptote.

💡 Research Summary

Gamma‑ray burst (GRB) afterglows arise from a relativistic jet that decelerates as it sweeps up the circumburst medium, producing broadband synchrotron emission. Traditionally, analytical treatments have relied on the ultra‑relativistic Blandford‑McKee (BM) solution for early times and the non‑relativistic Sedov‑Taylor (ST) solution for late times, while the intermediate trans‑relativistic phase has been approximated with simplified 1‑D models. This paper demonstrates that the full evolution of the afterglow—covering BM, the intermediate stage, and ST—can be captured by a single set of high‑resolution two‑dimensional relativistic hydrodynamic (RHD) simulations combined with a set of exact scaling relations.

The authors first identify that the fluid state at any radius r, angle θ and lab‑frame time tₑ can be expressed through a small number of dimensionless combinations (A = r/ctₑ, B = E_iso tₑ²/ρ₀ r⁵, the initial jet half‑opening angle θ₀, and the observer angle θ_obs). Because the governing equations are scale‑invariant, a transformation of the basic physical parameters—energy (E_iso), ambient density (ρ₀), radius and time—leads to a simple scaling: E′ = κ E_iso, ρ′ = λ ρ₀, r′ = (κ/λ)¹ᐟ³ r, t′ = (κ/λ)¹ᐟ³ t. Crucially, this scaling holds not only in the asymptotic BM and ST limits but also throughout the intermediate regime where the jet spreads laterally and becomes inhomogeneous.

The synchrotron radiation formalism (j_ν ∝ ξ_N n B f(ν,ν_m,ν_c)/γ²(1−βμ)²) together with the radiative transfer equation shows that the observed flux in any spectral segment scales as κ λ¹ᐟ³, provided the peak flux F_peak and the characteristic frequencies ν_m (synchrotron break) and ν_c (cooling break) are scaled accordingly. Table 1 in the paper lists the explicit scaling exponents for all four spectral regimes (self‑absorbed, ν_a < ν < ν_m, ν_m < ν < ν_c, ν > ν_c) in both BM and ST limits. Therefore, once a baseline simulation is performed for a given jet opening angle and observer angle, the entire family of light curves and spectra for arbitrary E_iso, n₀, θ₀, and θ_obs can be generated by simple algebraic transformations—no additional costly RHD runs are required.

Numerical verification (Figures 2–3) confirms that the scaling works at all times. Light curves for two isotropic energies (10⁴⁸ erg and 10⁵⁰ erg) overlay perfectly after applying the scaling, regardless of observer angle. The evolution of ν_m, ν_c, and F_peak for an on‑axis observer shows the expected asymptotic slopes: ν_m transitions from t⁻³ᐟ² (BM) to t⁻³ (ST); ν_c initially follows t⁻¹ᐟ², steepens temporarily to ~t⁻¹·² around the jet break, then settles to the ST slope t⁻³ᐟ⁵, even rising for a period after the break. This rise of ν_c matches observations of GRB 091127 and illustrates that the cooling break is strongly affected by jet dynamics. The temporal slope of ν_c is less steep at late times, and its ST asymptote lies higher than the BM asymptote, reflecting the change from isotropic to true jet energy (E_j ≈ E_iso θ₀²/2).

The study also shows that larger observer angles delay the jet break and accentuate the ν_c rise, emphasizing the importance of the full 2‑D jet structure for interpreting multi‑wavelength data. Because the transition to the ST regime is now known analytically for any set of parameters, late‑time radio calorimetry no longer depends on the precise moment the light curve reaches the ST asymptote; the full shape can be modeled accurately.

In summary, the paper establishes a powerful, physics‑based framework: a single high‑resolution 2‑D simulation per jet opening angle, combined with exact energy‑density‑time scalings, yields the complete afterglow evolution for any explosion parameters and viewing geometry. This dramatically reduces computational cost, enables robust fitting of jet‑break shapes in X‑ray and optical bands, and improves energy estimates from radio observations. The authors suggest future extensions to include evolving microphysical parameters (ε_B, ε_e) and non‑uniform external media, which would further broaden the applicability of this scaling‑based approach.