Accurate initial conditions for cosmological N-body simulations: Minimizing truncation and discreteness errors

Abstract

Inaccuracies in the initial conditions for cosmological N-body simulations could easily be the largest source of systematic error in predicting the non-linear large-scale structure. From the theory side, initial conditions are usually provided by using low-order truncations of the displacement field from Lagrangian perturbation theory, with the first and second-order approximations being the most common ones. Here we investigate the improvement brought by using initial conditions based on third-order Lagrangian perturbation theory (3LPT). We show that with 3LPT, truncation errors are vastly suppressed, thereby opening the portal to initializing simulations accurately as late as z=12 (for the resolution we consider). We analyse the competing effects of perturbative truncation and particle discreteness on various summary statistics. Discreteness errors are essentially decaying modes and thus get strongly amplified for earlier initialization times. We show that late starting times with 3LPT provide the most accurate configuration, which we find to coincide with the continuum fluid limit within 1 per cent for the power- and bispectrum at z=0 up to the particle Nyquist wave number of our simulations (k ~ 3h/Mpc). In conclusion, to suppress non-fluid artefacts, we recommend initializing simulations as late as possible with 3LPT. We make our 3LPT initial condition generator publicly available.

arXiv comment22 pages, 13 figures, minor revisions (matches published version), code available from https://bitbucket.org/ohahn/monofonic

References in corpus (90)

• The cosmological simulation code GADGET-2
• Efficient Computation of CMB anisotropies in closed FRW models
• Populating a cluster of galaxies - I. Results at z=0
• The Cosmic Linear Anisotropy Solving System (CLASS) II: Approximation schemes
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• Toward a halo mass function for precision cosmology: the limits of universality
• Mass function of dark matter halos
• Euclid Definition Study Report
• LSST Science Book, Version 2.0
• Multi-scale initial conditions for cosmological simulations
• Transients from Initial Conditions in Cosmological Simulations
• First Data Release of the Hyper Suprime-Cam Subaru Strategic Program
• Evolution of structure in cold dark matter universes
• Precision Determination of the Mass Function of Dark Matter Halos
• Scaling relations for galaxy clusters in the Millennium-XXL simulation
• Fundamental Physics with the Square Kilometre Array
• Transients from Initial Conditions: A Perturbative Analysis
• The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
• nbodykit: an open-source, massively parallel toolkit for large-scale structure
• Correcting for the alias effect when measuring the power spectrum using FFT
• PINOCCHIO: pinpointing orbit-crossing collapsed hierarchical objects in a linear density field
• Accurate Estimators of Correlation Functions in Fourier Space
• Dark Quest. I. Fast and Accurate Emulation of Halo Clustering Statistics and Its Application to Galaxy Clustering
• Generation of Vorticity and Velocity Dispersion by Orbit Crossing
• The BACCO Simulation Project: Exploiting the full power of large-scale structure for cosmology
• Cosmological $N$-body simulations with suppressed variance
• Euclid preparation: II. The EuclidEmulator -- A tool to compute the cosmology dependence of the nonlinear matter power spectrum
• The Bispectrum: From Theory to Observations
• Matter power spectrum and the challenge of percent accuracy
• The Bispectrum in the Effective Field Theory of Large Scale Structure
• Reconstruction of the early Universe as a convex optimization problem
• Lagrangian theory of gravitational instability of Friedman-Lemaitre cosmologies - generic third-order model for non-linear clustering
• The Mira-Titan Universe: Precision Predictions for Dark Energy Surveys
• Connection between Newtonian simulations and general relativity
• Second-order Lagrangian perturbation theory initial conditions for resimulations
• The Properties of Cosmic Velocity Fields
• The Outer Rim Simulation: A Path to Many-Core Supercomputers
• A new approach to simulating collisionless dark matter fluids
• Lagrangian Dynamics in Non-Flat Universes and Non-Linear Gravitational Evolution
• Numerical convergence of simulations of galaxy formation: the abundance and internal structure of cold dark matter haloes
• Recursive Solutions of Lagrangian Perturbation Theory
• An adaptively refined phase-space element method for cosmological simulations and collisionless dynamics
• Lagrangian perturbations and the matter bispectrum I: fourth-order model for non-linear clustering
• Non-Linear Approximations to Gravitational Instability: A Comparison in the Quasi-Linear Regime
• Time-analyticity of Lagrangian particle trajectories in ideal fluid flow
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Improving fast generation of halo catalogs with higher-order Lagrangian perturbation theory
• Higher-order initial conditions for mixed baryon-CDM simulations
• The recursion relation in Lagrangian perturbation theory
• Toward an accurate mass function for precision cosmology
• ColDICE: a parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation
• How closely do baryons follow dark matter on large scales?
• Linear perturbative theory of the discrete cosmological N-body problem
• GridSPT: Grid-based calculation for perturbation theory of large-scale structure
• Baryonic effects on the matter bispectrum
• Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions
• General relativistic weak-field limit and Newtonian N-body simulations
• Fast Estimation of Gravitational and Primordial Bispectra in Large Scale Structures
• Post-collapse perturbation theory in 1D cosmology -- beyond shell-crossing
• How smooth are particle trajectories in a $\Lambda$CDM Universe?
• Effects of the initial conditions on cosmological $N$-body simulations
• Gravitational evolution of a perturbed lattice and its fluid limit
• Large-scale structure perturbation theory without losing stream crossing
• Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM
• Vlasov-Poisson in 1D for initially cold systems: post-collapse Lagrangian perturbation theory
• Simulating the Complexity of the Dark Matter Sheet I: Numerical Algorithms
• Shell-crossing in quasi-one-dimensional flow
• Lagrangian cosmological perturbation theory at shell-crossing
• Structure formation beyond shell-crossing: nonperturbative expansions and late-time attractors
• More accurate simulations with separate initial conditions for baryons and dark matter
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Large-Scale Velocity Dispersion and the Cosmic Web
• Semiclassical path to cosmic large-scale structure
• Shell-crossing in a $\Lambda$CDM Universe
• On the starting redshift for cosmological simulations: Focusing on halo properties
• The Cauchy-Lagrangian method for numerical analysis of Euler flow
• General relativistic 'screening' in cosmological simulations
• Accurate initial conditions in mixed Dark Matter--Baryon simulations
• An alternative method to generate pre-initial conditions for cosmological $N$-body simulations
• An alternative to grids and glasses: Quaquaversal pre-initial conditions for N-body simulations
• Unveiling the singular dynamics in the cosmic large-scale structure
• Toward Accurate Modeling of the Nonlinear Matter Bispectrum: Standard Perturbation Theory and Transients from Initial Conditions
• Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations
• Efficient parallel algorithm for estimating higher-order polyspectra
• Quasi-spherical collapse of matter in $\Lambda$CDM
• Particle linear theory on a self-gravitating perturbed cubic Bravais lattice
• CUBE -- Towards an Optimal Scaling of Cosmological N-body Simulations
• Advancing the matter bispectrum estimation of large-scale structure: a comparison of dark matter codes
• Transients from initial conditions based on Lagrangian perturbation theory in $N$-body simulations II: the effect of the transverse mode
• Projection of the gravitational dynamics on a subspace of probability distributions: curl-free Gaussian ansatz

Cited by in corpus (64)

• The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys
• Large-scale dark matter simulations
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• The cosmology dependence of galaxy clustering and lensing from a hybrid $N$-body-perturbation theory model
• Higher-order initial conditions for mixed baryon-CDM simulations
• The BACCO simulation project: biased tracers in real space
• An $n$-th order Lagrangian Forward Model for Large-Scale Structure
• Non-universality of the mass function: dependence on the growth rate and power spectrum shape
• Aemulus $\nu$: Precise Predictions for Matter and Biased Tracer Power Spectra in the Presence of Neutrinos
• Higher-order initial conditions with massive neutrinos
• Measuring the Tidal Response of Structure Formation: Anisotropic Separate Universe Simulations using TreePM
• Sigma-Eight at the Percent Level: The EFT Likelihood in Real Space
• Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation
• SPICE: the connection between cosmic reionisation and stellar feedback in the first galaxies
• Cosmological perturbations for two cold fluids in $\Lambda$CDM
• Cosmology inference at the field level from biased tracers in redshift-space
• Precision Redshift-Space Galaxy Power Spectra using Zel'dovich Control Variates
• Shell-crossing in a $\Lambda$CDM Universe
• A unified model for the clustering of quasars and galaxies at $z\approx 6$
• Bayesian Inference of Initial Conditions from Non-Linear Cosmic Structures using Field-Level Emulators
• Cosmological constraints from the cross-correlation of DESI Luminous Red Galaxies with CMB lensing from Planck PR4 and ACT DR6
• The e-MANTIS emulator: fast predictions of the non-linear matter power spectrum in $f(R)$CDM cosmology
• Primordial non-Gaussianity and non-Gaussian Covariance
• Perturbation theory with dispersion and higher cumulants: non-linear regime
• Cosmological Vlasov-Poisson equations for dark matter: Recent developments and connections to selected plasma problems
• The AGORA High-resolution Galaxy Simulations Comparison Project. V: Satellite Galaxy Populations In A Cosmological Zoom-in Simulation of A Milky Way-mass Halo
• Differentiable Cosmological Simulation with Adjoint Method
• Geodesic motion and phase-space evolution of massive neutrinos
• The RayGalGroupSims cosmological simulation suite for the study of relativistic effects: an application to lensing-matter clustering statistics
• Indra: a Public Computationally Accessible Suite of Cosmological $N$-body Simulations
• The FLAMINGO project: Baryonic impact on weak gravitational lensing convergence peak counts
• Cold dark matter protohalo structure around collapse: Lagrangian cosmological perturbation theory versus Vlasov simulations
• Exploring the effects of primordial non-Gaussianity at galactic scales
• Analytical growth functions for cosmic structures in a $\Lambda$CDM Universe
• GEO-FPT: a model of the galaxy bispectrum at mildly non-linear scales
• Improving initialization and evolution accuracy of cosmological neutrino simulations
• BullFrog: Multi-step perturbation theory as a time integrator for cosmological simulations
• The effect of local universe constraints on halo abundance and clustering
• Radiative transfer of Lyman-$\alpha$ photons at cosmic dawn with realistic gas physics
• Renormalization group and UV completion of cosmological perturbations: Gravitational collapse as a critical phenomenon
• On general-relativistic Lagrangian perturbation theory and its non-perturbative generalization
• Halo mass function in scale invariant models
• Numerical Discreteness Errors in Multi-Species Cosmological N-body Simulations
• COLIBRE: calibrating subgrid feedback in cosmological simulations that include a cold gas phase
• Perturbation-theory informed integrators for cosmological simulations
• Neutrinos in N-body simulations
• Starting Cosmological Simulations from the Big Bang
• Relativistic matter bispectrum of cosmic structures on the light cone
• Relativistic second-order initial conditions for simulations of large-scale structure
• Fast and accurate collapse-time predictions for collisionless matter
• Effect of Separate Initial Conditions on the Lyman-$\alpha$ Forest in Simulations
• Field-level Reconstruction from Foreground-Contaminated 21-cm Maps
• PySCo: A fast Particle-Mesh $N$-body code for modified gravity simulations in Python
• On the Robustness of Cluster Clustering Covariance Calibration
• Aemulus $\nu$: Precision halo mass functions in w$\nu$CDM cosmologies
• On the origin of transient features in cosmological N-Body Simulations
• Simulation budgeting for hybrid effective field theories
• Parity Violation in Galaxy Shapes: Primordial Non-Gaussianity
• Euclid preparation. L. Calibration of the linear halo bias in $\Lambda (\nu )$CDM cosmologies
• Cosmic voids and BAO with relative baryon-CDM perturbations
• The Manticore Project II: Bayesian digital twins of cosmic structure across the SDSS and BOSS volumes
• Ginnungagap -- a massively parallel cosmological initial conditions generator
• Post-collapse Lagrangian perturbation theory in three dimensions
• Tidal adaptive softening and artificial fragmentation in cosmological simulations