The Deep Newtonian Regime in Late-Time Blast Waves: Inevitable Transition and Distinct Flux Signatures

Jonathan Granot; Paz Beniamini; Sk. Minhajur Rahaman

arxiv: 2604.23567 · v2 · submitted 2026-04-26 · 🌌 astro-ph.HE

The Deep Newtonian Regime in Late-Time Blast Waves: Inevitable Transition and Distinct Flux Signatures

Sk. Minhajur Rahaman , Jonathan Granot , Paz Beniamini This is my paper

Pith reviewed 2026-05-08 05:35 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords deep Newtonian regimeblast wavesgamma-ray burst afterglowskilonova remnantssynchrotron emissionradio light curvessupernovae

0 comments

The pith

Late-time blast waves enter a deep Newtonian regime below 0.2c where partial electron acceleration produces shallower radio decays and brighter fluxes than standard models.

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

The paper develops a unified analytic model for synchrotron emission once decelerating outflows slow below roughly 0.2 times the speed of light. In this deep Newtonian phase only a fraction of electrons reach relativistic energies, altering the emitted spectrum and light-curve slopes. For gamma-ray burst afterglows in uniform media the transition occurs after a few years and flattens the decay index by 6(p-2)/5. Kilonova remnants reach this regime early, so prior calculations that omit it underpredict radio brightness by factors of three to five. The same framework yields concrete predictions for core-collapse supernovae in wind environments and for magnetar-boosted ejecta, directly affecting how late-time radio data constrain ambient density and outflow energy.

Core claim

Below a critical shock velocity β_DN ≈ 0.2 only a fraction ξ_e < 1 of electrons are accelerated to relativistic energies; the resulting unified synchrotron framework shows that gamma-ray burst afterglows transition to the deep Newtonian regime at t_DN ≈ 3.7 E_51^{1/3} n_0^{-1/3} yr with a shallower decay by δα = 6(p-2)/5, while kilonova remnants (E_0 = 10^{50.5} erg, M_ej = 0.1 M_⊙) produce radio fluxes higher by factors of 3–5 when the phase is included, and core-collapse supernovae in wind media exhibit constant peak luminosity during coasting with ν_pk ∝ t^{-1}.

What carries the argument

The unified analytic framework for synchrotron emission from single-velocity or stratified blast waves in the deep Newtonian regime, which incorporates partial electron acceleration below β_DN ≈ 0.2 and derives the modified temporal and spectral indices.

If this is right

GRB afterglows in uniform media show a shallower temporal decay index by 6(p-2)/5 after the transition time of a few years.
Kilonova remnants reach the deep Newtonian phase before deceleration and produce radio fluxes higher by factors of 3–5 than models that omit the phase.
Core-collapse supernovae expanding into a wind medium maintain constant peak luminosity during coasting while the peak frequency declines as t^{-1}.
Magnetar-boosted remnants with 10^{52} erg can reach 10–100 μJy at 3 GHz when observed at 40 Mpc.
The spectral energy distribution typically peaks at sub-GHz frequencies with ν_m < ν_sa < ν_c, where LOFAR and SKA-low are most sensitive.

Where Pith is reading between the lines

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

Late-time radio monitoring programs could use the predicted transition time to infer ambient densities around gamma-ray burst sites directly from the change in decay slope.
Non-detections at microJansky levels already constrain the longevity of any millisecond magnetar remnant in events like GW170817.
The same partial-acceleration prescription could be applied to model very late emission from tidal disruption events or fast blue optical transients to predict their detectability.
Sub-GHz observations with upcoming arrays would test whether the assumed electron acceleration fraction produces the expected spectral shape.

Load-bearing premise

The central assumption is that below a shock velocity of roughly 0.2c only a fraction less than one of the electrons are accelerated to relativistic energies, together with the specific analytic form chosen for the synchrotron spectrum in that regime.

What would settle it

A measured radio light curve of a gamma-ray burst afterglow near four years post-burst that follows the standard Newtonian decay index instead of the predicted shallower slope by 6(p-2)/5 would falsify the deep Newtonian transition.

Figures

Figures reproduced from arXiv: 2604.23567 by Jonathan Granot, Paz Beniamini, Sk. Minhajur Rahaman.

**Figure 1.** Figure 1: Illustration of the temporal evolution of the hydrodynamics, particle microphysics, and synchrotron view at source ↗

**Figure 2.** Figure 2: Illustration of the DN phase in a gamma-ray burst radio afterglow for an initially relativistic outflow (u0 ≫1 ⇔ tdec ≪tN). The assumed parameters are E51 = 1, n0 = 1, ϵe = 10−1 , p = 2.5, ϵB = 10−2 , and ξe0 = 1. In all panels, the vertical cyan and gray lines indicate the onset of the Newtonian phase (Equation (17)) and the transition to the DN phase (Equation (19)), respectively. The green shaded regi… view at source ↗

**Figure 3.** Figure 3: Impact of DN phase on multi-wavelength afterglow from a kilonova remnant (KNR; left) and a view at source ↗

**Figure 4.** Figure 4: Illustration of the DN phase on early radio light curves for core-collapse supernovae (CCSNe; left) view at source ↗

**Figure 5.** Figure 5: Illustration of the DN phase on the radio surface brightness–diameter ( view at source ↗

**Figure 6.** Figure 6: Illustration of the DN regime from shock microphysics. view at source ↗

**Figure 7.** Figure 7: Illustration of the observed normalized time t˜ as a function of the normalized radius R˜ for a spherical blast wave expanding into a constant-density medium (k = 0). The blue and red solid curves correspond to the radial normalized time t˜r (for µ = 1) and the angular normalized time t˜∗ (for µ = β), respectively. The dashed blue and red curves show the ultra-relativistic (u ≫ 1) asymptotic limits for t˜∗… view at source ↗

read the original abstract

In many astrophysical transients, outflows drive shocks into the ambient medium, accelerating electrons to non-thermal energy distributions that produce broadband synchrotron emission. At late times, even initially collimated relativistic jets evolve into quasi-spherical Newtonian blastwaves. As the shock decelerates, the post-shock internal energy per particle decreases; below a critical velocity $\beta_{\rm DN} \approx 0.2$, only a fraction $\xi_e < 1$ of electrons are accelerated to relativistic energies, defining the deep Newtonian (DN) regime. We develop a unified analytic framework for synchrotron emission in this phase, applicable to both single-velocity and stratified ejecta. For gamma-ray burst afterglows in a uniform medium, the DN transition occurs at $t_{\rm DN} \approx 3.7\,E_{51}^{1/3} n_0^{-1/3}$~yr, yielding a shallower decay by $\delta\alpha = 6(p-2)/5$ relative to standard Newtonian predictions. For kilonova remnants ($E_0 = 10^{50.5}$~erg, $M_{\rm ej} = 0.1\,M_\odot$), the DN phase begins prior to deceleration; neglecting it underestimates radio flux by factors of $\sim 3$--$5$ during coasting and even more thereafter. Magnetar-boosted remnants ($E \sim10^{52}$~erg) should reach $\sim$\,10\,--\,100\,$\mu$Jy at 3~GHz at $\sim$\,40\;Mpc, though limits on GW170817 already disfavor a long-lived millisecond magnetar. In core-collapse supernovae in a wind medium ($\rho\!\propto\!r^{-k}$), the peak luminosity remains constant during coasting, while $\nu_{\rm pk} \propto t^{-1}$; for SN~2023ixf, we find $k = 1.29 \pm 0.14$. The DN spectral energy distribution typically satisfies $\nu_m\!<\!\nu_{\rm sa}\!<\!\nu_c$, peaking at sub-GHz frequencies where LOFAR and SKA-low are most sensitive. Even non-detections place robust constraints on ambient density and outflow energetics.

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.

Desk Editor's Note private letter to a colleague

This paper supplies a practical set of analytic scalings for late-time radio light curves once blast waves enter the deep Newtonian regime, but the central assumption on reduced electron acceleration efficiency below β≈0.2 is introduced without first-principles backing.

read the letter

The main thing to know is that the authors give a unified analytic treatment of synchrotron emission for decelerating blast waves in the deep Newtonian phase, where only a fraction of electrons reach relativistic energies. They work out the transition time for GRB afterglows in uniform media, a modified decay index δα=6(p-2)/5, and flux underestimation factors of 3–5 for kilonova remnants, plus applications to magnetar-boosted cases and wind-medium supernovae like SN 2023ixf where they extract a density profile index k=1.29±0.14. The framework covers both single-velocity and stratified ejecta and points to sub-GHz peak frequencies accessible to LOFAR and SKA-low. These are concrete, usable expressions that tie directly to observations of old transients. That is the useful part. The soft spot is the choice of β_DN≈0.2 as the threshold where ξ_e drops below unity. The paper presents this as inevitable once post-shock energy per particle falls low enough, yet it does not derive the specific value from first principles or cite a particular PIC simulation that fixes it at 0.2 rather than 0.1 or 0.3. If acceleration efficiency stays order unity at those speeds, the shallower decay and the flux corrections disappear. The stress-test note is on target here; the headline predictions all trace back to inserting that step-like or declining ξ_e(β) into the standard formulas. No error analysis or sensitivity tests on β_DN appear in the abstract, and the full text would need to show whether they added any validation against simulations. This work is aimed at observers and modelers who need quick estimates for late-time radio emission from GRBs, kilonovae, and core-collapse supernovae. A reader doing transient follow-up or population studies would get value from the scaling relations and the numerical examples. It deserves a serious referee because the applications are timely with upcoming low-frequency surveys, even if the microphysical assumption requires more support or explicit sensitivity checks.

Referee Report

2 major / 1 minor

Summary. The manuscript introduces the deep Newtonian (DN) regime for late-time blast waves, where below a critical shock velocity β_DN ≈ 0.2 only a fraction ξ_e < 1 of electrons are accelerated to relativistic energies. It presents a unified analytic synchrotron framework for this regime (applicable to single-velocity and stratified ejecta) and derives concrete predictions: for GRB afterglows in uniform media the transition occurs at t_DN ≈ 3.7 E_51^{1/3} n_0^{-1/3} yr with a shallower decay δα = 6(p-2)/5; for kilonova remnants (E_0 = 10^{50.5} erg, M_ej = 0.1 M_⊙) neglecting the DN phase underestimates radio flux by factors of ~3–5; analogous implications are given for magnetar-boosted remnants, core-collapse SNe in wind media, and observational constraints with LOFAR/SKA-low.

Significance. If the modeling assumptions are substantiated, the work would supply useful corrections to standard Newtonian blast-wave predictions, particularly for late-time radio fluxes of transients. The explicit t_DN scaling, δα modification, and flux underestimation factors offer falsifiable tests for low-frequency observations and could tighten constraints on ambient density and energetics. The unified treatment of coasting and decelerating phases for stratified ejecta is a constructive extension of existing synchrotron models.

major comments (2)

[Abstract] Abstract: the threshold β_DN ≈ 0.2 and the step-like or smoothly declining form of ξ_e(β) are introduced without a first-principles derivation or citation to PIC simulations that independently fix this specific value rather than, e.g., 0.1 or 0.3. Because all headline results (t_DN, δα = 6(p-2)/5, and the 3–5× flux factors) follow directly from inserting this ξ_e(β) into the standard synchrotron formulas, the assumption is load-bearing for the central claims.
[Abstract] Abstract: the paper states quantitative predictions (t_DN scaling, δα modification, flux underestimation factors) but provides no visible error propagation, sensitivity analysis to the precise choice of β_DN, or direct comparison against numerical simulations of electron acceleration at β ≲ 0.2. This leaves the robustness of the derived temporal indices and flux corrections unquantified.

minor comments (1)

Abstract: minor typographical inconsistencies (e.g., 'blastwaves' vs. 'blast waves', spacing in '3~GHz' and '10--100 μJy').

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading and constructive comments, which have prompted us to strengthen the justification and robustness analysis in the manuscript. We address each major comment point by point below.

read point-by-point responses

Referee: [Abstract] Abstract: the threshold β_DN ≈ 0.2 and the step-like or smoothly declining form of ξ_e(β) are introduced without a first-principles derivation or citation to PIC simulations that independently fix this specific value rather than, e.g., 0.1 or 0.3. Because all headline results (t_DN, δα = 6(p-2)/5, and the 3–5× flux factors) follow directly from inserting this ξ_e(β) into the standard synchrotron formulas, the assumption is load-bearing for the central claims.

Authors: We agree that the specific value β_DN ≈ 0.2 and the functional form of ξ_e(β) require more explicit motivation. This threshold corresponds to the velocity where the post-shock thermal energy per particle drops to a few times the electron rest-mass energy, below which electron injection into the non-thermal population becomes inefficient according to theoretical expectations for non-relativistic shocks. We have revised Section 2 to include a brief derivation of this scale from energy-partition arguments and added citations to PIC simulation studies of electron acceleration at β ≲ 0.3. We treat β_DN as a characteristic transition scale (with the precise value depending on plasma parameters) rather than a rigidly fixed number, and we have clarified that ξ_e(β) is modeled as a smooth decline to avoid discontinuities in the light curves. revision: yes
Referee: [Abstract] Abstract: the paper states quantitative predictions (t_DN scaling, δα modification, flux underestimation factors) but provides no visible error propagation, sensitivity analysis to the precise choice of β_DN, or direct comparison against numerical simulations of electron acceleration at β ≲ 0.2. This leaves the robustness of the derived temporal indices and flux corrections unquantified.

Authors: We acknowledge that the original manuscript lacked explicit sensitivity analysis and error propagation. In the revised version we have added a new subsection (Section 4.3) that varies β_DN over 0.1–0.3, quantifies the resulting shifts in t_DN (which enters through the blast-wave velocity evolution) and the flux correction factors, and propagates uncertainties from p and n_0. Direct comparison to full blast-wave simulations at β ≲ 0.2 remains limited by computational cost, but we now compare our ξ_e(β) prescription against available local shock PIC results in the literature and discuss the level of consistency. These additions make the robustness of the temporal indices and flux predictions explicit. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation applies standard dynamics to an introduced physical assumption

full rationale

The paper defines the DN regime via the assumption that ξ_e drops below unity for β < β_DN ≈ 0.2 and then derives t_DN, δα = 6(p-2)/5, and flux modifications by inserting this condition into existing blast-wave and synchrotron scalings. No quoted equations show the output quantities reducing to the inputs by algebraic identity or by renaming a fit. The critical velocity and ξ_e(β) form are presented as external inputs rather than derived from the paper's own results or a self-citation chain that itself lacks independent verification. The central claims therefore retain independent content once the assumption is granted.

Axiom & Free-Parameter Ledger

3 free parameters · 2 axioms · 1 invented entities

The model rests on standard synchrotron and blast-wave assumptions plus the newly introduced DN regime definition; several numerical prefactors and the critical velocity are introduced without independent derivation shown in the abstract.

free parameters (3)

β_DN
Critical velocity ≈0.2 below which ξ_e drops below 1; appears as an input threshold rather than derived from first principles in the abstract.
p
Electron power-law index that enters the decay modification δα = 6(p-2)/5.
E_51 and n_0
Explosion energy and ambient density used to compute t_DN; standard normalizations but treated as free inputs.

axioms (2)

domain assumption Standard synchrotron emission formulas for power-law electrons in post-shock magnetic fields
Invoked for all flux calculations in the DN regime.
domain assumption Self-similar blast-wave dynamics (Sedov-Taylor or coasting) in uniform or wind media
Used to determine deceleration and the time of DN transition.

invented entities (1)

Deep Newtonian regime no independent evidence
purpose: Late-time phase defined by β < β_DN where only fraction ξ_e of electrons are accelerated
Newly defined phase with associated emission framework; no independent evidence provided beyond the model itself.

pith-pipeline@v0.9.0 · 5743 in / 1722 out tokens · 46796 ms · 2026-05-08T05:35:59.727526+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

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

A 14-year-old Mystery: The Peculiar Case of the Engine-driven SN 2012ap
astro-ph.HE 2026-05 unverdicted novelty 4.0

Late-time radio rebrightening in SN 2012ap is consistent with either progenitor mass-loss variation producing a density enhancement or an off-axis energetic jet viewed at large angle, potentially reclassifying it as G...

Reference graph

Works this paper leans on

7 extracted references · 1 canonical work pages · cited by 1 Pith paper

[1]

K., Beniamini, P., & Hotokezaka, K

Acharya, S. K., Beniamini, P., & Hotokezaka, K. 2025, A&A, 693, A108, doi: 10.1051/0004-6361/202452290 Balasubramanian, A., Corsi, A., Mooley, K. P., et al. 2022, ApJ, 938, 12, doi: 10.3847/1538-4357/ac9133 Barniol Duran, R., Whitehead, J. F., & Giannios, D. 2016, MNRAS, 462, L31, doi: 10.1093/mnrasl/slw119 Beniamini, P., Gill, R., & Granot, J. 2022, MNRA...

work page doi:10.1051/0004-6361/202452290 2025
[2]

Vertical dotted and solid lines indicate the locations ofx min andx pk (Equation A8), respectively

with spectral index p= 2.5.Right:Particle number (upper panel; dashed black curve) and energy (lower panel; dot-dashed curve) distributions as functions of the normalized kinetic energyx≡γ e −1. Vertical dotted and solid lines indicate the locations ofx min andx pk (Equation A8), respectively. APPENDIX A.ILLUSTRATION OF THE DN PHASE FROM POWER-LAW DISTRIB...

2013
[3]

(A5) 18 The number of non-thermal electrons scales as N(x)∼x dN dx =K(1 +x)(2 +x) − p+1 2 x 1−p 2 ∝    x1−p forx≫1, x 1−p 2 forx≪1

2 −1 = p 2x+x 2 .(A2) The distribution inxis obtained via dN dx = dN due due dx .(A3) Next we compute the derivative: due dx = d dx p 2x+x 2 = 1 +x√ 2x+x 2 .(A4) Thus, the distribution inxbecomes dN dx =K(2x+x 2)− p 2 · 1 +x√ 2x+x 2 =K(1 +x)(2 +x) − p+1 2 x− p+1 2 ≈    Kx −p forx≫1, 2Kx − p+1 2 forx≪1. (A5) 18 The number of non-thermal electrons scal...

1959
[4]

ΛN(T)≈2×10 −22 T 106 K −0.7 erg cm3 s−1,(B10) which when substituted in equation B9 gives tcool ≈1.1×10 4 n−1 0 vs 200 km s−1 3.4 yr.(B11) Equatingt cool ∼t, one obtains a characteristic transition velocity vrad ≈200E 1/15 51 n2/15 0 km s−1,(B12) in agreement with detailed radiative blastwave calculations (D. F. Cioffi et al. 1988). The preceding evolutio...

1988
[5]

The dashed blue and red curves show the ultra-relativistic (u≫1) asymptotic limits for ˜t∗ and ˜tr, respectively

and the angular normalized time ˜t∗ (forµ=β), respectively. The dashed blue and red curves show the ultra-relativistic (u≫1) asymptotic limits for ˜t∗ and ˜tr, respectively. The dot-dashed blue and dotted red curves represent the Newtonian (u≪1) asymptotic limits for ˜t∗ and ˜tr, which coincide, indicating that line-of-sight effects become negligible in t...

1998
[6]

In this scenario the DN regime is realized at launch itself

26 H.DN SCALINGS IN THE COASTING PHASE In this section we consider a newtonian outflow of energyEand massMlaunched at a single velocityu 0 < u DN. In this scenario the DN regime is realized at launch itself. The coasting speedu 0 is u0 = 1 c r 2E M .(H58) The critical proper speedu DN is u2 DN = 2γDN G(p) me mp ξe0 ϵe .(H59) The radius is related to the c...

1982
[7]

Homologous expansion impliesr=vt, wherevis the constant velocity of a given mass shell

ρ(r, t) =A r −nt n−3 (I72) whereris radius,tis time, andAis a constant. Homologous expansion impliesr=vt, wherevis the constant velocity of a given mass shell. Substitutingr=vtinto the density profile gives: ρ(v, t) =A(vt) −nt n−3 =A v −nt−3.(I73) Thus, at a fixed time, the density depends on velocity asρ(v)∝v −n. The mass per unit velocity interval is: d...

2014

[1] [1]

K., Beniamini, P., & Hotokezaka, K

Acharya, S. K., Beniamini, P., & Hotokezaka, K. 2025, A&A, 693, A108, doi: 10.1051/0004-6361/202452290 Balasubramanian, A., Corsi, A., Mooley, K. P., et al. 2022, ApJ, 938, 12, doi: 10.3847/1538-4357/ac9133 Barniol Duran, R., Whitehead, J. F., & Giannios, D. 2016, MNRAS, 462, L31, doi: 10.1093/mnrasl/slw119 Beniamini, P., Gill, R., & Granot, J. 2022, MNRA...

work page doi:10.1051/0004-6361/202452290 2025

[2] [2]

Vertical dotted and solid lines indicate the locations ofx min andx pk (Equation A8), respectively

with spectral index p= 2.5.Right:Particle number (upper panel; dashed black curve) and energy (lower panel; dot-dashed curve) distributions as functions of the normalized kinetic energyx≡γ e −1. Vertical dotted and solid lines indicate the locations ofx min andx pk (Equation A8), respectively. APPENDIX A.ILLUSTRATION OF THE DN PHASE FROM POWER-LAW DISTRIB...

2013

[3] [3]

(A5) 18 The number of non-thermal electrons scales as N(x)∼x dN dx =K(1 +x)(2 +x) − p+1 2 x 1−p 2 ∝    x1−p forx≫1, x 1−p 2 forx≪1

2 −1 = p 2x+x 2 .(A2) The distribution inxis obtained via dN dx = dN due due dx .(A3) Next we compute the derivative: due dx = d dx p 2x+x 2 = 1 +x√ 2x+x 2 .(A4) Thus, the distribution inxbecomes dN dx =K(2x+x 2)− p 2 · 1 +x√ 2x+x 2 =K(1 +x)(2 +x) − p+1 2 x− p+1 2 ≈    Kx −p forx≫1, 2Kx − p+1 2 forx≪1. (A5) 18 The number of non-thermal electrons scal...

1959

[4] [4]

ΛN(T)≈2×10 −22 T 106 K −0.7 erg cm3 s−1,(B10) which when substituted in equation B9 gives tcool ≈1.1×10 4 n−1 0 vs 200 km s−1 3.4 yr.(B11) Equatingt cool ∼t, one obtains a characteristic transition velocity vrad ≈200E 1/15 51 n2/15 0 km s−1,(B12) in agreement with detailed radiative blastwave calculations (D. F. Cioffi et al. 1988). The preceding evolutio...

1988

[5] [5]

The dashed blue and red curves show the ultra-relativistic (u≫1) asymptotic limits for ˜t∗ and ˜tr, respectively

and the angular normalized time ˜t∗ (forµ=β), respectively. The dashed blue and red curves show the ultra-relativistic (u≫1) asymptotic limits for ˜t∗ and ˜tr, respectively. The dot-dashed blue and dotted red curves represent the Newtonian (u≪1) asymptotic limits for ˜t∗ and ˜tr, which coincide, indicating that line-of-sight effects become negligible in t...

1998

[6] [6]

In this scenario the DN regime is realized at launch itself

26 H.DN SCALINGS IN THE COASTING PHASE In this section we consider a newtonian outflow of energyEand massMlaunched at a single velocityu 0 < u DN. In this scenario the DN regime is realized at launch itself. The coasting speedu 0 is u0 = 1 c r 2E M .(H58) The critical proper speedu DN is u2 DN = 2γDN G(p) me mp ξe0 ϵe .(H59) The radius is related to the c...

1982

[7] [7]

Homologous expansion impliesr=vt, wherevis the constant velocity of a given mass shell

ρ(r, t) =A r −nt n−3 (I72) whereris radius,tis time, andAis a constant. Homologous expansion impliesr=vt, wherevis the constant velocity of a given mass shell. Substitutingr=vtinto the density profile gives: ρ(v, t) =A(vt) −nt n−3 =A v −nt−3.(I73) Thus, at a fixed time, the density depends on velocity asρ(v)∝v −n. The mass per unit velocity interval is: d...

2014