A Semi-analytic Formulation for Relativistic Blast Waves with a Long-lived Reverse Shock

This paper performs a semi-analytic study of relativistic blast waves in the context of gamma-ray bursts (GRBs). Although commonly used in a wide range of analytical and numerical studies, the equation of state (EOS) with a constant adiabatic index is a poor approximation for relativistic hydrodynamics. Adopting a more realistic EOS with a variable adiabatic index, we present a simple form of jump conditions for relativistic hydrodynamical shocks. Then we describe in detail our technique of modeling a very general class of GRB blast waves with a long-lived reverse shock. Our technique admits an arbitrary radial stratification of the ejecta and ambient medium. We use two different methods to find dynamics of the blast wave: (1) customary pressure balance across the blast wave and (2) the “mechanical model”. Using a simple example model, we demonstrate that the two methods yield significantly different dynamical evolutions of the blast wave. We show that the pressure balance does not satisfy the energy conservation for an adiabatic blast wave while the mechanical model does. We also compare two sets of afterglow light curves obtained with the two different methods.

💡 Research Summary

This paper presents a semi‑analytic framework for modeling relativistic blast waves in gamma‑ray bursts (GRBs), with particular emphasis on the dynamics of a long‑lived reverse shock (RS). The authors begin by pointing out a well‑known limitation of the commonly used equation of state (EOS) with a constant adiabatic index (γ = 4/3). Such an EOS accurately describes ultra‑relativistic gas but fails for gas that transitions from non‑relativistic to relativistic temperatures, a situation that naturally occurs in the RS region as the blast wave decelerates.

To overcome this, the authors adopt a more realistic EOS originally proposed by Mathews (1971) and later reformulated by Mignone et al. (2005). This EOS treats the post‑shock gas as a mono‑energetic relativistic fluid, i.e., all particles share the same Lorentz factor. In this formulation the effective adiabatic index κ varies smoothly between 2/3 (non‑relativistic limit) and 1/3 (ultra‑relativistic limit) as a function of the shock strength. The key advantage is that κ can be expressed analytically in terms of the shock Lorentz factor γ₁₂, eliminating the need for modified Bessel functions while retaining the correct limiting behavior.

Using this EOS, the authors re‑derive the relativistic shock jump conditions. By enforcing continuity of mass, energy, and momentum fluxes across the shock front, they obtain compact expressions for the downstream Lorentz factors (γ₁, γ₂), pressure (p₂), and compression ratio (a = ρ₂/ρ₁) that depend solely on the shock strength γ₁₂. Equations (22)–(25) in the paper give these relations in a remarkably simple form, and they reduce to the classic Blandford‑McKee results when κ is constant.

The blast‑wave structure is then described in the standard four‑region picture: (1) unshocked external medium, (2) shocked external medium (forward‑shock region), (3) shocked ejecta (reverse‑shock region), and (4) unshocked ejecta. The authors assume spherical symmetry and that the entire shocked region (regions 2 and 3) moves with a common Lorentz factor Γ, justified by sub‑sonic internal motions observed in numerical simulations. A novel aspect of the work is the treatment of a radially stratified ejecta. Each ejecta shell is labeled by its ejection time τ, which serves as a Lagrangian coordinate. The authors solve the continuity equation for the ejecta (Eq. 26) to obtain analytic expressions for the density profile ρ_ej(r, t) and the Lorentz factor distribution Γ_ej(τ). This allows them to follow the RS as it propagates through ejecta with arbitrary radial structure.

Two distinct methods are employed to determine the blast‑wave dynamics:

1. Pressure‑balance method – The traditional approach assumes that the pressure immediately behind the forward shock (FS) equals that behind the reverse shock. This yields a simple ordinary differential equation for Γ(t). However, when applied to the example model, the authors find that total energy is not conserved; the blast loses about a factor of five in energy, indicating a fundamental inconsistency.

2. Mechanical model – Building on Beloborodov & Uhm (2006), the authors write down conservation equations for mass, energy, and momentum fluxes across the FS and RS without imposing pressure equality. The resulting system of equations is solved simultaneously for Γ(t) and the RS radius r_RS(t). This “mechanical” approach respects global energy conservation and produces a markedly different evolution: the RS persists longer, the deceleration of the blast is slower, and the pressure behind the FS is higher at early times compared with the pressure‑balance solution.

To illustrate the observational consequences, the authors compute afterglow light curves in the X‑ray and optical bands for both dynamical prescriptions. The pressure‑balance model predicts a rapid decline once the RS weakens, while the mechanical model yields a flatter, brighter afterglow because the long‑lived RS continues to inject energy into the shocked region. These differences are significant for interpreting the diversity of observed GRB afterglows, especially those showing plateaus or re‑brightenings that are difficult to reconcile with a short‑lived RS.

In summary, the paper makes three principal contributions: (i) it introduces a variable‑κ EOS that accurately captures the thermodynamics of both non‑relativistic and relativistic post‑shock gas; (ii) it derives a set of compact jump conditions that depend only on the shock Lorentz factor, greatly simplifying analytical treatments; and (iii) it demonstrates, through a semi‑analytic mechanical model, that enforcing global conservation laws yields physically consistent blast‑wave dynamics and markedly different afterglow predictions compared with the traditional pressure‑balance assumption. The framework is flexible enough to accommodate arbitrary radial stratification of the ejecta and ambient medium, making it a valuable tool for both analytic studies and as a benchmark for numerical simulations of GRB blast waves.