We propose a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$\times$galaxy $-$ galaxy$\times$IM; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As a pure cross-correlation estimator, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. For cosmic variance dominated surveys the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of DES cross-correlated with HIRAX, the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

Future galaxy clustering surveys will probe small scales where non-linearities become important. Since the number of modes accessible on intermediate to small scales is very high, having a precise model at these scales is important especially in the context of discriminating alternative cosmological models from the standard one. In the mildly non-linear regime, such models typically differ from each other, and galaxy clustering data will become very precise on these scales in the near future. As the observable quantity is the angular power spectrum in redshift space, it is important to study the effects of non-linear density and redshift space distortion (RSD) in the angular power spectrum. We compute non-linear contributions to the angular power spectrum using a flat-sky approximation that we introduce in this work, and compare the results of different perturbative approaches with $N$-body simulations. We find that the TNS perturbative approach is significantly closer to the $N$-body result than Eulerian or Lagrangian 1-loop approximations, effective field theory of large scale structure or a halofit-inspired model. However, none of these prescriptions is accurate enough to model the angular power spectrum well into the non-linear regime. In addition, for narrow redshift bins, $\Delta z \lesssim 0.01$, the angular power spectrum acquires non-linear contributions on all scales, right down to $\ell=2$, and is hence not a reliable tool at this time. To overcome this problem, we need to model non-linear RSD terms, for example as TNS does, but for a matter power spectrum that remains reasonably accurate well into the deeply non-linear regime, such as halofit.

In this paper, we study lensing of 21cm intensity mapping (IM). Like in the cosmic microwave background (CMB), there is no first order lensing in intensity mapping. The first effects in the power spectrum are therefore of second and third order. Despite this, lensing of the CMB power spectrum is an important effect that needs to be taken into account, which motivates the study of the impact of lensing on the IM power spectrum. We derive a general formula up to third order in perturbation theory including all the terms with two derivatives of the gravitational potential, i.e. the dominant terms on sub-Hubble scales. We then show that in intensity mapping there is a new lensing term which is not present in the CMB. We obtain that the signal-to-noise of 21 cm lensing for futuristic surveys like SKA2 is about 10. We find that surveys probing only large scales, lmax < 700, can safely neglect the lensing of the intensity mapping power spectrum, but that otherwise this effect should be included.