The metallicity dependence of envelope inflation in massive stars

Recently it has been found that models of massive stars reach the Eddington limit in their interior, which leads to dilute extended envelopes. We perform a comparative study of the envelope properties of massive stars at different metallicities, with the aim to establish the impact of the stellar metallicity on the effect of envelope inflation. We analyse published grids of core-hydrogen burning massive star models computed with metallicities appropriate for massive stars in the Milky Way, the LMC and the SMC, the very metal poor dwarf galaxy I Zwicky 18, and for metal-free chemical composition. Stellar models of all the investigated metallicities reach and exceed the Eddington limit in their interior, aided by the opacity peaks of iron, helium and hydrogen, and consequently develop inflated envelopes. Envelope inflation leads to a redward bending of the zero-age main sequence and a broadening of the main sequence band in the upper part of the Hertzsprung-Russell diagram. We derive the limiting L/M-values as function of the stellar surface temperature above which inflation occurs, and find them to be larger for lower metallicity. While Galactic models show inflation above ~29 Msun, the corresponding mass limit for Population III stars is ~150 Msun. While the masses of the inflated envelopes are generally small, we find that they can reach 1-100 Msun in models with effective temperatures below ~8000 K, with higher masses reached by models of lower metallicity. Envelope inflation is expected to occur in sufficiently massive stars at all metallicities, and is expected to lead to rapidly growing pulsations, high macroturbulent velocities, and might well be related to the unexplained variability observed in Luminous Blue Variables like S Doradus and Eta Carina.Comment: 16 pages (with Appendix), accepted in A&

Asteroseismic test of rotational mixing in low-mass white dwarfs

We exploit the recent discovery of pulsations in mixed-atmosphere (He/H), extremely low-mass white dwarf precursors (ELM proto-WDs) to test the proposition that rotational mixing is a fundamental process in the formation and evolution of low-mass helium core white dwarfs. Rotational mixing has been shown to be a mechanism able to compete efficiently against gravitational settling, thus accounting naturally for the presence of He, as well as traces of metals such as Mg and Ca, typically found in the atmospheres of ELM proto-WDs. Here we investigate whether rotational mixing can maintain a sufficient amount of He in the deeper driving region of the star, such that it can fuel, through HeII-HeIII ionization, the observed pulsations in this type of stars. Using state-of-the-art evolutionary models computed with MESA, we show that rotational mixing can indeed explain qualitatively the very existence and general properties of the known pulsating, mixed-atmosphere ELM proto-WDs. Moreover, such objects are very likely to pulsate again during their final WD cooling phase.Comment: accepted for publication in A&A Letter

Massive main-sequence stars evolving at the Eddington limit

Context. Massive stars play a vital role in the Universe, however, their evolution even on the main-sequence is not yet well understood. Aims. Because of the steep mass-luminosity relation, massive main-sequence stars become extremely luminous. This brings their envelopes very close to the Eddington limit. We analyse stellar evolutionary models in which the Eddington limit is reached and exceeded, explore the rich diversity of physical phenomena that take place in their envelopes, and investigate their observational consequences. Methods. We use published grids of detailed stellar models, computed with a state-of-the-art, one-dimensional hydrodynamic stellar evolution code using LMC composition, to investigate the envelope properties of core hydrogen burning massive stars. Results. We find that the Eddington limit is almost never reached at the stellar surface, even for stars up to 500 M⊙. When we define an appropriate Eddington limit locally in the stellar envelope, we can show that most stars more massive than ~40 M⊙ actually exceed this limit, in particular, in the partial ionisation zones of iron, helium, or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions, and, possibly to, pulsations. We find that all models with luminosities higher than 4 × 105L⊙, i.e. stars above ~40 M⊙ show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates for producing major luminous blue variable eruptions. Conclusions. Our models show that the Eddington limit is expected to be reached in all stars above ~40 M⊙ in the LMC, even in lower mass stars in the Galaxy, or in close binaries or rapid rotators. While our results do not support the idea of a direct super-Eddington wind driven by continuum photons, the consequences of the Eddington limit in the form of inflation, pulsations and possibly eruptions may well give rise to a significant enhancement of the time averaged mass-loss rate

X-Shooting ULLYSES: Massive stars at low metallicity: I. Project description

Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational wave events involving spectacular black hole mergers indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observing ∼250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES programme. The complementary X-Shooting ULLYSES (XShootU) project provides an enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESOa's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates as a function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of astrophysics, the data and modelling of the XShootU project is expected to be a game changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope spectra of the first stellar generations, the individual spectra of low-Z stars need to be understood, which is exactly where XShootU can deliver

Properties of OB star−black hole systems derived from detailed binary evolution models

Context. The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. Aims. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples, such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. With this purpose in mind, we analysed a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 M⊙, and identified the model systems that potentially evolve into a binary consisting of a black hole and a massive main-sequence star. We then derived the observable properties of such systems, as well as peculiarities of the OB star component. Results. We find that ∼3% of the LMC late-O and early-B stars in binaries are expected to possess a black hole companion when stars with a final helium core mass above 6.6 M⊙ are assumed to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these black holes may be identified in spectroscopic binaries, either by large amplitude radial velocity variations (≳50 km s−1) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity (≳10 km s−1) and simultaneous rapid rotation of the OB star. The predicted mass ratios are such that main-sequence companions can be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be and X-ray binaries, and known massive black hole binaries supports our conclusion. Conclusions. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general and for the formation of double black hole mergers in particular

On the optically thick winds of Wolf-Rayet stars

Context. The classical Wolf-Rayet (WR) phase is believed to mark the end stage of the evolution of massive stars with initial masses higher than ~25M⊙. Stars in this phase expose their stripped cores with the products of H- or He-burning at their surface. They develop strong, optically thick stellar winds that are important for the mechanical and chemical feedback of massive stars, and that determine whether the most massive stars end their lives as neutron stars or black holes. The winds of WR stars are currently not well understood, and their inclusion in stellar evolution models relies on uncertain empirical mass-loss relations. Aims. We investigate theoretically the mass-loss properties of H-free WR stars of the nitrogen sequence (WN stars). Methods. We connected stellar structure models for He stars with wind models for optically thick winds and assessed the degree to which these two types of models can simultaneously fulfil their respective sonic-point conditions. Results. Fixing the outer wind law and terminal wind velocity ν∞, we obtain unique solutions for the mass-loss rates of optically thick, radiation-driven winds of WR stars in the phase of core He-burning. The resulting mass-loss relations as a function of stellar parameters agree well with previous empirical relations. Furthermore, we encounter stellar mass limits below which no continuous solutions exist. While these mass limits agree with observations of WR stars in the Galaxy, they contradict observations in the LMC. Conclusions. While our results in particular confirm the slope of often-used empirical mass-loss relations, they imply that only part of the observed WN population can be understood in the framework of the standard assumptions of a smooth transonic flow and compact stellar core. This means that alternative approaches such as a clumped and inflated wind structure or deviations from the diffusion limit at the sonic point may have to be invoked. Qualitatively, the existence of mass limits for the formation of WR-type winds may be relevant for the non-detection of low-mass WR stars in binary systems, which are believed to be progenitors of Type Ib/c supernovae. The sonic-point conditions derived in this work may provide a possibility to include optically thick winds in stellar evolution models in a more physically motivated form than in current models

Relating turbulent pressure and macroturbulence across the HR diagram with a possible link to γ Doradus stars

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational γ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the γ Dor phenomenon, calling for more detailed theoretical modeling in this direction

Relating turbulent pressure and macroturbulence across the HR diagram with a possible link to γ Doradus stars

A significant fraction of the envelope of low- and intermediate-mass stars is unstable to convection, leading to sub-surface turbulent motion. Here, we consider and include the effects of turbulence pressure in our stellar evolution calculations. In search of an observational signature, we compare the fractional contribution of turbulent pressure to the observed macroturbulent velocities in stars at different evolutionary stages. We find a strong correlation between the two quantities, similar to what was previously found for massive OB stars. We therefore argue that turbulent pressure fluctuations of finite amplitude may excite high-order, high-angular degree stellar oscillations, which manifest themselves at the surface an additional broadening of the spectral lines, i.e., macroturbulence, across most of the HR diagram. When considering the locations in the HR diagram where we expect high-order oscillations to be excited by stochastic turbulent pressure fluctuations, we find a close match with the observational γ Doradus instability strip, which indeed contains high-order, non-radial pulsators. We suggest that turbulent pressure fluctuations on a percentual level may contribute to the γ Dor phenomenon, calling for more detailed theoretical modeling in this direction