High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (f(esc)) of hydrogen-ionizing photons often assume stellar radiation based on single-stellar population synthesis models. However, strong evidence suggests the binary fraction of massive stars is greater than or similar to 70%. Moreover, simulations so far have yielded values of f(esc) falling only on the lower end of the similar to 10%-20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study how f(esc) changes when we include two different products of binary stellar evolution-stars stripped of their hydrogen envelopes and massive blue stragglers. Both produce significant amounts of ionizing photons 10-200 Myr after each starburst. We find the relative importance of these photons to be amplified with respect to escaped ionizing photons, because peaks in star formation rates (SFRs) and f(esc) are often out of phase by this 10-200 Myr. Additionally, low-mass, bursty galaxies emit Lyman continuum radiation primarily from binary products when SFRs are low. Observations of these galaxies by the James Webb Space Telescope could provide crucial information on the evolution of binary stars as a function of redshift. Overall, including stripped stars and massive blue stragglers increases our photon-weighted mean escape fraction (< f(esc)>) by similar to 13% and similar to 10%, respectively, resulting in a < f(esc)> = 17%. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.

Binary interaction can cause stellar envelopes to be stripped, which significantly reduces the radius of the star. The orbit of a binary composed of a stripped star and a compact object can therefore be so tight that the gravitational radiation the system produces reaches frequencies accessible to the Laser Interferometer Space Antenna (LISA). Two such stripped stars in tight orbits with white dwarfs are known so far (ZTF J2130+4420 and CD-30 degrees 11223), but many more are expected to exist. These binaries provide important constraints for binary evolution models and may be used as LISA verification sources. We develop a Monte Carlo code that uses detailed evolutionary models to simulate the Galactic population of stripped stars in tight orbits with either neutron star or white dwarf companions. We predict 0-100 stripped star + white dwarf binaries and 0-4 stripped star + neutron star binaries with a signal-to-noise ratio >5 after 10 yr of observations with LISA. More than 90% of these binaries are expected to show large radial velocity shifts of greater than or similar to 200 km s(-1), which are spectroscopically detectable. Photometric variability due to tidal deformation of the stripped star is also expected and has been observed in ZTF J2130+4420 and CD-30 degrees 11223. In addition, the stripped star + neutron star binaries are expected to be X-ray bright with L-X greater than or similar to 10(33)-10(36) erg s(-1). Our results show that stripped star binaries are promising multimessenger sources for the upcoming electromagnetic and gravitational wave facilities.

Massive binaries that merge as compact objects are the progenitors of gravitational-wave sources. Most of these binaries experience one or more phases of mass transfer, during which one of the stars loses all or part of its outer envelope and becomes a stripped-envelope star. The evolution of the size of these stripped stars is crucial in determining whether they experience further interactions and understanding their ultimate fate. We present new calculations of stripped-envelope stars based on binary evolution models computed with MESA. We use these to investigate their radius evolution as a function of mass and metallicity. We further discuss their pre-supernova observable characteristics and potential consequences of their evolution on the properties of supernovae from stripped stars. At high metallicity, we find that practically all of the hydrogen-rich envelope is removed, which is in agreement with earlier findings. Only progenitors with initial masses below 10 M-circle dot expand to large radii (up to 100 R-circle dot), while more massive progenitors remain compact. At low metallicity, a substantial amount of hydrogen remains and the progenitors can, in principle, expand to giant sizes (> 400 R-circle dot) for all masses we consider. This implies that they can fill their Roche lobe anew. We show that the prescriptions commonly used in population synthesis models underestimate the stellar radii by up to two orders of magnitude. We expect that this has consequences for the predictions for gravitational-wave sources from double neutron star mergers, particularly with regard to their metallicity dependence.

Most massive stars experience binary interactions in their lifetimes that can alter both the surface and core structure of the stripped star with significant effects on their ultimate fate as core-collapse supernovae. However, core-collapse supernovae simulations to date have focused almost exclusively on the evolution of single stars. We present a systematic simulation study of single and binary-stripped stars with the same initial mass as candidates for core-collapse supernovae (11-21 M (circle dot)). Generally, we find that binary-stripped stars core tend to have a smaller compactness parameter, with a more prominent, deeper silicon/oxygen interface, and explode preferentially to the corresponding single stars of the same initial mass. Such a dichotomy of behavior between these two modes of evolution would have important implications for supernovae statistics, including the final neutron star masses, explosion energies, and nucleosynthetic yields. Binary-stripped remnants are also well poised to populate the possible mass gap between the heaviest neutron stars and the lightest black holes. Our work presents an improvement along two fronts, as we self-consistently account for the pre-collapse stellar evolution and the subsequent explosion outcome. Even so, our results emphasize the need for more detailed stellar evolutionary models to capture the sensitive nature of explosion outcome.

Context. Observations of massive stars in open clusters younger than similar to 8 Myr have shown that a majority of them are in binary systems, most of which will interact during their life. While these can be used as a proxy of the initial multiplicity properties, studying populations of massive stars older than similar to 20 Myr allows us to probe the outcome of these interactions after a significant number of systems have experienced mass and angular momentum transfer and may even have merged.

Gravitational-wave detections are nowprobing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on time-scales shorter than the convective turnover time-scale. This invalidates the steady-state assumption on which the classic mixing length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code MESA: (i) using the default implementation of Paxton et al. (2018) and (ii) solving an additional equation accounting for the time-scale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to similar to 5 M-circle dot. The differences are much smaller for the progenitors that determine the maximum mass of BHs below the gap. This prediction is robust at M-BH,M- max similar or equal to 48 M-circle dot, at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational-wave detectors.

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.

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.

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.