pith. sign in

arxiv: 2511.08684 · v2 · pith:IXBL5TTJnew · submitted 2025-11-11 · 🌌 astro-ph.HE

Radiation-mediated shocks in gamma-ray bursts: spectral evolution

Filip Alamaa , Fr\'ed\'eric Daigne This is my paper

Pith reviewed 2026-05-21 18:19 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords gamma-ray burstsradiation-mediated shocksprompt emissionspectral evolutioninternal shocksrelativistic hydrodynamicsreverse shock
0
0 comments X

The pith

Radiation-mediated shocks from an internal collision in a GRB jet produce a short pulse that radiates 23 percent of the energy while the spectrum softens and develops a high-energy power-law tail.

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

The paper simulates a specific internal collision inside a gamma-ray burst jet with one-dimensional special relativistic hydrodynamics. It then tracks the photon distributions through the resulting forward and reverse radiation-mediated shocks and the shared downstream region all the way above the photosphere. For the chosen parameters and redshift of one, the calculation yields a single pulse lasting roughly 0.1 seconds that carries about 23 percent of the total burst energy into radiation. The spectrum settles into a generic shape featuring smooth curvature below a peak energy that steadily drops from 250 keV to 100 keV, plus a high-energy power law generated at low optical depth by the reverse shock.

Core claim

In the modeled internal collision, the forward and reverse radiation-mediated shocks together produce a light curve consisting of a short pulse of duration approximately 0.1 s at z=1. Roughly 23 percent of the total energy is radiated. After early times the spectrum develops a smooth curvature below the peak energy Ep together with a high-energy power-law tail that cuts off near 5 MeV in the observer frame; Ep decreases from about 250 keV to 100 keV during the decay while the high-energy index falls from roughly -2 to -2.4, the latter component arising from the reverse shock at optical depths below 30.

What carries the argument

Radiation-mediated shocks in the forward and reverse shocks of a relativistic internal collision, with photon distributions evolved by a dedicated simulation code from the shocks through the photosphere.

If this is right

  • Early radiation-mediated shocks can form short pulses that either constitute entire short GRBs or act as building blocks inside longer, highly variable light curves.
  • The radiative efficiency reaches approximately 23 percent for the modeled collision.
  • The time-resolved spectrum narrows and softens with time, with Ep decreasing steadily during the decay phase.
  • A high-energy power-law component with index between -2 and -2.4 is produced at low optical depth by the reverse shock.
  • The low-energy spectral index decreases from roughly -0.5 to -1 during the bright part of the pulse.

Where Pith is reading between the lines

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

  • Similar internal collisions at other Lorentz factors or collision radii could produce a range of pulse durations and spectral evolution patterns observed across the GRB population.
  • The model predicts that the high-energy tail should be strongest when the reverse shock is viewed at low optical depth, offering a possible observational signature to separate forward- and reverse-shock contributions.
  • Extending the calculation to include angular dependence or multiple successive collisions could test whether the reported spectral features persist in more complex light curves.

Load-bearing premise

The results rest on one particular set of internal-collision parameters together with the accuracy of the one-dimensional hydrodynamics and radiation-mediated-shock photon-tracking code.

What would settle it

High-time-resolution observations of a short GRB or an isolated pulse in a longer burst that show a pulse duration near 0.1 s, radiative efficiency near 23 percent, a peak energy falling from 250 keV to 100 keV, and a high-energy spectral index softening from -2 to -2.4 with a cutoff near 5 MeV would support the model; a clear mismatch in any of these quantities would falsify it for this setup.

Figures

Figures reproduced from arXiv: 2511.08684 by Filip Alamaa, Fr\'ed\'eric Daigne.

Figure 1
Figure 1. Figure 1: A flowchart showing the methodology of the paper. Each box shows one step of the chain, with the treated physics [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Top panel: Sketch of the Lorentz factor profile across the ejecta with the positions of the five different regions mentioned in subsection 2.4 marked. The dashed vertical lines indicate the position of each region relative to the KRA zones given in the bottom panel. Bottom panel: Ge￾ometry of the KRA with both reverse and forward shocks included. Green, red, and purple color indicate the up￾stream zones, t… view at source ↗
Figure 3
Figure 3. Figure 3: The Lorentz factor (left) and the comoving density (right) at different optical depths as obtained from the [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Comoving photon distribution in the five different zones as obtained by [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Observed spectral evolution (top left), light curve (top right), count spectrum of a Band fit to the time-integrated [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of a Band function fit (top) and a [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
read the original abstract

Radiation-mediated shocks (RMSs) occurring below the photosphere in a gamma-ray burst (GRB) jet could play a crucial role in shaping the prompt emission. In this paper, we study the time-resolved signal expected from such early shocks. We model an internal collision using a 1D special relativistic hydrodynamical simulation, and we follow the photon distributions in the resulting forward and reverse shocks as well as in the common downstream region to well above the photosphere using a designated RMS simulation code. We compute the light curve and time-resolved spectrum of the resulting single pulse taking into account the emission at different optical depths and angles to the line of sight. For the specific case considered, we find a light curve consisting of a short pulse lasting $\sim 0.1~$s for an assumed redshift of $z = 1$, which could constitute a whole short GRB or be a building block within a highly variable longer GRB light curve. The efficiency is large, with $\approx 23$% of the total burst energy being radiated. The spectrum has a complex shape at very early times, after which it settles into a more generic shape with a smooth curvature below the peak energy, $E_p$, and a clear high-energy power law that cuts off at $\sim 5~$MeV in the observer frame. The spectrum becomes narrower and softer at late times with $E_p$ steadily decreasing during the pulse decay from $E_p \approx 250~$keV to $E_p \approx 100~$keV. The low-energy index, $\alpha$, decreases during the bright part of the pulse from $\alpha \approx -0.5$ to $\alpha \approx -1$, although the low-energy part is better fit with a broken power law when the signal-to-noise ratio is high. The high-energy power law is generated by the reverse shock at low optical depths ($\tau < 30$) and has an index that decreases from $\beta \approx -2$ to $\beta \approx -2.4$.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript models an internal collision in a GRB jet using 1D special relativistic hydrodynamics to simulate forward and reverse shocks, followed by a designated RMS simulation code to evolve photon distributions from the shocks through the photosphere. For a specific case with assumed redshift z=1, it reports a short ~0.1 s pulse with ~23% radiative efficiency, and time-resolved spectra that evolve from complex early shapes to a generic form with smooth curvature below Ep (decreasing from ~250 keV to ~100 keV), low-energy index alpha evolving from ~-0.5 to ~-1, and a high-energy power-law component (beta from ~-2 to ~-2.4) generated by the reverse shock at tau <30, with a cutoff at ~5 MeV.

Significance. If validated, these results offer a detailed prediction for the spectral and temporal evolution of prompt emission arising from radiation-mediated shocks in GRB internal collisions. This could help explain observed GRB pulse characteristics such as time-dependent peak energies and high-energy cutoffs, providing a physical basis for sub-photospheric dissipation models. The approach yields concrete, falsifiable outputs like the 23% efficiency and specific beta indices that can be compared to observations.

major comments (1)
  1. [Numerical methods / RMS simulation description] The central spectral features, including the high-energy power law index beta evolving from -2 to -2.4 and the cutoff at ~5 MeV, depend on the accuracy of the RMS code in tracking photon distributions at low optical depths (tau < 30). However, no convergence tests, resolution studies, or benchmark comparisons for the RMS code are reported, which is necessary to rule out numerical artifacts influencing these results.
minor comments (2)
  1. [Abstract] The abstract would benefit from specifying the key internal collision parameters (e.g., relative Lorentz factors or densities) used in the simulation to allow better assessment of generality.
  2. [Results] Inclusion of error bars or uncertainty estimates on the reported quantities such as the 23% efficiency and Ep values would strengthen the presentation, especially given the numerical nature of the study.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their positive evaluation of the significance of our results and for the constructive comment on the numerical methods. We address the major comment below and will incorporate revisions to strengthen the manuscript.

read point-by-point responses

  1. Referee: [Numerical methods / RMS simulation description] The central spectral features, including the high-energy power law index beta evolving from -2 to -2.4 and the cutoff at ~5 MeV, depend on the accuracy of the RMS code in tracking photon distributions at low optical depths (tau < 30). However, no convergence tests, resolution studies, or benchmark comparisons for the RMS code are reported, which is necessary to rule out numerical artifacts influencing these results.

    Authors: We agree that explicit convergence tests, resolution studies, and benchmark comparisons for the RMS code are important to demonstrate the robustness of the spectral features at low optical depths. Although the RMS code builds on methods validated in our prior work on radiation-mediated shocks, these details were not included in the present manuscript. In the revised version, we will add an appendix or dedicated subsection presenting resolution studies (varying spatial and energy grid resolutions) focused on tau < 30, along with benchmark comparisons to analytic limits and previous simulations of photon spectra in shocks. This will confirm that the reported evolution of beta from approximately -2 to -2.4 and the ~5 MeV cutoff are not numerical artifacts. revision: yes

Circularity Check

0 steps flagged

Simulation outputs emerge directly from hydro and radiative-transfer runs

full rationale

The paper's central results—a ~0.1 s pulse, 23% radiative efficiency, time-evolving spectrum with Ep decreasing from ~250 keV to ~100 keV, alpha from -0.5 to -1, and beta from -2 to -2.4—are obtained by running a 1D special-relativistic hydro simulation of an internal collision followed by a designated RMS code that propagates photon distributions from the shocks through the downstream region to well above the photosphere. These quantities are direct numerical outputs for the chosen initial conditions and z=1; they are not obtained by fitting parameters to a subset of data and then predicting closely related observables, nor do they reduce via self-citation or definitional closure to the inputs. The derivation chain is therefore self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The simulation rests on standard relativistic hydrodynamics and radiative transfer assumptions plus chosen initial conditions for one specific collision; no new particles or forces are introduced.

free parameters (2)
  • redshift z = 1
    Assumed value of 1 used to convert to observer-frame quantities such as pulse duration and energies
  • internal collision parameters
    Specific initial Lorentz factors, densities, and radii for the colliding shells that define the 'specific case considered'
axioms (2)
  • domain assumption 1D special relativistic hydrodynamics sufficiently captures the shock structure and downstream flow for this internal collision
    Invoked to generate the forward and reverse shocks whose photon distributions are then tracked
  • domain assumption The designated RMS simulation code correctly evolves photon distributions through scattering and emission up to and above the photosphere
    Core assumption enabling the spectral and light-curve calculations at different optical depths and angles

pith-pipeline@v0.9.0 · 5913 in / 1874 out tokens · 96574 ms · 2026-05-21T18:19:41.940129+00:00 · methodology

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. From Internal Collision to Photon Escape: First-Principles Modeling of Radiation-Mediated Shocks in Gamma-Ray Burst Photospheres

    astro-ph.HE 2026-04 unverdicted novelty 8.0

    Self-consistent simulations reveal that reverse shocks in GRB photospheres stay radiation-mediated down to optical depths of a few tenths, with photons decoupling over a broad radial range and forming a quasi-thermal ...

  2. Magnetized Shocks Mediated by Radiation from Leptonic and Hadronic Processes

    astro-ph.HE 2025-11 unverdicted novelty 4.0

    Synchrotron self-absorption alters radiation-mediated shock profiles for magnetizations above 10^{-8}, subshocks appear above 0.1, and hadronic processes add a high-energy photon tail with negligible effect on overall...

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages · cited by 2 Pith papers

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....

    work page