Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation

Abstract

We study different variants of the Gibbs sampler algorithm from the perspective of their applicability to the estimation of power spectra of the cosmic microwave background (CMB) anisotropies. These include approaches studied earlier in the CMB literature as well as new ones which are proposed in this work. We demonstrate all these variants on full and cut sky simulations and compare their performance, assessing both their computational and statistical efficiency. For this we employ a consistent comparison metric, an effective sample size (ESS) per second, commonly used in this context in the statistical literature. We show that one of the proposed approaches, referred to as Centered overrelax, which capitalizes on additional, auxiliary variables to minimize computational time needed per sample, and uses overrelaxation to decorrelate subsequent samples, performs better than the standard Gibbs sampler by a factor between one and two orders of magnitude in the nearly full-sky, satellite-like cases. It therefore potentially provides an interesting alternative to the currently favored approaches.

References in corpus (16)

• Array Programming with NumPy
• The Simons Observatory: Science goals and forecasts
• A General Framework for the Parametrization of Hierarchical Models
• Maximum Likelihood algorithm for parametric component separation in CMB experiments
• Polarized CMB power spectrum estimation using the pure pseudo-cross-spectrum approach
• Estimation of cosmological parameters using adaptive importance sampling
• Estimation of Polarized Power Spectra by Gibbs sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Efficient Wiener filtering without preconditioning
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Efficient Gaussian Sampling for Solving Large-Scale Inverse Problems using MCMC Methods
• Properties and use of CMB power spectrum likelihoods
• Likelihood methods for CMB experiments
• Exact joint likelihood of pseudo-$C_{\ell }$ estimates from correlated Gaussian cosmological fields
• Solving linear equations with messenger-field and conjugate gradients techniques - an application to CMB data analysis
• Robust Likelihoods for Inflationary Gravitational Waves from Maps of Cosmic Microwave Background Polarization