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

Imagine a Gamma-Ray Burst (GRB) as the universe's most spectacular fireworks display. For decades, astronomers have been trying to figure out exactly how these fireworks work. We know the initial "boom" (the prompt gamma-ray emission) comes from a massive explosion. But then, something strange happens: instead of fading away quickly like a normal firework, the light stays bright for a long time, forming a flat "plateau" in the data before finally fading out.

This paper proposes a new, simpler explanation for that plateau. It suggests we've been looking at the fireworks wrong. Instead of a single, uniform shell of debris, the explosion is actually a layered cake of debris, with fast layers on the outside and slower layers on the inside.

Here is the breakdown of the paper's idea using everyday analogies:

1. The Old Problem: The "Magic Battery" Mystery

Previously, astronomers thought the plateau happened because the central engine (the "battery" of the explosion) kept pumping energy into the blast long after the initial boom.

The Analogy: Imagine a car crash. The car hits a wall and stops. But then, for 10 minutes, the car keeps glowing brightly. To explain this, scientists had to assume the car had a secret, infinite battery that kept recharging itself.
The Problem: This "secret battery" theory had issues. It required the explosion to be incredibly inefficient at first, then suddenly become super efficient later. It felt like a "magic trick" rather than physics.

2. The New Idea: The "Traffic Jam" of Debris

The authors, Sadeh, Hotokezaka, and Shibata, suggest the explosion isn't a single shell. It's a stratified outflow.

The Analogy: Imagine a convoy of cars leaving a city.
- The fastest cars (the ultra-relativistic ejecta) zoom ahead. These are the ones that create the initial gamma-ray boom.
- Behind them are slower cars, and behind those are even slower cars.
- The "road" they are driving on is the empty space of the universe (the external medium).

When the fastest cars hit the "wall" of the external medium, they create a shockwave (the Forward Shock). This is what we see in X-rays.

But here is the twist: The slower cars are still catching up. As they catch up to the fast cars that are already slowing down, they crash into them from behind. This creates a second shockwave moving backward through the debris (the Reverse Shock).

3. How the Plateau is Formed: The "Energy Refill"

Because the slower cars are constantly catching up and crashing into the front group, they keep dumping their kinetic energy into the main shockwave.

The Analogy: Think of a runner (the shockwave) who is getting tired and slowing down. But, a steady stream of fresh runners (the slower debris) keeps tapping them on the shoulder and pushing them forward.
The Result: The runner doesn't slow down as fast as they should. They maintain a steady, slow pace for a long time. In the data, this looks like a flat plateau.
Why it ends: Eventually, the slowest car in the convoy catches up. Once the last bit of debris has crashed into the front, the "pushing" stops. The runner finally slows down naturally, and the light curve drops off steeply, just like a normal explosion.

4. The "Millimeter" Surprise

The paper also predicts something cool about what we should see in other colors of light.

The Analogy: The front shock (the fast cars) is hot and energetic, so it glows brightly in high-energy light (X-rays). The backward shock (the slower cars catching up) is cooler and less energetic.
The Prediction: While the X-rays come from the front, the backward shock should be glowing very brightly in millimeter waves (a type of radio light).
The Takeaway: If you point a radio telescope at a GRB with a plateau, you should see a massive, long-lasting glow in the millimeter band that is much brighter than the X-rays. This is a "smoking gun" test for their theory.

5. Why This Matters

This model is elegant because it doesn't need "magic batteries" or weird geometry.

Physical Reality: It makes sense that an explosion wouldn't launch everything at the exact same speed. Just like a real explosion, some bits go faster than others.
Unified Picture: The same debris that makes the initial gamma-ray flash is also responsible for the long X-ray plateau. It's all one continuous process of a layered explosion interacting with space.

Summary

The paper argues that the mysterious "X-ray plateau" isn't a sign of a mysterious engine keeping the light on. Instead, it's a traffic jam in space. The explosion sends out debris at different speeds; the slow stuff catches up to the fast stuff, keeping the shockwave alive and glowing for hours. Once the slowest debris catches up, the light finally fades away, and we can see the "millimeter glow" of the backward shock as proof of the theory.

Here is a detailed technical summary of the paper "The hydrodynamics of stratified ultra-relativistic outflows and the origin of GRB X-ray plateaus" by Sadeh, Hotokezaka, and Shibata.

1. Problem Statement

The origin of the X-ray plateau phase observed in a significant fraction of Gamma-Ray Burst (GRB) afterglows remains a major unsolved problem in high-energy astrophysics.

Observational Context: Swift satellite observations revealed that many GRB afterglows exhibit a prolonged phase of shallow flux decay ( $F_X \propto t^{-\alpha}$ , with $0.3 \lesssim \alpha \lesssim 0.6 $) lasting$ 10^3 $–$ 10^5$ seconds, followed by a transition to a steeper decay.
Limitations of Existing Models:
- Energy Injection: The standard "refreshed shock" model requires late-time energy injection from a central engine (e.g., magnetar spin-down or fallback accretion). This often leads to energy budget inconsistencies, implying radiative efficiencies far below canonical values ( $\eta \sim 10\%$ ) or requiring unphysically high prompt efficiencies.
- Geometric Effects: Models invoking high-latitude emission or off-axis viewing of structured jets struggle to explain the energetics and the smooth transition to standard afterglow decay without fine-tuning.
- Coasting Phase: Models relying on low initial Lorentz factors ( $\Gamma \sim 40$ ) in wind environments conflict with constraints from prompt $\gamma$ -ray emission, which requires ultra-relativistic outflows ( $\Gamma \gtrsim 100$ ).

The authors propose that the plateau arises not from new energy injection or geometry, but from the hydrodynamic evolution of a radially stratified ultra-relativistic (UR) outflow interacting with the external medium.

2. Methodology

The authors developed a novel analytic framework to describe the interaction between a stratified ejecta shell and an external medium, moving beyond the standard "single Lorentz factor" shell approximation.

Stratified Ejecta Model: The ejecta is modeled as a sequence of shells with a continuous distribution of Lorentz factors ( $\Gamma$ ). The cumulative mass profile is defined as $M_{iso}(>\Gamma) \propto \Gamma^{-s}$ , where $s$ is the stratification index. This accounts for the physical expectation that GRB engines launch material with a range of velocities due to internal variability.
Hydrodynamic Derivation:
- They derived closed-form analytic expressions for the Forward Shock (FS) and Reverse Shock (RS) dynamics.
- Crucially, they accounted for a continuous distribution of Lorentz factors within the ejecta, allowing the RS to be long-lived and mildly relativistic rather than Newtonian or short-lived.
- They solved the shock jump conditions and energy conservation equations to determine the evolution of the shocked plasma Lorentz factors ( $\gamma_2$ for FS, $\bar{\gamma}_3$ for RS) and the RS crossing time ( $t_{cross}$ ).
Synchrotron Emission Calculation: Using standard microphysical assumptions (power-law electron distribution, constant fractions of energy in electrons $\epsilon_e$ and magnetic fields $\epsilon_B$ ), they computed the synchrotron break frequencies ( $\nu_m, \nu_c, \nu_a$ ) and fluxes for both the FS and RS across different spectral regimes.
External Medium: The analysis considers both a uniform Interstellar Medium (ISM, $k=0$ ) and a wind-like profile ( $k=2$ ).

3. Key Contributions

Analytic Solution for Stratified Flows: The paper provides the first compact analytic expressions describing the hydrodynamics of UR outflows with a broad Lorentz factor distribution, recovering standard Blandford-McKee (BM) solutions in appropriate limits.
Unified Origin for Prompt and Plateau: The model demonstrates that the same outflow responsible for the prompt $\gamma$ -ray emission (requiring $\Gamma \gtrsim 100$ ) naturally produces the X-ray plateau through its internal stratification, eliminating the need for prolonged central-engine activity.
Prediction of Millimeter Emission: A novel prediction is that the long-lived RS in stratified outflows generates a bright, long-lived emission component in the millimeter (mm) band, which outshines the FS at these wavelengths during the plateau phase.

4. Key Results

A. Hydrodynamic Evolution

RS Crossing Time: The time required for the RS to cross the ejecta is significantly extended compared to thin-shell approximations. For moderate stratification indices ($2 \lesssim s \lesssim 4 $), the crossing time is increased by an order of magnitude ($ t_{cross} \propto h(s)$).
Deceleration Law: As the RS propagates through the stratified ejecta, it gradually transfers energy to the blast wave. This results in a smooth, shallow deceleration of the FS ( $\gamma_2 \propto t^{-\alpha}$ ) rather than the abrupt transition seen in single-shell models.
Transition to BM: Once the slowest ejecta are processed, the system smoothly transitions to the standard adiabatic Blandford-McKee self-similar evolution.

B. Reproducing X-ray Plateaus

Slope Consistency: For a uniform ISM ( $k=0$ ) and electron index $p \approx 2.0-2.4$ , a stratification index of $s \approx 2.5 - 4$ naturally yields a temporal decay index $\alpha_X \approx 0.3 - 0.6$ , matching observed plateau slopes.
Duration and $\Gamma_{min}$ : The plateau duration ( $t_{break} \sim 10^3-10^4$ s) constrains the minimum Lorentz factor of the ejecta to $\Gamma_{min} \gtrsim 100$ . This is consistent with the UR outflow required for prompt emission, resolving the tension with low- $\Gamma$ coasting models.
Dainotti Relation: The model naturally reproduces the observed anti-correlation between plateau luminosity and duration ( $L_X \propto t_{break}^{-1}$ ) without fine-tuning.
Spectral Regime: The model is consistent with X-ray spectra lying above both the peak ( $\nu_m$ ) and cooling ( $\nu_c$ ) frequencies, requiring $\epsilon_B \gtrsim 10^{-3}$ .

C. Reverse Shock Emission

Millimeter Brightness: Because the RS processes a larger mass of electrons (sweeping the ejecta) and has lower characteristic electron energies than the FS, its synchrotron peak shifts to lower frequencies.
Prediction: During the plateau, the RS emission in the millimeter band is predicted to be >10 times brighter than the FS emission. The RS contribution is negligible in the X-ray band.
Light Curve Behavior: The mm emission persists throughout the plateau and fades rapidly once the RS crossing is complete, providing a distinct temporal signature.

5. Significance

Physical Economy: This model offers a unified, physically motivated explanation for the prompt emission, X-ray plateau, and subsequent afterglow using a single outflow structure. It removes the need for ad-hoc energy injection mechanisms or extreme geometric fine-tuning.
Energy Budget Resolution: By attributing the plateau energy to the kinetic energy already present in the stratified ejecta, the model avoids the radiative efficiency paradoxes associated with late-time energy injection models.
Testable Predictions: The paper provides a clear, falsifiable prediction: GRBs with X-ray plateaus should exhibit a bright, long-lived millimeter component during the plateau phase. This can be tested with current and future facilities like ALMA or AtLAST. The absence of such mm emission would constrain the stratification index or microphysical parameters.
Theoretical Advancement: The analytic framework bridges the gap between idealized single-shell models and complex numerical simulations, offering a robust tool for interpreting GRB afterglow data.

In conclusion, Sadeh et al. demonstrate that radial stratification of ultra-relativistic ejecta is a sufficient and natural mechanism to explain the X-ray plateau phenomenon, linking the hydrodynamics of the outflow directly to the observed light curve morphology.