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 similar to 10(8)-10(12) M-circle dot (M* similar to 10(4)-10(10) M-circle dot) at z = 5-12. We find the sample average < f(esc)> increases with halo mass for M-vir similar to 10(8)-10(9.5) M-circle dot, becomes nearly constant for 109.5-1011 M-circle dot, and decreases at greater than or similar to 10(11) M-circle dot. Equivalently, < f(esc)> increases with stellar mass up to M* similar to 10(8) M-circle dot and decreases at higher masses. Even applying single-star stellar population synthesis models, we find a moderate < f(esc)> similar to 0.2 for galaxies at M* similar to 10(8) M-circle dot. 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 similar to 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.