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>.

An Approximation to the Likelihood Function for Band-Power Estimates of CMB Anisotropies

Matching Bayesian and frequentist coverage probabilities when using an approximate data covariance matrix

Exploring the likelihood of the 21-cm power spectrum with simulation-based inference

Assumptions of the primordial spectrum and cosmological parameter estimation

Galaxy 2-Point Covariance Matrix Estimation for Next Generation Surveys

A foreground-marginalized 'BK-lite' likelihood for the tensor-to-scalar ratio

A Partial-Sky Gibbs ILC Approach for the Estimation of CMB Posterior over Large Angular Scales of the Sky

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

A Minkowski Functional Analysis of the Cosmic Microwave Background Weak Lensing Convergence

candl: cosmic microwave background analysis with a differentiable likelihood

Frequentist Estimation of Cosmological Parameters from the MAXIMA-1 Cosmic Microwave Background Anisotropy Data

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

Non-Parametric Approximations for Anisotropy Estimation in Two-dimensional Differentiable Gaussian Random Fields

Compressed 'CMB-lite' Likelihoods Using Automatic Differentiation

Parameter constraints for flat cosmologies from CMB and 2dFGRS power spectra

Cosmological Parameter Estimation with a Joint-Likelihood Analysis of the Cosmic Microwave Background and Big Bang Nucleosynthesis

Flat-sky Angular Power Spectra Revisited

Beyond Gaussian Noise: A Generalized Approach to Likelihood Analysis with non-Gaussian Noise

Cosmological measurements from the CMB and BAO are insensitive to the tail probability in the assumed likelihood

Estimating the gravitational wave background anisotropy: a Bayesian approach boosted by cross-correlation angular power spectrum