X-ray flares from propagation instabilities in long Gamma-Ray Burst jets

We present a numerical simulation of a gamma-ray burst jet from a long-lasting engine in the core of a 16 solar mass Wolf-Rayet star. The engine is kept active for 6000 s with a luminosity that decays in time as a power-law with index -5/3. Even though there is no short time-scale variability in the injected engine luminosity, we find that the jet’s kinetic luminosity outside the progenitor star is characterized by fluctuations with relatively short time scale. We analyze the temporal characteristics of those fluctuations and we find that they are consistent with the properties of observed flares in X-ray afterglows. The peak to continuum flux ratio of the flares in the simulation is consistent with some, but not all, the observed flares. We propose that propagation instabilities, rather than variability in the engine luminosity, are responsible for the X-ray flares with moderate contrast. Strong flares such as the one detected in GRB 050502B, instead, cannot be reproduced by this model and require strong variability in the engine activity.

💡 Research Summary

The authors present a high‑resolution relativistic hydrodynamic simulation of a long‑duration gamma‑ray burst (GRB) jet launched from the core of a 16 M☉ Wolf‑Rayet star. The central engine is assumed to operate continuously for 6000 s, with its total power declining as a power‑law L(t) ∝ t⁻⁵ᐟ³, a functional form motivated by fallback accretion models. Crucially, no short‑timescale variability is imposed on the injected luminosity; the engine is deliberately kept smooth to test whether the observed X‑ray flares in GRB afterglows can arise from jet dynamics alone.

The simulation follows the jet as it drills through the dense stellar envelope, emerges at the stellar surface, and propagates into the low‑density circumstellar medium. During this transition the jet experiences strong pressure mismatches, shear layers, and internal shocks that give rise to what the authors term “propagation instabilities.” These instabilities modulate the kinetic luminosity of the jet on timescales far shorter than the engine’s secular decay—typically from a few seconds up to a few tens of seconds. The authors extract the kinetic luminosity at a radius well outside the progenitor and identify a series of distinct spikes superimposed on the smooth power‑law decline.

A statistical analysis of the spikes shows that (i) they appear preferentially between a few hundred and a few thousand seconds after engine ignition, (ii) their duration Δt scales roughly as 0.1–0.3 × t (where t is the time of occurrence), and (iii) the peak‑to‑continuum flux ratio lies in the range 2–10. These properties match the bulk of the X‑ray flares observed by Swift/XRT, which typically show Δt/t ≈ 0.1 and moderate contrast with the underlying afterglow. The simulated flares therefore provide a natural explanation for “moderate‑contrast” flares without invoking any intrinsic variability in the central engine.

However, the model fails to reproduce the most extreme flares, such as the bright, high‑contrast event seen in GRB 050502B, where the peak flux exceeds the underlying afterglow by a factor > 20. The authors argue that such events must involve genuine engine variability—perhaps episodic accretion or magnetic reconnection—on top of the propagation‑induced modulation.

The paper discusses several caveats. The calculations are performed primarily in two‑dimensional axisymmetry, which may suppress fully three‑dimensional turbulent modes and underestimate the amplitude of instabilities. Magnetic fields are omitted, so magnetohydrodynamic effects (e.g., kink instabilities) that could enhance variability are not captured. Radiative transfer is not modeled; the conversion from kinetic to X‑ray luminosity is assumed rather than computed self‑consistently. Finally, numerical resolution limits the smallest scales of shear‑driven Kelvin‑Helmholtz turbulence that can be resolved.

In summary, the study demonstrates that propagation instabilities inherent to a relativistic jet breaking out of a massive star can generate X‑ray flares with timing and contrast comparable to the majority of observed afterglow flares. This mechanism offers a compelling alternative to engine‑driven variability for moderate‑contrast flares, while still allowing for a hybrid scenario where the most energetic flares require genuine central‑engine outbursts. The work thus advances our understanding of the diversity of GRB afterglow phenomenology and highlights the need for future three‑dimensional, magnetized, radiation‑hydrodynamic simulations to fully capture the flare‑generation process.