Ylva

Götberg

Magnetic fields can drastically change predictions of evolutionary models of massive stars via mass-loss quenching, magnetic braking, and efficient angular momentum transport, which we aim to quantify in this work. We use the mesa software instrument to compute an extensive main-sequence grid of stellar structure and evolution models, as well as isochrones, accounting for the effects attributed to a surface fossil magnetic field. The grid is densely populated in initial mass (3-60 M-circle dot), surface equatorial magnetic field strength (0-50 kG), and metallicity (representative of the Solar neighbourhood and the Magellanic Clouds). We use two magnetic braking and two chemical mixing schemes and compare the model predictions for slowly rotating, nitrogen-enriched ('Group 2') stars with observations in the Large Magellanic Cloud. We quantify a range of initial field strengths that allow for producing Group 2 stars and find that typical values (up to a few kG) lead to solutions. Between the subgrids, we find notable departures in surface abundances and evolutionary paths. In our magnetic models, chemical mixing is always less efficient compared to non-magnetic models due to the rapid spin-down. We identify that quasi-chemically homogeneous main sequence evolution by efficient mixing could be prevented by fossil magnetic fields. We recommend comparing this grid of evolutionary models with spectropolarimetric and spectroscopic observations with the goals of (i) revisiting the derived stellar parameters of known magnetic stars, and (ii) observationally constraining the uncertain magnetic braking and chemical mixing schemes.

The majority of massive stars live in binary or multiple systems and will interact with a companion during their lifetimes, which helps to explain the observed diversity of core-collapse supernovae. Donor stars in binary systems can lose most of their hydrogen-rich envelopes through mass transfer. As a result, not only are the surface properties affected, but so is the core structure. However, most calculations of the core-collapse properties of massive stars rely on single-star models. We present a systematic study of the difference between the pre-supernova structures of single stars and stars of the same initial mass (11-21M(circle dot)) that have been stripped due to stable post-main-sequence mass transfer at solar metallicity. We present the pre-supernova core composition with novel diagrams that give an intuitive representation of the isotope distribution. As shown in previous studies, at the edge of the carbon-oxygen core, the binary-stripped star models contain an extended gradient of carbon, oxygen, and neon. This layer remains until core collapse and is more extended in mass for higher initial stellar masses. It originates from the receding of the convective helium core during core helium burning in binary-stripped stars, which does not occur in single-star models. We find that this same evolutionary phase leads to systematic differences in the final density and nuclear energy generation profiles. Binary-stripped star models have systematically higher total masses of carbon at the moment of core collapse compared to single-star models, which likely results in systematically different supernova yields. In about half of our models, the silicon-burning and oxygen-rich layers merge after core silicon burning. We discuss the implications of our findings for the "explodability", supernova observations, and nucleosynthesis of these stars. Our models are publicly available and can be readily used as input for detailed supernova simulations.

Helium star-carbon-oxygen white dwarf (CO WD) binaries are potential single-degenerate progenitor systems of thermonuclear supernovae. Revisiting a set of binary evolution calculations using the stellar evolution code MESA, we refine our previous predictions about which systems can lead to a thermonuclear supernova and then characterize the properties of the helium star donor at the time of explosion. We convert these model properties to near-UV/optical magnitudes assuming a blackbody spectrum and support this approach using a matched stellar atmosphere model. These models will be valuable to compare with pre-explosion imaging for future supernovae, though we emphasize the observational difficulty of detecting extremely blue companions. The pre-explosion source detected in association with SN 2012Z has been interpreted as a helium star binary containing an initially ultra-massive WD in a multiday orbit. However, extending our binary models to initial CO WD masses of up to 1.2 M (circle dot), we find that these systems undergo off-center carbon ignitions and thus are not expected to produce thermonuclear supernovae. This tension suggests that, if SN 2012Z is associated with a helium star-WD binary, then the pre-explosion optical light from the system must be significantly modified by the binary environment and/or the WD does not have a carbon-rich interior composition.

Most massive stars are born in binaries close enough for mass transfer episodes. These modify the appearance, structure, and future evolution of both stars. We compute the evolution of a 100-day-period binary, consisting initially of a 25 M (circle dot) star and a 17 M (circle dot) star, which experiences stable mass transfer. We focus on the impact of mass accretion on the surface composition, internal rotation, and structure of the accretor. To anchor our models, we show that our accretor broadly reproduces the properties of zeta Ophiuchi, which has long been proposed to have accreted mass before being ejected as a runaway star when the companion exploded. We compare our accretor to models of single rotating stars and find that the later and stronger spin-up provided by mass accretion produces significant differences. Specifically, the core of the accretor retains higher spin at the end of the main sequence, and a convective layer develops that changes its density profile. Moreover, the surface of the accretor star is polluted by CNO-processed material donated by the companion. Our models show effects of mass accretion in binaries that are not captured in single rotating stellar models. This possibly impacts the further evolution (either in a binary or as single stars), the final collapse, and the resulting spin of the compact object.

While most simulations of the epoch of reionization have focused on single-stellar populations in star-forming dwarf galaxies, products of binary evolution are expected to significantly contribute to emissions of hydrogen-ionizing photons. Among these products are stripped stars (or helium stars), which have their envelopes stripped from interactions with binary companions, leaving an exposed helium core. Previous work has suggested these stripped stars can dominate the Lyman Continuum (LyC) photon output of high-redshift, low-luminosity galaxies post-starburst. Other sources of hard radiation in the early universe include zero-metallicity Population III stars, which may have similar spectral energy distribution (SED) properties to galaxies with radiation dominated by stripped-star emissions. Here, we use four metrics (the power-law exponent over wavelength intervals 240-500 angstrom, 600-900 angstrom, and 1200-2000 angstrom, and the ratio of total luminosity in FUV wavelengths to LyC wavelengths) to compare the SEDs of simulated galaxies with only single-stellar evolution, galaxies containing stripped stars, and galaxies containing Population iii stars, with four different initial mass functions (IMFs). We find that stripped stars significantly alter SEDs in the LyC range of galaxies at the epoch of reionization. SEDs in galaxies with stripped stars have lower power-law indices in the LyC range and lower FUV to LyC luminosity ratios. These differences in SEDs are present at all considered luminosities (M-UV > -15, AB system), and are most pronounced for lower-luminosity galaxies. Intrinsic SEDs as well as those with interstellar medium absorption of galaxies with stripped stars and Population iii stars are found to be distinct for all tested Population III IMFs.