Power spectrum estimation from high-resolution maps by Gibbs sampling

Abstract

We revisit a recently introduced power spectrum estimation technique based on Gibbs sampling, with the goal of applying it to the high-resolution WMAP data. In order to facilitate this analysis, a number of sophistications have to be introduced, each of which is discussed in detail. We have implemented two independent versions of the algorithm to cross-check the computer codes, and to verify that a particular solution to any given problem does not affect the scientific results. We then apply these programs to simulated data with known properties at intermediate (N_side = 128) and high (N_side = 512) resolutions, to study effects such as incomplete sky coverage and white vs. correlated noise. From these simulations we also establish the Markov chain correlation length as a function of signal-to-noise ratio, and give a few comments on the properties of the correlation matrices involved. Parallelization issues are also discussed, with emphasis on real-world limitations imposed by current super-computer facilities. The scientific results from the analysis of the first-year WMAP data are presented in a companion letter.

15 pages, 11 figures, accepted for publication in ApJS

References in corpus (20)

• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Foreground Emission
• MASTER of the CMB Anisotropy Power Spectrum: A Fast Method for Statistical Analysis of Large and Complex CMB Data Sets
• The significance of the largest scale CMB fluctuations in WMAP
• Estimating the Power Spectrum of the Cosmic Microwave Background
• On Foreground Removal from the Wilkinson Microwave Anisotropy Probe Data by an Internal Linear Combination Method: Limitations and Implications
• Is the Low CMB Quadrupole a Signature of Spatial Curvature?
• The Pseudo-$C_l$ method: Cosmic microwave background anisotropy power spectrum statistics for high precision cosmology
• Cosmic microwave background anisotropy power spectrum statistics for high precision cosmology
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Beam Profiles and Window Functions
• 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
• The Statistical Significance of the Low CMB Mulitipoles
• Gabor Transforms on the Sphere with Applications to CMB Power Spectrum Estimation
• Fast, exact CMB power spectrum estimation for a certain class of observational strategies
• Fast CMB Power Spectrum Estimation of Temperature and Polarisation with Gabor Transforms
• Harmonic analysis of cosmic microwave background data I: ring reductions and point-source catalogue
• Harmonic analysis of cosmic microwave background data II: from ring-sets to the sky
• CMB Power Spectrum Estimation via Hierarchical Decomposition

Cited by in corpus (169)

• 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
• Three-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Temperature Analysis
• Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters
• Planck 2013 results. XV. CMB power spectra and likelihood
• Planck 2018 results. IV. Diffuse component separation
• A unified pseudo-$C_{\ell }$ framework
• Detection of Gravitational Lensing in the Cosmic Microwave Background
• Bayesian physical reconstruction of initial conditions from large scale structure surveys
• 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
• Planck 2015 results. IX. Diffuse component separation: CMB maps
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Physical Bayesian modelling of the non-linear matter distribution: new insights into the Nearby Universe
• Bayesian analysis of sparse anisotropic universe models and application to the 5-yr WMAP data
• CMB component separation by parameter estimation
• CMB multipole measurements in the presence of foregrounds
• Past and present cosmic structure in the SDSS DR7 main sample
• Bayesian power-spectrum inference for Large Scale Structure data
• Bayesian analysis of an anisotropic universe model: systematics and polarization
• First Season QUIET Observations: Measurements of CMB Polarization Power Spectra at 43 GHz in the Multipole Range 25 <= ell <= 475
• Detection of Cosmological 21 cm Emission with the Canadian Hydrogen Intensity Mapping Experiment
• Sensitivity and foreground modelling for large-scale CMB B-mode polarization satellite missions
• 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
• Bayesian optimisation for likelihood-free cosmological inference
• Hierarchical Cosmic Shear Power Spectrum Inference
• Anomalous parity asymmetry of WMAP power spectrum data at low multpoles: is it cosmological or systematics?
• Bayesian Power Spectrum Analysis of the First-Year WMAP data
• Exploring Cosmic Origins with CORE: B-mode Component Separation
• 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
• Results from the Wilkinson Microwave Anisotropy Probe
• Needlet Detection of Features in WMAP CMB Sky and the Impact on Anisotropies and Hemispherical Asymmetries
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• A Case Against Spinning PAHs as the Source of the Anomalous Microwave Emission
• Cosmological Parameter Constraints as Derived from the Wilkinson Microwave Anisotropy Probe Data via Gibbs Sampling and the Blackwell-Rao Estimator
• CMB power spectrum estimation using wavelets
• Large-Angle Correlations in the Cosmic Microwave Background
• Methods for Bayesian power spectrum inference with galaxy surveys
• Reconstructing cosmic growth with kSZ observations in the era of Stage IV experiments
• Bayesian inference from photometric redshift surveys
• Multiscale Inference of Matter Fields and Baryon Acoustic Oscillations from the Ly-alpha Forest
• Approximate Bayesian Computation in Large Scale Structure: constraining the galaxy-halo connection
• 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
• Filling in CMB map missing data using constrained Gaussian realizations
• The Simons Observatory: Constraining inflationary gravitational waves with multi-tracer B-mode delensing
• 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
• A minimal power-spectrum-based moment expansion for CMB B-mode searches
• Impact of Galactic polarized emission on B-mode detection at low multipoles
• Cross-Correlation of Planck CMB Lensing with DESI-Like LRGs
• 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
• Performance forecasts for the primordial gravitational wave detection pipelines for AliCPT-1
• Bayesian inference of CMB gravitational lensing
• 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
• The integrated bispectrum as a test of CMB non-Gaussianity: detection power and limits on f_NL with WMAP data
• Matrix-free Large Scale Bayesian inference in cosmology
• Bayesian analysis of the low-resolution polarized 3-year WMAP sky maps
• Alignments of parity even/odd-only multipoles in CMB
• Search for anomalous alignments of structures in Planck data using Minkowski Tensors
• The VIMOS Public Extragalactic Redshift Survey (VIPERS): Reconstruction of the redshift-space galaxy density field
• Likelihood methods for CMB experiments
• Efficient Power Spectrum Estimation for High Resolution CMB Maps
• Methods for pixel domain correction of EB leakage
• How to make a clean separation between CMB E and B modes with proper foreground masking
• 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
• CMB likelihood approximation by a Gaussianized Blackwell-Rao estimator
• BeyondPlanck XIII. Intensity foreground sampling, degeneracies, and priors
• Multi-resolution Bayesian CMB component separation through Wiener-filtering with a pseudo-inverse preconditioner
• Tensor-to-scalar ratio forecasts for extended LiteBIRD frequency configurations
• Wiener filter reloaded: fast signal reconstruction without preconditioning
• The BINGO Project V: Further steps in Component Separation and Bispectrum Analysis
• Harmonic in-painting of CMB sky by constrained Gaussian realization
• LiteBIRD Science Goals and Forecasts: Improving Sensitivity to Inflationary Gravitational Waves with Multitracer Delensing
• The bispectra of galactic CMB foregrounds and their impact on primordial non-Gaussianity estimation
• A comparison of CMB Angular Power Spectrum Estimators at Large Scales: the TT case
• Wiener filtering and pure E/B decomposition of CMB maps with anisotropic correlated noise
• Symmetry and anti-symmetry of the CMB anisotropy pattern
• On the impact of large angle CMB polarization data on cosmological parameters
• BeyondPlanck X. Bandpass and beam leakage corrections
• Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis
• BeyondPlanck XII. Cosmological parameter constraints with end-to-end error propagation
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• Bayesian power-spectrum inference with foreground and target contamination treatment
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• A multi-level solver for Gaussian constrained CMB realizations
• Hybrid-basis inference for large-scale galaxy clustering: combining spherical and Cartesian Fourier analyses
• Bayesian angular power spectrum analysis of interferometric data
• Data analysis methods for the cosmic microwave background
• Analysis of WMAP 7-year Temperature Data: Astrophysics of the Galactic Haze
• Joint Bayesian analysis of large angular scale CMB temperature anomalies
• Cosmo++: An Object-Oriented C++ Library for Cosmology
• Effects of Destriping Errors on CMB Polarisation Power Spectra and and Pixel Noise Covariances
• The nature of non-Gaussianity and statistical isotropy of the 408 MHz Haslam synchrotron map
• Hierarchical Bayesian CMB Component Separation with the No-U-Turn Sampler
• Reassessment of the dipole in the distribution of quasars on the sky
• LiteBIRD Science Goals and Forecasts: A full-sky measurement of gravitational lensing of the CMB
• Testing for spectral index variations in polarised CMB foregrounds
• A Novel CMB Component Separation Method: Hierarchical Generalized Morphological Component Analysis
• Needlet thresholding methods in component separation
• Marginal distributions for cosmic variance limited CMB polarization data
• Recovering the Tidal Field in the Projected Galaxy Distribution
• Bayesian analysis of white noise levels in the 5-year WMAP data
• Parametric tension between even and odd multipole data of WMAP power spectrum: unaccounted contamination or missing parameters?
• Testing the ABS method with the simulated Planck temperature maps
• Cosmological Parameters from CMB Maps without Likelihood Approximation
• Optimal and fast E/B separation with a dual messenger field
• Impact of Galactic dust non-Gaussianity on searches for B-modes from inflation
• Large-Scale Polarized Foreground Component Separation for Planck
• Optimized Large-Scale CMB Likelihood And Quadratic Maximum Likelihood Power Spectrum Estimation
• Optimal analysis of azimuthal features in the CMB
• Joint Bayesian estimation of tensor and lensing B-modes in the power spectrum of CMB polarization data
• An Unbiased Estimator of the Full-sky CMB Angular Power Spectrum at Large Scales using Neural Networks
• Estimation of Full Sky Power Spectrum between Intermediate to Large Angular Scales from Partial Sky CMB Anisotropies using Artificial Neural Network
• BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• Cosmic dipole tensions: confronting the cosmic microwave background with infrared and radio populations of cosmological sources
• Bayesian foreground analysis with CMB data
• Cosmological Parameters from a re-analysis of the WMAP-7 low resolution maps
• 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
• QUIJOTE scientific results -- VI. The Haze as seen by QUIJOTE
• Mitigating Bias in CMB B-modes from Foreground Cleaning Using a Moment Expansion
• Preconditioner-free Wiener filtering with a dense noise matrix
• A Bayesian estimate of the CMB-large-scale structure cross-correlation
• Enhancing CMB map reconstruction and power spectrum estimation with convolutional neural networks
• An Accurate Reconstruction of CMB E Mode Signal over Large Angular Scales using Prior Information of CMB Covariance Matrix in ILC Algorithm
• A CMB Gibbs sampler for localized secondary anisotropies
• A Monte Carlo comparison between template-based and Wiener-filter CMB dipole estimators
• A new look at the integrated radio/microwave continuum spectrum of Galactic supernova remnant IC443
• Bayesian analysis of cosmic structures
• Examining statistical isotropy of CMB low multipoles from Planck PR4 data
• CMB likelihood approximation for banded probability distributions
• Statistical properties of Galactic synchrotron temperature and polarization maps -- a multi-frequency comparison
• The two-and three-point correlation functions of the polarized five-year WMAP sky maps
• Improved Gibbs samplers for Cosmic Microwave Background power spectrum estimation
• Limits on Polarized Dust Spectral Index Variations for CMB Foreground Analysis
• ParamANN: A Neural Network to Estimate Cosmological Parameters for $\Lambda$CDM Universe Using Hubble Measurements
• Application of Genetic Algorithm to Estimate the Large Angular Scale Features of Cosmic Microwave Background
• Enhancing Primordial B-mode Detection: Comprehensive Delensing Pipelines for Improved Sensitivity to $r$
• An Optimal and Model-Independent Measurement of the Intracluster Pressure Profile I: Methodology and First Applications
• On the Impact of Random Residual Calibration Error on the Gibbs ILC Estimates over Large Angular Scales
• Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations
• Halo detection via large-scale Bayesian inference
• SymPix: A spherical grid for efficient sampling of rotationally invariant operators
• The CMB angular power spectrum via component separation: a study on Planck data
• Shrinking the Quadratic Estimator
• Making Maps from Planck LFI 30GHz Data with Asymmetric Beams and Cooler Noise
• LiteBIRD science goals and forecasts: improved full-sky reconstruction of the gravitational lensing potential through the combination of Planck and LiteBIRD data
• Pixel domain multi-resolution minimum variance painting of CMB maps
• Mitigating Cosmic Microwave Background Shadow Degradation of Tensor-to-scalar Ratio Measurements through Map-based Studies
• Non-Gaussianity in Cosmology: from Inflation to the CMB
• B-mode polarization forecasts for GreenPol
• Do We Have Sufficient Knowledge of the Galactic Foreground Emission in Cosmic Microwave Background Science?
• Accurate and efficient likelihood modeling for large-scale CMB data