Cosmological Analysis with Calibrated Neural Quantile Estimation and Approximate Simulators

Abstract

A major challenge in extracting information from current and upcoming surveys of cosmological Large-Scale Structure (LSS) is the limited availability of computationally expensive high-fidelity simulations. We introduce calibrated Neural Quantile Estimation (NQE), a new Simulation-Based Inference (SBI) method that leverages a large number of approximate simulations for training and a small number of high-fidelity simulations for calibration. This approach guarantees an unbiased posterior regardless of approximate simulation accuracy, while achieving near-optimal constraining power when the approximate simulations are reasonably accurate. As a proof of concept, we demonstrate that cosmological parameters can be inferred at field level from projected 2-dim dark matter density maps up to $k_{\max }\sim 1.5h$/Mpc at $z=0$ by training on $\sim 10^4$ Particle-Mesh (PM) simulations with transfer function correction and calibrating with $\sim 10^2$ Particle-Particle (PP) simulations. The calibrated posteriors closely match those obtained by directly training on $\sim 10^4$ expensive PP simulations, but at a fraction of the computational cost. Our method offers a practical and scalable framework for SBI of cosmological LSS, enabling precise inference across vast volumes and down to small scales.

5+5 pages, 5+4 figures, published in PRL

References in corpus (22)

• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• Euclid. I. Overview of the Euclid mission
• Overview of the Instrumentation for the Dark Energy Spectroscopic Instrument
• The Arepo public code release
• Simulating cosmic structure formation with the GADGET-4 code
• The Quijote simulations
• Power Spectrum Super-Sample Covariance
• Learning to Predict the Cosmological Structure Formation
• FastPM: a new scheme for fast simulations of dark matter and halos
• A new approach to observational cosmology using the scattering transform
• The $\text{Abacus}$ Cosmological $N$-body Code
• Towards an Optimal Estimation of Cosmological Parameters with the Wavelet Scattering Transform
• How much information can be extracted from galaxy clustering at the field level?
• Robust Field-level Likelihood-free Inference with Galaxies
• A gradient based method for modeling baryons and matter in halos of fast simulations
• A Parameter-Masked Mock Data Challenge for Beyond-Two-Point Galaxy Clustering Statistics
• A comparative study of cosmological constraints from weak lensing using Convolutional Neural Networks
• Super-sample covariance of the power spectrum, bispectrum, halos, voids, and their cross covariances
• A field-level emulator for modeling baryonic effects across hydrodynamic simulations
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• Inferring Cosmological Parameters on SDSS via Domain-Generalized Neural Networks and Lightcone Simulations