Joint Bayesian component separation and CMB power spectrum estimation

Abstract

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are 1) conditional sampling of foreground spectral parameters, and 2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground-CMB posterior distribution, and therefore all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multi-resolution observations. To verify the method, we analyse simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3-yr WMAP data, downgraded to a common resolution of 3 degree FWHM. The results from the actual 3-yr WMAP temperature analysis are presented in a companion Letter.

23 pages, 16 figures; version accepted for publication in ApJ -- only minor changes, all clarifications. More information about the WMAP3 analysis available at http://www.astro.uio.no/~hke under the Research tab

References in corpus (10)

• Wilkinson Microwave Anisotropy Probe (WMAP) Three Year Results: Implications for Cosmology
• Three Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Polarization Analysis
• Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Temperature Analysis
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Improved Measurements of the CMB Power Spectrum with ACBAR
• A re-analysis of the three-year WMAP temperature power spectrum and likelihood
• Estimation of Polarized Power Spectra by Gibbs sampling
• First season QUaD CMB temperature and polarization power spectra
• The joint large-scale foreground-CMB posteriors of the 3-year WMAP data
• Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps

Cited by in corpus (211)

• Planck 2013 results. XVI. Cosmological parameters
• Planck 2013 results. I. Overview of products and scientific results
• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Data Processing, Sky Maps, and Basic Results
• Planck 2018 results. I. Overview and the cosmological legacy of Planck
• Planck 2018 results. V. CMB power spectra and likelihoods
• Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
• Planck 2013 results. XV. CMB power spectra and likelihood
• Planck 2013 Results. XXIV. Constraints on primordial non-Gaussianity
• Planck 2018 results. IV. Diffuse component separation
• Planck intermediate results. LVII. Joint Planck LFI and HFI data processing
• Planck 2013 results. XII. Component separation
• Planck 2015 results. X. Diffuse component separation: Foreground maps
• Planck 2015 results. IX. Diffuse component separation: CMB maps
• Component separation methods for the Planck mission
• Planck 2013 results. XIII. Galactic CO emission
• Cosmological parameters derived from the final (PR4) Planck data release
• A refined model for spinning dust radiation
• MILCA, a Modified Internal Linear Combination Algorithm to extract astrophysical emissions from multi-frequency sky maps
• Planck constraints on the tensor-to-scalar ratio
• Bayesian analysis of sparse anisotropic universe models and application to the 5-yr WMAP data
• Maximum Likelihood algorithm for parametric component separation in CMB experiments
• Planck 2015 results. XXV. Diffuse low-frequency Galactic foregrounds
• Bayesian analysis of an anisotropic universe model: systematics and polarization
• Multi-resolution internal template cleaning: An application to the Wilkinson Microwave Anisotropy Probe 7-yr polarization data
• Planck 2018 results. II. Low Frequency Instrument data processing
• Planck 2013 results. II. Low Frequency Instrument data processing
• Sensitivity and foreground modelling for large-scale CMB B-mode polarization satellite missions
• Estimating the large-scale angular power spectrum in the presence of systematics: a case study of Sloan Digital Sky Survey quasars
• Libsharp - spherical harmonic transforms revisited
• Updated design of the CMB polarization experiment satellite LiteBIRD
• Exploring Cosmic Origins with CORE: B-mode Component Separation
• Peeling off foregrounds with the constrained moment ILC method to unveil primordial CMB $B$-modes
• 21cm foregrounds and polarization leakage: a user's guide on cleaning and mitigation strategies
• Planck Intermediate Results. IX. Detection of the Galactic haze with Planck
• Planck intermediate results. XLVIII. Disentangling Galactic dust emission and cosmic infrared background anisotropies
• Spatial variations in the spectral index of polarized synchrotron emission in the 9 yr WMAP sky maps
• Joint Planck and WMAP CMB Map Reconstruction
• Planck intermediate results. L. Evidence for spatial variation of the polarized thermal dust spectral energy distribution and implications for CMB $B$-mode analysis
• A Case Against Spinning PAHs as the Source of the Anomalous Microwave Emission
• The C-Band All-Sky Survey (C-BASS): Design and capabilities
• The Galactic Faraday depth sky revisited
• QUIJOTE scientific results -- IV. A northern sky survey in intensity and polarization at 10-20GHz with the Multi-Frequency Instrument
• Large-scale Maps of the Cosmic Infrared Background from Planck
• Portrait of a Dark Horse: a Photometric and Spectroscopic Study of the Ultra-faint Milky Way Satellite Pegasus III
• Galactic foreground contributions to the WMAP5 maps
• Searching for Modified Gravity: Scale and Redshift Dependent Constraints from Galaxy Peculiar Velocities
• Sparse component separation for accurate CMB map estimation
• Full-sky Cosmic Microwave Background Foreground Cleaning Using Machine Learning
• Mitigating Complex Dust Foregrounds in Future CMB Polarization Experiments
• Filling in CMB map missing data using constrained Gaussian realizations
• The joint large-scale foreground-CMB posteriors of the 3-year WMAP data
• Bayesian component separation and CMB estimation for the 5-year WMAP temperature data
• A minimal power-spectrum-based moment expansion for CMB B-mode searches
• Reconstruction of high-resolution SZ maps from heterogeneous datasets using needlets
• Impact of Galactic polarized emission on B-mode detection at low multipoles
• The bispectrum of polarized galactic foregrounds
• A Markov Chain Monte Carlo Algorithm for analysis of low signal-to-noise CMB data
• Stellar Proper Motions in the Orion Nebula Cluster
• CMB reconstruction from the WMAP and Planck PR2 data
• BeyondPlanck III. Commander3
• Direct detection of the kinetic Sunyaev-Zel'dovich effect in galaxy clusters
• Simulated forecasts for primordial B-mode searches in ground-based experiments
• Combining ILC and moment expansion techniques for extracting average-sky signals and CMB anisotropies
• Alignments of parity even/odd-only multipoles in CMB
• Search for anomalous alignments of structures in Planck data using Minkowski Tensors
• The ISW-tSZ cross correlation: ISW extraction out of pure CMB data
• Likelihood methods for CMB experiments
• $B$-mode forecast of CMB-Bh$\overline{a}$rat
• Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression
• Extracting foreground-obscured $\mu$-distortion anisotropies to constrain primordial non-Gaussianity
• 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
• Recovering the CMB Signal with Machine Learning
• Detection of Spectral Variations of Anomalous Microwave Emission with QUIJOTE and C-BASS
• CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator
• Maps of the number of HI clouds along the line of sight at high galactic latitude
• A Millimeter-Wave Galactic Plane Survey With The BICEP Polarimeter
• BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
• An improved Compton parameter map of thermal Sunyaev-Zeldovich effect from Planck PR4 data
• Impact on the tensor-to-scalar ratio of incorrect Galactic foreground modelling
• Framework for analysis of next generation, polarised CMB data sets in the presence of galactic foregrounds and systematic effects
• Multi-Clustering Needlet-ILC for CMB B-modes component separation
• Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• The BINGO Project V: Further steps in Component Separation and Bispectrum Analysis
• Statistical recovery of 21cm visibilities and their power spectra with Gaussian constrained realisations and Gibbs sampling
• Determining thermal dust emission from Planck HFI data using a sparse, parametric technique
• The Atacama Cosmology Telescope: Map-Based Noise Simulations for DR6
• Disentangling the Faraday rotation sky
• The bispectra of galactic CMB foregrounds and their impact on primordial non-Gaussianity estimation
• An alternative validation strategy for the Planck cluster catalog and $y$-distortion maps
• SILC: a new Planck Internal Linear Combination CMB temperature map using directional wavelets
• LiteBIRD Science Goals and Forecasts: Improving Sensitivity to Inflationary Gravitational Waves with Multitracer Delensing
• Local non-Gaussianity in the Cosmic Microwave Background the Bayesian way
• QUIJOTE scientific results -- VIII. Diffuse polarized foregrounds from component separation with QUIJOTE-MFI
• An Improved Diffuse Foreground Subtraction by ILC method: CMB Map and Angular Power Spectrum using Planck and WMAP Observations
• Tracing the cosmological evolution of stars and cold gas with CMB spectral surveys
• A Hint of a Blue Axion Isocurvature Spectrum?
• BeyondPlanck X. Bandpass and beam leakage corrections
• Exploring two-spin internal linear combinations for the recovery of the CMB polarization
• On the impact of large angle CMB polarization data on cosmological parameters
• CENN: A fully convolutional neural network for CMB recovery in realistic microwave sky simulations
• The C-Band All-Sky Survey (C-BASS): Template Fitting of Diffuse Galactic Microwave Emission in the Northern Sky
• Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis
• The C-Band All-Sky Survey (C-BASS): Constraining diffuse Galactic radio emission in the North Celestial Pole region
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• Bayesian power-spectrum inference with foreground and target contamination treatment
• BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
• Spin-SILC: CMB polarisation component separation with spin wavelets
• The Effects of Bandpass Variations on Foreground Removal Forecasts for Future CMB Experiments
• New insights into foreground analysis of the WMAP five-year data using FASTICA
• A multi-level solver for Gaussian constrained CMB realizations
• Improved Galactic Foreground Removal for B-Modes Detection with Clustering Methods
• Full-sky beam convolution for cosmic microwave background applications
• Semi-blind Bayesian inference of CMB map and power spectrum
• A Global ILC Approach in Pixel Space over Large Angular Scales of the Sky using CMB Covariance Matrix
• Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction
• Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze
• Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
• Spectral imaging of galaxy clusters with Planck
• Analysis of NILC performance on B-modes data of sub-orbital experiments
• A hybrid map-$C_{\ell }$ component separation method for primordial CMB $B$-mode searches
• CMB Polarization and Temperature Power Spectra Estimation using Linear Combination of WMAP 5-year Maps
• Hemispherical Power Asymmetry in intensity and polarization for Planck PR4 data
• Joint Bayesian analysis of large angular scale CMB temperature anomalies
• Effect of Foregrounds on the CMBR Multipole Alignment
• Cosmo++: An Object-Oriented C++ Library for Cosmology
• BeyondPlanck XVI. Limits on Large-Scale Polarized Anomalous Microwave Emission from Planck LFI and WMAP
• QUIJOTE scientific results -- X. Spatial variations of Anomalous Microwave Emission along the Galactic plane
• Recovering Cosmic Microwave Background Polarization Signals with Machine Learning
• A LEKID-based CMB instrument design for large-scale observations in Greenland
• Characterizing line-of-sight variability of polarized dust emission with future CMB experiments
• Linearized iterative least-squares (LIL): A parameter fitting algorithm for component separation in multifrequency CMB experiments such as Planck
• Modeling parity-violating spectra in Galactic dust polarization with filaments and its applications to cosmic birefringence searches
• Robust constraints on tensor perturbations from cosmological data: a comparative analysis from Bayesian and frequentist perspectives
• Testing for spectral index variations in polarised CMB foregrounds
• Impact of SZ cluster residuals in CMB maps and CMB-LSS cross-correlations
• Impact of modelling foreground uncertainties on future CMB polarization satellite experiments
• Marginal distributions for cosmic variance limited CMB polarization data
• X-ray Clusters of Galaxies as Cosmological Tools
• A Perturbative Analysis of Synchrotron Spectral Index Variation over Microwave Sky
• Needlet thresholding methods in component separation
• Bayesian analysis of white noise levels in the 5-year WMAP data
• Large-Scale Polarized Foreground Component Separation for Planck
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Statistical properties of polarized CMB foreground maps
• Impact of Galactic dust non-Gaussianity on searches for B-modes from inflation
• Revisiting the large-scale CMB anomalies: The impact of the SZ signal from the Local Universe
• Foreground removal with ILC methods for AliCPT-1
• Cross-correlation analysis of CMB with foregrounds for residuals
• Planck intermediate results. LV. Reliability and thermal properties of high-frequency sources in the Second Planck Catalogue of Compact Sources
• Impact of half-wave plate systematics on the measurement of CMB $B$-mode polarization
• The integrated angular bispectrum
• Testing for foreground residuals in the Planck foreground cleaned maps: A new method for designing confidence masks
• Joint Bayesian estimation of tensor and lensing B-modes in the power spectrum of CMB polarization data
• The C-Band All-Sky Survey (C-BASS): New Constraints on the Integrated Radio Spectrum of M 31
• An Unbiased Estimator of the Full-sky CMB Angular Power Spectrum at Large Scales using Neural Networks
• Forecast of foreground cleaning strategies for AliCPT-1
• Non-parametric maximum likelihood component separation for CMB polarization data
• Cosmic dipole tensions: confronting the cosmic microwave background with infrared and radio populations of cosmological sources
• Sensitivity of Bayesian 21 cm power spectrum estimation to foreground model errors
• Bayesian noise estimation for non-ideal CMB experiments
• A Differentiable, End-to-End Forward Model for 21 cm Cosmology: Estimating the Foreground, Instrument, and Signal Joint Posterior
• BeyondPlanck II. CMB map-making through Gibbs sampling
• Testing synchrotron models and frequency resolution in BINGO 21 cm simulated maps using GNILC
• QUIJOTE scientific results -- VI. The Haze as seen by QUIJOTE
• Observations of Galactic star-forming regions with the Cosmic Background Imager at 31 GHz
• Field-level multiprobe analysis of the CMB, integrated Sachs-Wolfe effect, and the galaxy density maps
• The Detector System for the Stratospheric Kinetic Inductance Polarimeter (SKIP)
• New probe of inflationary gravitational waves: cross-correlations of lensed primary CMB B-modes with large-scale structure
• Identifying frequency decorrelated dust residuals in B-mode maps by exploiting the spectral capability of bolometric interferometry
• Cosmoglobe DR1. III. First full-sky model of polarized synchrotron emission from all WMAP and Planck LFI data
• Bayesian decomposition of the Galactic multi-frequency sky using probabilistic autoencoders
• Unexplained correlation between the Cosmic Microwave Background temperature and the local matter density distribution
• Statistical estimation of full-sky radio maps from 21cm array visibility data using Gaussian constrained realizations
• A CMB Gibbs sampler for localized secondary anisotropies
• Examining statistical isotropy of CMB low multipoles from Planck PR4 data
• A new look at the integrated radio/microwave continuum spectrum of Galactic supernova remnant IC443
• Reconstructing the epoch of reionisation with Planck PR4
• An Accurate Reconstruction of CMB E Mode Signal over Large Angular Scales using Prior Information of CMB Covariance Matrix in ILC Algorithm
• Level correlations of the CMB temperature angular power spectrum
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Accelerating linear system solvers for time domain component separation of cosmic microwave background data
• Application of Genetic Algorithm to Estimate the Large Angular Scale Features of Cosmic Microwave Background
• Extended Delta-map: a map-based foreground removal method for CMB polarization observations
• Statistical properties of Galactic synchrotron temperature and polarization maps -- a multi-frequency comparison
• Self-supervised component separation for the extragalactic submillimeter sky
• CMB likelihood approximation for banded probability distributions
• CMB-PAInT: An inpainting tool for the cosmic microwave background
• Deep Needlet: A CNN based full sky component separation method in Needlet space
• The two-and three-point correlation functions of the polarized five-year WMAP sky maps
• On the Impact of Random Residual Calibration Error on the Gibbs ILC Estimates over Large Angular Scales
• Polarized CMB recovery with sparse component separation
• Foreground removal and angular power spectrum estimation of 21 cm signal using harmonic space ILC method
• Delensing for Precision Cosmology: Optimizing Future CMB B-mode Surveys to Constrain r
• Thermal Sunyaev-Zeldovich measurements and cosmic infrared background leakage mitigation combining upcoming ground-based telescopes
• Unbiased Primordial Gravitational Wave Inference from the CMB with SMICA
• Bayesian recalibration of flux scale factors in diffuse radio maps using low-resolution absolute radiometers
• A Foreground Model Independent Bayesian CMB Temperature and Polarization Signal Reconstruction and Cosmological Parameter Estimation over Large Angular Scales
• Implications of scattering for CMB foreground emission modelling
• Cosmic backgrounds from the radio to the far-infrared: recent results and perspectives from cosmological and astrophysical surveys
• The CMB angular power spectrum via component separation: a study on Planck data
• Pixel domain multi-resolution minimum variance painting of CMB maps
• Accurate and efficient likelihood modeling for large-scale CMB data
• Accidental deep field bias in CMB T and SNe $z$ correlation
• Mitigating Cosmic Microwave Background Shadow Degradation of Tensor-to-scalar Ratio Measurements through Map-based Studies
• Joint Bayesian separation and restoration of CMB from convolutional mixtures
• Do We Have Sufficient Knowledge of the Galactic Foreground Emission in Cosmic Microwave Background Science?
• Forecasts of constraining isotropic cosmic birefringence on AliCPT-1
• A new approach of estimating the Galactic thermal dust and synchrotron polarized emission template in the microwave bands
• B-mode polarization forecasts for GreenPol