Energy-Stratified Shock Model Overview

• Energy-Stratified Shock Model is a framework where shocks partition energy by region, species, or scale, setting the basis for analyzing different radiative and kinetic processes.
• The model underlies diverse applications including GRB external shocks, CME-driven shocks, plasma stratification, and energetic material hotspots, each emphasizing distinct energy channels.
• By distinguishing early vs. late emission, spatial hierarchies, and spectral cutoffs, the model refines traditional shock theory and enhances predictive capabilities across disciplines.

“Energy-Stratified Shock Model” (Editor's term) denotes a family of shock descriptions in which energy is not treated as a single homogeneous downstream quantity, but is partitioned by shock region, species, wave mode, spatial scale, or epoch. In the literature summarized here, the term covers component-stratified gamma-ray-burst external shocks with distinct reverse-shock and forward-shock synchrotron self-Compton components (Fraija et al., 2012), diffusive-shock systems in which higher-energy particles occupy broader precursors while lower-energy particles remain near the shock (Wang et al., 2015), multiscale hotspot localization behind shocks in heterogeneous energetic materials (Nguyen et al., 2022), ion-species and temperature stratification in multi-ion plasma shocks (Keenan et al., 2017), and radially stratified relativistic ejecta whose forward and reverse shocks are evolved by global energy conservation (Wang et al., 21 Jul 2025). The common feature is not a single canonical mechanism, but a recurrent structure: different parts of the shock system dominate different energies, times, or transport channels.

1. Conceptual basis and major variants

In the gamma-ray-burst external-shock formulation, the stratification is explicitly “component-stratified” rather than a radial energy profile inside one zone. A relativistic forward shock propagating into the circumburst medium naturally produces long-lasting GeV emission through synchrotron self-Compton, while a reverse shock crossing the ejecta produces a short-duration prompt-like high-energy component with a timescale tied to $T_{90}$; in the thick-shell regime, the reverse shock becomes relativistic before crossing the shell, and its SSC emission can be strong (Fraija et al., 2012).

In isolated CME-driven shocks, the stratification is spatial and transport-based. Lower-energy particles remain closer to the shock, whereas higher-energy particles sample a thicker precursor because the mean free path scales as $\lambda=v_L\tau$. Yet the resulting spectrum is not a broken power law: the downstream spectrum is a Maxwellian followed by an approximately single power-law tail and then a high-energy cutoff near a few MeV (Wang et al., 2015).

In multiscale shock-to-detonation modeling of heterogeneous energetic materials, stratification is defined by spatially and temporally non-uniform internal energy localization. Shock passage through pressed HMX produces localized hotspots associated with void collapse, crack collapse, interfaces, and hotspot growth; the mesoscale state is summarized by hotspot area $A_{\text{hs}}(t)$, average hotspot temperature $T_{\text{hs}}(t)$, and their rates of change (Nguyen et al., 2022).

In strong collisional shocks of two-ion plasmas, stratification is species-selective. The shock front is enriched with the lighter ion species, the enrichment scales as $M^4$ for $M\gg 1$, and the temperature separation has a nearly universal structure in which the lighter species is hotter in the electron pre-heat layer and the heavier species tends to be slightly hotter in the ion compression and equilibration layers (Keenan et al., 2017).

These uses establish that the phrase denotes a cross-domain modeling pattern rather than a single formal theory. This suggests that “energy stratification” is best understood as a shock-induced hierarchy of dominant channels—radiative, kinetic, thermal, magnetic, or chemical—resolved by region, species, or scale.

2. Relativistic emission models and shock-component stratification

A particularly clear realization appears in leptonic GRB external-shock models. Electrons accelerated in both forward and reverse shocks are assumed to follow a power law,

$$N(\gamma_e)\,d\gamma_e \propto \gamma_e^{-p}\,d\gamma_e,\qquad \gamma_e>\gamma_m,$$

with minimum Lorentz factors such as

$$\gamma_{m,f}=\epsilon_{e,f}\frac{p-2}{p-1}\frac{m_p}{m_e}\Gamma_f, \qquad \gamma_{m,r}=\epsilon_{e,r}\frac{p-2}{p-1}\frac{m_p}{m_e}\Gamma_r.$$

In this framework, reverse-shock SSC gives a short-lived prompt-like component, while forward-shock SSC gives a longer GeV–TeV component. For GRB 940217, forward-shock SSC with $E^{IC}_{m,f}\approx 33.1$ GeV and duration $\sim 1000$ s is consistent with the observed extended GeV emission and the late 18 GeV photon; for GRB 941017, the high-energy tail is modeled as a superposition of reverse-shock SSC peaking at $\lambda=v_L\tau$0 MeV and forward-shock SSC peaking at $\lambda=v_L\tau$1 GeV; for GRB 970417A, forward-shock SSC reaches $\lambda=v_L\tau$2 GeV, while reverse-shock SSC contributes at $\lambda=v_L\tau$3 MeV (Fraija et al., 2012).

A related but distinct GRB realization occurs in sub-photospheric radiation-mediated shocks. A train of shocks with different upstream four-velocities $\lambda=v_L\tau$4 and formation depths $\lambda=v_L\tau$5 produces a superposition of thermal components whose time-integrated spectral energy distribution below the peak behaves as

$\lambda=v_L\tau$6

provided the shock strengths span a sufficiently broad distribution. In this model, a single stronger shock followed by a sequence of weaker shocks produces a double-peak SED, and constant-velocity versus decelerating shocks change both the low-energy slope and the amount of spectral substructure (Keren et al., 2014).

The most explicit contemporary formulation of radial stratification in relativistic ejecta treats the jet as axisymmetric and two-dimensional, with arbitrary angular profiles $\lambda=v_L\tau$7, $\lambda=v_L\tau$8, and a radial distribution $\lambda=v_L\tau$9, where $A_{\text{hs}}(t)$0. The forward–reverse shock system is evolved with an energy conservation prescription rather than pressure balance, and this changes the predicted reverse-shock emission substantially. In particular, reverse-shock emission in the thin-shell case is significantly overestimated in analytic pressure-balance models; the code also shows that off-axis observers can see a thin-to-thick transition in structured jets, although the resulting light-curve morphology is hard to distinguish from pure thin- or thick-shell cases. Because the same radial structure also provides hydrodynamic energy injection, the framework naturally extends to refreshed shocks and kilonova afterglows (Wang et al., 21 Jul 2025).

3. Particle acceleration, spectral breaks, and cutoff-limited stratification

In heliospheric and cosmic-ray applications, energy stratification is often controlled by diffusion, escape, and age. In the isolated CME-driven Monte Carlo model, pitch-angle scattering is prescribed by a constant scattering time $A_{\text{hs}}(t)$1 with

$A_{\text{hs}}(t)$2

and the downstream nonthermal tail has an average spectral index $A_{\text{hs}}(t)$3, close to the observed low-energy index $A_{\text{hs}}(t)$4 for the 14 December 2006 event. Varying $A_{\text{hs}}(t)$5 changes the extent of the tail but not its basic shape. The largest simulated maximum energy is $A_{\text{hs}}(t)$6 MeV at $A_{\text{hs}}(t)$7, and no case exceeds the upper end of the observed $A_{\text{hs}}(t)$8–$A_{\text{hs}}(t)$9 MeV break range. The model therefore yields a cutoff near a few MeV rather than a true broken power law, supporting the conclusion that an isolated shock does not by itself explain the observed spectral breaks (Wang et al., 2015).

A complementary dynamical Monte Carlo study of non-relativistic planar shocks shows that the scattering law itself can stratify the shock energetics. With Gaussian angular distributions of width $T_{\text{hs}}(t)$0, smaller $T_{\text{hs}}(t)$1 gives more anisotropic scattering, higher escape losses, softer subshock spectra, and lower cutoff energies, whereas isotropic scattering minimizes energy losses and hardens the subshock spectrum. The reported cutoff energies increase from $T_{\text{hs}}(t)$2 MeV in Case A to $T_{\text{hs}}(t)$3 MeV in Case D, while the escaped-particle counts decrease from $T_{\text{hs}}(t)$4 to $T_{\text{hs}}(t)$5; the main conclusion is that smaller energy losses in the prescribed scattering law produce harder spectra (Wang et al., 2011).

For plane shocks accelerating relativistic electrons, energy stratification becomes explicitly spatial. With Bohm-type diffusion,

$T_{\text{hs}}(t)$6

and synchrotron/IC losses $T_{\text{hs}}(t)$7, the electron spectrum at the shock reaches a steady cutoff at the equilibrium momentum $T_{\text{hs}}(t)$8, but downstream the cutoff momentum decreases with distance from the shock because of radiative losses. The thickness of the spatial distribution scales as $T_{\text{hs}}(t)$9, the downstream integrated electron spectrum steepens by one power for $M^4$0, and the break momentum decreases as $M^4$1. In CR-modified shocks, both proton and electron spectra become concave, and near the cutoff the upstream integrated electron spectrum can dominate the downstream one (Kang, 2011).

At larger astrophysical scales, the maximum proton energy of a non-relativistic shock can be estimated empirically as

$M^4$2

with $M^4$3 when diffusing CRs drive magnetic-field amplification and $M^4$4 when escaping CRs drive it. The multi-zone model shows that SNRs reach PeV energies only under a restricted set of conditions, notably if $M^4$5 and escaping particles drive magnetic amplification; it also emphasizes that older remnants may still display signatures of PeV particles accelerated at earlier epochs (Diesing, 2023).

These results separate two recurring spectral outcomes. One is a true stratified transport hierarchy, in which different energies populate different spatial zones and epochs. The other is a cutoff-limited hierarchy, in which stratification ends at a terminal $M^4$6 without generating a broken power law. The isolated-CME case is an explicit example of the second outcome.

4. Plasma, multi-fluid, and partially ionised shock stratification

In multi-ion collisional plasmas, strong shocks generate both composition and temperature stratification. The local light-species mass fraction

$M^4$7

departs from its upstream value $M^4$8, with the lighter species enriched at the shock front and into the electron pre-heat layer. The integrated enrichment scales as $M^4$9 for strong shocks, and the relative enrichment is largest when the light species is initially rare. The ion temperature separation satisfies a balance involving $M\gg 1$0 work, viscous heating, ion heat flux, and inter-species energy exchange; in the pre-heat layer the heat-flux term dominates, giving

$M\gg 1$1

so the lighter species is hotter there, while the heavier species tends to be slightly hotter in the compression region (Keenan et al., 2017).

A single-fluid discrete Boltzmann plasma-shock model resolves a different form of stratification, one mediated by charge separation and thermodynamic non-equilibrium. Electrons are taken as inertialess and Boltzmann-distributed,

$M\gg 1$2

while ions are kinetic and obey a BGK-based discrete Boltzmann equation coupled to Poisson’s equation. The resulting plasma shock differs from a neutral-fluid shock by developing a peak in density, temperature, pressure, and flow speed at the front, a negative charge layer at the wave front, a positive charge layer at the wave rear, and strong thermodynamic non-equilibrium that intensifies with Mach number. The paper further states that the characteristics of thermodynamic non-equilibrium can distinguish plasma shock waves from detonation waves (Liu et al., 2021).

In partially ionised hydrogen shocks relevant to the chromosphere, energy stratification appears through ion–neutral decoupling, collisional ionisation, and radiative ionisation/recombination. A multi-level hydrogen two-fluid model shows that collisional ionisation removes macroscopic thermal energy, so both the post-shock temperature and the temperature inside the finite-width shock are significantly cooler than in single-fluid MHD, while the compression is greater, especially in the mid to lower chromosphere. The study concludes that such shocks are not accurately described by Rankine–Hugoniot relations, so it may be incorrect to infer lower-atmospheric shock properties from single-fluid jump conditions (Snow et al., 2023).

A complementary two-fluid study of an isothermal, stratified atmosphere shows that an upward acoustic shock encounters a height where sound and Alfvén speeds are equal and then separates into fast and slow components, with the energy distribution depending on the magnetic-field angle. In the finitely coupled regime, the fast-mode shock produces the largest ion–neutral drift and the broadest structure; the paper finds that fast-mode shock widths can exceed the pressure scale height, suggesting a direct observable of two-fluid effects in the lower solar atmosphere (Snow et al., 2020).

Taken together, these plasma and two-fluid models show that shock energy can be stratified not only by position or time, but also by species, wave mode, and collisional state. This suggests that single-temperature or single-fluid closures suppress physically important internal layers of the shock.

5. Multiscale materials, thermoacoustic cascades, and implosive shocks

In heterogeneous energetic materials, the shock-to-detonation transition is governed by mesoscale localization rather than by a single homogenized shock state. A hotspot is defined by the threshold

$M\gg 1$3

and the mesoscale state is reduced to hotspot area $M\gg 1$4 and temperature $M\gg 1$5. The macroscale reaction progress variable is linked to hotspot growth through

$M\gg 1$6

A physics-aware recurrent convolutional neural network learns the evolution operator

$M\gg 1$7

from DNS of pressed HMX microstructures, and the resulting multiscale model reproduces hotspot metrics and the James criticality envelope while reducing mesoscale evaluation time from hours-to-days on HPC clusters to less than a second on a commodity laptop (Nguyen et al., 2022).

In thermoacoustic shock waves, energy stratification is spectral. Fully compressible Navier–Stokes simulations of a 2.58 m looped resonator show three regimes: monochromatic harmonic growth, hierarchical spectral broadening, and a shock-wave-dominated limit cycle. The acoustic energy density is

$M\gg 1$8

with balance

$M\gg 1$9

At the limit cycle, the spectral energy density obeys

$$N(\gamma_e)\,d\gamma_e \propto \gamma_e^{-p}\,d\gamma_e,\qquad \gamma_e>\gamma_m,$$0

and energy production is balanced by dissipation at the captured shock-thickness scale (Gupta et al., 2017).

In cylindrical and spherical shock implosion in water, the stratification is expressed in the balance of internal and kinetic energy densities and their fluxes behind the front. Using a polytropic EOS for water and strong-shock Rankine–Hugoniot relations, the analysis shows that throughout the entire range of implosion radii considered, the internal and kinetic energy density fluxes are equal. The same framework derives threshold radii at which the shock velocity begins a rapid quasi self-similar increase and yields power-law forms $$N(\gamma_e)\,d\gamma_e \propto \gamma_e^{-p}\,d\gamma_e,\qquad \gamma_e>\gamma_m,$$1 whose exponents are close to known self-similar converging-shock scalings (Chefranov et al., 2021).

These examples broaden the subject beyond astrophysical radiation and particle acceleration. In one case, the relevant energy layers are hotspot populations inside a shocked microstructure; in another, they are harmonic bands in a spectral cascade; in another, they are internal and kinetic flux channels during implosion. The unifying element is again the existence of a constrained hierarchy of shock-mediated energy reservoirs.

6. Recurring themes, misconceptions, and limitations

Several misconceptions are directly contradicted by the cited literature. First, energy stratification does not necessarily imply a broken power law. In isolated CME-driven shocks, the model produces a single power-law tail with a high-energy cutoff near a few MeV and specifically argues that, without interaction with another shock, an isolated shock does not form the observed $$N(\gamma_e)\,d\gamma_e \propto \gamma_e^{-p}\,d\gamma_e,\qquad \gamma_e>\gamma_m,$$2–$$N(\gamma_e)\,d\gamma_e \propto \gamma_e^{-p}\,d\gamma_e,\qquad \gamma_e>\gamma_m,$$3 MeV spectral break (Wang et al., 2015).

Second, energy stratification is not always a radial energy profile internal to one emission zone. In the GRB external-shock model, the authors explicitly distinguish their interpretation from that usage: the stratification is component-based, with reverse-shock SSC dominating early and forward-shock SSC dominating later (Fraija et al., 2012).

Third, commonly used analytic closures can overestimate reverse-shock emission. The energy-conserving structured-jet calculation finds that pressure-balance analytic models significantly overestimate thin-shell reverse-shock emission, in part because they ignore the velocity gradient across a spreading shell and in part because the blast decelerates during reverse-shock crossing (Wang et al., 21 Jul 2025).

Fourth, single-fluid jump conditions can fail in partially ionised shocks. In the chromospheric hydrogen model, collisional ionisation and excitation act as macroscopic thermal-energy sinks, so the post-shock state can be cooler and more compressed than its MHD analogue; the paper therefore warns that Rankine–Hugoniot relations may be incorrect for inferring lower-atmospheric shock properties (Snow et al., 2023).

The limitations are equally recurrent. Several models are effectively one-dimensional or single-zone: the GRB external-shock treatment assumes constant microphysics and does not include radial stratification within the ejecta or multi-zone radiative transfer (Fraija et al., 2012); the isolated CME shock calculations prescribe a single diffusion regime and a finite free-escape boundary (Wang et al., 2015); the energetic-material surrogate is two-dimensional and reduces chemistry to a hotspot threshold and averaged closure (Nguyen et al., 2022). This suggests that “energy-stratified shock model” should be read as a structural descriptor of how shocks distribute energy, not as a guarantee of full microphysical completeness.

Across these domains, the most stable common result is that shocks often organize energy into non-equivalent layers—early and late, forward and reverse, light and heavy species, hotspot and bulk, inertial range and dissipation range, or plasma and neutral components. The detailed realization changes with the governing physics, but the hierarchy itself is the defining feature.