Abstract:N-body simulations are the most powerful method to study the non-linear evolution of large-scale structure. However, they require large amounts of computational resources, making unfeasible their direct adoption in scenarios that require broad explorations of parameter spaces. In this work, we show that it is possible to perform fast dark matter density field emulations with competitive accuracy using simple machine-learning approaches. We build an emulator based on dimensionality reduction and machine learning regression combining simple Principal Component Analysis and supervised learning methods. For the estimations with a single free parameter, we train on the dark matter density parameter, $\Omega_m$, while for emulations with two free parameters, we train on a range of $\Omega_m$ and redshift. The method first adopts a projection of a grid of simulations on a given basis; then, a machine learning regression is trained on this projected grid. Finally, new density cubes for different cosmological parameters can be estimated without relying directly on new N-body simulations by predicting and de-projecting the basis coefficients. We show that the proposed emulator can generate density cubes at non-linear cosmological scales with density distributions within a few percent compared to the corresponding N-body simulations. The method enables gains of three orders of magnitude in CPU run times compared to performing a full N-body simulation while reproducing the power spectrum and bispectrum within $\sim 1\%$ and $\sim 3\%$, respectively, for the single free parameter emulation and $\sim 5\%$ and $\sim 15\%$ for two free parameters. This can significantly accelerate the generation of density cubes for a wide variety of cosmological models, opening the doors to previously unfeasible applications, such as parameter and model inferences at full survey scales as the ESA/NASA Euclid mission.

Mean-Field Simulation-Based Inference for Cosmological Initial Conditions

Bayesian Inference of Initial Conditions from Non-Linear Cosmic Structures using Field-Level Emulators

Field-level simulation-based inference of galaxy clustering with convolutional neural networks

COSMIC BIRTH: Efficient Bayesian Inference of the Evolving Cosmic Web from Galaxy Surveys

SimBIG: Field-level Simulation-Based Inference of Galaxy Clustering

Reconstructing galaxy star formation histories from COSMOS2020 photometry using simulation-based inference

Simulation-Based Inference Benchmark for LSST Weak Lensing Cosmology

Posterior Sampling of the Initial Conditions of the Universe from Non-linear Large Scale Structures using Score-Based Generative Models

Fast emulation of cosmological density fields based on dimensionality reduction and supervised machine-learning

Field-Level Inference with Microcanonical Langevin Monte Carlo

Exploring the likelihood of the 21-cm power spectrum with simulation-based inference

A simulation-based inference pipeline for cosmic shear with the Kilo-Degree Survey

Cosmological Analysis with Calibrated Neural Quantile Estimation and Approximate Simulators

Cosmological inference from Bayesian forward modelling of deep galaxy redshift surveys

Cosmological Field Emulation and Parameter Inference with Diffusion Models

Simulation-based inference of deep fields: galaxy population model and redshift distributions

Cosmological Inflation in N-Dimensional Gaussian Random Fields with Algorithmic Data Compression

Constraining Cosmology with Simulation-based inference and Optical Galaxy Cluster Abundance

Tracking the Multifield Dynamics with Cosmological Data: A Monte Carlo approach

Inpainting Galaxy Counts onto N-Body Simulations over Multiple Cosmologies and Astrophysics

The future of cosmological likelihood-based inference: accelerated high-dimensional parameter estimation and model comparison