The DESI $N$-body Simulation Project I: Testing the Robustness of Simulations for the DESI Dark Time Survey

Abstract

Analysis of large galaxy surveys requires confidence in the robustness of numerical simulation methods. The simulations are used to construct mock galaxy catalogs to validate data analysis pipelines and identify potential systematics. We compare three $N$-body simulation codes, ABACUS, GADGET, and SWIFT, to investigate the regimes in which their results agree. We run $N$-body simulations at three different mass resolutions, $6.25\times 10^8$, $2.11\times 10^9$, and $5.00\times 10^9{h}^{-1}$M${}_{\odot }$, matching phases to reduce the noise within the comparisons. We find systematic errors in the halo clustering between different codes are smaller than the DESI statistical error for $s>20h^{-1}$Mpc in the correlation function in redshift space. Through the resolution comparison we find that simulations run with a mass resolution of $2.1\times 10^9{h}^{-1}$M${}_{\odot }$ are sufficiently converged for systematic effects in the halo clustering to be smaller than the DESI statistical error at scales larger than $20h^{-1}$Mpc. These findings show that the simulations are robust for extracting cosmological information from large scales which is the key goal of the DESI survey. Comparing matter power spectra, we find the codes agree to within 1% for $k\le 10h$Mpc${}^{-1}$. We also run a comparison of three initial condition generation codes and find good agreement. In addition, we include a quasi-$N$-body code, FastPM, since we plan use it for certain DESI analyses. The impact of the halo definition and galaxy-halo relation will be presented in a follow up study.

arXiv comment18 pages, 17 figures. Submitted to MNRAS

References in corpus (37)

• A Universal Density Profile from Hierarchical Clustering
• The cosmological simulation code GADGET-2
• Simulating the joint evolution of quasars, galaxies and their large-scale distribution
• Statistical Properties of X-ray Clusters: Analytic and Numerical Comparisons
• An analytic model for the spatial clustering of dark matter haloes
• The Rockstar Phase-Space Temporal Halo Finder and the Velocity Offsets of Cluster Cores
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Using the Evolution of Clusters to Constrain Omega
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Halo occupation numbers and galaxy bias
• The Connection between Galaxies and their Dark Matter Halos
• Transients from Initial Conditions in Cosmological Simulations
• Connecting Galaxies, Halos, and Star Formation Rates Across Cosmic Time
• The Nature of Galaxy Bias and Clustering
• Scaling relations for galaxy clusters in the Millennium-XXL simulation
• The Kilo-Degree Survey
• Statistical Analysis of Galaxy Surveys - I. Robust error estimation for 2-point clustering statistics
• The DESI Experiment, a whitepaper for Snowmass 2013
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• AbacusSummit: A Massive Set of High-Accuracy, High-Resolution $N$-Body Simulations
• FastPM: a new scheme for fast simulations of dark matter and halos
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Modeling the clustering and halo occupation distribution of BOSS-CMASS galaxies in the Final Data Release
• Matter power spectrum and the challenge of percent accuracy
• Modified Gravity N-body Code Comparison Project
• The Aemulus Project I: Numerical Simulations for Precision Cosmology
• The $\text{Abacus}$ Cosmological $N$-body Code
• A new way of setting the phases for cosmological multi-scale Gaussian initial conditions
• Second-order Lagrangian perturbation theory initial conditions for resimulations
• A High-Fidelity Realization of the Euclid Code Comparison $N$-body Simulation with Abacus
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Linear perturbative theory of the discrete cosmological N-body problem
• SWIFT: Using task-based parallelism, fully asynchronous communication, and graph partition-based domain decomposition for strong scaling on more than 100,000 cores
• Dark Energy Survey Year 3 results: cosmology from combined galaxy clustering and lensing -- validation on cosmological simulations
• Ringing the initial Universe: the response of overdensity and transformed-density power spectra to initial spikes
• Reducing the Variance of Redshift Space Distortion Measurements from Mock Galaxy Catalogues with Different Lines of Sight
• Good and Proper: Self-similarity of N-body Simulations with Proper Force Softening
• SWIFT: task-based hydrodynamics and gravity for cosmological simulations

Cited by in corpus (13)

• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• The MillenniumTNG Project: High-precision predictions for matter clustering and halo statistics
• Baryon Acoustic Oscillation Theory and Modelling Systematics for the DESI 2024 results
• Constraints on $S_8$ from a full-scale and full-shape analysis of redshift-space clustering and galaxy-galaxy lensing in BOSS
• Higher-order initial conditions with massive neutrinos
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• Analysis of Unified Galaxy Power Spectrum Multipole Measurements
• A comparison between Shapefit compression and Full-Modelling method with PyBird for DESI 2024 and beyond
• Galaxy clustering from the bottom up: A Streaming Model emulator I
• The HalfDome Multi-Survey Cosmological Simulations: N-body Simulations
• DESI Mock Challenge: Constructing DESI galaxy catalogues based on FastPM simulations
• Field-Level Comparison and Robustness Analysis of Cosmological N-body Simulations
• A lightcone catalogue from the Millennium-XXL simulation: improved spatial interpolation and colour distributions for the DESI BGS