Zero-age Main-sequence Mass Research Articles

Enhanced mass-loss rate evolution of stars with ≳18 M⊙ and missing optically observed type II core-collapse supernovae

ABSTRACT We evolve stellar models with zero-age main-sequence (ZAMS) mass of MZAMS ≳ 18 M⊙ under the assumption that they experience an enhanced mass-loss rate when crossing the instability strip at high luminosities and conclude that most of them end as type Ibc supernovae (SNe Ibc) or dust-obscured SNe II. We explore what level of enhanced mass-loss rate during the instability strip would be necessary to explain the ‘red supergiant problem’. This problem refers to the dearth of observed core-collapse supernovae progenitors with MZAMS ≳ 18 M⊙. Namely, we examine what enhanced mass-loss rate could make it possible for all these stars actually to explode as core-collapse supernovae (CCSNe). We find that the mass-loss rate should increase by a factor of at least about 10. We reach this conclusion by analysing the hydrogen mass in the stellar envelope and the optical depth of the dusty wind at the explosion, and crudely estimate that under our assumptions only about a fifth of these stars explode as unobscured SNe II and SNe IIb. About 10–15 per cent end as obscured SNe II that are infrared-bright but visibly very faint, and the rest, about 65–70 per cent, end as SNe Ibc. However, the statistical uncertainties are still too significant to decide whether many stars with MZAMS ≳ 18 M⊙ do not explode as expected in the neutrino driven explosion mechanism, or whether all of them explode as CCSNe, as expected by the jittering jets explosion mechanism.

Magnetorotational core collapse of possible GRB progenitors – I. Explosion mechanisms

ABSTRACT We investigate the explosion of stars with zero-age main-sequence masses between 20 and 35 M⊙ and varying degrees of rotation and magnetic fields including ones commonly considered progenitors of gamma-ray bursts (GRBs). The simulations, combining special relativistic magnetohydrodynamics, a general relativistic approximate gravitational potential, and two-moment neutrino transport, demonstrate the viability of different scenarios for the post-bounce evolution. Having formed a highly massive proto-neutron star (PNS), several models launch successful explosions, either by the standard supernova mechanism based on neutrino heating and hydrodynamic instabilities or by magnetorotational processes. It is, however, quite common for the PNS to collapse to a black hole (BH) within a few seconds. Others might produce proto-magnetar-driven explosions. We explore several ways to describe the different explosion mechanisms. The competition between the time-scales for advection of gas through the gain layer and heating by neutrinos provides an approximate explanation for models with insignificant magnetic fields. The fidelity of this explosion criterion in the case of rapid rotation can be improved by accounting for the strong deviations from spherical symmetry and mixing between pole and equator. We furthermore study an alternative description including the ram pressure of the gas falling through the shock. Magnetically driven explosions tend to arise from a strongly magnetized region around the polar axis. In these cases, the onset of the explosion corresponds to the equality between the advection time-scale and the time-scale for the propagation of Alfvén waves through the gain layer.

Impact of the Rotation and Compactness of Progenitors on the Mass of Black Holes

Abstract We investigate the impact of stellar rotation on the formation of black holes (BHs) by means of our population synthesis code sevn. Rotation affects the mass function of BHs in several ways. In massive metal-poor stars, fast rotation reduces the minimum zero-age main sequence (ZAMS) mass for a star to undergo pair instability and pulsational pair instability. Moreover, stellar winds are enhanced by rotation, peeling off the entire hydrogen envelope. As a consequence of these two effects, the maximum BH mass we expect from the collapse of a rotating metal-poor star is only ∼45 M ⊙, while the maximum mass of a BH born from a nonrotating star is ∼60 M ⊙. Furthermore, stellar rotation reduces the minimum ZAMS mass for a star to collapse into a BH from ∼18–25 M ⊙ to ∼13–18 M ⊙. Finally, we have investigated the impact of different core-collapse supernova (CCSN) prescriptions on our results. While the threshold value of compactness for direct collapse and the fallback efficiency strongly affect the minimum ZAMS mass for a star to collapse into a BH, the fraction of the hydrogen envelope that can be accreted onto the final BH is the most important ingredient in determining the maximum BH mass. Our results confirm that the interplay between stellar rotation, CCSNe and pair instability plays a major role in shaping the BH mass spectrum.

The 50–100 pc scale parent stellar populations of Type II supernovae and limitations of single star evolution models

ABSTRACT There is observational evidence of a dearth in core-collapse supernova (ccSN) explosions from stars with zero-age main-sequence (ZAMS) mass M0 ≈ 17–30M⊙, referred to as the ‘red supergiant problem’. However, simulations now predict that above 20 M⊙ we should indeed only expect stars within certain pockets of M0 to produce a visible SN explosion. Validating these predictions requires large numbers of ccSNe of different types with measured M0, which is challenging. In this paper, we explore the reliability of using host galaxy emission lines and the H α equivalent width to constrain the age, and thus the M0 of ccSNe progenitors. We use Binary Population and Spectral Synthesis models to infer a stellar population age from MUSE observations of the ionized gas properties and H α EW at the location of eleven ccSNe with reliable M0 measurements. Comparing our results to published M0 values, we find that models that do not consider binary systems yield stellar ages that are systematically too young (thus M0 too large), whereas accounting for binary system interactions typically overpredict the stellar age (thus underpredict M0). Taking into account the effects of photon leakage bring our M0 estimates in much closer agreement with expectations. These results highlight the need for careful modelling of diffuse environments, such as are present in the vicinity of Type II SNe, before ionized emission line spectra can be used as reliable tracers of progenitor stellar age.

Excavating the Explosion and Progenitor Properties of Type IIP Supernovae via Modeling of their Optical Light Curves.

The progenitors of Type IIP supernovae (SNe) are known to be red supergiants, but their properties are not well determined. We employ hydrodynamical modeling to investigate the explosion characteristics of eight Type IIP SNe and the properties of their progenitor stars. We create evolutionary models using the mesa stellar evolution code, explode these models, and simulate the optical light curves using the stella code. We fit the optical light curves, Fe ii 5169 Å velocity, and photospheric velocity to the observational data. Recent research has suggested that the progenitors of Type IIP SNe have a zero-age main-sequence (ZAMS) mass not exceeding ~18 M ⊙. Our fits give a progenitor ZAMS mass of ⩽18 M ⊙ for seven of the SNe. Where previous progenitor mass estimates exist from various sources, such as hydrodynamical modeling, multiwavelength observations, or semi-analytic calculations, our modeling generally tends toward the lower-mass values. This result is in contrast to results from previous hydrodynamical modeling but consistent with those obtained using general-relativistic radiation-hydrodynamical codes. We do find that one event, SN 2015ba, has a progenitor whose mass is closer to 24 M ⊙, although we are unable to fit it well. We also derive the amount of 56Ni required to reproduce the tail of the light curve and find values generally larger than previous estimates. Overall, we find that it is difficult to characterize the explosion by a single parameter, and that a range of parameters is needed.

ZTF18aalrxas: A Type IIb Supernova from a Very Extended Low-mass Progenitor

Abstract We investigate ZTF18aalrxas, a double-peaked Type IIb core-collapse supernova (SN) discovered during science validation of the Zwicky Transient Facility. ZTF18aalrxas was discovered while the optical emission was still rising toward the initial cooling peak (0.7 mag over 2 days). Our observations consist of multi-band (ultraviolet and optical) light curves (LCs), and optical spectra spanning from ≈0.7 to ≈180 days past the explosion. We use a Monte-Carlo based non-local thermodynamic equilibrium model that simultaneously reproduces both the 56Ni-powered bolometric LC and our nebular spectrum. This model is used to constrain the synthesized radioactive nickel mass (0.17 M ☉) and the total ejecta mass (1.7 M ☉) of the SN. The cooling emission is modeled using semi-analytical extended envelope models to constrain the progenitor radius (790–1050 R ⊙) at the time of explosion. Our nebular spectrum shows signs of interaction with a dense circumstellar medium (CSM), and this spectrum is modeled and analyzed to constrain the amount of ejected oxygen (0.3–0.5 M ☉) and the total hydrogen mass (≈0.15 M ☉) in the envelope of the progenitor. The oxygen mass of ZTF18aalrxas is consistent with a low (12–13 M ☉) zero-age main-sequence mass progenitor. The LCs and spectra of ZTF18aalrxas are not consistent with massive single-star SN Type IIb progenitor models. The presence of an extended hydrogen envelope of low mass, the presence of a dense CSM, the derived ejecta mass, and the late-time oxygen emission can all be explained in a binary model scenario.

The difficulty of inferring progenitor masses from type-II-Plateau supernova light curves

Much controversy surrounds the inferred progenitor masses of type-II-Plateau (II-P) supernovae (SNe). The debate is nourished by the discrepant results from radiation-hydrodynamics simulations, pre-explosion imaging, and studies of host stellar populations. Here, we present a controlled experiment using four solar-metallicity models with zero-age main sequence masses of 12, 15, 20, and 25 M⊙. Because of the effects of core burning and surface mass loss, these models reach core collapse as red-supergiant (RSG) stars with a similar H-rich envelope mass of 8 to 9 M⊙ but with final masses in the range 11 to 16 M⊙. We explode the progenitors using a thermal bomb, adjusting the energy deposition to yield an asymptotic ejecta kinetic energy of 1.25 × 1051 erg and an initial 56Ni mass of 0.04 M⊙. The resulting SNe produce similar photometric and spectroscopic properties from 10 to 200 d. The spectral characteristics are degenerate. The scatter in early-time color results from the range in progenitor radii, while the differences in late-time spectra reflect the larger oxygen yields in more massive progenitors. Because the progenitors have a comparable H-rich envelope mass, the photospheric phase duration is comparable for all models; the difference in He-core mass is invisible. As different main sequence masses can produce progenitors with a similar H-rich envelope mass, light-curve modeling cannot provide a robust and unique solution for the ejecta mass of type-II-P SNe. The numerous uncertainties in massive-star evolution and wind-mass loss also prevent a robust association with a main sequence star mass. Light-curve modeling can at best propose compatibility.

A hybrid envelope-stripping mechanism for massive stars from supernova nebular spectroscopy

When nuclear fuel in the core of a massive star with a zero-age main-sequence mass $M_{\rm ZAMS} \gtrsim 8M_\odot$ is exhausted, the central part of the iron or magnesium core collapses and forms a neutron star or a black hole. At the same time, the material above the collapsing core is rapidly ejected, leading to a stripped-envelope supernova (SESN) explosion if the outer hydrogen envelope of the star was removed before its explosion. The envelope is presumably stripped either via strong stellar winds or due to mass transfer to a companion star in a close binary orbit. It is not clear which process is dominant, and whether different mechanisms are at work for different classes of SESNe; type IIb, Ib, and Ic SNe in order of increasing degree of envelope stripping. In this work, a new analysis of late-time nebular spectra of SESNe is presented, which is more sensitive to differences in the core structure than early-phase spectral analysis. The results show that the progenitors of SNe IIb and Ib are indistinguishable except for the residual amount of hydrogen envelope while the progenitors of SNe Ic are not only deficient in hydrogen and helium, but are also distinctly more massive than SNe IIb and Ib. These findings strongly suggest that more than one mechanism is responsible for the removal of the outer hydrogen envelope and the deeper helium layer, with the former most likely due to binary interaction, and the latter involving a mass-dependent process such as strong stellar winds or episodic pre-explosion mass ejection.

A progenitor candidate for the type II-P supernova SN 2018aoq in NGC 4151

We present our findings based on pre- and post-explosion data of the type II-Plateau SN 2018aoq that exploded in NGC 4151. As distance estimates to NGC 4151 vary by an order of magnitude, we utilised the well-known correlation between ejecta velocity and plateau brightness, i.e. the standard candle method, to obtain a distance of 18.2 ± 1.2 Mpc, which is in very good agreement with measurements based on geometric methods. The above distance implies a mid-plateau absolute magnitude of MV50 = − 15.76 ± 0.14 suggesting that it is of intermediate brightness when compared to IIP SNe such as SN 2005cs at the faint end, and more typical events such as SN 1999em. This is further supported by relatively low expansion velocities (Fe IIλ5169 ∼ 3000 km s−1 at +42 d). Using archival HST/WFC3 imaging data, we find a point source coincident with the supernova position in the F350LP, F555W, F814W, and F160W filters. This source shows no significant variability over the ∼2 month time span of the data. From fits to the spectral energy distribution of the candidate progenitor, we find log(L/L⊙) ∼ 4.7 and Teff ∼ 3.5 kK, implying an M-type red supergiant progenitor. From comparisons to single and binary star models, we find that both favour the explosion of a star with a zero-age main sequence mass of ∼10 M⊙.

Three-dimensional simulations of neutrino-driven core-collapse supernovae from low-mass single and binary star progenitors

We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.

Quark deconfinement as a supernova explosion engine for massive blue supergiant stars

Blue-supergiant stars develop into core-collapse supernovae --- one of the most energetic outbursts in the universe --- when all nuclear burning fuel is exhausted in the stellar core. Previous attempts failed to explain observed explosions of such stars which have a zero-age main sequence mass of 50~M$_\odot$ or more. Here we exploit the largely uncertain state of matter at high density, and connect the modeling of such stellar explosions with a first-order phase transition from nuclear matter to the quark-gluon plasma. The resulting energetic supernova explosions can account for a large variety of lightcurves, from peculiar type II to super-luminous events. The remnants are neutron stars with quark matter core, known as hybrid stars, of about 2~M$_\odot$ at birth. A galactic event of this kind could be observable due to the release of a second neutrino burst. Its observation would confirm such a first-order phase transition at densities relevant for astrophysics.

Progenitor Mass Distribution of Core-collapse Supernova Remnants in Our Galaxy and Magellanic Clouds Based on Elemental Abundances

Abstract We investigate a progenitor mass distribution of core-collapse supernova remnants (CCSNRs) in our Galaxy and Large and Small Magellanic Clouds for the first time. We count the number of CCSNRs in three mass ranges divided by the zero-age main-sequence (ZAMS) mass, M ZAMS; A: M ZAMS < 15 M ⊙, B: 15 M ⊙ < M ZAMS < 22.5 M ⊙, C: M ZAMS > 22.5 M ⊙. A simple compilation of progenitor masses in the literature yields a progenitor mass distribution of f A : f B : f C = 0.27:0.27:0.46, where f is the number fraction of the progenitors. The distribution is inconsistent with any standard initial mass functions (IMFs). We notice, however, that previous mass estimates are subject to large systematic uncertainties because most of the relative abundances (X/Si) are not good probes for the progenitor masses. Instead, we rely only on the Fe/Si ratio, which is sensitive to the CO core mass (M COcore) and M ZAMS. Comparing Fe/Si ratios in SN remnants in the literature with the newest theoretical model, we estimate 33 M COcore and M ZAMS, leading to a revised progenitor mass distribution of f A : f B : f C = 0.47:0.32:0.21. This is consistent with the standard Salpeter IMF. However, the relation between M COcore and M ZAMS could be affected by binary evolution, which is not taken into account in this study and should be considered in future work to derive a better progenitor mass distribution estimate.

Bayesian hierarchical modelling of initial–final mass relations acrossstar clusters

The initial-final mass relation (IFMR) of white dwarfs (WDs) plays an important role in stellar evolution. To derive precise estimates of IFMRs and explore how they may vary among star clusters, we propose a Bayesian hierarchical model that pools photo- metric data from multiple star clusters. After performing a simulation study to show the benefits of the Bayesian hierarchical model, we apply this model to five star clus- ters: the Hyades, M67, NGC 188, NGC 2168, and NGC 2477, leading to reasonable and consistent estimates of IFMRs for these clusters. We illustrate how a cluster-specific analysis of NGC 188 using its own photometric data can produce an unreasonable IFMR since its WDs have a narrow range of zero-age main sequence (ZAMS) masses. However, the Bayesian hierarchical model corrects the cluster-specific analysis by bor- rowing strength from other clusters, thus generating more reliable estimates of IFMR parameters. The data analysis presents the benefits of Bayesian hierarchical modelling over conventional cluster-specific methods, which motivates us to elaborate the pow- erful statistical techniques in this article.

SN 2015ba: a Type IIP supernova with a long plateau

We present optical photometry and spectroscopy from about a week after explosion to $\sim$272 d of an atypical Type IIP supernova, SN 2015ba, which exploded in the edge-on galaxy IC 1029. SN 2015ba is a luminous event with an absolute V-band magnitude of -17.1$\pm$0.2 mag at 50 d since explosion and has a long plateau lasting for $\sim$123 d. The distance to the SN is estimated to be 34.8$\pm$0.7 Mpc using the expanding photosphere and standard candle methods. High-velocity H-Balmer components constant with time are observed in the late-plateau phase spectra of SN 2015ba, which suggests a possible role of circumstellar interaction at these phases. Both hydrodynamical and analytical modelling suggest a massive progenitor of SN 2015ba with a pre-explosion mass of 24-26 M$_\odot$. However, the nebular spectra of SN 2015ba exhibit insignificant levels of oxygen, which is otherwise expected from a massive progenitor. This might be suggestive of the non-monotonical link between O-core masses and the zero-age main-sequence mass of pre-supernova stars and/or uncertainties in the mixing scenario in the ejecta of supernovae.

The lowest-metallicity type II supernova from the highest-mass red supergiant progenitor

Red supergiants have been confirmed as the progenitor stars of the majority of hydrogen-rich type II supernovae1. However, while such stars are observed with masses >25 M⊙ (ref. 2), detections of >18 M⊙ progenitors remain elusive1. Red supergiants are also expected to form at all metallicities, but discoveries of explosions from low-metallicity progenitors are scarce. Here, we report observations of the type II supernova, SN 2015bs, for which we infer a progenitor metallicity of ≤0.1 Z⊙ from comparison to photospheric-phase spectral models3, and a zero-age main-sequence mass of 17–25 M⊙ through comparison to nebular-phase spectral models4,5. SN 2015bs displays a normal ‘plateau’ light-curve morphology, and typical spectral properties, implying a red supergiant progenitor. This is the first example of such a high-mass progenitor for a ‘normal’ type II supernova, suggesting a link between high-mass red supergiant explosions and low-metallicity progenitors.

Measuring the Progenitor Masses and Dense Circumstellar Material of Type II Supernovae

Abstract Recent modeling of hydrogen-rich Type II supernova (SN II) light curves suggests the presence of dense circumstellar material (CSM) surrounding the exploding progenitor stars. This has important implications for the activity and structure of massive stars near the end of their lives. Since previous work focused on just a few events, here we expand to a larger sample of 20 well-observed SNe II. For each event we are able to constrain the progenitor zero-age main-sequence (ZAMS) mass, explosion energy, and the mass and radial extent of the dense CSM. We then study the distribution of each of these properties across the full sample of SNe. The inferred ZAMS masses are found to be largely consistent with a Salpeter distribution with minimum and maximum masses of 10.4 and 22.9 M ⊙, respectively. We also compare the individual ZAMS masses we measure with specific SNe II that have pre-explosion imaging to check their consistency. Our masses are generally comparable to or higher than the pre-explosion imaging masses, potentially helping ease the red supergiant problem. The explosion energies vary from (0.1–1.3) × 1051 erg, and for ∼70% of the SNe we obtain CSM masses in the range between 0.18 and 0.83 M ⊙. We see a potential correlation between the CSM mass and explosion energy, which suggests that pre-explosion activity has a strong impact on the structure of the star. This may be important to take into account in future studies of the ability of the neutrino mechanism to explode stars. We also see a possible correlation between the CSM radial extent and ZAMS mass, which could be related to the time with respect to explosion when the CSM is first generated.

Subphotospheric fluctuations in magnetized radiative envelopes: contribution from unstable magnetosonic waves

We examine the excitation of unstable magnetosonic waves in the radiative envelopes of intermediate- and high-mass stars with a magnetic field of ~kG strength. Wind clumping close to the star and microturbulence can often be accounted for when including small-scale, sub-photospheric density or velocity perturbations. Compressional waves - with wavelengths comparable to or shorter than the gas pressure scale height - can be destabilized by the radiative flux in optically-thick media when a magnetic field is present, in a process called the Radiation-Driven Magneto-Acoustic Instability (RMI). The instability does not require radiation or magnetic pressure to dominate over gas pressure, and acts independently of sub-surface convection zones. Here we evaluate the conditions for the RMI to operate on a grid of stellar models covering a mass range $3-40M_\odot$ at solar metallicity. For a uniform 1kG magnetic field, fast magnetosonic modes are unstable down to an optical depth of a few tens, while unstable slow modes extend beyond the depth of the iron convection zone. The qualitative behavior is robust to magnetic field strength variations by a factor of a few. When combining our findings with previous results for the saturation amplitude of the RMI, we predict velocity fluctuations in the range ~0.1-10 km/s. These amplitudes are a monotonically increasing function of the ratio of radiation to gas pressure, or alternatively, of the zero-age main sequence mass.

The Progenitor Dependence of Core-collapse Supernovae from Three-dimensional Simulations with Progenitor Models of 12–40 M⊙

Abstract We present a first study of the progenitor star dependence of the three-dimensional (3D) neutrino mechanism of core-collapse supernovae. We employ full 3D general-relativistic multi-group neutrino radiation-hydrodynamics and simulate the postbounce evolutions of progenitors with zero-age main sequence masses of 12, 15, 20, 27, and 40 M ⊙. All progenitors, with the exception of the 12 M ⊙ star, experience shock runaway by the end of their simulations. In most cases, a strongly asymmetric explosion will result. We find three qualitatively distinct evolutions that suggest a complex dependence of explosion dynamics on progenitor density structure, neutrino heating, and 3D flow. (1) Progenitors with massive cores, shallow density profiles, and high post-core-bounce accretion rates experience very strong neutrino heating and neutrino-driven turbulent convection, leading to early shock runaway. Accretion continues at a high rate, likely leading to black hole formation. (2) Intermediate progenitors experience neutrino-driven, turbulence-aided explosions triggered by the arrival of density discontinuities at the shock. These occur typically at the silicon/silicon–oxygen shell boundary. (3) Progenitors with small cores and density profiles without strong discontinuities experience shock recession and develop the 3D standing-accretion shock instability (SASI). Shock runaway ensues late, once declining accretion rate, SASI, and neutrino-driven convection create favorable conditions. These differences in explosion times and dynamics result in a non-monotonic relationship between progenitor and compact remnant mass.

Constraining the Final Fates of Massive Stars by Oxygen and Iron Enrichment History in the Galaxy

Abstract Recent observational studies of core-collapse supernovae suggest that only stars with zero-age main-sequence masses smaller than 16– explode when they are red supergiants, producing Type IIP supernovae. This may imply that more massive stars produce other types of supernovae or they simply collapse to black holes without giving rise to bright supernovae. This failed supernova hypothesis can lead to significantly inefficient oxygen production because oxygen abundantly produced in inner layers of massive stars with zero-age main-sequence masses around 20– might not be ejected into the surrounding interstellar space. We first assume an unspecified population of oxygen injection events related to massive stars and obtain a model-independent constraint on how much oxygen should be released in a single event and how frequently such events should happen. We further carry out one-box galactic chemical enrichment calculations with different mass ranges of massive stars exploding as core-collapse supernovae. Our results suggest that the model assuming that all massive stars with 9– explode as core-collapse supernovae is still most appropriate in explaining the solar abundances of oxygen and iron and their enrichment history in the Galaxy. The oxygen mass in the Galaxy is not explained when assuming that only massive stars with zero-age main-sequence masses in the range of 9–17 contribute to the galactic oxygen enrichment. This finding implies that a good fraction of stars more massive than should eject their oxygen layers in either supernova explosions or some other mass-loss processes.

Systematic survey of the effects of wind mass loss algorithms on the evolution of single massive stars

Mass loss is a key uncertainty in the evolution of massive stars. Stellar evolution calculations must employ parametric algorithms for mass loss, and usually only include stellar winds. We carry out a parameter study of the effects of wind mass loss on massive star evolution using the open-source stellar evolution code MESA. We provide a systematic comparison of wind mass loss algorithms for solar-metallicity, nonrotating, single stars in the initial mass range of $15-35\,M_\odot$. We consider combinations drawn from two hot phase algorithms, three cool phase algorithms, and two Wolf-Rayet algorithms. We consider linear wind efficiency scale factors of $1$, $0.33$, and $0.1$ to account for reductions in mass loss rates due to wind inhomogeneities. We find that the initial to final mass mapping for each zero-age main-sequence (ZAMS) mass has a $\sim 50\%$ uncertainty if all algorithm combinations and wind efficiencies are considered. The ad-hoc efficiency scale factor dominates this uncertainty. While the final total mass and internal structure of our models vary tremendously with mass loss treatment, final observable parameters are much less sensitive for ZAMS mass $\lesssim 30\,M_\odot$. This indicates that uncertainty in wind mass loss does not negatively affect estimates of the ZAMS mass of most single-star supernova progenitors from pre-explosion observations. Furthermore, we show that the internal structure of presupernova stars is sensitive to variations in both main sequence and post main-sequence mass loss. We find that the compactness parameter $\xi\propto M/R(M)$ varies by as much as $30\%$ for a given ZAMS mass evolved with different wind efficiencies and mass loss algorithm combinations. [abridged]