Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published July 11, 2016 | Submitted + Published
Journal Article Open

Binary Stars Can Provide the "Missing Photons" Needed for Reionization

Abstract

Empirical constraints on reionization require galactic ionizing photon escape fractions f_(esc) ≳ 20 per cent, but recent high-resolution radiation-hydrodynamic calculations have consistently found much lower values ∼1–5 per cent. While these models include strong stellar feedback and additional processes such as runaway stars, they almost exclusively consider stellar evolution models based on single (isolated) stars, despite the fact that most massive stars are in binaries. We re-visit these calculations, combining radiative transfer and high-resolution cosmological simulations with detailed models for stellar feedback from the Feedback in Realistic Environments project. For the first time, we use a stellar evolution model that includes a physically and observationally motivated treatment of binaries (the Binary Population and Spectral Synthesis model). Binary mass transfer and mergers enhance the population of massive stars at late times (≳3 Myr) after star formation, which in turn strongly enhances the late-time ionizing photon production (especially at low metallicities). These photons are produced after feedback from massive stars has carved escape channels in the interstellar medium, and so efficiently leak out of galaxies. As a result, the time-averaged 'effective' escape fraction (ratio of escaped ionizing photons to observed 1500 Å photons) increases by factors ∼4–10, sufficient to explain reionization. While important uncertainties remain, we conclude that binary evolution may be critical for understanding the ionization of the Universe.

Additional Information

© 2016 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. Accepted 2016 April 19. Received 2016 April 19; in original form 2016 January 27. We thank Chuck Steidel for helpful discussions and the referee for useful comments. We also thank John Beacom and Mike Shull for helpful suggestions after the paper was submitted to arXiv. The simulations used in this paper were run on XSEDE computational resources (allocations TG-AST120025, TG-AST130039, and TG-AST140023). The analysis was performed on the Caltech compute cluster 'Zwicky' (NSF MRI award #PHY-0960291). Support for PFH was provided by an Alfred P. Sloan Research Fellowship, NASA ATP Grant NNX14AH35G, and NSF Collaborative Research Grant #1411920 and CAREER grant #1455342. D. Kasen is supported in part by a Department of Energy Office of Nuclear Physics Early Career Award, and by the Director, Office of Energy Research, Office of High Energy and Nuclear Physics, Divisions of Nuclear Physics, of the US Department of Energy under Contract No. DE-AC02-05CH11231 and by the NSF through grant AST-1109896. D. Kereš was supported by NSF grant AST-1412153 and funds from the University of California, San Diego. CAFG was supported by NSF through grants AST-1412836 and AST-1517491, by NASA through grant NNX15AB22G, and by STScI through grant HST-AR-14293.001-A. EQ was supported by NASA ATP grant 12-APT12-0183, a Simons Investigator award from the Simons Foundation, and the David and Lucile Packard Foundation.

Attached Files

Published - MNRAS-2016-Ma-3614-9.pdf

Submitted - 1601.07559v1.pdf

Files

1601.07559v1.pdf
Files (678.8 kB)
Name Size Download all
md5:84928e85c7d7ec314b49027992a2e261
334.8 kB Preview Download
md5:0b068df34e68decf8ec9897b79e4ba6a
344.0 kB Preview Download

Additional details

Created:
August 20, 2023
Modified:
October 17, 2023