The N5K Challenge: Non-Limber Integration for LSST Cosmology

Abstract

The rapidly increasing statistical power of cosmological imaging surveys requires us to reassess the regime of validity for various approximations that accelerate the calculation of relevant theoretical predictions. In this paper, we present the results of the 'N5K non-Limber integration challenge', the goal of which was to quantify the performance of different approaches to calculating the angular power spectrum of galaxy number counts and cosmic shear data without invoking the so-called 'Limber approximation', in the context of the Rubin Observatory Legacy Survey of Space and Time (LSST). We quantify the performance, in terms of accuracy and speed, of three non-Limber implementations: $\mathtt{FKEM}(\mathtt{CosmoLike})$, $\mathtt{Levin}$, and $\mathtt{matter}$, themselves based on different integration schemes and approximations. We find that in the challenge's fiducial 3x2pt LSST Year 10 scenario, $\mathtt{FKEM}(\mathtt{CosmoLike})$ produces the fastest run time within the required accuracy by a considerable margin, positioning it favourably for use in Bayesian parameter inference. This method, however, requires further development and testing to extend its use to certain analysis scenarios, particularly those involving a scale-dependent growth rate. For this and other reasons discussed herein, alternative approaches such as $\mathtt{matter}$ and $\mathtt{Levin}$ may be necessary for a full exploration of parameter space. We also find that the usual first-order Limber approximation is insufficiently accurate for LSST Year 10 3x2pt analysis on $\ell =200-1000$, whereas invoking the second-order Limber approximation on these scales (with a full non-Limber method at smaller $\ell$) does suffice.

15 pages, 8 figures. Journal-accepted version

References in corpus (5)

• The NumPy array: a structure for efficient numerical computation
• Extended Limber Approximation
• Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program
• The Wide Field Infrared Survey Telescope: 100 Hubbles for the 2020s
• Dark Energy Survey Year 3 Results: Multi-Probe Modeling Strategy and Validation

Cited by in corpus (16)

• KiDS-SBI: Simulation-based inference analysis of KiDS-1000 cosmic shear
• On the asymptotic connection between full- and flat-sky angular correlators
• Accuracy requirements on intrinsic alignments for Stage-IV cosmic shear
• LimberJack.jl: auto-differentiable methods for angular power spectra analyses
• KiDS-Legacy: Covariance validation and the unified OneCovariance framework for projected large-scale structure observables
• Window convolution of the galaxy clustering bispectrum
• Constraints on primordial non-Gaussianity from Quaia
• ${\text{SwiftC}}_{\ell }$: fast differentiable angular power spectra beyond Limber
• Calibrating baryonic effects in cosmic shear with external data in the LSST era
• Cross-correlating radial peculiar velocities and CMB lensing convergence
• Linear Relativistic Corrections in the Spherical Fourier-Bessel Power Spectrum
• Cosmological forecasts from the combination of Stage-IV photometric galaxy surveys and the magnification from forthcoming GW observatories
• Reconstructing redshift distributions with photometric galaxy clustering
• 6x2pt: Forecasting gains from joint weak lensing and galaxy clustering analyses with spectroscopic-photometric galaxy cross-correlations
• Probing baryonic feedback with fast radio bursts: joint analyses with cosmic shear and galaxy clustering
• pylevin: Efficient numerical integration of integrals containing up to three Bessel functions