Three Dimensional Radiation Hydrodynamic Simulations of Massive Star Envelopes

Abstract

(Abridged) Stars more massive than $20-30M_{\odot }$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff}\approx 2-4\times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff}\approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5}{M}_{\odot }/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.

This is the authors' original version. The paper is published in the September 27th issue of Nature: https://rdcu.be/7PgW

References in corpus (8)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Convective Boundaries, Element Diffusion, and Massive Star Explosions
• A Global Three Dimensional Radiation Magneto-hydrodynamic Simulation of Super-Eddington Accretion Disks
• Luminous Blue Variables are Antisocial: Their Isolation Implies that they are Kicked Mass Gainers in Binary Evolution
• An Algorithm for Radiation Magnetohydrodynamics Based on Solving the Time-dependent Transfer Equation
• Instabilities in the Envelopes and Winds of Very Massive Stars
• Spectroscopic and photometric oscillatory envelope variability during the S Doradus outburst of the Luminous Blue Variable R71
• Pulsations as a Driver for LBV Variability

Cited by in corpus (63)

• The Athena++ Adaptive Mesh Refinement Framework: Design and Magnetohydrodynamic Solvers
• Different to the core: the pre-supernova structures of massive single and binary-stripped stars
• Formation of GW190521 from stellar evolution: the impact of the hydrogen-rich envelope, dredge-up and ${}^{12}$C($\alpha$, $\gamma$)${}^{16}$O rate on the pair-instability black hole mass gap
• It has to be cool: on supergiant progenitors of binary black hole mergers from common-envelope evolution
• The redshift evolution of the binary black hole merger rate: a weighty matter
• Magnetohydrodynamic Simulations of Active Galactic Nucleus Disks and Jets
• Stellar Evolution in AGN Disks
• Maximum Black Hole mass across Cosmic Time
• Is GW190521 the merger of black holes from the first stellar generations?
• Sensitivity of the lower-edge of the pair instability black hole mass gap to the treatment of time dependent convection
• Catalogue of new Herbig Ae/Be and classical Be stars. A machine learning approach to Gaia DR2
• Long-term evolution of a magnetic massive merger product
• Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes
• JWST Imaging of Earendel, the Extremely Magnified Star at Redshift $z=6.2$
• The Stellar Merger Scenario for Black Holes in the Pair-instability Gap
• Accretion onto Stars in the Disks of Active Galactic Nuclei
• First 3D Radiation-Hydrodynamic Simulations of Wolf-Rayet Winds
• On the Origin of Stochastic, Low-Frequency Photometric Variability in Massive Stars
• Opacity Driven Convection and Variability in Accretion Disks around Supermassive Black Holes
• The fates of massive stars: exploring uncertainties in stellar evolution with METISSE
• Pre-supernova evolution and final fate of stellar mergers and accretors of binary mass transfer
• Wind-envelope interaction as the origin of the slow cyclic brightness variations of luminous blue variables
• On the Gaia DR2 distances for Galactic Luminous Blue Variables
• All at once: transient pulsations, spin down and a glitch from the Pulsating Ultraluminous X-ray Source M82 X-2
• Stochastic Low Frequency Variability in 3-Dimensional Radiation Hydrodynamical Models of Massive Star Envelopes
• It's written in the massive stars: The role of stellar physics in the formation of black holes
• Supernova Precursor Emission and the Origin of Pre-Explosion Stellar Mass-Loss
• An Atlas of Convection in Main-Sequence Stars
• Convectively Driven Three Dimensional Turbulence in Massive Star Envelopes: I. A 1D Implementation of Diffusive Radiative Transport
• Short Term Variability of Evolved Massive Stars with TESS
• A new and unusual LBV-like outburst from a Wolf-Rayet star in the outskirts of M33
• The Role of Inhomogeneities in Supernova Shock Breakout Emission
• MOCMC: Method of Characteristics Moment Closure, a Numerical Method for Covariant Radiation Magnetohydrodynamics
• X-Shooting ULLYSES: Massive Stars at low metallicity IX: Empirical constraints on mass-loss rates and clumping parameters for OB supergiants in the Large Magellanic Cloud
• Photometric detection of internal gravity waves in upper main-sequence stars. III. Comparison of amplitude spectrum fitting and Gaussian process regression using celerite2
• Gap Transients Interacting with Circumstellar Medium
• A systematic study of super-Eddington layers in the envelopes of massive stars
• Centrifugally Driven Mass Loss and Outbursts of Massive Stars
• Common Envelope Evolution of Massive Stars
• On the apparent absence of WR+NS systems: the curious case of WR124
• Dissecting the microphysics behind the metallicity-dependence of massive stars radii
• Method and new tabulations for flux-weighted line-opacity and radiation line-force in supersonic media
• Empirical mass-loss rates and clumping properties of O-type stars in the LMC
• Turbulence Supported Massive Star Envelopes
• Theoretical wind clumping predictions from 2D LDI models of O-star winds at different metallicities
• SN 1961V: A Pulsational Pair-Instability Supernova
• Repeating periodic eruptions of the supernova impostor SN 2000ch
• The landscape of binary core-collapse supernova progenitors and the late emergence of Wolf-Rayet winds
• Synthesizing Spectra from 3D Radiation Hydrodynamic Models of Massive Stars Using Monte Carlo Radiation Transport
• Three Dimensional Natures of Massive Star Envelopes
• Fast Ejecta Resulted from Jet-Wind Interaction in the Great Eruption of Eta Carinae
• Variability of massive stars in M31 from the Palomar Transient Factory
• The Influence of Disk Composition on the Evolution of Stars in the Disks of Active Galactic Nuclei
• Luminous Blue Variables
• 13 Years of P Cygni Spectropolarimetry: Investigating Mass-loss Through H$\alpha$, Periodicity, and Ellipticity
• Five years of BRITE-Constellation photometry of the luminous blue variable P Cygni: properties of the stochastic low-frequency variability
• How long can LBVs sleep? A long-term photometric vaiability and spectral study of the Galactic candidate luminous blue variable MN112
• The drastic impact of Eddington-limit induced mass ejections on massive star populations
• A constant upper luminosity limit of cool supergiant stars down to the extremely low metallicity of I Zw 18
• A multi-frequency, multi-epoch radio continuum study of the Arches cluster with the VLA
• Spherically symmetric model atmospheres using approximate lambda operators V. Static inhomogeneous atmospheres of hot dwarf stars
• Viscous and centrifugal instabilities of massive stars
• Radiation-Driven Stellar Eruptions