Delayed photons from binary evolution help reionize the universe

High-resolution numerical simulations including feedback and aimed at calculating the escape fraction (fesc) 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 ≳70%. Moreover, simulations so far have yielded values of fesc falling only on the lower end of the ∼10%–20% range, the amount presumed necessary to reionize the universe. Analyzing a high-resolution (4 pc) cosmological radiation-hydrodynamic simulation, we study how fesc 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 fesc 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 ($\langle {f}_{\mathrm{esc}}\rangle $) by ∼13% and ∼10%, respectively, resulting in $\langle {f}_{\mathrm{esc}}\rangle =17 \% $. Our results emphasize that using updated stellar population synthesis models with binary stellar evolution provides a more sound physical basis for stellar reionization.