Simulation-Based Inference Benchmark for Weak Lensing Cosmology

Abstract

Standard cosmological analysis, which relies on two-point statistics, fails to extract the full information of the data. This limits our ability to constrain with precision cosmological parameters. Thus, recent years have seen a paradigm shift from analytical likelihood-based to simulation-based inference. However, such methods require a large number of costly simulations. We focus on full-field inference, considered the optimal form of inference. Our objective is to benchmark several ways of conducting full-field inference to gain insight into the number of simulations required for each method. We make a distinction between explicit and implicit full-field inference. Moreover, as it is crucial for explicit full-field inference to use a differentiable forward model, we aim to discuss the advantages of having this property for the implicit approach. We use the sbi_lens package which provides a fast and differentiable log-normal forward model. This forward model enables us to compare explicit and implicit full-field inference with and without gradient. The former is achieved by sampling the forward model through the No U-Turns sampler. The latter starts by compressing the data into sufficient statistics and uses the Neural Likelihood Estimation algorithm and the one augmented with gradient. We perform a full-field analysis on LSST Y10 like weak lensing simulated mass maps. We show that explicit and implicit full-field inference yield consistent constraints. Explicit inference requires 630 000 simulations with our particular sampler corresponding to 400 independent samples. Implicit inference requires a maximum of 101 000 simulations split into 100 000 simulations to build sufficient statistics (this number is not fine tuned) and 1 000 simulations to perform inference. Additionally, we show that our way of exploiting the gradients does not significantly help implicit inference.

20 pages, 16 figures, accepted to A&A

References in corpus (39)

• Array Programming with NumPy
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• The Hyper Suprime-Cam SSP Survey: Overview and Survey Design
• Cosmology with cosmic shear observations: a review
• Non-linear regression models for Approximate Bayesian Computation
• The Kilo-Degree Survey
• CFHTLenS: The Canada-France-Hawaii Telescope Lensing Survey - Imaging Data and Catalogue Products
• CFHTLenS: Cosmological constraints from a combination of cosmic shear two-point and three-point correlations
• Improving lognormal models for cosmological fields
• Cosmology Constraints from the Weak Lensing Peak Counts and the Power Spectrum in CFHTLenS
• Mining gold from implicit models to improve likelihood-free inference
• Cosmology constraints from shear peak statistics in Dark Energy Survey Science Verification data
• A new approach to observational cosmology using the scattering transform
• KiDS-450: Cosmological Constraints from Weak Lensing Peak Statistics - II: Inference from Shear Peaks using N-body Simulations
• Non-Gaussian information from weak lensing data via deep learning
• Likelihood-free inference with neural compression of DES SV weak lensing map statistics
• KiDS-450: Cosmological Constraints from Weak Lensing Peak Statistics-I: Inference from Analytical Prediction of High Signal-to-Noise Ratio Convergence Peaks
• Automatic physical inference with information maximising neural networks
• Cosmological Constraints From Weak Lensing Peak Statistics With CFHT Stripe 82 Survey
• Weak lensing cosmology with convolutional neural networks on noisy data
• Weak lensing from space: first cosmological constraints from three-point shear statistics
• Cosmological constraints from noisy convergence maps through deep learning
• Hierarchical Cosmic Shear Power Spectrum Inference
• A new model to predict weak-lensing peak counts II. Parameter constraint strategies
• Testing the lognormality of the galaxy and weak lensing convergence distributions from Dark Energy Survey maps
• Cosmological parameters, shear maps and power spectra from CFHTLenS using Bayesian hierarchical inference
• Primordial non-Gaussianity without tails -- how to measure fNL with the bulk of the density PDF
• JAX-COSMO: An End-to-End Differentiable and GPU Accelerated Cosmology Library
• A Full $w$CDM Analysis of KiDS-1000 Weak Lensing Maps using Deep Learning
• Cosmological constraints with weak lensing peak counts and second-order statistics in a large-field survey
• The integrated 3-point correlation function of cosmic shear
• Bayesian forward modelling of cosmic shear data
• MUSE: Marginal Unbiased Score Expansion and Application to CMB Lensing
• KiDS-SBI: Simulation-based inference analysis of KiDS-1000 cosmic shear
• Hierarchical probabilistic inference of cosmic shear
• Bayesian weak lensing tomography: Reconstructing the 3D large-scale distribution of matter with a lognormal prior
• A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey
• Data Compression and Inference in Cosmology with Self-Supervised Machine Learning
• Accurate field-level weak lensing inference for precision cosmology

Cited by in corpus (6)

• Observable Optimization for Precision Theory: Machine Learning Energy Correlators
• Quantifying Weighted Morphological Content of Large-Scale Structures via Simulation-Based Inference
• How to embed any likelihood into SBI: Application to Planck + Stage IV galaxy surveys and Dynamical Dark Energy
• GalSBI-SPS: a stellar population synthesis-based galaxy population model for cosmology and galaxy evolution applications
• Benchmarking field-level cosmological inference from galaxy redshift surveys
• Diagnosing Systematic Effects Using the Inferred Initial Power Spectrum