Abstract:Standard cosmic microwave background (CMB) analyses constrain cosmological and astrophysical parameters by fitting parametric models to multifrequency power spectra (MFPS). However, such methods do not optimally weight maps in power spectrum (PS) measurements for non-Gaussian CMB foregrounds. We propose needlet internal linear combination (NILC), operating on wavelets with compact support in pixel and harmonic space, as a weighting scheme to yield more optimal parameter constraints. In a companion paper, we derived an analytic formula for NILC map PS, which is physically insightful but computationally difficult to use in parameter inference pipelines. In this work, we analytically show that fitting parametric templates to MFPS and harmonic ILC PS yields identical parameter constraints when the number of sky components equals or exceeds the number of frequency channels. We numerically show that, for Gaussian random fields, the same holds for NILC PS. This suggests that NILC can reduce parameter error bars in the presence of non-Gaussian fields since it uses non-Gaussian information. As Gaussian likelihoods may be inaccurate, we use likelihood-free inference (LFI) with neural posterior estimation. We show that performing inference with auto- and cross-PS of NILC component maps as summary statistics yields smaller parameter error bars than inference with MFPS. For a model with CMB, an amplified thermal Sunyaev--Zel'dovich (tSZ) signal, and noise, we find a 60% reduction in the area of the 2D 68% confidence region for component amplitude parameters inferred from NILC PS, as compared to inference from MFPS. Primordial $B$-mode searches are a promising application for our new method, as the amplitude of the non-Gaussian dust foreground is known to be larger than a potential signal. Our code is available in <a class="link-external link-https" href="https://github.com/kmsurrao/NILC-Inference-Pipeline" rel="external noopener nofollow">this https URL</a>.

Catalog-based pseudo-$C_\ell$s

Pseudo-$C_\ell$s for spin-$s$ fields with component-wise weighting

BLAST: Beyond Limber Angular power Spectra Toolkit. A fast and efficient algorithm for 3x2 pt analysis

Galaxy 2-Point Covariance Matrix Estimation for Next Generation Surveys

Almanac: MCMC-based signal extraction of power spectra and maps on the sphere

Euclid preparation. LIX. Angular power spectra from discrete observations

Avoiding lensing bias in cosmic shear analysis

Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars

How to measure CMB power spectra without losing information

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

Modeling the 3-point correlation function of projected scalar fields on the sphere

Studying the Multi-frequency Angular Power Spectrum of the Cosmic Dawn 21-cm Signal

Fast estimation of aperture mass statistics II: Detectability of higher order statistics in current and future surveys

Hierarchical Cosmic Shear Power Spectrum Inference

Accurate power spectrum estimation toward Nyquist limit

Constraining Cosmological Parameters with Needlet Internal Linear Combination Maps II: Likelihood-Free Inference on NILC Power Spectra

The Angular Power Spectra of Photometric Sloan Digital Sky Survey Luminous Red Galaxies

The Probability Distribution Function of the Sunyaev‐Zel’dovich Power Spectrum: an Analytical Approach

I. Apples to apples $A^2$: realistic galaxy simulated catalogs and photometric redshift predictions for next-generation surveys

Noise power spectrum estimation and fast map making for CMB experiments

Bayesian inference methodology for Primordial Power Spectrum reconstructions from Large Scale Structure