Abstract:We revisit the flat-sky approximation for evaluating the angular power spectra of projected random fields by retaining information about the correlations along the line of sight. For the case of projections with broad, overlapping radial window functions, these line-of-sight correlations are suppressed and are ignored in the commonly adopted Limber approximation. However, retaining the correlations is important for narrow window functions or unequal-time spectra but introduces significant computational difficulties due to the highly oscillatory nature of the integrands involved. We deal with the integral over line-of-sight wave-modes in the flat-sky approximation analytically, using the FFTlog expansion of the 3D power spectrum. This results in an efficient computational method, which is a substantial improvement compared to any full-sky approaches. We apply our results to galaxy clustering (with and without redshift-space distortions), CMB lensing and galaxy lensing observables in a flat ΛCDM universe. In the case of galaxy clustering, we find excellent agreement with the full-sky results on large (percent-level agreement) and intermediate or small (subpercent agreement) scales, dramatically out-performing the Limber approximation for both wide and narrow window functions, and in equal- and unequal-time cases. In the cases of lensing, we show on the full-sky that the angular power spectrum of the lensing convergence can be very well approximated by projecting the 3D Laplacian (rather than the correct angular Laplacian) of the gravitational potential, even on large scales. Combining this approximation with our flat-sky techniques provides an efficient and accurate evaluation of the CMB lensing angular power spectrum on all scales. We further analyse the clustering and lensing angular power spectra by isolating the projection effects due to the observable- and survey-specific window functions, separating them from the effects due to integration along the line of sight and unequal-time mixing in the 3D power spectrum. All of the angular power spectrum results presented in this paper are obtained using a Python code implementation, which we make publicly available.

Gabor Transforms on the Sphere with Applications to CMB Power Spectrum Estimation

CMB Power Spectrum Estimation Using Wavelets

Minimum-variance multitaper spectral estimation on the sphere

Optimal Estimation of the Binned Mask-Free Power Spectrum, Bispectrum, and Trispectrum on the Full Sky: Scalar Edition

The Galaxy Bispectrum in the Spherical Fourier-Bessel Basis

Exact Modeling of Power Spectrum Multipole through Spherical Fourier-Bessel Basis

Optimal Estimation of the Binned Mask-Free Power Spectrum, Bispectrum, and Trispectrum on the Full Sky: Tensor Edition

Spherical bispectrum expansion and quadratic estimators

Galaxy bispectrum in the spherical Fourier-Bessel basis

Flat-sky angular power spectra revisited

Galaxy 2-Point Covariance Matrix Estimation for Next Generation Surveys

CMB polarization analysis on circular scans

The visibility based Tapered Gridded Estimator (TGE) for the redshifted 21-cm power spectrum

How to measure CMB power spectra without losing information

Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations

Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling

Exactly Azimuthal Pixelizations of the Sky

Minkowski Functionals of CMB polarisation intensity with Pynkowski: theory and application to Planck and future data