The DESI $N$-body Simulation Project -- II. Suppressing sample variance with fast simulations

Abstract

Dark Energy Spectroscopic Instrument (DESI) will construct a large and precise three-dimensional map of our Universe. The survey effective volume reaches $\sim20\Gpchcube$. It is a great challenge to prepare high-resolution simulations with a much larger volume for validating the DESI analysis pipelines. \textsc{AbacusSummit} is a suite of high-resolution dark-matter-only simulations designed for this purpose, with $200\Gpchcube$ (10 times DESI volume) for the base cosmology. However, further efforts need to be done to provide a more precise analysis of the data and to cover also other cosmologies. Recently, the CARPool method was proposed to use paired accurate and approximate simulations to achieve high statistical precision with a limited number of high-resolution simulations. Relying on this technique, we propose to use fast quasi-$N$-body solvers combined with accurate simulations to produce accurate summary statistics. This enables us to obtain 100 times smaller variance than the expected DESI statistical variance at the scales we are interested in, e.g. $k < 0.3\hMpc$ for the halo power spectrum. In addition, it can significantly suppress the sample variance of the halo bispectrum. We further generalize the method for other cosmologies with only one realization in \textsc{AbacusSummit} suite to extend the effective volume $\sim 20$ times. In summary, our proposed strategy of combining high-fidelity simulations with fast approximate gravity solvers and a series of variance suppression techniques sets the path for a robust cosmological analysis of galaxy survey data.

Matched version accepted by MNRAS, should be clearer

References in corpus (12)

• Dynamics of dark energy
• Massive neutrinos and cosmology
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• Corrfunc --- A Suite of Blazing Fast Correlation Functions on the CPU
• The $\text{Abacus}$ Cosmological $N$-body Code
• Constraining $M_{\nu }$ with the Bispectrum II: The Total Information Content of the Galaxy Bispectrum
• \textsc{CompaSO}: A new halo finder for competitive assignment to spherical overdensities
• nIFTy Cosmology: Galaxy/halo mock catalogue comparison project on clustering statistics
• The host dark matter halos of [OII] emitters at 0.5< z< 1.5
• Constructing high-fidelity halo merger trees in AbacusSummit
• Statistics of Fourier Modes in Non-Gaussian Fields
• Numerical investigation of non-Gaussianities in the phase and modulus of density Fourier modes

Cited by in corpus (14)

• Baryon Acoustic Oscillation Theory and Modelling Systematics for the DESI 2024 results
• Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation
• The MillenniumTNG Project: The large-scale clustering of galaxies
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• PNG-UNITsims: Halo clustering response to primordial non-Gaussianities as a function of mass
• Fast and realistic large-scale structure from machine-learning-augmented random field simulations
• Mitigating the noise of DESI mocks using analytic control variates
• Bayesian Control Variates for optimal covariance estimation with pairs of simulations and surrogates
• Quijote-PNG: Quasi-maximum likelihood estimation of Primordial Non-Gaussianity in the non-linear halo density field
• The HalfDome Multi-Survey Cosmological Simulations: N-body Simulations
• Validating galaxy clustering models with Fixed & Paired and Matched-ICs simulations: application to Primordial Non-Gaussianities
• Fast simulation mapping: from standard to modified gravity cosmologies using the bias assignment method
• Lognormal semi-numerical simulations of the Lyman-$\alpha$ forest: comparison with full hydrodynamic simulations
• Simulation budgeting for hybrid effective field theories