No missing photons for reionization: moderate ionizing photon escape fractions from the FIRE-2 simulations

We present the escape fraction of hydrogen ionizing photons (f_(esc)) from a sample of 34 high-resolution cosmological zoom-in simulations of galaxies at z ≥ 5 in the Feedback in Realistic Environments project, post-processed with a Monte Carlo radiative transfer code for ionizing radiation. Our sample consists of 8500 haloes in M_(vir) ∼ 10⁸–10¹²M⊙ (M* ∼ 10⁴–10¹⁰M⊙⁠) at z = 5–12. We find the sample average〈f_(esc)〉increases with halo mass for M_(vir) ∼ 10⁸–10^(9.5) M⊙⁠, becomes nearly constant for 10^(9.5)–10¹¹M⊙, and decreases at ≳10¹¹M⊙. Equivalently,〈f_(esc)〉 increases with stellar mass up to M∗ ∼ 10⁸M⊙ and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate〈f_(esc)〉∼ 0.2 for galaxies at M∗ ∼ 10⁸M⊙⁠. Nearly half of the escaped ionizing photons come from stars 1–3 Myr old and the rest from stars 3–10 Myr old. Binaries only have a modest effect, boosting〈f_(esc)〉by ∼25–35 per cent and the number of escaped photons by 60–80 per cent. Most leaked ionizing photons are from vigorously star-forming regions that usually contain a feedback-driven kpc-scale superbubble surrounded by a dense shell. The shell is forming stars while accelerated, so new stars formed earlier in the shell are already inside the shell. Young stars in the bubble and near the edge of the shell can fully ionize some low-column-density paths pre-cleared by feedback, allowing a large fraction of their ionizing photons to escape. The decrease of〈f_(esc)〉at the high-mass end is due to dust attenuation, while at the low-mass end,〈f_(esc)〉decreases owing to inefficient star formation and hence feedback. At fixed mass,〈f_(esc)〉tends to increase with redshift. Although the absolute〈f_(esc)〉does not fully converge with resolution in our simulations, the mass- and redshift-dependence of〈f_(esc)〉 is likely robust. Our simulations produce sufficient ionizing photons for cosmic reionization