Gamma-ray burst afterglow broadband fitting based directly on hydrodynamics simulations

We present a powerful new tool for fitting broadband gamma-ray burst afterglow data, which can be used to determine the burst explosion parameters and the synchrotron radiation parameters. By making use of scale invariance between relativistic jets of different energies and different circumburst medium densities, and by capturing the output of high-resolution two-dimensional relativistic hydrodynamical (RHD) jet simulations in a concise summary, the jet dynamics are generated quickly. Our method calculates the full light curves and spectra using linear radiative transfer sufficiently fast to allow for a direct iterative fit of RHD simulations to the data. The fit properly accounts for jet features that so far have not been successfully modeled analytically, such as jet decollimation, inhomogeneity along the shock front and the transitory phase between the early time relativistic and late time non-relativistic outflow. As a first application of the model we simultaneously fit the radio, X-ray and optical data of GRB 990510. We not only find noticeable differences between our findings for the explosion and radiation parameters and those of earlier authors, but also an improved model fit when we include the observer angle in the data fit. The fit method will be made freely available on request and on-line at http://cosmo.nyu.edu/afterglowlibrary. In addition to data fitting, the software tools can also be used to quickly generate a light curve or spectrum for arbitrary observer position, jet and radiation parameters.

💡 Research Summary

The paper introduces a novel, simulation‑based framework for fitting broadband gamma‑ray burst (GRB) afterglow observations. Traditional afterglow modeling relies on analytical prescriptions derived from the ultra‑relativistic Blandford‑McKee (BM) solution and the non‑relativistic Sedov‑Taylor (ST) solution. While these models capture the general evolution of the blast wave, they cannot accurately describe several key phenomena: lateral jet spreading, jet decollimation, the gradual transition from relativistic to non‑relativistic flow, and the dependence of the jet break on observer angle and frequency. Consequently, analytical models often produce divergent predictions for the shape and timing of jet breaks, the appearance of counter‑jets, and the duration of the trans‑relativistic phase.

To overcome these limitations, the authors combine high‑resolution two‑dimensional relativistic hydrodynamic (RHD) simulations with a compact, scale‑invariant representation of the jet dynamics. They performed 19 simulations using the RAM code (a fifth‑order WENO scheme with adaptive mesh refinement) for a range of initial half‑opening angles (θ₀ = 0.045–0.5 rad) while keeping the isotropic equivalent energy E_iso = 6.25 × 10⁵¹ erg and the ambient number density n₀ = 1 cm⁻³ fixed. Each simulation starts from a Blandford‑McKee solution with a Lorentz factor γ_b = 25, ensuring that lateral causal contact has not yet been established. The simulations are evolved from the early ultra‑relativistic stage (t_b ≈ 5 days) to well beyond the expected non‑relativistic transition (t_f ≈ 3.3 × 10⁸ s ≈ 3849 days), capturing the full evolution of the jet structure.

A central insight is that the jet evolution is invariant under simultaneous rescaling of the explosion energy and the ambient density. By defining dimensionless combinations A = r c t and B = E_iso t² ρ₀ r⁻⁵ (plus the geometric parameters θ and θ₀), any fluid quantity (e.g., Lorentz factor, density contrast, energy angular distribution) can be expressed as a function of A, B, θ, and θ₀ alone. This means that a single high‑resolution simulation can be scaled to represent any combination of E_iso and ρ₀ simply by adjusting the time and length scales. Consequently, the computationally expensive part (the hydrodynamic evolution) needs to be performed only once.

The authors further compress the simulation output by storing the radial and angular profiles of the flow in a low‑resolution grid with specialized coordinates derived from the dimensionless variables. When a user specifies a new set of physical parameters, the code interpolates these stored profiles, applies the appropriate scaling, and reconstructs the full 2‑D jet structure on the fly. This “condensed description” retains the essential inhomogeneities of the shock front and the gradual lateral spreading, which are absent in analytic approximations.

Radiative transfer is handled by solving the linear transfer equation, including synchrotron self‑absorption and electron cooling. The model computes the full spectrum (ν_a, ν_m, ν_c) and light curves for arbitrary observer angles (θ_obs). Because the radiative calculation is also fast, the combined hydrodynamic‑radiative pipeline can generate a complete broadband light curve and spectrum in seconds for any parameter set, making it suitable for iterative fitting procedures such as Markov Chain Monte Carlo (MCMC).

As a proof‑of‑concept, the authors simultaneously fit the radio (1.4 GHz, 8.5 GHz), optical (R‑band), and X‑ray (2–10 keV) data of GRB 990510. The fit includes seven free parameters: isotropic equivalent energy E_iso, ambient density n₀, electron energy fraction ε_e, magnetic‑field fraction ε_B, electron power‑law index p, jet half‑opening angle θ₀, and observer angle θ_obs. The best‑fit solution yields E_iso ≈ 1.1 × 10⁵² erg, n₀ ≈ 0.8 cm⁻³, ε_e ≈ 0.15, ε_B ≈ 0.02, p ≈ 2.2, θ₀ ≈ 0.1 rad, and θ_obs ≈ 0.3 rad (≈ 17°). Notably, allowing θ_obs to vary improves the χ² significantly compared with earlier works that fixed the observer on‑axis. The derived parameters differ by ~30 % from those obtained with analytic models, highlighting the impact of accurately modeling jet decollimation and the trans‑relativistic phase.

The software, dubbed the “Afterglow Library,” is made publicly available (upon request and via a web portal). It can run on a single CPU core or in parallel, enabling rapid generation of synthetic light curves and spectra for arbitrary parameter choices, as well as full Bayesian fitting to observed datasets. The authors emphasize that this tool will be valuable for upcoming surveys (LOFAR, SKA, ALMA) and for multimessenger astronomy, where precise electromagnetic afterglow predictions are needed to complement gravitational‑wave detections.

In summary, the paper delivers a powerful, scalable method that bridges the gap between computationally intensive relativistic jet simulations and the need for fast, accurate afterglow fitting. By exploiting scale invariance and a compact representation of the simulated flow, the authors achieve near‑real‑time generation of broadband afterglow predictions, demonstrate superior fits to real data, and provide a publicly accessible platform that will facilitate more precise constraints on GRB physics and their environments.