Abstract:
21 cm radiation from neutral hydrogen serves as a direct probe of the Epoch of Reionization. However, both its detection and physical interpretation are severely hindered by contamination from astrophysical foreground emission and instrumental systematics that are several orders of magnitude brighter than the signal of interest. A promising way to tackle these challenges is to cross-correlate the 21 cm signal with other independent tracers of large-scale structure, most notably high-redshift galaxies. Besides validating putative 21 cm detections, such joint analyses are expected to provide independent insights into the properties of ionizing sources and the evolving morphology of ionized regions during reionization. The 21 cm signal, however, is intrinsically highly non-Gaussian, limiting the effectiveness of conventional two-point cross-correlation statistics, which capture information only up to the second order. In this work, we therefore investigate the utility of k-nearest-neighbour cumulative distribution functions (kNN CDF), which encode information from the joint clustering at all orders, as an alternative framework for probing 21 cm-galaxy cross-correlations. Using self-consistently simulated mock 21 cm fields and a catalog of line-emitting galaxies at z = 7 that satisfy all available observations, we conduct a proof-of-concept study comparing the kNN CDF formalism and the two-point cross-correlation approach. We find that the kNN CDF statistics outperform the two-point statistics in detecting 21 cm-galaxy cross-correlations, even in the presence of instrumental noise and aggressive foreground filtering. Moreover, at a fixed global ionized fraction, it is even able to differentiate between reionization models that remain indistinguishable using two-point statistics. These results demonstrate the power and relatively unexplored potential of exploiting higher order statistics for extracting maximal information from 21 cm-galaxy synergies.