Fast optimal CMB power spectrum estimation with Hamiltonian sampling

Abstract

We present a method for fast optimal estimation of the temperature angular power spectrum from observations of the cosmic microwave background. We employ a Hamiltonian Monte Carlo (HMC) sampler to obtain samples from the posterior probability distribution of all the power spectrum coefficients given a set of observations. We compare the properties of the HMC and the related Gibbs sampling approach on low-resolution simulations and find that the HMC method performs favourably even in the regime of relatively low signal-to-noise. We also demonstrate the method on high-resolution data by applying it to simulated WMAP data. Analysis of a WMAP-sized data set is possible in a around eighty hours on a high-end desktop computer. HMC imposes few conditions on the distribution to be sampled and provides us with an extremely flexible approach upon which to build.

9 pages, 5 figures, submitted to Monthly Notices of the Royal Astronomical Society. Content changed in response to reviewers comments and new developments

References in corpus (20)

• Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results
• An All-Sky Analysis of Polarization in the Microwave Background
• MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets
• Estimating the Power Spectrum of the Cosmic Microwave Background
• Fast and reliable MCMC for cosmological parameter estimation
• Myths and Truths Concerning Estimation of Power Spectra
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Power spectrum estimation from high-resolution maps by Gibbs sampling
• An Efficient Technique to Determine the Power Spectrum from Cosmic Microwave Background Sky Maps
• Application of Monte Carlo Algorithms to the Bayesian Analysis of the Cosmic Microwave Background
• On Determining the Spectrum of Primordial Inhomogeneity from the COBE DMR Sky Maps: I. Method
• The Statistical Significance of the Low CMB Mulitipoles
• Bayesian Power Spectrum Analysis of the First-Year WMAP data
• Maximum-likelihood estimation of the CMB power spectrum from interferometer observations
• Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• Fast, exact CMB power spectrum estimation for a certain class of observational strategies
• Harmonic analysis of cosmic microwave background data II: from ring-sets to the sky

Cited by in corpus (52)

• A unified pseudo-$C_{\ell }$ framework
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• Markov Chain Monte Carlo Methods for Bayesian Data Analysis in Astronomy
• ELUCID - Exploring the Local Universe with reConstructed Initial Density field I: Hamiltonian Markov Chain Monte Carlo Method with Particle Mesh Dynamics
• Hyper-efficient model-independent Bayesian method for the analysis of pulsar timing data
• Fast Hamiltonian sampling for large scale structure inference
• Reconstructing the Initial Density Field of the Local Universe: Method and Test with Mock Catalogs
• Estimation of cosmological parameters using adaptive importance sampling
• Bayesian inference of cosmic density fields from non-linear, scale-dependent, and stochastic biased tracers
• Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys
• Methods for Bayesian power spectrum inference with galaxy surveys
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• COSMIC BIRTH: Efficient Bayesian Inference of the Evolving Cosmic Web from Galaxy Surveys
• Statistical methods for cosmological parameter selection and estimation
• 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
• Cross-Correlation of Planck CMB Lensing with DESI-Like LRGs
• Bayesian inference of CMB gravitational lensing
• CosmoPower-JAX: high-dimensional Bayesian inference with differentiable cosmological emulators
• Bayesian estimation of our local motion from the Planck-2018 CMB temperature map
• Wide-band Profile Domain Pulsar Timing Analysis
• KaRMMa -- Kappa Reconstruction for Mass Mapping
• Cosmology with the largest galaxy cluster surveys: Going beyond Fisher matrix forecasts
• The Clustering of Galaxies in the Completed SDSS-III Baryon Oscillation Spectroscopic Survey: Cosmic Flows and Cosmic Web from Luminous Red Galaxies
• A comparison of CMB Angular Power Spectrum Estimators at Large Scales: the TT case
• Bayesian cosmic density field inference from redshift space dark matter maps
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• ELUCID -- Exploring the Local Universe with the reConstructed Initial Density field. II: Reconstruction diagnostics, applied to numerical halo catalogs
• Reconstructing dark matter distribution with peculiar velocities: Bayesian forward modelling with corrections for inhomogeneous Malmquist bias
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• Full-sky beam convolution for cosmic microwave background applications
• Bayesian angular power spectrum analysis of interferometric data
• A practical guide to Basic Statistical Techniques for Data Analysis in Cosmology
• Joint Bayesian analysis of large angular scale CMB temperature anomalies
• Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
• Radio Weak Lensing Shear Measurement in the Visibility Domain - I. Methodology
• Estimating the Integrated Bispectrum from Weak Lensing Maps
• Marginal distributions for cosmic variance limited CMB polarization data
• Probing local non-Gaussianities within a Bayesian framework
• Radio Galaxy Shape Measurement with Hamiltonian Monte Carlo in the Visibility Domain
• Higher Order Hamiltonian Monte Carlo Sampling for Cosmological Large-Scale Structure Analysis
• Sparse Inpainting and Isotropy
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• Bayesian inference on the sphere beyond statistical isotropy
• Bayesian large-scale structure inference and cosmic web analysis
• Bayesian analysis of cosmic structures
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• SIToolBox : A package for Bayesian estimation of the isotropy violation in the CMB sky
• MCMC generation of cosmological fields far beyond Gaussianity
• Bayesian inference methodology to characterize the dust emissivity at far-infrared and submillimeter frequencies
• Accurate and efficient likelihood modeling for large-scale CMB data
• A Bayesian approach to high fidelity interferometric calibration I: mathematical formalism