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

Abstract

We derive optimal estimators for the two-, three-, and four-point correlators of statistically isotropic scalar fields defined on the sphere, such as the Cosmic Microwave Background temperature fluctuations, allowing for arbitrary (linear) masking and inpainting schemes. In each case, we give the optimal unwindowed estimator (obtained via a maximum-likelihood prescription, with an associated Fisher deconvolution matrix), and an idealized form, and pay close attention to their efficient computation. For the trispectrum, we include both parity-even and parity-odd contributions, as allowed by symmetry. The estimators can include arbitrary weighting of the data (and remain unbiased), but are shown to be optimal in the limit of inverse-covariance weighting and Gaussian statistics. The normalization of the estimators is computed via Monte Carlo methods, with the rate-limiting steps (involving spherical harmonic transforms) scaling linearly with the number of bins. An accompanying code package, PolyBin, implements these estimators in Python, and we demonstrate the estimators' efficacy via a suite of validation tests.

33 pages, 10 figures, code available at https://github.com/oliverphilcox/PolyBin. Accepted by Phys. Rev. D

References in corpus (10)

• The shape of primordial non-Gaussianity and the CMB bispectrum
• CMB Constraints on Primordial non-Gaussianity from the Bispectrum (f_{NL}) and Trispectrum (g_{NL} and τ_{NL}) and a New Consistency Test of Single-Field Inflation
• Probing Parity-Violation with the Four-Point Correlation Function of BOSS Galaxies
• Colliding Ghosts: Constraining Inflation with the Parity-Odd Galaxy Four-Point Function
• Cosmology Without Window Functions: Cubic Estimators for the Galaxy Bispectrum
• The CMB Bispectrum, Trispectrum, non-Gaussianity, and the Cramer-Rao Bound
• CMB lensing reconstruction biases from masking extragalactic sources
• Efficient Computation of $N$-point Correlation Functions in $D$ Dimensions
• Constraints on ${\tau }_{\mathrm{NL}}$ from Planck temperature and polarization
• Fast scalar quadratic maximum likelihood estimators for the CMB B-mode power spectrum

Cited by in corpus (2)

• The size and shape dependence of the SDSS galaxy bispectrum
• Searching for Unparticles with the Cosmic Microwave Background