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

TL;DR

This study uses advanced radiation-hydrodynamic simulations to model the atmospheres and winds of luminous blue variable stars during their quiescent and eruptive phases, revealing key differences in their properties.

Contribution

It provides the first unified 2D and 1D RHD models that simulate the transition from quiescent to eruptive LBV phases, capturing the complex wind and atmospheric behaviors.

Findings

Models reproduce observed LBV wind speeds and mass-loss rates.

Simulations show a natural transition from turbulent to super-Eddington outflows.

Results suggest different stellar surface conditions during quiescent and eruptive phases.

Abstract

Mass loss from massive stars located in the part of the Hertzsprung-Russell diagram (HRD) where we find luminous blue variables (LBVs) is profoundly important for stellar evolution yet poorly understood. We use time-dependent radiation-hydrodynamic (RHD) simulations to examine the atmosphere and wind properties of such massive stars, computing 2D and 1D RHD models of the coupled envelopes, atmospheres, and wind outflows, tuned to this region in the HRD. Our unified simulations start deep in the stellar envelope (well below T ~ 200 kK) and include the outflowing wind, accounting for line-driving, radiative enthalpy, and photon tiring. Mass-loss rates, wind speeds, and the radiative luminosity at the photosphere are emergent properties in the simulations. A grid of models is created by slightly increasing the stellar energy at the lower boundary. This results in a natural transition from very turbulent atmospheres with line-driven winds to effectively stationary super-Eddington massive outflows. Our sub-Eddington models are essentially blue hypergiant stars with very variable surfaces, effective mass-loss rates $\dot{M} \sim 2 - 5 \times 10^{-5}$ $M_{\odot}$/year, and wind speeds $v_{\infty} \sim 200 - 300$ km/s, resembling quiescent LBVs like P Cygni. The super-Eddington models have optically thick wind envelopes and extremely inflated yellow surfaces (Teff ~ 5000 K), $\dot{M} \sim 0.1 - 1$ $M_{\odot}$/year, and $v_{\infty} \sim 400 - 500$ km/s, resembling a massive star during a great eruption like eta Carinae's. Our models naturally reproduce the overall characteristic stellar and wind parameters inferred for massive stars in their quiescent LBV and yellow giant eruptive phases. It remains an open question whether the energy increase needed to trigger a giant eruption can be obtained solely by the internal evolution of the star itself or if it requires an external energy source.