Collapsar Accretion and the Gamma-Ray Burst X-Ray Light Curve

We present axisymmetric hydrodynamical simulations of the long-term accretion of a rotating GRB progenitor star, a “collapsar,” onto the central compact object. The simulations were carried out with the adaptive mesh refinement code FLASH in two spatial dimensions and with an explicit shear viscosity. The evolution of the central accretion rate exhibits phases reminiscent of the long GRB gamma-ray and X-ray light curve, which lends support to the proposal that the luminosity is modulated by the central accretion rate. After a few tens of seconds, an accretion shock sweeps outward through the star. The formation and outward expansion of the accretion shock is accompanied with a sudden and rapid power-law decline in the central accretion rate Mdot ~ t^{-2.8}, which resembles the L_X ~ t^{-3} decline observed in the X-ray light curves. The collapsed, shock-heated stellar envelope settles into a thick, low-mass equatorial disk embedded within a massive, pressure-supported atmosphere. After a few hundred seconds, the inflow of low-angular-momentum material in the axial funnel reverses into an outflow from the surface of the thick disk. Meanwhile, the rapid decline of the accretion rate slows down, or even settles a in steady state with Mdot ~ 5x10^{-5} Msun/s, which resembles the “plateau” phase in the X-ray light curve. While the duration of the “prompt” phase depends on the resolution in our simulations, we provide an analytical model taking into account neutrino losses that estimates the duration to be ~20 s. The model suggests that the steep decline in GRB X-ray light curves is triggered by the circularization of the infalling stellar envelope at radii where the virial temperature is below ~10^{10} K, such that neutrino cooling shuts off and an outward expansion of the accretion shock becomes imminent.

💡 Research Summary

The paper presents a comprehensive numerical study of the long‑term accretion of a rotating massive star (a “collapsar”) onto a central compact object, using the adaptive‑mesh‑refinement (AMR) code FLASH in two‑dimensional axisymmetric geometry. An explicit shear viscosity term is added to model angular‑momentum transport and viscous heating within the forming accretion flow. The authors follow the evolution from the initial free‑fall of the stellar envelope through the formation of a centrifugally supported disk, the emergence of an accretion shock, and the subsequent development of distinct phases in the central mass‑accretion rate (Ṁ).

Key results and physical interpretation

1. Prompt‑phase behavior (first few tens of seconds). As the high‑angular‑momentum material circularizes, a strong accretion shock propagates outward through the star. When the shock reaches radii where the post‑shock virial temperature drops below ~10¹⁰ K, neutrino cooling becomes inefficient. The loss of this cooling channel causes the shocked gas to become pressure‑supported, driving rapid outward expansion of the shock. During this interval the central accretion rate declines steeply, following an approximate power law Ṁ ∝ t⁻²·⁸. This decline mirrors the observed X‑ray light‑curve steep‑decay segment (L_X ∝ t⁻³) seen in many long GRBs.

2. Transition to a plateau (hundreds of seconds). After the shock has traversed most of the stellar envelope, the low‑angular‑momentum material that continues to flow down the polar funnel reverses direction and is expelled from the surface of the thick, low‑mass equatorial disk that remains embedded in a massive, pressure‑supported atmosphere. The accretion rate then settles to a quasi‑steady value of order 5 × 10⁻⁵ M_⊙ s⁻¹. This plateau in Ṁ reproduces the “plateau” phase of the X‑ray afterglow, during which the luminosity decays only slowly or remains roughly constant for 10³–10⁴ s.

3. Analytical model for the prompt‑phase duration. The authors construct a simple analytic estimate that includes neutrino cooling losses. By equating the cooling time to the dynamical time at the circularization radius, they find that the prompt phase should last ≈20 s for the progenitor model used. This timescale is consistent with the typical duration of the γ‑ray emission in long GRBs.

4. Resolution dependence. The length of the prompt phase in the simulations is found to depend on numerical resolution: higher resolution yields a shorter, sharper decline. The authors argue that the physical prompt duration is set by the neutrino‑cooling cutoff rather than by numerical diffusion, and that their analytic estimate provides a resolution‑independent benchmark.

Methodological strengths

• Use of explicit viscosity allows a realistic treatment of angular‑momentum transport, which is often approximated by an α‑prescription or omitted entirely in earlier collapsar simulations.
• The AMR approach resolves the shock front and the thin equatorial disk while still covering the entire stellar radius, a demanding computational task.
• The study directly connects central engine physics (Ṁ(t)) to observable X‑ray light‑curve features, providing a testable prediction that the steep‑decay slope should be close to –3.

Limitations and future directions

• The calculations are axisymmetric; three‑dimensional effects such as non‑axisymmetric instabilities, magnetic fields, and jet launching are not captured.
• Magnetic stresses, which could dominate angular‑momentum transport in the inner disk, are omitted. Incorporating magnetohydrodynamics (MHD) would be essential to assess jet collimation and the role of magnetic braking.
• The outer envelope is treated with a simplified equation of state; detailed nuclear burning and composition changes could affect the shock propagation.

Overall significance
The work provides a physically motivated, quantitative bridge between the dynamics of a collapsing massive star and the multi‑phase structure of long‑GRB X‑ray afterglows. By identifying the shutdown of neutrino cooling at a critical temperature as the trigger for the steep‑decay phase, the authors offer a concrete mechanism that can be incorporated into broader GRB central‑engine models. The agreement between simulated Ṁ(t) and observed light‑curve segments (prompt, steep decline, plateau) strengthens the collapsar paradigm and suggests that future high‑resolution, three‑dimensional MHD simulations will be able to reproduce the full diversity of GRB afterglow phenomenology.