VegasAfterglow: Modular GRB Afterglow Framework

• VegasAfterglow is a high-performance, modular C++ framework designed to model gamma-ray burst afterglows using integrated shock dynamics, radiation calculations, and arbitrary ambient density profiles.
• It balances computational efficiency and physical comprehensiveness, supporting MCMC/Bayesian inference with forward and reverse shock treatments, inverse Compton emission, and self-absorption effects.
• Its modular design and support for varied jet structures, energy injection, and environmental profiles make it ideal for interpreting complex events like off-axis GRBs and multi-messenger observations.

Searching arXiv for the VegasAfterglow paper and closely related GRB afterglow modeling references. VegasAfterglow is a high-performance, modular C++ framework for gamma-ray burst afterglow modeling, developed to bridge the gap between very fast but simplified semi-analytic codes and very detailed but computationally expensive hydrodynamic simulations. It is designed to remain physically comprehensive while being fast enough for broad parameter inference such as MCMC and Bayesian fitting. The framework self-consistently solves forward and reverse shock dynamics, computes synchrotron radiation across all standard spectral regimes including self-absorption, and includes inverse Compton emission with Klein–Nishina corrections. It also supports arbitrary ambient density profiles, central-engine activity histories, viewing angles, and jet structures in energy, Lorentz factor, and magnetization, while covering relativistic, trans-relativistic, and deep-Newtonian evolution and including lateral jet spreading effects (Wang et al., 14 Jul 2025).

1. Scientific setting and motivation

Gamma-ray burst afterglows are the long-lived emission that follows the prompt gamma-ray flash. Since their discovery by BeppoSAX in 1997, afterglows have been studied from the GeV band down to radio and from seconds to months after trigger. The pre-Swift standard model attributed most afterglow emission to synchrotron radiation from shock-accelerated electrons in a forward shock, with the observable law commonly written as $F_\nu \propto t^{-\alpha}\nu^{-\beta}$. Swift and Fermi then exposed major complexities in the early afterglow phase, including steep early decays, plateaus, flares, chromatic breaks, and long-lived high-energy emission beyond the prompt phase (Godet et al., 2012).

VegasAfterglow is motivated by precisely this post-Swift and post-Fermi landscape. The framework is built around the problem of modeling GRB afterglows accurately across all relevant physical regimes without making the code too slow for inference. Its stated rationale is that simple top-hat models are often adequate only for on-axis cases, whereas off-axis events such as GW170817/GRB 170817A demand structured jets and accurate geometry. Existing public tools are described as often lacking some combination of reverse-shock physics, thick-shell treatment, self-absorption, inverse Compton cooling and emission, Klein–Nishina corrections, arbitrary density or injection histories, and full non-relativistic to relativistic evolution (Wang et al., 14 Jul 2025).

This positioning places VegasAfterglow within the broader transition of afterglow theory from a relatively compact forward-shock paradigm toward a multi-component description in which reverse shocks, off-axis geometry, long-lived engine activity, and evolving radiative regimes can all become observationally relevant. A plausible implication is that the framework is intended less as a single-scenario light-curve generator than as a general inference engine for modern broadband and multi-messenger afterglow datasets.

2. Physical scope and model architecture

VegasAfterglow models the outflow on a spherical grid in $(t,\theta,\phi)$, while exploiting symmetries to reduce computational cost in the radiation step. The jet is specified through user-defined angular profiles of the initial Lorentz factor $\Gamma_0(\phi,\theta)$, the energy per solid angle $dE/d\Omega(\phi,\theta)$, and the magnetization $\sigma(\phi,\theta)$. The framework includes top-hat, Gaussian, and power-law jet profiles as explicit examples, and it permits arbitrary user-defined angular structures (Wang et al., 14 Jul 2025).

The ambient medium is likewise user-defined through density and mass profiles. Standard forms include a homogeneous interstellar medium,

$\rho(\phi,\theta,r)=n_0 m_{\rm p},$

and a wind-like medium,

$\rho(\phi,\theta,r)=\frac{A}{r^k}.$

This allows the code to represent both merger-like uniform environments and collapsar-like winds. Central-engine activity can also be imposed through arbitrary injection histories $L_{\rm inj}(\phi,\theta,t)$, including a default magnetar spin-down example with $\sigma \to \infty$ to mimic a Poynting-flux-dominated wind (Wang et al., 14 Jul 2025).

The architectural claim of the framework is therefore not merely that it contains many physical options, but that these options are integrated into one computational scheme. Forward and reverse shocks, off-axis viewing, arbitrary jet geometry, energy injection, lateral spreading, synchrotron self-absorption, inverse Compton emission, and deep-Newtonian evolution are all treated within the same package. This differs from the more limited configuration space typical of simplified public afterglow tools.

3. Shock dynamics and hydrodynamical treatment

A central component of VegasAfterglow is its self-consistent treatment of both forward and reverse shocks, including the transition between thin-shell and thick-shell ejecta. For arbitrary upstream magnetization $\sigma_u$ and arbitrary relative Lorentz factor $(t,\theta,\phi)$0, the framework uses an analytical shock-jump solution based on the shock-frame continuity condition

$(t,\theta,\phi)$1

The downstream four-velocity $(t,\theta,\phi)$2 is obtained from a cubic equation, and the framework emphasizes that only one root is physical, with the others violating energy conservation. In the unmagnetized limit, the solution reduces to

$(t,\theta,\phi)$3

This places the dynamical core in a regime that remains valid when the shock is mildly relativistic or magnetized, where simpler formulae can fail (Wang et al., 14 Jul 2025).

The reverse shock is treated separately in the crossing and post-crossing phases. During crossing, forward and reverse shocks can be linked by pressure balance,

$(t,\theta,\phi)$4

but the framework notes that this approximation can violate energy conservation for thick shells. VegasAfterglow therefore adopts an energy-conservation-based method,

$(t,\theta,\phi)$5

tracks the shocked-proton number in region 3, and defines reverse-shock crossing by the condition $(t,\theta,\phi)$6 (Wang et al., 14 Jul 2025).

One of the framework’s explicit physical conclusions is that in the thin-shell regime the reverse shock is weaker than often assumed, with a typical peak $(t,\theta,\phi)$7 rather than the often-quoted $(t,\theta,\phi)$8. The stated significance is that reverse-shock optical flashes may otherwise be overpredicted. After the crossing phase, the reverse-shocked region expands adiabatically, while the forward shock continues to decelerate into the external medium and is evolved through the Blandford–McKee and eventually deep-Newtonian regimes. Lateral spreading is included through a smoother transition prescription rather than the simplest $(t,\theta,\phi)$9 prescription, which is described as better matching hydrodynamic expectations (Wang et al., 14 Jul 2025).

4. Radiation physics and observer mapping

The radiation module is computed in the fluid comoving frame and then Doppler-transformed to the observer frame. For axisymmetric jets, the intrinsic emission calculation is reduced to a 2D grid in $\Gamma_0(\phi,\theta)$0, which the framework identifies as a major reason for speed (Wang et al., 14 Jul 2025).

Electron injection is handled with an explicit deep-Newtonian correction. The minimum electron Lorentz factor $\Gamma_0(\phi,\theta)$1 is written with a terminal “$\Gamma_0(\phi,\theta)$2” term, and the framework states that this ensures $\Gamma_0(\phi,\theta)$3 and correctly handles low-energy, deep-Newtonian electrons. To correct synchrotron suppression when electrons become only mildly relativistic, it introduces

$\Gamma_0(\phi,\theta)$4

Cooling is computed with synchrotron and inverse Compton losses through

$\Gamma_0(\phi,\theta)$5

with corresponding expressions for $\Gamma_0(\phi,\theta)$6 and $\Gamma_0(\phi,\theta)$7 (Wang et al., 14 Jul 2025).

Synchrotron self-absorption is modeled explicitly through the intersection of optically thin synchrotron emission with a blackbody-like low-frequency limit,

$\Gamma_0(\phi,\theta)$8

with $\Gamma_0(\phi,\theta)$9. The framework treats the three standard absorption cases $dE/d\Omega(\phi,\theta)$0, $dE/d\Omega(\phi,\theta)$1, and $dE/d\Omega(\phi,\theta)$2, and it lists six broken-power-law electron distributions covering slow cooling or fast cooling under weak or strong absorption. The synchrotron intensity itself is represented as piecewise broken power laws with an exponential cutoff $dE/d\Omega(\phi,\theta)$3 (Wang et al., 14 Jul 2025).

Inverse Compton treatment is one of the framework’s defining extensions beyond simplified afterglow models. It uses a $dE/d\Omega(\phi,\theta)$4-dependent Compton parameter $dE/d\Omega(\phi,\theta)$5, modifies cooled electron and synchrotron distributions accordingly, and computes the full synchrotron self-Compton spectrum numerically with the Klein–Nishina cross section. This means that the inverse Compton component is not restricted to approximate SSC power laws (Wang et al., 14 Jul 2025).

Observer geometry is handled through the Doppler factor

$dE/d\Omega(\phi,\theta)$6

with

$dE/d\Omega(\phi,\theta)$7

The observed flux is then obtained by integration over the equal-arrival-time surface, with

$dE/d\Omega(\phi,\theta)$8

This formulation is essential for off-axis structured jets, because it accounts simultaneously for relativistic beaming and light-travel-time effects (Wang et al., 14 Jul 2025).

5. Implementation, performance, and benchmarking

VegasAfterglow is implemented in C++ with Python interfaces. The reported performance result is that, on a 2022 Apple M2 laptop, a single-frequency light curve can be computed in roughly $dE/d\Omega(\phi,\theta)$9 for a forward-shock synchrotron-only calculation at sufficiently low resolution, while still converging at higher resolution. The framework attributes this speed to log-log interpolation for power-law-like functions, reduction of expensive repeated power-law evaluations, sequential memory access to improve cache behavior, and 2D intrinsic emission grids for axisymmetric jets with observer geometry applied afterward (Wang et al., 14 Jul 2025).

Validation is carried out through comparison with several public codes: AfterglowPy, JetsimPy, PyFRS, and PyBlastAfterglow. In top-hat cases, the reported outcome is good agreement. In structured cases, the framework notes that differences arise from how $\sigma(\phi,\theta)$0 is modeled and whether the $\sigma(\phi,\theta)$1 correction is used. VegasAfterglow applies the deep-Newtonian correction by default and predicts a late-time flattening consistent with its electron treatment. In radio, VegasAfterglow and PyFRS suppress flux relative to codes without self-absorption. In radiative fireball scenarios, the code produces fainter light curves and earlier jet breaks when radiative losses are strong. For reverse shocks, it reproduces both thin-shell and thick-shell behavior and finds that thick-shell reverse-shock light curves are delayed and weaker. It can also generate plateaus for magnetar-like injection and supports arbitrary angular jet profiles evaluated at multiple viewing angles (Wang et al., 14 Jul 2025).

These benchmarks clarify the framework’s intended identity. “High-performance” does not denote a stripped-down implementation, and “physically comprehensive” does not denote a numerically prohibitive one. The stated contribution is precisely the combination of computational tractability and unusually broad physical coverage.

6. Observational relevance and relation to afterglow phenomenology

The framework is intended for broadband fitting of GRB afterglows from radio to X-rays, interpretation of off-axis jets, reverse-shock diagnostics from early optical and radio data, modeling of plateaus from continued engine activity, studies of magnetized jets and structured outflows, fitting of GW+GRB events such as GW170817-like mergers, and population studies and Bayesian parameter estimation (Wang et al., 14 Jul 2025).

This use-case profile corresponds closely to major observational problems in the afterglow literature. The possible Fermi-GBM event associated with GW150914 was modeled as either a weak on-axis short GRB or a normal short GRB seen far off axis, with detectability depending strongly on jet geometry, density, and timing; in the brightest off-axis $\sigma(\phi,\theta)$2 case, the predicted radio afterglow reached $\sigma(\phi,\theta)$3–$\sigma(\phi,\theta)$4 at $\sigma(\phi,\theta)$5 and $\sigma(\phi,\theta)$6–$\sigma(\phi,\theta)$7 at $\sigma(\phi,\theta)$8, peaking about $\sigma(\phi,\theta)$9–$\rho(\phi,\theta,r)=n_0 m_{\rm p},$0 months after the merger, with a hard-to-soft low-frequency spectral evolution as self-absorption lifted (Morsony et al., 2016). A plausible implication is that a framework combining off-axis geometry, self-absorption, arbitrary media, and late-time dynamics is directly aligned with such “late-rising radio source” scenarios.

The same is true for reverse-shock and high-energy problems. GRB 130427A displayed an optical flash peaking at $\rho(\phi,\theta,r)=n_0 m_{\rm p},$1 between $\rho(\phi,\theta,r)=n_0 m_{\rm p},$2 and $\rho(\phi,\theta,r)=n_0 m_{\rm p},$3, followed by an optical afterglow whose post-$\rho(\phi,\theta,r)=n_0 m_{\rm p},$4 evolution closely tracked the $\rho(\phi,\theta,r)=n_0 m_{\rm p},$5 LAT light curve. That event was interpreted in an external-shock framework with reverse-shock emission dominating the optical flash and a forward shock emerging later in a wind-like medium, with multiple episodes of energy injection required to fit the data (Vestrand et al., 2013). VegasAfterglow’s explicit support for reverse shocks, energy injection, wind-like media, synchrotron emission across spectral regimes, and inverse Compton cooling with Klein–Nishina corrections makes it technically suited to this class of problem.

Within the broader development of afterglow theory, VegasAfterglow therefore occupies the methodological space opened by Swift and Fermi: a space in which late afterglows may still be treated successfully within external-shock physics, but early and off-axis observations often demand reverse shocks, structured jets, injection histories, spectral self-absorption, and multi-band geometric effects (Godet et al., 2012). Its distinguishing feature is not a new afterglow paradigm, but a computational framework broad enough to encode many of the physically motivated extensions that the post-Swift literature has made difficult to ignore.