High-Resolution Numerical Calculations of GRB Afterglow Plateaus arising from Stratified Ejecta

Paul C. Duffell; Ranadeep Ghosh Dastidar

Stratified GRB ejecta alone can produce afterglow plateaus in wind and ISM media, needing only modest Lorentz factors γ₀ ~ 10–50.

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 · grok-4.5

2026-07-15 01:56 UTC pith:WCOA4MLJ

load-bearing objection High-res numerical check that stratified ejecta can make GRB plateaus in both media, with useful γ₀ scalings; idea is old, contribution is computational verification. the 3 major comments →

arxiv 2607.12967 v1 pith:WCOA4MLJ submitted 2026-07-14 astro-ph.HE

High-Resolution Numerical Calculations of GRB Afterglow Plateaus arising from Stratified Ejecta

Paul C. Duffell , Ranadeep Ghosh Dastidar This is my paper

classification astro-ph.HE

keywords gamma-ray burstsafterglow plateausstratified ejectareverse shockrelativistic hydrodynamicsLorentz factorwind mediumISM

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 long plateaus seen in gamma-ray burst afterglows do not require continuous late-time energy injection. Instead, a stratified ejecta structure—already expected in GRB jets—is enough. In that picture the plateau is simply the early phase before the reverse shock finishes crossing the outflow, during which the blast wave either coasts or decelerates only very slowly. High-resolution relativistic hydrodynamics (Δr/r ≲ 10^{-5}) plus radiative post-processing show that the same mechanism works in both wind-like and constant-density environments. The duration of a typical plateau then fixes a surprisingly low characteristic Lorentz factor for the ejecta, γ₀ ~ 10–50, and the observed slope of the plateau directly constrains how the fastest-moving material is stratified. If correct, the result supplies a clean diagnostic of the outer ejecta that may also be responsible for the prompt emission.

Core claim

A stratified ejecta profile is sufficient to generate GRB afterglow plateaus in both wind (k=2) and ISM (k=0) media; the long observed plateau durations require only a modest characteristic Lorentz factor γ₀ ~ 10–50, and the plateau slopes themselves can be used to constrain the ejecta stratification.

What carries the argument

Extreme-resolution relativistic hydrodynamics of a stratified outflow (JET code, Δr/r ≲ 10^{-5}) followed by afterglow light-curve post-processing (Firefly). The resolution is required to evolve the tiny ultra-relativistic outer mass that sets the reverse-shock crossing time and therefore the plateau duration and slope.

Load-bearing premise

The adopted stratified density profiles and the extreme but finite numerical resolution are assumed to be adequate to capture the dynamics of the minute ultra-relativistic outer shell, so that the measured plateau durations and slopes are physical rather than numerical artifacts.

What would settle it

A high-resolution calculation with the same stratified profiles that fails to produce a plateau of the observed duration and slope, or an independent afterglow observation whose plateau duration and slope cannot be matched by any γ₀ in the 10–50 range under the reported scalings.

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

If this is right

Stratified ejecta become a viable default explanation for afterglow plateaus, reducing the need for continuous late-time energy injection.
Observed plateau slopes furnish direct constraints on the radial stratification of the fastest GRB ejecta.
Typical plateau lengths imply characteristic Lorentz factors γ₀ ~ 10–50 for the ejecta.
The same mechanism operates in both wind-like and constant-density circumburst media.
Measured scalings of plateau duration versus γ₀ can be used as quantitative diagnostics for future afterglow modeling.

Where Pith is reading between the lines

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

If the outer stratified layers also power the prompt emission, the same γ₀ ~ 10–50 range would apply to the prompt-emitting material, tightening multi-messenger constraints.
Plateau slope versus duration diagrams for large GRB samples could be inverted to map the distribution of ejecta stratification indices.
The required resolution Δr/r ≲ 10^{-5} sets a practical benchmark for any future numerical study that claims to resolve reverse-shock crossing in stratified jets.

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

High-res numerical check that stratified ejecta can make GRB plateaus in both media, with useful γ₀ scalings; idea is old, contribution is computational verification.

read the letter

The one thing to know: this is a high-resolution numerical verification that stratified ejecta alone can produce afterglow plateaus in both ISM (k=0) and wind (k=2) media, plus explicit measured scalings of plateau duration with characteristic Lorentz factor γ₀ ~ 10–50. The mechanism itself is not new—the abstract says so and cites prior suggestions of pre-reverse-shock-crossing plateaus—but the computational product is.

What they do well is state the resolution problem clearly (Δr/r ≲ 10^{-5} needed for the tiny ultra-relativistic outer mass) and then actually run it with JET + Firefly. That is real work. Getting clean break-time scalings versus γ₀ in both ambient media is useful; if those numbers hold up, people can use plateau slopes as a diagnostic of ejecta stratification without automatically invoking late continuous energy injection. Circular fitting burden looks low: γ₀ is a physical input whose range is compared to observations, not a free parameter tuned to force a plateau.

Soft spots are mostly documentation, not logic. We only have the abstract, so we cannot audit convergence, the exact stratification functional forms, microphysical parameters, or post-processing assumptions. The extreme-resolution claim is load-bearing for the tiny outer mass that sets reverse-shock crossing time; without resolution studies or code/data, that remains unchecked. That is a real gap for reproducibility, not an internal contradiction. On the text we have, the central claim is coherent and the free parameters (stratification profile, γ₀, k) are the expected physical ones.

This is for people who model GRB afterglows and care whether plateaus require ongoing energy injection. A serious referee should see the full paper, figures, and resolution tests. I would accept it for peer review; the numerical verification and scalings are worth the referee time even if the idea is already in the literature. Bring it to reading group only if someone is actively working on reverse-shock or plateau diagnostics.

Referee Report

3 major / 3 minor

Summary. The manuscript claims that high-resolution relativistic hydrodynamics (Δr/r ≲ 10^{-5}) of stratified GRB ejecta, evolved with the JET code and post-processed with Firefly, produce afterglow plateaus prior to reverse-shock crossing in both ISM (k=0) and wind (k=2) media. Plateau durations are reported to require only a modest characteristic Lorentz factor γ₀ ∼ 10–50, and plateau slopes are said to constrain the ejecta stratification. The work is presented as numerical verification that stratified ejecta—expected to be generic in GRB jets—are sufficient to explain observed plateaus without continuous late-time energy injection.

Significance. If the numerical results and scalings hold, the paper would supply concrete computational support for a long-standing alternative to continuous energy injection as the origin of Swift-era afterglow plateaus. Explicit break-time measurements and γ₀–duration scalings would be useful for population interpretation and for linking afterglow structure to the prompt-emitting ejecta. The extreme-resolution hydrodynamics and dual-code (JET + Firefly) pipeline are methodological strengths, provided they are accompanied by documented convergence tests and reproducible setups.

major comments (3)

Only the abstract is available for this review. The load-bearing numerical claim—that Δr/r ≲ 10^{-5} is both necessary and sufficient to evolve the ultra-relativistic outer mass without artificial diffusion that would erase or distort the plateau—cannot be audited without resolution studies, convergence plots, and the measured plateau durations/slopes. Until those appear in the full manuscript, the central result remains unverified.
The abstract asserts that plateau slopes constrain ejecta stratification and that stratified structure is generically sufficient in both k=0 and k=2 media. Without the functional forms of the adopted density/velocity profiles, the mapping from slope to stratification index, and the explicit break-time measurements, it is not possible to assess whether those conclusions are robust or profile-dependent.
Microphysical parameters (ε_e, ε_B, p, electron participation fraction) and the precise definition of the characteristic Lorentz factor γ₀ are not stated in the abstract. Because afterglow light-curve slopes and break times depend on these choices, their values and any sensitivity tests are load-bearing for the claim that the plateaus match observed phenomenology.

minor comments (3)

The abstract would be clearer if it named the stratification functional form(s) used (e.g., power-law index in mass or energy vs. Lorentz factor) rather than referring only generically to “stratified ejecta structure.”
A brief statement of whether the reported γ₀ ∼ 10–50 scalings are analytic expectations confirmed numerically, or purely numerical fits, would help readers place the result relative to prior analytic work on reverse-shock crossing plateaus.
Code and data availability (JET/Firefly versions, input decks, resolution-study outputs) should be stated explicitly once the full manuscript is provided, given the extreme resolution claimed.

Circularity Check

0 steps flagged

No significant circularity: numerical verification of a previously suggested stratified-ejecta mechanism against external afterglow phenomenology.

full rationale

The abstract-only paper reports high-resolution relativistic hydrodynamics (JET) plus afterglow post-processing (Firefly) showing that stratified ejecta generically produce plateaus before reverse-shock crossing in both k=0 and k=2 media. The characteristic Lorentz factor γ₀ is treated as a physical input whose modest range (~10–50) is required to match observed plateau durations; it is not fitted so that the plateau is forced by construction. Plateau slopes are presented as diagnostics of ejecta stratification rather than as tautological restatements of the input profiles. No equations, uniqueness theorems, or load-bearing self-citations appear that would reduce the claimed result to its inputs. Minor use of the authors’ own codes is ordinary numerical practice and does not constitute circularity under the stated criteria. The derivation chain is therefore self-contained numerical confirmation of an independent physical mechanism.

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 0 invented entities

The central claim rests on standard relativistic hydrodynamics and afterglow radiation theory plus an assumed stratified ejecta profile whose functional form and characteristic Lorentz factor are inputs. No new particles or forces are introduced. Free parameters are the stratification law and γ₀ (and standard microphysical ε_e, ε_B, p if used in Firefly), which are not fitted to force the plateau but are scanned to match observed durations.

free parameters (3)

characteristic Lorentz factor γ₀ = ~10–50
Scanned to match observed long plateau durations; abstract reports the required range γ₀ ~ 10–50.
ejecta stratification profile
The radial/velocity structure of the fastest ejecta is an assumed input that sets the reverse-shock crossing history and therefore the plateau slope; functional form not specified in the abstract.
ambient density index k = 0 or 2
Environment is set to either ISM (k=0) or wind (k=2); discrete choice, not a continuous fit, but still an external modeling assumption.

axioms (3)

domain assumption Relativistic hydrodynamics of a cold baryonic jet interacting with a power-law ambient medium correctly describes the reverse-shock crossing phase.
Invoked throughout; the plateau is identified with the pre-crossing or slow-crossing phase of the stratified ejecta.
domain assumption Standard synchrotron afterglow emission (Firefly post-processing) maps the hydrodynamical evolution to the observed light curve without additional late energy injection.
Required to convert JET hydro output into the claimed afterglow plateau.
ad hoc to paper Δr/r ≲ 10^{-5} is sufficient numerical resolution to evolve the ultra-relativistic outer mass without artificial diffusion that would erase the plateau.
The abstract’s central technical claim; convergence is not demonstrated in the abstract itself.

pith-pipeline@v1.1.0-grok45 · 6212 in / 2608 out tokens · 31608 ms · 2026-07-15T01:56:49.029086+00:00 · methodology

0 comments

read the original abstract

Since the discovery of plateaus in GRB afterglows by Swift, they have been modeled predominantly by late-time energy injection. However, many studies have suggested that the plateau may be modeled by an early phase before reverse shock crossing (either coasting with constant Lorentz factor or decelerating very slowly as the reverse shock crosses the ejecta). The slope of the early plateau provides some constraints the stratification of the fastest-moving ejecta, which could be the ejecta responsible for the prompt emission. However, numerical studies typically do not model the jet as a stratified outflow; the reason being the extremely high resolution required in order to accurately evolve this tiny amount of highly relativistic material. In this study, we perform high-resolution numerical calculations ($\Delta r/r \lesssim 10^{-5}$) verifying that a stratified ejecta structure can indeed produce an afterglow plateau in both wind ($k=2$) and ISM ($k=0$) environments, and computing break times explicitly. We evolve the relativistic hydrodynamics using the JET code, and post-process this to compute the afterglow using the Firefly code. Our results show that a stratified ejecta structure (which should generically be present in GRB jets) is sufficient to explain GRB afterglows, and the plateau slopes can be used to constrain the ejecta stratification. We additionally provide precise measured scalings for the plateau duration as a function of the characteristic Lorentz factor of the ejecta. The long duration of typical plateaus requires a very modest characteristic Lorentz factor for the ejecta ($\gamma_0 \sim 10-50$), in agreement with other afterglow plateau models.

High-Resolution Numerical Calculations of GRB Afterglow Plateaus arising from Stratified Ejecta

Core claim

What carries the argument

Load-bearing premise

What would settle it

If this is right

Where Pith is reading between the lines

discussion (0)