🔭 astrophysics

Simulations of massive star atmospheres and winds during giant eruptive and quiescent luminous blue variable phases

This paper presents time-dependent radiation-hydrodynamic simulations that successfully reproduce the distinct atmospheric and wind properties of massive stars in both quiescent luminous blue variable phases and giant eruptive phases, demonstrating how a gradual increase in stellar energy naturally transitions the star from a turbulent, line-driven wind state to a super-Eddington, optically thick outflow resembling a great eruption.

Original authors: P. Schillemans, J. O. Sundqvist, D. Debnath, L. Delbroek, N. Moens, C. Van der Sijpt

Published 2026-03-24

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: P. Schillemans, J. O. Sundqvist, D. Debnath, L. Delbroek, N. Moens, C. Van der Sijpt

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine a massive star not as a static, glowing ball of fire, but as a living, breathing giant that is constantly struggling to keep its cool. This paper is about simulating what happens when these giants get "feverish"—when they gain a little too much energy and start to erupt.

The researchers used powerful computer simulations to watch these stars from deep inside their cores all the way out into the winds they blow into space. They discovered that depending on how much extra energy the star has, it behaves in two very different ways: like a turbulent, churning ocean or like a giant, slow-motion explosion.

Here is the breakdown of their findings using simple analogies:

1. The Setup: The Star's "Thermostat"

The team started with a model of a massive star (about 60 times heavier than our Sun). They then played with a "thermostat" at the bottom of the star's atmosphere.

The "Cool" Setting: They kept the energy relatively low.
The "Hot" Setting: They cranked the energy up significantly.

They wanted to see how the star's "skin" (atmosphere) and its "breath" (stellar wind) reacted to these changes.

2. Scenario A: The Turbulent Ocean (The "Cool" Star)

When the star has a moderate amount of energy, it behaves like a boiling pot of water or a stormy sea.

What happens: Deep inside the star, there is a layer where the gas gets "sticky" (due to iron opacity). This causes the gas to churn violently. It tries to shoot upward like a geyser, but it doesn't have enough speed to escape the star's gravity.
The Cycle: The gas shoots up, slows down, and then crashes back down like a failed fountain. This crash heats up the layer below, causing another geyser to shoot up. It's a constant cycle of "failed winds."
The Result: The star's surface is a mess of turbulence. However, the outer layers are still hot enough that light pressure can push some gas away, creating a steady, fast wind (like a strong jet stream).
Real-life match: This looks like a "quiet" Luminous Blue Variable star, such as P Cygni. It's active and variable, but not exploding.

3. Scenario B: The Giant Inflating Balloon (The "Hot" Star)

When the researchers turned up the energy, the star didn't just get hotter; it changed its entire personality. It stopped churning and started inflating like a giant balloon.

The "Photon Tiring" Effect: Imagine the star is trying to push a heavy load (its own atmosphere) out into space. If the load is too heavy, the energy gets "tired" and used up just trying to lift the gas, rather than speeding it up.
The Transformation: Because the star is so energetic, it pushes a massive amount of gas outward. This gas is so thick and dense that it traps the light, creating a huge, puffy envelope. The star's surface cools down (turning yellow/orange) because it has expanded so much.
The Wind: Instead of a fast jet, the star blows a slow, thick, massive wind. It's like a slow-moving fog rolling off a mountain, but on a cosmic scale.
Real-life match: This perfectly mimics a star during a "Great Eruption," like the famous Eta Carinae in the 1800s, which blew off several suns' worth of mass in a single event.

4. The Big Discovery: It's All One Spectrum

The most exciting part of the paper is that they didn't have to program two different types of stars. By simply turning up the "energy dial" on the same star, it naturally transitioned from the Turbulent Ocean (Scenario A) to the Inflating Balloon (Scenario B).

This suggests that stars like P Cygni and Eta Carinae aren't fundamentally different species. They are likely the same type of star, just at different stages of an energy crisis.

Low Energy: Turbulent, fast wind, hot surface.
High Energy: Smooth, slow wind, cool surface, massive mass loss.

5. The Mystery: Where Does the Extra Energy Come From?

The simulations show what happens when a star gets too much energy, but they can't explain why it gets that energy in the first place.

The Question: Does the star just naturally evolve to this point? Or does it need a "push," like crashing into another star (a merger) or having a violent interaction with a binary partner?
The Analogy: Think of a car. The simulation shows exactly what happens when you floor the gas pedal (the star erupts). But the simulation doesn't tell us if the driver (the star's internal evolution) pressed the pedal, or if someone else kicked the car from behind (an external merger).

Summary

In short, this paper uses super-computers to show that massive stars are incredibly sensitive to their internal energy. A small increase can turn a churning, hot star into a puffy, erupting giant. This helps astronomers understand why some stars are just "moody" (variable) while others go "berserk" (eruption), and suggests they might all be part of the same dramatic family.

1. Problem Statement

Massive stars in the Luminous Blue Variable (LBV) region of the Hertzsprung-Russell Diagram (HRD) exhibit extreme mass loss and instability, ranging from semi-periodic variability (S Doradus type) to violent giant eruptions (e.g., $\eta$ Carinae). The physical mechanisms governing the transition between these states, particularly the coupled dynamics of the stellar envelope, atmosphere, and wind, remain poorly understood.

Key Gaps: Previous 1D models often neglected critical physics like radiative enthalpy and photon tiring, while previous multi-dimensional (2D/3D) simulations were typically tuned to lower luminosities or different evolutionary stages.
Core Question: Can the transition from a "quiescent" blue hypergiant state to a "giant eruption" state occur naturally through internal stellar evolution (energy increase), or is an external trigger (e.g., merger) required? Furthermore, how do radiative enthalpy and line-driving interact to drive these distinct wind regimes?

2. Methodology

The authors employed time-dependent Radiation-Hydrodynamics (RHD) simulations using the MPI-AMRVAC code with a "box-in-a-star" approach.

Physics Included:
- Coupled Envelope & Wind: Simulations start deep within the stellar envelope (below the iron opacity bump) and extend well into the outflowing wind.
- Radiative Transfer: Solved in the co-moving frame (CMF) using the Flux-Limited Diffusion (FLD) approximation.
- Key Mechanisms: Included gas and radiative enthalpy, photon tiring, line-driving (using the Sobolev approximation with fixed line-force parameters $\bar{Q}=2000, \alpha=2/3$ ), and turbulent pressure.
- Opacity: Used a "hybrid opacity" combining Rosseland mean opacities (OPAL + Ferguson) and line opacities.
Initial Conditions:
- Based on a $60 M_\odot$ galactic star evolved with MESA to the LBV instability strip ( $57 M_\odot$ remaining).
- Lower boundary set at $R_c = 38.6 R_\odot$ with a fixed density and a variable temperature ( $T_c$ ) to control the emergent luminosity.
Model Grid:
- Created a grid of 7 models by varying the lower boundary temperature ( $T_c$ ) from 360 kK to 460 kK (relative to a baseline of 400 kK).
- 2D Models: Used for sub-Eddington and lower super-Eddington cases to capture turbulence.
- 1D Models: Used for high super-Eddington cases to extend the computational domain to resolve the photosphere and wind terminal velocity (up to $31,000 R_\odot$ ).

3. Key Contributions

Unified Simulation Framework: The first set of unified RHD models that naturally transition from turbulent, line-driven winds to stationary, super-Eddington outflows solely by increasing the stellar energy input.
Role of Radiative Enthalpy: Demonstrated that radiative enthalpy is the primary energy transport mechanism in the deep envelope of super-Eddington stars, keeping the deep layers sub-Eddington while the surface becomes highly super-Eddington due to the inefficiency of enthalpy transport near the surface.
Natural Transition: Showed that the distinct observational classes of LBVs (quiescent vs. eruptive) may not be separate categories but rather different points on a continuum driven by total stellar energy.

4. Results

A. Sub-Eddington Regime ("Coolplate" Models)

Characteristics: Represent quiescent LBVs (e.g., P Cygni).
Dynamics: Highly turbulent atmospheres driven by the iron opacity bump.
Wind Mechanism: "Failed winds" occur where gas is accelerated but falls back due to insufficient momentum. This triggers semi-periodic cycles.
Mass Loss: Driven by efficient line-driving in the turbulent photosphere.
Parameters:
- $T_{\text{eff}} \sim 10,000 - 20,000$ K.
- $\dot{M} \sim 2 - 5 \times 10^{-5} M_\odot/\text{yr}$ .
- $v_\infty \sim 200 - 300$ km/s.
- Photospheric radius oscillates by $\sim 20 - 200 R_\odot$ .

B. Super-Eddington Regime ("Hotplate" Models)

Characteristics: Represent giant eruptive phases (e.g., $\eta$ Carinae).
Dynamics: The atmosphere becomes effectively stationary and quasi-steady; turbulence is suppressed.
Wind Mechanism: Driven by radiative enthalpy launching a flow from deep within the envelope. As the flow moves outward, enthalpy transport becomes inefficient, causing the radiative flux (and local Eddington factor) to spike, driving a massive outflow.
Photosphere: The star becomes highly inflated. Hydrogen recombination occurs at $T \sim 5000 - 5500$ K, causing a sharp drop in opacity. This defines the photosphere, resulting in a "yellow" surface temperature regardless of the deep luminosity.
Parameters:
- $T_{\text{eff}} \sim 5,000$ K (recombined surface).
- $\dot{M} \sim 0.1 - 1 M_\odot/\text{yr}$ (extreme mass loss).
- $v_\infty \sim 400 - 500$ km/s.
- Photospheric radius inflates to $\sim 10,000 R_\odot$ .
Oscillations: Giant oscillations in the 1D models are triggered by the reionization of hydrogen in partial ionization zones, causing the photosphere to jump outward by a factor of $\sim 4$ in days.

C. The Transition

The transition between these two states is abrupt. As total stellar energy increases, the system does not gradually shift; it undergoes a sudden switch where the Eddington ratio at the photosphere jumps from $<1$ to $>10$ . This suggests that the "Great Eruption" is a natural consequence of the star exceeding a critical energy threshold where enthalpy-driven outflows take over line-driving.

5. Significance and Conclusions

Reproduction of Observations: The models successfully reproduce the observed parameters of both quiescent LBVs (P Cygni) and eruptive giants ( $\eta$ Carinae) without tuning to specific stars, validating the underlying physics of RHD in this regime.
Mass Loss Interpretation: The study highlights that empirical mass-loss rates derived from stationary models may be inaccurate for LBVs, as the simulations show significant "failed winds" where material is launched but falls back, leading to high instantaneous rates that do not contribute to net mass loss.
Origin of Eruptions: While the simulations prove that an energy increase can trigger a giant eruption naturally, the paper concludes that it remains an open question whether this energy increase arises solely from internal stellar evolution or requires external triggers (e.g., mergers, binary interactions).
Future Work: The authors note limitations regarding fixed line-force parameters and the assumption of a constant adiabatic index ( $\gamma_g = 5/3$ ), suggesting future work should include non-LTE corrections and feedback mechanisms between the envelope and the core.

In summary, this paper provides a robust physical framework linking the internal energy of massive stars to their external eruptive behavior, suggesting that the diversity of LBV phenomena is a continuum driven by the efficiency of radiative enthalpy transport versus line-driving.