LIMFAST. III. Timing Cosmic Reionization with the 21 cm and Near-Infrared Backgrounds

TL;DR

This paper proposes a novel method to determine the timing of cosmic reionization by cross-correlating 21 cm signals with near-infrared backgrounds, overcoming foreground contamination issues.

Contribution

The study introduces a squared 21 cm signal cross-correlation technique to probe reionization history, validated through simulations and applicable to upcoming surveys.

Findings

The squared 21 cm signal correlates with NIRB during early reionization.

The correlation switches from positive to negative as reionization progresses.

The method can be detected with high significance by future surveys like SKA and SPHEREx.

Abstract

The timeline of cosmic reionization remains uncertain despite sustained efforts to study how the ionizing output of early galaxies shaped the intergalactic medium (IGM). Using the semi-numerical code LIMFAST, we investigate the prospects for timing the reionization process by cross-correlating the 21 cm signal with the cosmic near-infrared background (NIRB) contributed by galaxies at $z>5$. Tracing opposite phases of the IGM on large scales during reionization, the two signals together serve as a powerful probe for the reionization history. However, because long-wavelength, line-of-sight Fourier modes -- the only modes probed by NIRB fluctuations -- are contaminated by 21 cm foregrounds and thus inevitably lost to foreground cleaning or avoidance, a direct cross-correlation of the two signals vanishes. We show that this problem can be circumvented by squaring the foreground-filtered 21 cm signal and cross-correlating the squared field with the NIRB. This statistic is related to the 21 cm--21 cm--NIRB cross-bispectrum and encodes valuable information regarding the reionization timeline. Particularly, the 21 cm$^2$ and NIRB signals are positively correlated during the early phases of reionization and negatively correlated at later stages. We demonstrate that this behavior is generic across several different reionization models and compare our simulated results with perturbative calculations. We show that this cross-correlation can be detected at high significance by forthcoming 21 cm and NIRB surveys such as SKA and SPHEREx. Our methodology is more broadly applicable to cross-correlations between line intensity mapping data and 2D tracers of the large-scale structure, including photometric galaxy surveys and CMB lensing mass maps, among others.