Bayesian inference of CMB gravitational lensing

Abstract

The Planck satellite, along with several ground based telescopes, have mapped the cosmic microwave background (CMB) at sufficient resolution and signal-to-noise so as to allow a detection of the subtle distortions due to the gravitational influence of the intervening matter distribution. A natural modeling approach is to write a Bayesian hierarchical model for the lensed CMB in terms of the unlensed CMB and the lensing potential. So far there has been no feasible algorithm for inferring the posterior distribution of the lensing potential from the lensed CMB map. We propose a solution that allows efficient Markov Chain Monte Carlo sampling from the joint posterior of the lensing potential and the unlensed CMB map using the Hamiltonian Monte Carlo technique. The main conceptual step in the solution is a re-parameterization of CMB lensing in terms of the lensed CMB and the "inverse lensing" potential. We demonstrate a fast implementation on simulated data including noise and a sky cut, that uses a further acceleration based on a very mild approximation of the inverse lensing potential. We find that the resulting Markov Chain has short correlation lengths and excellent convergence properties, making it promising for application to high resolution CMB data sets of the future.

References in corpus (9)

• A General Framework for the Parametrization of Hierarchical Models
• A Measurement of the Cosmic Microwave Background B-Mode Polarization Power Spectrum at Sub-Degree Scales with POLARBEAR
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Efficient Wiener filtering without preconditioning
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Methods for Bayesian power spectrum inference with galaxy surveys
• Matrix-free Large Scale Bayesian inference in cosmology
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• A multi-level solver for Gaussian constrained CMB realizations

Cited by in corpus (30)

• Towards a cosmological neutrino mass detection
• Maximum a posteriori CMB lensing reconstruction
• CMB Delensing Beyond the B Modes
• Internal delensing of Planck CMB temperature and polarization
• Exploring cosmic origins with CORE: gravitational lensing of the CMB
• Bayesian delensing delight: sampling-based inference of the primordial CMB and gravitational lensing
• Bayesian delensing of CMB temperature and polarization
• Optimal CMB Lensing Reconstruction and Parameter Estimation with SPTpol Data
• A Demonstration of Improved Constraints on Primordial Gravitational Waves with Delensing
• CMB-S4: Iterative Internal Delensing and $r$ Constraints
• MUSE: Marginal Unbiased Score Expansion and Application to CMB Lensing
• Internal Delensing of Cosmic Microwave Background Acoustic Peaks
• CMB lensing power spectrum estimation without instrument noise bias
• Bayesian estimation of our local motion from the Planck-2018 CMB temperature map
• Minimizing gravitational lensing contributions to the primordial bispectrum covariance
• Limitations of CMB B-mode template delensing
• Wiener filter reloaded: fast signal reconstruction without preconditioning
• Comparison of delensing methodologies and assessment of the delensing capabilities of future experiments
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Joint Bayesian analysis of large angular scale CMB temperature anomalies
• Optimal and fast E/B separation with a dual messenger field
• Cosmic Microwave Background Mapmaking with a Messenger Field
• The primordial information content of Rayleigh Anisotropies
• Field-level multiprobe analysis of the CMB, integrated Sachs-Wolfe effect, and the galaxy density maps
• Preconditioner-free Wiener filtering with a dense noise matrix
• Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
• Approximate Likelihood Approaches for Detecting the Influence of Primordial Gravitational Waves in Cosmic Microwave Background Polarization
• LiteBIRD science goals and forecasts: improved full-sky reconstruction of the gravitational lensing potential through the combination of Planck and LiteBIRD data
• A Hamiltonian, post-Born, three-dimensional, on-the-fly ray tracing algorithm for gravitational lensing
• Bayesian inference methodology to characterize the dust emissivity at far-infrared and submillimeter frequencies