Cosmology with Intensity Mapping via Statistics Beyond the Power Spectrum in the SKAO Era

The cosmological distribution of neutral hydrogen (HI) during the post-reionization era is highly non-Gaussian due to the underlying non-linear structure formation, complex galaxy biasing, and potential primordial non-Gaussianity. One needs higher-order (beyond two-point) statistics to maximally extract the non-Gaussian information out of the 21-cm intensity maps. This chapter summarizes the potential of several higher-order statistics, including voxel intensity distribution, emission line stacking, probability density functions, $\ell_1$-norm, bispectrum, and various marked statistics. Additionally, image-based morphological descriptors, such as the Largest Cluster Statistic, local dimensions, and Minkowski functionals, etc., can potentially characterize the morphology and geometry of the cosmic web encoded in the 21-cm intensity maps. This chapter presents forecasts of the detectability of these higher-order statistics in the context of the future SKAO observations. These forecasts incorporate instrumental noise, observational effects, and, in some cases, foreground removal in their analyses. With its unprecedented sensitivity, the future SKAO 21-cm observations will enable us to measure these higher-order statistics more precisely, possibly helping to break degeneracies between astrophysical and cosmological parameters, and maximizing the science outcome from these surveys.