Resolving the ionizing photon budget crisis with JWST/NIRCam HII clumping constraints at z=6

We present a comprehensive study of the ionizing properties of 1721 galaxies at $5.6<z<6.5$ using deep JWST/NIRCam photometric imaging from the NEP, JADES, and PRIMER surveys spanning an unmasked area $\sim550$arcmin$^2$ across UV magnitudes $-22\lesssim M_{\rm UV}\lesssim-17.5$. Our $90%$ stellar mass complete sample suggests little relation of UV slope with magnitude, $β_{\rm UV}=(-0.040\pm0.022)M_{\rm UV}-2.88^{+0.43}{-0.44}$, implying $f{\rm esc}^{\rm LyC}\simeq5%$ based on calibrations from the Low-redshift Lyman Continuum Survey (LzLCS). We measure a constant ionizing photon production efficiency with UV magnitude, $\log_{10}(ξ_{\rm ion, 0}/\rm Hz,erg^{-1}) = -0.006^{+0.019}{-0.017}~M{\rm UV} + 25.05^{+0.39}{-0.34}$, consistent with HST canonical values. The total production rate of photons escaping into the IGM is computed as $\log{10}(\dot{n}{\rm ion}/\rm s^{-1}Mpc^{-3})=50.31^{+0.07}{-0.06}$ for $M_{\rm UV}<-17$ galaxies from our star forming and smouldering UV luminosity functions (UVLFs), which differ in the faint-end slope ($α_{\rm SFG}=-2.2\pm0.2$; $α_{\rm sm}=-1.7\pm0.2$). Extrapolating to the latest UVLF turnover limits from the massive lensing galaxy cluster Abell S1063 ($M_{\rm UV, lim}=-13.5$) implies that a recombination-weighted HII clumping factor $C_{\rm HII, rec}=6.2^{+4.1}_{-2.1}$ is required to produce fully stably reionized at $z\simeq6$. A clumping factor of this magnitude resolves the ionizing photon budget crisis. Our methodology paves the way for indirect clumping measurements from galaxies which will provide insight into earlier stages of the EoR when the Ly$α$-forest becomes saturated and more direct quasar measurements become impossible.

💡 Research Summary

This paper presents a comprehensive analysis of the ionizing properties of galaxies at the tail end of the Epoch of Reionization (z ≈ 6), utilizing deep imaging data from the JWST/NIRCam instrument. By combining data from the NEP, JADES, and PRIMER surveys over an unmasked area of ~550 arcmin², the study constructs a robust sample of 1721 galaxies within the redshift range 5.6 < z < 6.5. The sample is further divided into “star-forming galaxies” (SFGs) and “smouldering galaxies” based on their star-formation burstiness parameter (Φ = SFR_10Myr / SFR_100Myr).

The core analysis involves deriving key galactic properties directly from the photometry. The study finds a negligible relationship between the UV continuum slope (β_UV) and absolute UV magnitude (M_UV). Using calibrations from the low-redshift Lyman Continuum Survey (LzLCS), this implies an average Lyman continuum escape fraction (f_esc^LyC) of approximately 5%. Crucially, the ionizing photon production efficiency (ξ_ion) is found to be constant with M_UV, with a value of log₁₀(ξ_ion / Hz erg⁻¹) = 25.05^{+0.39}_{-0.34}. This result aligns with canonical values from the HST era and contrasts with some recent JWST studies suggesting very high, magnitude-dependent ξ_ion, potentially highlighting selection biases in those works.

The researchers then construct separate UV luminosity functions (UVLFs) for the SFG and smouldering subsamples, which exhibit different faint-end slopes (α_SFG = -2.2 ± 0.2; α_sm = -1.7 ± 0.2). Integrating these UVLFs down to M_UV < -17, they compute the total production rate of ionizing photons escaping into the intergalactic medium (IGM) as log₁₀(ṅ_ion / s⁻¹ Mpc⁻³) = 50.31^{+0.07}_{-0.06}.

The pivotal finding comes from extrapolating this photon production rate to the faintest galaxies currently observable, down to a turnover limit of M_UV,lim = -13.5 as probed by strong lensing clusters like Abell S1063. The study performs a reverse calculation to determine the conditions required for the universe to be fully and stably reionized by z ≈ 6. They find that a recombination-weighted HII clumping factor of C_HII,rec = 6.2^{+4.1}_{-2.1} is necessary. A clumping factor of this magnitude increases the recombination rate within ionized bubbles in the IGM, thereby requiring a higher injection rate of ionizing photons to maintain ionization. This result directly resolves the so-called “ionizing photon budget crisis,” where high ξ_ion measurements from JWST seemed to imply an overproduction of photons, forcing reionization to end earlier than observed. The required clumping factor of ~6 brings the budget back into balance with observational constraints on the reionization timeline.

The methodology established in this work paves the way for using galaxy observables to indirectly constrain the clumpiness of the IGM. This approach will be particularly valuable for probing the earlier, more opaque stages of reionization (z > 6) when the Lyα forest is saturated and direct measurements using background quasars become impossible.