The hydrodynamics of stratified ultra-relativistic outflows and the origin of GRB X-ray plateaus

Gilad Sadeh; Kenta Hotokezaka; Masaru Shibata

REVIEW 3 major objections 4 minor 2 references

A radially stratified relativistic outflow with Lorentz factors down to ~70–100 can explain both the observed X-ray plateau phase of gamma-ray bursts and the subsequent transition to normal afterglow.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-02 20:02 UTC pith:NQDSWNNP

load-bearing objection Careful analytic framework for stratified ejecta; the plateau application is plausible but not yet a quantitative fit, and the mm prediction is the part worth taking seriously. the 3 major comments →

arxiv 2602.24085 v4 pith:NQDSWNNP submitted 2026-02-27 astro-ph.HE

The hydrodynamics of stratified ultra-relativistic outflows and the origin of GRB X-ray plateaus

Gilad Sadeh , Kenta Hotokezaka , Masaru Shibata This is my paper

classification astro-ph.HE

keywords gamma-ray burstsX-ray plateauforward shockreverse shockstratified outflowLorentz factor distributionafterglowmillimeter emission

verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper argues that the puzzling X-ray plateau phase seen in many gamma-ray burst afterglows is not a sign of prolonged engine activity or exotic geometry, but the natural hydrodynamics of an outflow whose Lorentz factor is spread over a continuous range. It models the ejecta as a power-law mass distribution and shows that the forward shock driven into the external medium then decays slowly, matching observed plateau slopes, durations, spectra, and the luminosity–duration anti-correlation. The reverse shock crossing time sets the plateau end, and once the slowest material is processed the blast wave smoothly joins the standard self-similar afterglow. If right, the same outflow that made the prompt gamma-rays powers the entire afterglow, with a bright millimeter reverse-shock component as a distinctive signature.

Core claim

The central claim is that a radially stratified ultra-relativistic outflow whose cumulative mass above Lorentz factor γ falls as M(>γ) ∝ γ^{-s}, with s ≈ 2.5–4 and a minimum Lorentz factor γ_min ≈ 70–100, produces the observed X-ray plateau through its forward-shock emission. The plateau slope fixes s, the plateau duration fixes γ_min, and both inferred values put the electron energy fraction, magnetic fraction, and energetics in standard ranges. The reverse shock in the same ejecta is long-lived and mildly relativistic, and its synchrotron emission peaks in the millimeter band, outshining the forward shock there. When the reverse shock finishes crossing the slowest ejecta, both components f

What carries the argument

The central object is the cumulative mass profile M_iso(>γ_4) = M_iso (γ_4/γ_min)^{-s} for the ejecta, combined with an analytic forward–reverse shock solution in which the ratio x = γ_4/γ_2 of unshocked-ejecta to shocked-plasma Lorentz factor depends only on s. This ratio controls the mildly relativistic reverse shock, its crossing time (roughly an order of magnitude longer than the thin-shell estimate), and the resulting synchrotron break frequencies and fluxes of both shocks. The crossing-time expression sets the plateau duration, and the analytic closure relations link observed temporal and spectral slopes to s and the electron power-law index p.

Load-bearing premise

The load-bearing premise is that the cumulative ejecta mass follows a single power law in Lorentz factor, M(>γ) ∝ γ^{-s}, with the same exponent over the whole range the reverse shock processes; if the true profile bends or s varies, the inferred γ_min and the interpretation of the plateau break as the reverse-shock crossing time lose their foundation.

What would settle it

A GRB with a canonical X-ray plateau but no bright millimeter excess during the plateau, or a millimeter signal that dies well before the plateau break, would falsify the model. Alternatively, direct measurement of the ejecta Lorentz-factor distribution through early radio scintillation or the deceleration onset could test the required γ_min ~ 70–100.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

X-ray plateaus require ejecta with Lorentz factors extending down to ~70–100; steeper or shallower profiles would produce either impulsive or too-slow decay.
The plateau break marks the end of reverse-shock energy transfer, so the post-plateau decay should be a smooth transition to the standard afterglow, with little spectral evolution.
The same outflow accounts for prompt emission, plateau, and afterglow, avoiding the energy-budget tensions of refreshed-shock models.
Plateaus should be accompanied by a long-lived, bright millimeter component that outshines the forward shock at those wavelengths and fades only after the reverse shock crosses the slowest ejecta.
The model reproduces the observed anti-correlation between plateau luminosity and duration, with scatter arising from variations in electron index and stratification.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the identification of the plateau break with reverse-shock crossing is correct, multi-band observations during the plateau should show simultaneous timing between the X-ray break and the turnover of the millimeter component; a measurable lag would point to a different mechanism.
The inferred γ_min ≈ 70–100 sits near the minimum Lorentz factor needed for gamma-ray escape in long bursts, possibly connecting the plateau to the prompt-emission mechanism; the paper does not pursue this link.
The power-law profile is assumed constant, but a broken power-law or lognormal ejecta distribution could change the inferred parameters; the analytic structure should still apply locally, which numerical simulations could test.
Millimeter monitoring of nearby plateau bursts within the first hour could distinguish stratified-outflow plateaus from magnetar or geometry models, which do not predict a sustained reverse-shock millimeter excess.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Desk Editor's Note

Careful analytic framework for stratified ejecta; the plateau application is plausible but not yet a quantitative fit, and the mm prediction is the part worth taking seriously.

read the letter

The thing to know: this paper gives you a self-contained analytic treatment of a radially stratified ultra-relativistic outflow interacting with a uniform or wind-like medium, and it ends with a sharp, testable prediction—bright, long-lived millimeter emission from the reverse shock during X-ray plateaus. The hydrodynamics is the real contribution. The observational application is illustrative, not a rigorous fit.

What is genuinely new: the closed-form relation x(s) for the Lorentz-factor ratio across the shock system (Eq. 22), which depends only on the stratification index s, and the crossing-time correction factors h_gamma and h_a that can stretch the reverse-shock crossing time by an order of magnitude relative to thin-shell estimates. The emission scalings overlap with Sari & Meszaros 2000, and the authors say so. That is honest. The framework is internally consistent, recovers the expected limits, and the derivation shows care—the variable adiabatic index, the self-absorption treatment, and the closure relations in Table 1 are all useful. The paper also openly flags in Section 5 that real ejecta profiles are not expected to be single power laws, which is more transparency than many papers offer.

Where it gets softer: the claim of consistency with plateau observations is obtained by reading s off the observed slope and gamma_min off the observed duration, then showing the same quantities are reproduced. That is parameter inference, not a fit. No full sample, no error budget, no model comparison against energy injection or magnetar alternatives. The Dainotti-relation agreement reduces mostly to F_X ∝ E_iso/t_cross, a scaling rather than a distinctive prediction. And the single-power-law assumption is load-bearing for both gamma_min ~ 70–100 and the identification of the plateau break with the reverse-shock crossing time. If the ejecta profile is broken, those numbers shift. The zeta = 2 choice also propagates into every normalization. These are not fatal flaws, but they keep the paper at the level of a plausible framework rather than a demonstrated explanation.

The mm prediction is the part I would actually bet on being important. If real observations show a bright, long-lived mm excess during plateaus, this model becomes a leading contender. The note about GRB 260127A is suggestive but not yet a quantitative test.

Who is this for: GRB afterglow theorists and observers planning mm follow-up. It deserves a serious referee. Send it to review, with the expectation of revision to add caveats and ideally a real sample comparison.

Referee Report

3 major / 4 minor

Summary. The paper develops an analytic framework for the forward–reverse shock system in ultra-relativistic, radially stratified outflows with a cumulative ejecta mass profile M(>γ) ∝ γ^{-s}. It derives closed-form expressions for the shock Lorentz factors, the reverse-shock crossing time, and the synchrotron emission from both shocks. The authors argue that for s ≈ 2.5–4 and γ_min ≈ 70–100 the forward-shock X-ray emission reproduces the observed plateau decay slope, duration, flux, spectral ordering, and the Dainotti anti-correlation, while the reverse shock produces a long-lived millimeter excess. The paper identifies the end of the plateau with the reverse-shock crossing time, after which the flow relaxes to the Blandford–McKee self-similar solution.

Significance. If the central claims hold, the paper offers a physically economical, single-outflow explanation for X-ray plateaus that avoids late-time central-engine activity and preserves a reasonable prompt radiative efficiency. The analytic framework is a genuine contribution: it provides compact, closed-form expressions that recover the known single-shell and BM limits, and it was previously benchmarked against relativistic hydrodynamic simulations by the same group. The model also makes a falsifiable prediction—a bright, long-lived millimeter counterpart during the plateau—which is particularly valuable given the growing availability of fast mm follow-up. However, the observational 'consistency' is currently demonstrated through parameter calibration rather than a rigorous statistical comparison, and the inference of γ_min and s rests on a single-power-law ejecta profile that the paper itself concedes is an idealization.

major comments (3)

[§4.1–4.2 and §4.4] The consistency claim is partly circular. The observed plateau slope is used to select s (Eq. 56), and the observed plateau duration is used to solve for γ_min via Eq. (27)/(57). The same set of equations is then used in §4.4 to state that the model reproduces the plateau flux and the Dainotti relation (Eq. 63). The agreement is therefore a calibration, not an independent test. To support the 'origin of GRB X-ray plateaus' claim, the paper needs a forward-model fit to a sample of plateaus with error propagation and, ideally, a model comparison against the standard energy-injection scenario.
[§2.2, Eq. (9), §4.2, §5] The single-power-law profile M(>γ) ∝ γ^{-s} with constant s is load-bearing for the inferred γ_min ≈ 70–100 and for identifying the plateau break with t_cross (Eq. 27). If the true profile is a broken power law whose local slope varies as the reverse shock progresses, the slope inferred from the plateau phase (via the closure relations of Table 1) need not equal the effective s that enters t_cross. The paper itself states in §5 that the profile 'is not expected to follow a single power law over the entire LF range,' yet the analytic framework and all parameter inferences assume exactly that. The authors should quantify the sensitivity of γ_min and t_cross to a varying s, e.g., by considering a broken power law or by testing reconstruction against simulated stratified ejecta.
[§4.4, Eq. (63)] The derivation of the Dainotti relation as F_X ∝ E_iso / t_cross is a useful scaling, but the claimed scatter and normalization are not compared directly to the observed Dainotti data. The observed relation involves luminosities and break times with substantial scatter and selection effects; showing that a few fiducial parameter choices produce an anti-correlation is not yet a quantitative consistency test. A direct comparison to the published Dainotti samples, including the best-fit slope and scatter, would strengthen the claim.

minor comments (4)

[General] There are multiple typographical issues: 'Y amazaki2009' in the Introduction, inconsistent use of γ_f/γ_2 for the same quantity, and a garbled label in Fig. 7 ('FX/t^{-0.46}'). These should be corrected in a final pass.
[Eq. (20)] The parentheses in the third term of Eq. (20) appear unbalanced, making the expression hard to follow.
[Table 1] The table is dense and would benefit from a derivation guide or a short worked example of how s is mapped to α_X for a given p and k; currently the reader must trust the closure relations without a checkable example.
[§4.2, wind case] The wind case is dismissed quickly after deriving γ_min ~ 50. Given that many GRB environments are wind-like, the authors should either provide a more detailed discussion of why the ISM case is preferred observationally, or present the wind-case predictions alongside the ISM case.

Circularity Check

0 steps flagged

Parameter inference is explicitly stated; the model's independent prediction (mm reverse-shock emission) is not fitted. No circular reduction identified.

full rationale

The derivation chain is not circular in a way that invalidates the central claim. The free parameters s, gamma_min, and epsilon_e are explicitly inferred from the observed plateau slope (Sec. 4.1, Eq. 56), duration (Sec. 4.2, Eqs. 57-58), and flux (Sec. 4.4, Eq. 62), so the statements that the model 'naturally produces' a shallow X-ray decay are consistency checks after parameter inversion, not independent predictions smuggled in as outputs. The genuinely independent, falsifiable prediction is the long-lived millimeter reverse-shock component (Sec. 4.6, Figs. 6-7), which is not used to fit any parameter and is only compared to external data in the note added during review. The claimed consistency with the Dainotti relation (Sec. 4.4, Eq. 63) reduces to F_X ~ E_iso/t_cross; this is a relatively generic consequence of a finite energy reservoir released over the crossing time and is not fed back as an input, since Dainotti et al. is cited only as an external empirical relation. The paper's own caveat that the ejecta profile 'is not expected to follow a single power law over the entire LF range' (Sec. 5) is a limitation of the single-power-law ansatz in Eq. (9), but the inference of gamma_min ~ 70-100 from Eq. (27) is an honest, stated inversion of the model, not a hidden equivalence. Self-citations (Sadeh et al. 2023; Sadeh 2024; Sadeh & Waxman 2025) are used as numerical verification and as the starting point for the stratified treatment; they are not invoked as an unverified uniqueness theorem and do not carry the argument alone. The central quantitative claims—plateau duration scaling, spectral ordering, and mm-band brightness—are derived algebraically from stated assumptions rather than imported from prior work. No step was found where an equation equals its own input by definition. Score 2 reflects minor self-citation dependence and the presentation of post-fit consistency as 'natural,' but no substantive circularity.

Axiom & Free-Parameter Ledger

8 free parameters · 7 axioms · 0 invented entities

The central model rests on a power-law ejecta mass distribution and standard afterglow microphysics. The two key numbers s and gamma_min are fitted to plateau phenomenology. No new particles, fields, or dimensions are introduced.

free parameters (8)

s (ejecta stratification index) = ~2.5-4 (ISM case)
Power-law index of cumulative ejecta mass M(>gamma) ~ gamma^{-s}; inferred by matching observed plateau slopes alpha_X ~ 0.3-0.6 in Section 4.1; controls both slope and crossing time.
gamma_min (minimum Lorentz factor) = ~70-100 (ISM typical parameters)
Inferred from matching observed plateau durations t_b ~ 10^3-10^5 s via Eq. (57) with assumed E_iso ~ 10^53 erg, n ~ 10-100 cm^-3, z ~ 1.5-2.5; a central parameter of the model.
epsilon_e (electron energy fraction) = 0.05-0.5
Chosen to match observed plateau flux levels of ~0.1-10 microJy in Section 4.4.
epsilon_B (magnetic energy fraction) = >=10^-3; fiducial 10^-2
Lower bound set by requiring nu_X > max(nu_m, nu_c) in Section 4.3; value influences all synchrotron normalizations.
p (electron power-law index) = 2.0-2.4
Taken from observed X-ray spectral indices beta_X ~ 1.0-1.2; used in closure relations and flux normalization.
n (ISM number density) = n_-2 ~ 10-100
Assumed 'typical' long-GRB values; directly affects inferred gamma_min, flux, and frequencies.
E_iso (isotropic-equivalent energy) = 10^52.5-10^53.5 erg
Assumed typical GRB energetics; enters crossing time and flux scalings.
zeta (time-radius relation factor) = 2
Adopted as an intermediate value between coasting (zeta=1) and BM (zeta=4); affects all normalizations and inferred gamma_min, with no sensitivity analysis.

axioms (7)

standard math Relativistic shock jump conditions and energy-momentum conservation (Eqs. 1, 11-14) govern the forward/reverse shock system.
Standard relativistic hydrodynamics/afterglow framework; not in dispute.
domain assumption External medium density profile is a power law rho_1 = A r^{-k} with k < 3; the paper focuses on ISM k=0.
Used in Eqs. 24-28 and all temporal scalings; standard but an assumption about the environment.
ad hoc to paper Ejecta cumulative mass above Lorentz factor gamma_4 follows M_iso(>gamma_4) = M_iso (gamma_4/gamma_min)^{-s} with constant s > 1 and a sharp lower cutoff.
This is the central model ansatz (Section 2.2, Eq. 9); the plateau slope is used to infer s and the duration to infer gamma_min, so the claim and the ansatz are partly circular.
domain assumption The shocked regions have approximately uniform flow profiles and the reverse shock evolves quasi-statically.
Underlies Eqs. (4)-(22); the paper cites numerical validation for single shells (Sadeh et al. 2023; Sadeh & Waxman 2025) but does not re-verify the stratified case with new simulations.
domain assumption Standard synchrotron afterglow microphysics: fixed epsilon_e and epsilon_B, power-law electron distribution, adiabatic blast wave, and the same microphysical fractions in forward and reverse shocks.
Assumed throughout Section 3; no alternative radiation mechanism is considered.
domain assumption The plateau break time equals the reverse-shock crossing time, after which the flow joins the Blandford-McKee self-similar solution.
Makes t_cross the observable t_b in Section 4.2; this identification is central and not directly measured.
domain assumption The observer-time / radius relation uses zeta = 2 as a constant intermediate deceleration factor.
The paper states zeta=2 is adopted because the deceleration is intermediate; all normalization constants and inferred gamma_min depend on it.

pith-pipeline@v1.3.0-alltime-deepseek · 25205 in / 15983 out tokens · 142842 ms · 2026-08-02T20:02:52.676961+00:00 · methodology

0 comments

read the original abstract

The origin of the X-ray plateau phase observed in a large fraction of gamma-ray burst afterglows remains debated. We present a novel analytic framework for the hydrodynamics of ultra-relativistic, radially stratified outflows interacting with an external medium. By explicitly accounting for a continuous distribution of Lorentz factors within the ejecta, we derive analytic expressions describing the evolution of a long-lived, mildly relativistic reverse shock and determine its crossing time. Then, we compute the resulting synchrotron emission from both the forward and reverse shocks. The forward shock naturally produces a shallow, long-lasting X-ray decay consistent with the observed properties of X-ray plateaus, including the Dainotti relation, without requiring prolonged central-engine activity or an additional high-energy emission component. We further show that reproducing the observed plateau durations requires a broad distribution of ejecta Lorentz factors, extending down to $\gamma_\text{min}\sim70-100$, consistent with the ultra-relativistic outflow that powers the prompt $\gamma$-ray emission. The reverse shock generates a long-lived millimeter emission component that outshines the forward shock emission at these wavelengths. Both the plateau and reverse shock emission terminate smoothly once the slowest ejecta are processed, marking a transition to the standard Blandford-McKee self-similar evolution. Such stratified outflows are expected on physical grounds, as the ultra-relativistic ejecta responsible for the prompt $\gamma$-ray emission are unlikely to be launched with a single Lorentz factor. This model provides a unified picture in which the same outflow powers the prompt emission, the X-ray plateau, and the subsequent afterglow evolution.

Figures

Figures reproduced from arXiv: 2602.24085 by Gilad Sadeh, Kenta Hotokezaka, Masaru Shibata.

**Figure 1.** Figure 1: A schematic illustration of the forward-reverse shock structure, showing the four dynamical regions: (1) unshocked external medium, (2) shocked external medium, (3) shocked ejecta, and (4) unshocked ejecta. The pressure in the unshocked regions is negligible compared to that in the shocked regions, 𝑝1, 𝑝4 ≪ 𝑝2, 𝑝3. In our analytic treatment, the pressure and velocity are approximated as uniform throughout … view at source ↗

**Figure 2.** Figure 2: Exact analytic solutions for the LF ratio, 𝛾4/𝛾2, and for the shocked-ejecta LF, 𝛾¯3 (measured in the unshocked ejecta frame), shown as functions of the single dimensionless parameter 𝛾4/ p 𝑓 (with 𝑓 ≡ 𝜌4/𝜌1). Also shown are the asymptotic limits: the UR RS regime (𝛾 2 4 ≫ 𝑓 , Eq. (7)) and the Newtonian RS regime (𝛾 2 4 ≪ 𝑓 , Eq. (8)). leads to 𝛾2 ≈ 𝛾4 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 4.** Figure 4: In Blue: numerical solution of Eq. (22) showing the dependence of the LF ratio 𝑥 ≡ 𝛾4/𝛾2 on the ejecta power-law index 𝑠. In red: the reverse shock LF in the frame of the unshocked ejecta, 𝛾¯3 = 1 2 𝑥 + 1 𝑥 . In yellow: the ratio of the swept-up external medium mass, 𝑀ext, to the canonical deceleration mass, 𝑀iso (> 𝛾4 )/𝛾4 (Rees & Meszaros 1992). In green: The ratio between the shocked external mediu… view at source ↗

**Figure 5.** Figure 5: The corresponding correction factors for the RS crossing time, ℎ𝛾 (uniform ISM; 𝑘 = 0) and ℎ𝑎 (wind environment; 𝑘 = 2), are shown. For moderate values 2 ≲ 𝑠 ≲ 4, relevant when the RS propagates through the shallow part of the ejecta, the RS is mildly relativistic and the crossing time is significantly extended relative to the standard thin-shell estimates. first derived by Sari & Mészáros (2000). For brev… view at source ↗

**Figure 6.** Figure 6: illustrates a representative broadband spectrum of the FS and RS at 𝑡 ≃ 1 hr, while [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗

**Figure 7.** Figure 7: Representative synchrotron X-ray, optical and mm light curves for fiducial parameters 𝑝 = 2.1, 𝜀𝑒 = 0.1, 𝜀𝐵 = 10−2 , 𝑛 = 10−1 cm−3 , 𝐸 = 1053 erg, 𝛾min = 100, 𝑠 = 3, and 𝑧 = 2. As shown in [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗

Review history (2 revisions) →

discussion (0)

Reference graph

Works this paper leans on

2 extracted references · 1 linked inside Pith

[1]

D., McKee C

Beniamini P ., Nakar E., 2019, Monthly Notices of the Royal Astronomical Society, 482, 5430 Beniamini P ., Duque R., Daigne F., Mochkovitch R., 2020a,Monthly Notices of the Royal Astronomical Society, 492, 2847 Beniamini P ., Granot J., Gill R., 2020b, Monthly Notices of the Royal Astro- nomical Society, 493, 3521 Berger E., 2014, Annual Review of Astrono...

Pith/arXiv arXiv 2019
[2]

Thus, even though the photons emitted from the shocked ejecta go through the shocked external medium, they are not affected by its self-absorption

For the relevant parameter space, 𝜈𝑟 𝑎 ≫ 𝜈 𝑓 𝑎 . Thus, even though the photons emitted from the shocked ejecta go through the shocked external medium, they are not affected by its self-absorption. APPENDIX B: GENERAL SYNCHROTRON SCALINGS Here, we summarize the temporal scalings of the synchrotron observ- ables for a general external-density profile 𝜌1 = 𝐴...

2026

[1] [1]

D., McKee C

Beniamini P ., Nakar E., 2019, Monthly Notices of the Royal Astronomical Society, 482, 5430 Beniamini P ., Duque R., Daigne F., Mochkovitch R., 2020a,Monthly Notices of the Royal Astronomical Society, 492, 2847 Beniamini P ., Granot J., Gill R., 2020b, Monthly Notices of the Royal Astro- nomical Society, 493, 3521 Berger E., 2014, Annual Review of Astrono...

Pith/arXiv arXiv 2019

[2] [2]

Thus, even though the photons emitted from the shocked ejecta go through the shocked external medium, they are not affected by its self-absorption

For the relevant parameter space, 𝜈𝑟 𝑎 ≫ 𝜈 𝑓 𝑎 . Thus, even though the photons emitted from the shocked ejecta go through the shocked external medium, they are not affected by its self-absorption. APPENDIX B: GENERAL SYNCHROTRON SCALINGS Here, we summarize the temporal scalings of the synchrotron observ- ables for a general external-density profile 𝜌1 = 𝐴...

2026