A cross-bispectrum estimator for CMB-HI intensity mapping correlations
Intensity mapping of 21cm emission from neutral hydrogen (HI) promises to be a powerful probe of large-scale structure in the post-reionisation epoch. However, HI intensity mapping (IM) experiments will suffer the loss of long-wavelength line-of-sight HI modes in the galactic foreground subtraction process. The loss of these modes is particularly problematic for detecting HI IM cross-correlations with projected large-scale structure tracers, such as CMB secondary anisotropies. Here, we propose a cross-bispectrum estimator to recover the cross-correlation of the HI IM field, $\delta T_{21},$ with the CMB lensing field, $\kappa,$ constructed by correlating the position-dependent HI power spectrum with the mean overdensity traced by CMB lensing. We study the cross-bispectrum estimator, $B^{\bar \kappa \delta T_{21} \delta T_{21}},$ in the squeezed limit and forecast its detectability based on HI IM measurements from HIRAX and CMB lensing measurements from AdvACT. The cross-bispectrum improves constraints on cosmological parameters; in particular, the constraint on the dark energy equation-of-state parameter, $w_0,$ improves on the HI IM auto-power spectra constraint by 44\% (to 0.014), while the constraint on $w_a$ improves by 33\% (to 0.08), assuming Planck priors in each case. These results are robust to HI IM foreground removal because they largely derive from small-scale HI modes. The HI-HI-$\kappa$ cross-bispectrum thus provides a novel way to recover HI correlations with CMB lensing and constrain cosmological parameters at a level that is competitive with next-generation galaxy redshift surveys. As a striking example of this, we find that the combined constraint on the sum of the neutrino masses, while varying all redshift and standard cosmological parameters within a $w_0w_a\Omega_K$CDM model, is 5.5 meV.