Ionizing Photon Escape Fraction

• Ionizing Photon Escape Fraction is the proportion of hydrogen-ionizing photons escaping from galaxies, defined by the ratio of emergent to intrinsic stellar outputs.
• It depends on halo mass, gas distribution, and star formation geometry, with simulations showing values from <1% in massive halos to 10–40% in low-mass systems.
• Temporal variability and anisotropic escape, driven by feedback and ISM clumping, play a crucial role in cosmic reionization modeling.

The ionizing photon escape fraction, typically denoted $f_{\rm esc}$, quantifies the proportion of hydrogen-ionizing (Lyman continuum; $\lambda < 912$ Å) photons produced by stellar populations in galaxies that escape into the intergalactic medium (IGM), thereby contributing to the reionization of cosmic hydrogen. The determination and modeling of $f_{\rm esc}$ is pivotal for understanding whether observed galaxies provide sufficient ionizing photons to explain the timing and morphology of reionization, as constrained by observations such as the Thomson optical depth, Lyman alpha emitters, and high-redshift quasar damping wings.

1. Formal Definition and Physical Basis

$f_{\rm esc}$ is formally defined for a galaxy or halo as: $$f_{\rm esc} = \frac{N_{\rm phot}(r \geq r_{200})}{N_{\rm emitted}}$$

where $N_{\rm emitted}$ denotes the total number of produced ionizing photons and $N_{\rm phot}(r \geq r_{200})$ counts those emerging beyond the virial radius (Paardekooper et al., 2015). This can equivalently be expressed in terms of luminosity,

$f_{\rm esc} = \frac{L_{\rm esc}}{L_{\rm int}}$

with $L_{\rm esc}$ the emergent ionizing luminosity and $L_{\rm int}$ the intrinsic value.

At the population level: $$\dot n_{\rm ion,gal}(z) = \rho_{\rm UV}(z)\; \xi_{\rm ion}(z,M_{\rm UV})\; f_{\rm esc}(M_h,z)$$ where $\rho_{\rm UV}$ is the UV luminosity density, $\xi_{\rm ion}$ is the production efficiency, and $f_{\rm esc}$ encodes halo mass and redshift dependence (Finkelstein et al., 2019).

2. Physical Drivers and Dependencies

Halo Mass Dependence

Radiation hydrodynamics simulations and radiative transfer post-processing demonstrate a strong, nonlinear halo mass dependence (Paardekooper et al., 2015, Yajima et al., 2010). For halos with virial mass below $10^8$ M$_\odot$, shallow potential wells allow SNe and massive-star feedback to clear low column-density channels near young clusters, yielding $\langle f_{\rm esc} \rangle \sim 10$–40%. By contrast, in $M_h \gtrsim 10^9$–$10^{10}$ M$_\odot$ systems, dense gas and deep central embedding of star-forming regions reduce $f_{\rm esc}$ to below 1%.

Summarized scaling from simulations:

$M_h$ [$M_\odot$]	$\langle f_{\rm esc} \rangle$
$\sim 10^9$	$\sim 0.4$
$\sim 10^{10}$	$\sim 0.15$
$\sim 10^{11}$	$\sim 0.07$

In the “First Billion Years” project, the escape fraction probability density function $P(f_{\rm esc}\mid M_h)$ is broad in each mass bin, with up to $\sim$1 dex scatter (Paardekooper et al., 2015, Yajima et al., 2010).

Gas Distribution and Star Formation Geometry

The local column density of neutral gas within $\sim$10 pc of star clusters, $N_{\rm H}$, is the principal constraint in high-$z$ halos (Paardekooper et al., 2015). Porosity of the ISM, number and density of clumps, and the offset of young clusters from the neutral gas centroid modulate $f_{\rm esc}$. Analytical and simulation models (Fernandez et al., 2010) show that fewer, denser clumps yield higher escape fractions due to the decreased probability of intercepting a sightline through each clump.

3. Temporal and Angular Anisotropy

Due to bursty and clustered star formation, and delayed feedback effects, $f_{\rm esc}(t)$ is highly time-variable on $\sim$Myr timescales (Kimm et al., 2014). SNe typically create low-density escape channels some $\sim$10 Myr after the star formation peak, producing brief intervals when the instantaneous escape fraction can exceed 20%, though the photon-weighted average over time is lower (typically $\sim$11–14% depending on feedback and IMF) (Kimm et al., 2014).

The angular escape fraction distribution is strongly anisotropic with most photons escaping through narrow “beams” (few steradian cones) (Paardekooper et al., 2015). In halos with $f_{\rm esc} \gtrsim 50\%$, the escape solid angle is $\sim 2\pi$ sr, but if $f_{\rm esc} \sim 1\%$, typical escape cones cover $0.3\pi$ sr.

4. Population-Averaged, Luminosity, and Redshift Dependence

To reconcile low mean escape fractions with global reionization constraints, Finkelstein et al. (Finkelstein et al., 2019) propose a halo-mass–dependent $f_{\rm esc}(M_h)$, scaled by a global factor ($f_{\rm esc,scale} = 5.2$, 68% CI 3.3–7.5), allowing galaxies down to the atomic cooling limit to contribute. Posterior-averaged population escape fractions span:

• $\bar f_{\rm esc} \approx 1\%$ at $z=4$
• $\bar f_{\rm esc} < 5\%$ at $z<9$
• $\bar f_{\rm esc} \sim 10\%$ by $z \sim 15$

Fainter galaxies ($M_{\rm UV} > -15$) reach $\bar f_{\rm esc} \sim 6$–12% at $z \gtrsim 6$, while brighter ($-20 < M_{\rm UV} < -16$) systems yield only $\sim$1–3%. These trends are crucial, as the steep faint-end slope of galaxy LFs at high $z$ ensures that faint objects dominate the ionizing photon budget (Finkelstein et al., 2019).

5. Observational Constraints and Tension

Population-averaged escape fractions are matched to key observables:

• Becker & Bolton (2013) emissivity at $z=4$–4.75
• Planck optical depth $\tau_{\rm es} = 0.055 \pm 0.009$, with models achieving $\sim 0.071 \pm 0.005$ (a $1.6\sigma$ offset)
• McGreer et al. (2015) $Q_{\rm HII}(z)$: reionization completes by $z=5.6\pm0.5$, midpoint $z \simeq 8.6\pm0.7$, $Q_{\rm HII}=0.78\pm0.08$ at $z=7$, mildly ($\sim$1–2$\sigma$) in tension with Ly$\alpha$ emitter and quasar damping wing constraints which favor $Q_{\rm HII,z=7} \sim 0.4$–0.6 (Finkelstein et al., 2019).

Models with a single $f_{\rm esc}$ at all redshifts/luminosities generally underproduce the observed ionizing emissivity, reinforcing the necessity of mass- and redshift-dependent escape fractions.

6. Role of AGN and Other Ionizing Sources

In Finkelstein et al. (Finkelstein et al., 2019), AGN contribute a non-negligible ($\sim$30% at $z=6$) but subdominant component to the overall ionizing budget, never dominating before $z\sim4.6$. Parameter posterior constraints yield an AGN scale factor of $0.77 (>0.47)$ and a redshift slope of $-0.32$ ($-0.84$ to $-0.14$) for the AGN term.

7. Implications for Reionization Modeling and Physical Interpretation

The dominance of ultra-faint dwarfs in the ionizing budget, broad stochasticity in $f_{\rm esc}(M_h)$, and strong anisotropy and time variability present significant modeling challenges. Semi-analytic and cosmological-volume reionization codes are advised to adopt mass-dependent and probabilistic $f_{\rm esc}$ distributions rather than deterministic or uniform prescriptions (Paardekooper et al., 2015). In practical terms, power-law or broken-power-law $f_{\rm esc}(M_h)$ relations with log-normal scatter are supported by simulation results (see schematic scaling in Section 2).

The requirement for low mean escape fractions ($\lesssim 5\%$) to suffice for reionization is satisfied only if galaxies form stars to the atomic cooling limit pre-reionization and photosuppression mass $\log M_{h,\rm supp}/M_\odot \approx 9$ is invoked after (Finkelstein et al., 2019).

Finally, physical processes—including supernova feedback, ISM clumpiness, geometry of star cluster embedding, and rising ionizing photon production efficiency at higher $z$/fainter $M_{UV}$—are crucial to a realistic understanding of $f_{\rm esc}$.

---

In summary, contemporary models and simulations indicate that a low average escape fraction ($\lesssim 5\%$) can reionize the universe, but only via a sharply mass-dependent $f_{\rm esc}$ scaling favoring low-mass, ultra-faint galaxies, with significant contributions from small-scale ISM structure, strong anisotropy, temporal stochasticity, and secondary AGN emission, matching most observational constraints to within $\sim2\sigma$ (Finkelstein et al., 2019, Paardekooper et al., 2015, Yajima et al., 2010).