On the circumstellar medium of massive stars and how it may appear in GRB observations

Massive stars lose a large fraction of their original mass over the course of their evolution. These stellar winds shape the surrounding medium according to parameters that are the result of the characteristics of the stars, varying over time as the stars evolve, leading to both permanent and temporary features that can be used to constrain the evolution of the progenitor star. Because long Gamma-Ray Bursts (GRBs) are thought to originate from massive stars, the characteristics of the circumstellar medium (CSM) should be observable in the signal of GRBs. This can occur directly, as the characteristics of the GRB-jet are changed by the medium it collides with, and indirectly because the GRB can only be observed through the extended circumstellar bubble that surrounds each massive star. We use computer simulations to describe the circumstellar features that can be found in the vicinity of massive stars and discuss if, and how, they may appear in GRB observations. Specifically, we make hydrodynamical models of the circumstellar environment of a rapidly rotating, chemically near-homogeneous star, which is represents a GRB progenitor candidate model. The simulations show that the star creates a large scale bubble of shocked wind material, which sweeps up the interstellar medium in an expanding shell. Within this bubble, temporary circumstellar shells, clumps and voids are created as a result of changes in the stellar wind. Most of these temporary features have disappeared by the time the star reaches the end of its life, leaving a highly turbulent circumstellar bubble behind. Placing the same star in a high density environment simplifies the evolution of the CSM as the more confined bubble prohibits the formation of some of the temporary structures.

💡 Research Summary

This paper investigates how the circumstellar medium (CSM) surrounding massive stars—particularly those that are candidate progenitors of long gamma‑ray bursts (GRBs)—affects observable GRB signatures. The authors focus on a rapidly rotating, chemically near‑homogeneous 16 M⊙ stellar model taken from Yoon et al. (2006). Using the MPI‑AMRVAC hydrodynamics code with optically thin radiative cooling, they perform a two‑stage numerical experiment: an initial one‑dimensional (1‑D) spherical simulation to capture the large‑scale wind‑blown bubble, followed by a two‑dimensional (2‑D) simulation that resolves latitude‑dependent wind anisotropies during the star’s critical‑rotation phase.

In the 1‑D phase, the stellar wind creates a hot (∼10⁷ K) shocked‑wind bubble of radius ≈140 pc, bounded by a thin swept‑up interstellar shell. Two ambient density environments are explored: a low‑density ISM (2 cm⁻³, “Simulation A”) and a high‑density molecular cloud (2000 cm⁻³, “Simulation B”). The high‑density case experiences stronger radiative cooling, leading to a more compact bubble and a denser outer shell.

When the star approaches critical rotation (Ω≈0.99), its mass‑loss rate spikes while wind velocity drops dramatically. The authors map the 1‑D results onto a 2‑D r‑θ grid and apply the gravity‑darkening prescription of Dwarkadas & Owocki (2002) to impose latitude‑dependent mass loss and wind speed (higher at the poles, lower at the equator). This anisotropy drives the formation of a bipolar shell that expands preferentially along the polar axis. After the critical‑rotation episode ends, the wind speed recovers, generating a third, faster shell that eventually overtakes the bipolar shell. In the low‑density simulation, both shells survive for tens of thousands of years, interact, and finally dissolve into a highly turbulent shocked‑wind bubble. In the high‑density simulation, the smaller bubble forces the bipolar shell to collide with the outer swept‑up shell early, preventing the formation of the third shell and leaving a comparatively smooth interior.

The authors discuss two observational consequences. First, the density structure of the CSM directly influences the dynamics of the GRB jet. In low‑density environments, the jet traverses regions of varying density, potentially producing observable changes in afterglow light curves and spectra as kinetic energy is transferred to the shocked medium. In high‑density environments, the more uniform CSM leads to a steadier jet evolution and less pronounced afterglow variability. Second, the absorption spectrum of the afterglow can carry signatures of the CSM. By the time of core collapse, the transient shells have largely dissipated, leaving only the free‑streaming wind close to the star. This wind can imprint blue‑shifted, high‑ionisation absorption lines (e.g., C IV, Si IV) that are formed very near the GRB source, where only highly ionised species survive. The outer swept‑up shell moves at only ~25 km s⁻¹, making it indistinguishable from the ambient ISM in absorption. Consequently, GRB afterglow spectra are expected to show at most a single high‑ionisation component, with little evidence of the earlier, more massive shells.

In summary, the paper demonstrates that (1) the evolution of massive‑star winds, especially during episodes of near‑critical rotation, creates temporary circumstellar structures (bipolar shells, clumps, voids); (2) the ambient ISM density critically controls whether these structures survive to the moment of core collapse; and (3) the surviving CSM imprint on GRB observations is limited to modest density variations affecting jet dynamics and a narrow set of high‑ionisation absorption features. These findings provide a framework for interpreting GRB afterglow data in terms of progenitor wind histories and surrounding environmental conditions.