Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations

Abstract

We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. The net effective cost per independent sample is $\sim$2 000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for Planck LFI and WMAP in Cosmoglobe Data Release 1. Thus, although technically possible to run already in its current state, future work should aim to reduce the effective cost per independent sample by one order of magnitude to avoid excessive runtimes, for instance through multi-grid preconditioners and/or derivative-based Markov chain sampling schemes. This work demonstrates the computational feasibility of true Bayesian cosmological parameter estimation with end-to-end error propagation for high-precision CMB experiments without likelihood approximations, but it also highlights the need for additional optimizations before it is ready for full production-level analysis.

arXiv comment10 pages, 8 figures. Accepted in A&A

References in corpus (29)

• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
• Efficient Computation of CMB anisotropies in closed FRW models
• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results
• Planck 2018 results. V. CMB power spectra and likelihoods
• Cobaya: Code for Bayesian Analysis of hierarchical physical models
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Parameter Estimation Methodology
• Planck 2018 results. III. High Frequency Instrument data processing and frequency maps
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Likelihood Analysis of CMB Temperature and Polarization Power Spectra
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Power spectrum estimation from high-resolution maps by Gibbs sampling
• Application of Monte Carlo Algorithms to the Bayesian Analysis of the Cosmic Microwave Background
• Planck 2018 results. II. Low Frequency Instrument data processing
• Estimation of Polarized Power Spectra by Gibbs sampling
• Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing
• Optimal CMB Lensing Reconstruction and Parameter Estimation with SPTpol Data
• A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data
• MUSE: Marginal Unbiased Score Expansion and Application to CMB Lensing
• Bayesian inference of CMB gravitational lensing
• BeyondPlanck III. Commander3
• BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• Optimal constraints on primordial gravitational waves from the lensed CMB
• Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis
• A multi-level solver for Gaussian constrained CMB realizations
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• CMB likelihood approximation for banded probability distributions