Almanac: MCMC-based signal extraction of power spectra and maps on the sphere

Abstract

Inference in cosmology often starts with noisy observations of random fields on the celestial sphere, such as maps of the microwave background radiation, continuous maps of cosmic structure in different wavelengths, or maps of point tracers of the cosmological fields. Almanac uses Hamiltonian Monte Carlo sampling to infer the underlying all-sky noiseless maps of cosmic structures, in multiple redshift bins, together with their auto- and cross-power spectra. It can sample many millions of parameters, handling the highly variable signal-to-noise of typical cosmological signals, and it provides science-ready posterior data products. In the case of spin-weight 2 fields, Almanac infers $E$- and $B$-mode power spectra and parity-violating $EB$ power, and, by sampling the full posteriors rather than point estimates, it avoids the problem of $EB$-leakage. For theories with no $B$-mode signal, inferred non-zero $B$-mode power may be a useful diagnostic of systematic errors or an indication of new physics. Almanac's aim is to characterise the statistical properties of the maps, with outputs that are completely independent of the cosmological model, beyond an assumption of statistical isotropy. Inference of parameters of any particular cosmological model follows in a separate analysis stage. We demonstrate our signal extraction on a CMB-like experiment.

27 pages, 18 figures. v2 accepted for publication by The Open Journal of Astrophysics with minor changes. v3 no changes, missing acknowledgement added

References in corpus (18)

• Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Angular Power Spectra
• Likelihood-free inference with neural compression of DES SV weak lensing map statistics
• Cosmic Shear Cosmology Beyond 2-Point Statistics: A Combined Peak Count and Correlation Function Analysis of DES-Y1
• Estimation of Polarized Power Spectra by Gibbs sampling
• Fast optimal CMB power spectrum estimation with Hamiltonian sampling
• Efficient Cosmological Parameter Estimation with Hamiltonian Monte Carlo
• KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography
• Bayesian forward modelling of cosmic shear data
• CMB E/B decomposition of incomplete sky: a pixel space approach
• Redshift evolution of cosmic birefringence in CMB anisotropies
• On the accuracy and precision of correlation functions and field-level inference in cosmology
• BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation
• Almanac: Weak Lensing power spectra and map inference on the masked sphere
• Field-level inference of cosmic shear with intrinsic alignments and baryons
• The Information Content of Projected Galaxy Fields
• Euclid: Covariance of weak lensing pseudo-$C_{\ell }$ estimates. Calculation, comparison to simulations, and dependence on survey geometry
• KaRMMa 2.0 -- Kappa Reconstruction for Mass Mapping
• Map-based cosmology inference with lognormal cosmic shear maps

Cited by in corpus (4)

• Euclid preparation. LIX. Angular power spectra from discrete observations
• Stage IV baryonic feedback correction for non-Gaussianity inference
• Hybrid summary statistics: neural weak lensing inference beyond the power spectrum
• Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum