Dark-ages reionization and galaxy formation simulation -- XXI. Constraining the evolution of the ionizing escape fraction

Abstract: The fraction of ionizing photons that escape their host galaxies to ionize hydrogen in the inter-galactic medium (IGM) is a critical parameter in analyses of the reionization era. In this paper we use the Meraxes semi-analytic galaxy formation model to infer the mean ionizing photon escape fraction and its dependence on galaxy properties through joint modelling of the observed high redshift galaxy population and existing constraints on the reionization history. Using a Bayesian framework, and under the assumption that escape fraction is primarily related to halo mass, we find that the joint constraints of the UV luminosity function, CMB optical depth, and the Ly$\alpha$ forest require an escape fraction of $(18\pm5)\%$ for galaxies within haloes of $M\lesssim10^{9}$M$_\odot$ and $(5\pm2)\%$ for more massive haloes. In terms of galaxy properties, this transition in escape fraction occurs at stellar masses of $M_\star\sim10^7$M$_\odot$, nearly independent of redshift. As a function of redshift, reionization is dominated by the smaller $M_\star\lesssim10^7$M$_\odot$ galaxies with high escape fractions at $z\gtrsim6$ and by the larger $M_\star\gtrsim10^7$M$_\odot$ galaxies with lower escape fractions at $z\lesssim6$. Galaxies with star formation rates of $10^{-2.5}$M$_\odot$yr$^{-1}$ to $10^{-1.5}$M$_\odot$yr$^{-1}$ provide the dominant source of ionizing photons throughout reionization. Our results are consistent with recent direct measurements of a $\sim5\%$ escape fraction from massive galaxies at the end of reionization and support the picture of low mass galaxies being the dominant sources of ionizing photons during reionization.