Application of Monte Carlo Algorithms to the Bayesian Analysis of the Cosmic Microwave Background

Abstract

Power spectrum estimation and evaluation of associated errors in the presence of incomplete sky coverage; non-homogeneous, correlated instrumental noise; and foreground emission is a problem of central importance for the extraction of cosmological information from the cosmic microwave background. We develop a Monte Carlo approach for the maximum likelihood estimation of the power spectrum. The method is based on an identity for the Bayesian posterior as a marginalization over unknowns. Maximization of the posterior involves the computation of expectation values as a sample average from maps of the cosmic microwave background and foregrounds given some current estimate of the power spectrum or cosmological model, and some assumed statistical characterization of the foregrounds. Maps of the CMB are sampled by a linear transform of a Gaussian white noise process, implemented numerically with conjugate gradient descent. For time series data with N_{t} samples, and N pixels on the sphere, the method has a computational expense $KO[N^{2} +- N_{t} +AFw-log N_{t}], where K is a prefactor determined by the convergence rate of conjugate gradient descent. Preconditioners for conjugate gradient descent are given for scans close to great circle paths, and the method allows partial sky coverage for these cases by numerically marginalizing over the unobserved, or removed, region.

arXiv commentsubmitted to ApJ

References in corpus (13)

• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• A Line of Sight Approach to Cosmic Microwave Background Anisotropies
• Estimating the Power Spectrum of the Cosmic Microwave Background
• How to measure CMB power spectra without losing information
• A method for subtracting foregrounds from multi-frequency CMB sky maps
• Foreground separation methods for satellite observations of the cosmic microwave background
• II: Bayesian Methods for Cosmological Parameter Estimation from Cosmic Microwave Background Measurements
• An Efficient Technique to Determine the Power Spectrum from Cosmic Microwave Background Sky Maps
• The Age of the Universe and the Cosmological Constant Determined from Cosmic Microwave Background Anisotropy Measurements
• Making Maps Of The Cosmic Microwave Background: The MAXIMA Example
• First results from the Very Small Array -- IV. Cosmological parameter estimation
• Fast, exact CMB power spectrum estimation for a certain class of observational strategies
• Power Spectrum Estimators For Large CMB Datasets

Cited by in corpus (108)

• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP data
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Power Spectra and WMAP-Derived Parameters
• Planck 2013 results. XV. CMB power spectra and likelihood
• Planck 2018 results. IV. Diffuse component separation
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Joint Bayesian component separation and CMB power spectrum estimation
• Planck 2013 results. XII. Component separation
• Planck 2015 results. X. Diffuse component separation: Foreground maps
• Global, Exact Cosmic Microwave Background Data Analysis Using Gibbs Sampling
• Power spectrum estimation from high-resolution maps by Gibbs sampling
• The Atacama Cosmology Telescope: likelihood for small-scale CMB data
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Bayesian analysis of sparse anisotropic universe models and application to the 5-yr WMAP data
• Bayesian reconstruction of the cosmological large-scale structure: methodology, inverse algorithms and numerical optimization
• CMB component separation by parameter estimation
• Planck pre-launch status: the Planck-LFI programme
• Algorithms for bispectra: forecasting, optimal analysis, and simulation
• Bayesian power-spectrum inference for Large Scale Structure data
• Bayesian analysis of an anisotropic universe model: systematics and polarization
• Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
• Unmasking the Masked Universe: the 2M++ catalogue through Bayesian eyes
• Fast and Efficient Template Fitting of Deterministic Anisotropic Cosmological Models Applied to WMAP Data
• A re-analysis of the three-year WMAP temperature power spectrum and likelihood
• MADmap: A Massively Parallel Maximum-Likelihood Cosmic Microwave Background Map-Maker
• Cosmic Cartography of the Large-Scale Structure with Sloan Digital Sky Survey Data Release 6
• Bayesian Power Spectrum Analysis of the First-Year WMAP data
• Estimation of Polarized Power Spectra by Gibbs sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Reconstruction of signals with unknown spectra in information field theory with parameter uncertainty
• Planck Intermediate Results. IX. Detection of the Galactic haze with Planck
• Hybrid Estimation of CMB Polarization Power Spectra
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator
• Cosmic shear power spectra in practice
• Inflation and Early Dark Energy with a Stage II Hydrogen Intensity Mapping Experiment
• Methods for Bayesian power spectrum inference with galaxy surveys
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• Bayesian inference from photometric redshift surveys
• Inference with minimal Gibbs free energy in information field theory
• New estimates of the CMB angular power spectra from the WMAP 5 yrs low resolution data
• Searching for a dipole modulation in the large-scale structure of the Universe
• Bayesian component separation and CMB estimation for the 5-year WMAP temperature data
• TEASING: a fast and accurate approximation for the low multipole likelihood of the Cosmic Microwave Background temperature
• The joint large-scale foreground-CMB posteriors of the 3-year WMAP data
• Point-Source Power in 3 Year Wilkinson Microwave Anisotropy Probe Data
• A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data
• Bayesian semi-blind component separation for foreground removal in interferometric 21-cm observations
• The integrated bispectrum as a test of CMB non-Gaussianity: detection power and limits on f_NL with WMAP data
• Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps
• Matrix-free Large Scale Bayesian inference in cosmology
• The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
• Likelihood methods for CMB experiments
• Methods for pixel domain correction of EB leakage
• Probabilistic image reconstruction for radio interferometers
• Foreground Separation and Constraints on Primordial Gravitational Waves with the PICO Space Mission
• BeyondPlanck I. Global Bayesian analysis of the Planck Low Frequency Instrument data
• Efficient Optimal Reconstruction of Linear Fields and Band-powers from Cosmological Data
• CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator
• BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
• CMB Tomography: Reconstruction of Adiabatic Primordial Scalar Potential Using Temperature and Polarization Maps
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• Cosmology with the largest galaxy cluster surveys: Going beyond Fisher matrix forecasts
• Wiener filter reloaded: fast signal reconstruction without preconditioning
• A comparison of CMB Angular Power Spectrum Estimators at Large Scales: the TT case
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• Pure E and B polarization maps via Wiener filtering
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
• Bayesian power-spectrum inference with foreground and target contamination treatment
• Hybrid-basis inference for large-scale galaxy clustering: combining spherical and Cartesian Fourier analyses
• A multi-level solver for Gaussian constrained CMB realizations
• Bayesian angular power spectrum analysis of interferometric data
• Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze
• Data analysis methods for the cosmic microwave background
• Angular power spectrum of CMB anisotropy from WMAP
• Marginal distributions for cosmic variance limited CMB polarization data
• Bayesian analysis of white noise levels in the 5-year WMAP data
• Probing local non-Gaussianities within a Bayesian framework
• CMB power spectrum estimation with non-circular beam and incomplete sky coverage
• Non-Circular beam correction to the CMB power spectrum
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Large-Scale Polarized Foreground Component Separation for Planck
• Optimal and fast E/B separation with a dual messenger field
• Testing the ABS method with the simulated Planck temperature maps
• Optimized Large-Scale CMB Likelihood And Quadratic Maximum Likelihood Power Spectrum Estimation
• Bayesian Inference of Polarized CMB Power Spectra from Interferometric Data
• BeyondPlanck X. Planck LFI frequency maps with sample-based error propagation
• Fast scalar quadratic maximum likelihood estimators for the CMB B-mode power spectrum
• Forecast of foreground cleaning strategies for AliCPT-1
• Bayesian foreground analysis with CMB data
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• Bayesian noise estimation for non-ideal CMB experiments
• Cosmological Parameters from a re-analysis of the WMAP-7 low resolution maps
• Testing synchrotron models and frequency resolution in BINGO 21 cm simulated maps using GNILC
• A Bayesian estimate of the CMB-large-scale structure cross-correlation
• A CMB Gibbs sampler for localized secondary anisotropies
• PoPE: A population-based approach to model spatial structure of astronomical systems
• Bayesian analysis of cosmic structures
• A Monte Carlo comparison between template-based and Wiener-filter CMB dipole estimators
• CMB likelihood approximation for banded probability distributions
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
• SymPix: A spherical grid for efficient sampling of rotationally invariant operators
• Cosmology behind the mask: Constraining the parameters of $\Lambda$CDM with the unmasked galaxy density field from VIPERS
• Pixel domain multi-resolution minimum variance painting of CMB maps