On the Red Supergiant Wind Kink. A Universal mass-loss concept for massive stars

Red supergiants (RSG) are key objects for the evolution of massive stars and their endpoints, but uncertainties in their underlying mass-loss mechanism have thus far prevented an appropriate framework for massive star evolution. We analyse an empirical mass loss"kink" feature uncovered by Yang et al., and we highlight its similarity to hot star radiation-driven wind models and observations at the optically thin/thick transition point. We motivate a new RSG mass-loss prescription that depends on the Eddington factor Gamma (including both a steep L dependence and an inverse steep M dependence). We subsequently implement this new RSG mass-loss prescription in the stellar evolution code MESA. We find that our physically motivated mass-loss behaviour naturally reproduces the Humphreys-Davidson limit without a need for any ad-hoc tweaks. It also resolves the RSG supernova "problem". We argue that a universal behaviour of radiation-driven winds across the HR diagram, independent of the exact source of opacity, is a key feature of the evolution of the most massive stars.Comment: Accepted Letter in Astronomy & Astrophysics (A&A). 4 pages. 3 figure

The hydrogen clock to infer the upper stellar mass

The most massive stars dominate the chemical enrichment, mechanical and radiative feedback, and energy budget of their host environments. Yet how massive stars initially form and how they evolve throughout their lives is ambiguous. The mass loss of the most massive stars remains a key unknown in stellar physics, with consequences for stellar feedback and populations. In this work, we compare grids of very massive star (VMS) models with masses ranging from 80-1000Msun, for a range of input physics. We include enhanced winds close to the Eddington limit as a comparison to standard O-star winds, with consequences for present-day observations of ~50-100Msun stars. We probe the relevant surface H abundances (Xs) to determine the key traits of VMS evolution compared to O stars. We find fundamental differences in the behaviour of our models with the enhanced-wind prescription, with a convergence on the stellar mass at 1.6 Myr, regardless of the initial mass. It turns out that Xs is an important tool in deciphering the initial mass due to the chemically homogeneous nature of VMS above a mass threshold. We use Xs to break the degeneracy of the initial masses of both components of a detached binary, and a sample of WNh stars in the Tarantula nebula. We find that for some objects, the initial masses are unrestricted and, as such, even initial masses of the order 1000Msun are not excluded. Coupled with the mass turnover at 1.6 Myr, Xs can be used as a 'clock' to determine the upper stellar mass.Comment: Accepted for publication in MNRAS, 14 figure

Predicting the Heaviest Black Holes below the Pair Instability Gap

Traditionally, the pair instability (PI) mass gap is located between 50\,and 130\,$M_{\odot}$, with stellar mass black holes (BHs) expected to "pile up" towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an "impossibly" heavy BH of 85\,$M_{\odot}$ as part of GW\,190521 located inside the traditional PI gap, we realised that blue supergiant (BSG) progenitors with small cores but large Hydrogen envelopes at low metallicity ($Z$) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed MESA stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semi-convection, and $Z$. We evolve low $Z$ stars in the range $10^{-3} < Z / Z_{\odot} < Z_{\rm SMC}$, making no prior assumption regarding the mass of an envelope, but instead employing a wind mass loss recipe to calculate it. We compute critical Carbon-Oxygen and Helium core masses to determine our lower limit to PI physics, and we provide two equations for $M_{\text{core}}$ and $M_{\text{final}}$ that can also be of use for binary population synthesis. Assuming the H envelope falls into the BH, we confirm the maximum BH mass below PI is $M_{\text{BH}} \simeq 93.3$ $M_{\odot}$. Our grid allows us to populate the traditional PI gap, and we conclude that the distribution of BHs above the traditional boundary is not solely due to the shape of the initial mass function (IMF), but also to the same stellar interior physics (i.e. mixing) that which sets the BH maximum.Comment: 24 pages, 14 figures. Accepted in MNRA

Stellar Wind Yields of Very Massive Stars

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMS, M>100M) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50 - 500M, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne-Na and Mg-Al cycles. VMS with enhanced winds eject 5-10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anti-correlations, such as C-N and Na-O, in globular clusters. We find that for VMS 95% of the total wind yields is produced on the main sequence, while only ~5% is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf-Rayet winds had been suggested to be responsible for Galactic 26Al enrichment. Finally, 200M stars eject 100 times more of each heavy element in their winds than 50M stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50M stars.Comment: Accepted for publication in MNRAS. 14 pages, 10 figure

Stellar wind yields of very massive stars

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMSs, M>100) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50-500, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne-Na and Mg-Al cycles. VMS with enhanced winds eject 5-10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anticorrelations, such as C-N and Na-O, in globular clusters. We find that for VMS 95 per cent of the total wind yields is produced on the main sequence, while only ∼5 per cent is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf-Rayet winds had been suggested to be responsible for galactic 26Al enrichment. Finally, 200 stars eject 100 times more of each heavy element in their winds than 50 stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50 stars

Stellar wind yields of very massive stars

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMS, M&amp;gt;100 $\rm {\rm M}_{\odot }$) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50 – 500 $\rm {\rm M}_{\odot }$, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne-Na and Mg-Al cycles. VMS with enhanced winds eject 5-10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anti-correlations, such as C-N and Na-O, in globular clusters. We find that for VMS 95% of the total wind yields is produced on the main sequence, while only ∼ 5% is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf-Rayet winds had been suggested to be responsible for Galactic 26Al enrichment. Finally, 200 $\rm {\rm M}_{\odot }$ stars eject 100 times more of each heavy element in their winds than 50 $\rm {\rm M}_{\odot }$ stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50 $\rm {\rm M}_{\odot }$ stars

Predicting the heaviest black holes below the pair instability gap

Traditionally, the pair instability (PI) mass gap is located between 50 and 130 M⊙, with stellar mass black holes (BHs) expected to ”pile up” towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an ”impossibly” heavy BH of 85 M⊙ as part of GW 190521 located inside the traditional PI gap, we realised that blue supergiant (BSG) progenitors with small cores but large Hydrogen envelopes at low metallicity (Z) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed MESA stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semi-convection, and Z. We evolve low Z stars in the range 10−3 &amp;lt; Z/Z⊙ &amp;lt; ZSMC, making no prior assumption regarding the mass of an envelope, but instead employing a wind mass loss recipe to calculate it. We compute critical Carbon-Oxygen and Helium core masses to determine our lower limit to PI physics, and we provide two equations for Mcore and Mfinal that can also be of use for binary population synthesis. Assuming the H envelope falls into the BH, we confirm the maximum BH mass below PI is MBH ≃ 93.3 M⊙. Our grid allows us to populate the traditional PI gap, and we conclude that the distribution of BHs above the gap is not solely due to the shape of the initial mass function (IMF), but also to the same stellar interior physics (i.e. mixing) that which sets the BH maximum