Bayesian weak lensing tomography: Reconstructing the 3D large-scale distribution of matter with a lognormal prior

Abstract

We present a Bayesian reconstruction algorithm that infers the three-dimensional large-scale matter distribution from the weak gravitational lensing effects measured in the image shapes of galaxies. The algorithm is designed to also work with non-Gaussian posterior distributions which arise, for example, from a non-Gaussian prior distribution. In this work, we use a lognormal prior and compare the reconstruction results to a Gaussian prior in a suite of increasingly realistic tests on mock data. We find that in cases of high noise levels (i.e. for low source galaxy densities and/or high shape measurement uncertainties), both normal and lognormal priors lead to reconstructions of comparable quality, but with the lognormal reconstruction being prone to mass-sheet degeneracy. In the low-noise regime and on small scales, the lognormal model produces better reconstructions than the normal model: The lognormal model 1) enforces non-negative densities, while negative densities are present when a normal prior is employed, 2) better traces the extremal values and the skewness of the true underlying distribution, and 3) yields a higher pixel-wise correlation between the reconstruction and the true density.

23 pages, 12 figures; updated to match version accepted for publication in PRD

References in corpus (5)

• COSMOS: 3D weak lensing and the growth of structure
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys -I. Methodology and a Fast Fitting Algorithm
• Bayesian Galaxy Shape Measurement for Weak Lensing Surveys -II. Application to Simulations
• Measuring primordial non-Gaussianity with weak-lensing surveys
• Hierarchical probabilistic inference of cosmic shear

Cited by in corpus (15)

• Two- and three-dimensional wide-field weak lensing mass maps from the Hyper Suprime-Cam Subaru Strategic Program S16A data
• Deep learning dark matter map reconstructions from DES SV weak lensing data
• Lifting weak lensing degeneracies with a field-based likelihood
• Improving Weak Lensing Mass Map Reconstructions using Gaussian and Sparsity Priors: Application to DES SV
• Bayesian forward modelling of cosmic shear data
• KiDS-SBI: Simulation-based inference analysis of KiDS-1000 cosmic shear
• GLASS: Generator for Large Scale Structure
• Efficient Optimal Reconstruction of Linear Fields and Band-powers from Cosmological Data
• KaRMMa -- Kappa Reconstruction for Mass Mapping
• Accurate field-level weak lensing inference for precision cosmology
• Simulation-Based Inference Benchmark for Weak Lensing Cosmology
• Denoising, deconvolving and decomposing multi-domain photon observations- The D4PO algorithm
• Posterior sampling for inverse imaging problems on the sphere in seismology and cosmology
• Almanac: MCMC-based signal extraction of power spectra and maps on the sphere
• Generative modeling of convergence maps based on predicted one-point statistics