Accurate and efficient likelihood modeling for large-scale CMB data

Abstract

Accurate parameter estimation from cosmic microwave background data requires reliable likelihood modeling, particularly at large angular scales where angular power spectrum estimators exhibit non-Gaussian statistics. We present a novel approach, based on the Hamimeche-Lewis formalism, that marginalizes over auto-spectra, thus reducing residual biases from noise misestimation and partial sky coverage. We validate our approach by simulating three independent CMB channels, or data splits, in a multi-field setting, comparing to the pixel-based likelihood ground truth estimates for the optical depth $\tau$ and the tensor-to-scalar ratio $r$. We benchmark our method against the main power spectrum based alternatives available in the literature, showing that it outperforms all of them in terms of accuracy, while remaining fast and computationally efficient.

23 pages, 10 figures. One appendix moved to main body, minor modifications. Published on JCAP

References in corpus (39)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Detection of the Baryon Acoustic Peak in the Large-Scale Correlation Function of SDSS Luminous Red Galaxies
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
• BICEP2 I: Detection Of B-mode Polarization at Degree Angular Scales
• The Simons Observatory: Science goals and forecasts
• Planck 2018 results. V. CMB power spectra and likelihoods
• A Joint Analysis of BICEP2/Keck Array and Planck Data
• BICEP / Keck XIII: Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season
• Planck intermediate results. XLVI. Reduction of large-scale systematic effects in HFI polarization maps and estimation of the reionization optical depth
• The Quest for B Modes from Inflationary Gravitational Waves
• Joint Bayesian component separation and CMB power spectrum estimation
• A full sky, low foreground, high resolution CMB map from WMAP
• Improved limits on the tensor-to-scalar ratio using BICEP and Planck
• Likelihood Analysis of CMB Temperature and Polarization Power Spectra
• CMB power spectra and cosmological parameters from Planck PR4 with CamSpec
• Power spectrum estimation from high-resolution maps by Gibbs sampling
• Cosmological parameters derived from the final (PR4) Planck data release
• Planck constraints on the tensor-to-scalar ratio
• Maximum Likelihood algorithm for parametric component separation in CMB experiments
• Reionization optical depth determination from Planck HFI data with ten percent accuracy
• Updated constraints on amplitude and tilt of the tensor primordial spectrum
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• SRoll2: an improved mapmaking approach to reduce large-scale systematic effects in the Planck High Frequency Instrument legacy maps
• The impact of cosmic neutrinos on the gravitational-wave background
• Inference of the optical depth to reionization from low multipole temperature and polarisation Planck data
• Large-scale CMB temperature and polarization cross-spectra likelihoods
• Likelihood methods for the combined analysis of CMB temperature and polarisation power spectra
• TEASING: a fast and accurate approximation for the low multipole likelihood of the Cosmic Microwave Background temperature
• A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data
• Likelihood methods for CMB experiments
• Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations
• Multi-Clustering Needlet-ILC for CMB B-modes component separation
• On the impact of large angle CMB polarization data on cosmological parameters
• A novel CMB polarization likelihood package for large angular scales built from combined WMAP and Planck LFI legacy maps
• Analysis of NILC performance on B-modes data of sub-orbital experiments
• Robust constraints on tensor perturbations from cosmological data: a comparative analysis from Bayesian and frequentist perspectives
• B-mode constraints from Planck low multipole polarisation data