Probing Cosmology through Higher-Order CMB Lensing Statistics

TL;DR

This paper demonstrates that non-Gaussian statistics of CMB lensing, especially Minkowski functionals and peak counts, significantly improve constraints on cosmological parameters beyond traditional power spectrum analysis for upcoming experiments like SO.

Contribution

It introduces a comprehensive pipeline for modeling higher-order CMB lensing statistics and quantifies their added value in constraining cosmology, highlighting the importance of morphology-based measures.

Findings

Minkowski functionals and peak counts reduce uncertainties on $\

Combining all non-Gaussian statistics with the power spectrum yields up to 70% reduction in neutrino mass uncertainty.

Non-Gaussian statistics provide significant complementary information even when combined with extended power spectrum and other data.

Abstract

We investigate the cosmological information in higher-order statistics of the cosmic microwave background (CMB) lensing convergence field for a near-term experiment with noise properties similar to the Simons Observatory (SO). Using a fully field-level forward-modeling pipeline based on ray-traced simulations from the MassiveNuS suite and realistic SO-like CMB lensing reconstruction, we naturally include nonlinear structure formation, post-Born effects, and higher-order reconstruction noise. We measure several non-Gaussian statistics, including Minkowski functionals, peak and minima counts, moments, and wavelet-scattering coefficients. We train Gaussian-process emulators to model each statistic's dependence on the matter density fraction $\Omega_m$, the scalar power spectrum amplitude $A_s$, and the neutrino mass sum $M_\nu$. We quantify the relative information gain these statistics provide beyond the lensing power spectrum and identify which are most robust to reconstruction noise. We find that morphology-based statistics, particularly Minkowski functionals and peak/minima counts, offer significant complementary constraining power: combining all non-Gaussian statistics with the power spectrum yields reductions of 40% and 38% in the marginalized uncertainties on $\Omega_m$ and $A_s$, respectively, and a 70% reduction in the one-sided uncertainty on $M_\nu$. These gains remain non-negligible even when the power spectrum is extended to larger scales and combined with primary CMB and BAO data, with Minkowski functionals providing an additional 11% improvement in $\sigma(M_\nu)$ and 35% in $\sigma(\Omega_m)$ beyond the extended power spectrum. By contrast, moments and wavelet-scattering coefficients provide more limited gains at SO noise levels. Our results highlight the potential of non-Gaussian statistics to enhance cosmological constraints from SO and future CMB surveys.