A systematic study of super-Eddington layers in the envelopes of massive stars

Abstract

The proximity to the Eddington luminosity has been attributed as the cause of several observed effects in massive stars. Computationally, if the luminosity carried through radiation exceeds the local Eddington luminosity in the low-density envelopes of massive stars, it can result in numerical difficulties, inhibiting further computation of stellar models. This problem is exacerbated by the fact that very few massive stars are observed beyond the Humphreys-Davidson limit, the same region in the Hertzsprung-Russell diagram where the aforementioned numerical issues relating to the Eddington luminosity occur in stellar models. One-dimensional stellar evolution codes have to use pragmatic solutions to evolve massive stars through this computationally difficult phase. In this work, we quantify the impact of these solutions on the evolutionary properties of massive stars. We used the stellar evolution code MESA with commonly used input parameters for massive stellar models to compute the evolution of stars in the initial mass range of 10-110 M${}_{\odot }$ at one-tenth of solar metallicity. We find that numerical difficulties in stellar models with initial masses greater than or equal to 30 M${}_{\odot }$ cause these models to fail before the end of core helium burning. Recomputing these models using the same physical inputs but three different pragmatic solutions to treat the numerical instability, we find that the maximum radial expansion achieved by stars can vary by up to 2000 R${}_{\odot }$, while the remnant mass of the stars can vary by up to 14 M${}_{\odot }$ between the sets. These differences can have implications on studies such as binary population synthesis.

Revised version; Accepted for publication in A&A

References in corpus (30)

• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• Origin of the heavy elements in binary neutron-star mergers from a gravitational wave event
• Pre-Supernova Evolution of Massive Single and Binary Stars
• Mass loss from late-type WN stars and its Z-dependence: very massive stars approaching the Eddington limit
• The evolution of runaway stellar collision products
• Luminous Blue Variables are Antisocial: Their Isolation Implies that they are Kicked Mass Gainers in Binary Evolution
• A new mass-loss rate prescription for red supergiants
• The evolution of rotating very massive stars with LMC composition
• The VLT-FLAMES Tarantula Survey XVII. Physical and wind properties of massive stars at the top of the main sequence
• The IACOB project: III. New observational clues to understand macroturbulent broadening in massive O- and B-type stars
• Constraining mixing in massive stars in the Small Magellanic Cloud
• New predictions for radiation-driven, steady-state mass-loss and wind-momentum from hot, massive stars II. A grid of O-type stars in the Galaxy and the Magellanic Clouds
• On Variations Of Pre-Supernova Model Properties
• Maximum Black Hole mass across Cosmic Time
• Metallicity-dependent wind parameter predictions for OB stars
• Diverse Variability of O and B Stars Revealed from 2-minute Cadence Light Curves in Sectors 1 and 2 of the TESS Mission: Selection of an Asteroseismic Sample
• The metallicity dependence of envelope inflation in massive stars
• The fates of massive stars: exploring uncertainties in stellar evolution with METISSE
• Empirical mass-loss rates and clumping properties of Galactic early-type O supergiants
• The excess of cool supergiants from contemporary stellar evolution models defies the metallicity-independent Humphreys-Davidson limit
• Evolutionary models for R Coronae Borealis stars
• A theoretical investigation of the Humphreys-Davidson limit at high and low metallicity
• Instabilities in the Envelopes and Winds of Very Massive Stars
• Explaining the differences in massive star models from various simulations
• Physics and evolution of the most massive stars in 30 Dor. Mass loss, envelope inflation, and a variable upper stellar mass limit
• Convectively Driven Three Dimensional Turbulence in Massive Star Envelopes: I. A 1D Implementation of Diffusive Radiative Transport
• Superadiabaticity and the metallicity independence of the Humphreys-Davidson limit
• Neutrino cooling rates due to nickel isotopes for presupernova evolution of massive stars

Cited by in corpus (10)

• The boring history of Gaia BH3 from isolated binary evolution
• Three Dimensional Natures of Massive Star Envelopes
• Modelling stellar evolution in mass-transferring binaries and gravitational-wave progenitors with METISSE
• Gaia BH1 and BH2 -- evolutionary models with overshooting of the black hole progenitors within the present-day binary separation
• Evolution of rotating massive stars adopting a newer, self-consistent wind prescription at SMC metallicity
• Collision-induced mass loss and mass gain on an extremely massive star. An analytical approach and a static proto-globular cluster test-case
• Luminous Blue Variables
• Development of convective envelopes in massive stars: Implications for gravitational wave sources
• Hydrogen-free Wolf-Rayet stars: Helium stars with envelope-inflation structure and rotation
• Evolution of stars with 60 and 200 Msun: predictions for WNh stars in the Milky Way