The Humphreys-Davidson (HD) limit empirically defines a region of high luminosities (log L > 5.5) and low effective temperatures (T < 20kK) on the Hertzsprung-Russell Diagram in which hardly any supergiant stars are observed. Attempts to explain this
limit through instabilities arising in near- or super-Eddington winds have been largely unsuccessful. Using modern stellar evolution we aim to re-examine the HD limit, investigating the impact of enhanced mixing on massive stars. We construct grids of stellar evolution models appropriate for the Small and Large Magellanic Clouds (SMC, LMC), as well as for the Galaxy, spanning various initial rotation rates and convective overshooting parameters. Significantly enhanced mixing apparently steers stellar evolution tracks away from the region of the HD limit. To quantify the excess of over-luminous stars in stellar evolution simulations we generate synthetic populations of massive stars, and make detailed comparisons with catalogues of cool (T < 12.5kK) and luminous (log L > 4.7) stars in the SMC and LMC. We find that adjustments to the mixing parameters can lead to agreement between the observed and simulated red supergiant populations, but for hotter supergiants the simulations always over-predict the number of very luminous (log L > 5.4) stars compared to observations. The excess of luminous supergiants decreases for enhanced mixing, possibly hinting at an important role mixing has in explaining the HD limit. Still, the HD limit remains unexplained for hotter supergiants.
We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we us
e an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds.
We find significant fluctuations of angular momentum within the convective helium shell of a pre-collapse massive star - a core-collapse supernova progenitor - which may facilitate the formation of accretion disks and jets that can explode the star.
The convective flow in our model of an evolved M_ZAMS=15Msun star, computed with the sub-sonic hydrodynamic solver MAESTRO, contains entire shells with net angular momentum in different directions. This phenomenon may have important implications for the late evolutionary stages of massive stars, and for the dynamics of core-collapse.
