Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations

Abstract

The effects of discreteness arising from the use of the N-body method on the accuracy of simulations of cosmological structure formation are not currently well understood. After a discussion of how the relevant discretisation parameters introduced should be extrapolated to recover the Vlasov-Poisson limit, we study numerically, and with analytical methods we have developed recently, the central issue of how finite particle density affects the precision of results. In particular we focus on the power spectrum at wavenumbers around and above the Nyquist wavenumber, in simulations in which the force resolution is taken smaller than the initial interparticle spacing. Using simulations of identical theoretical initial conditions sampled on four different "pre-initial" configurations (three different Bravais lattices, and a glass) we obtain a {\it lower bound} on the real discreteness error. With the guidance of our analytical results, we establish with confidence that the measured dispersion is not contaminated either by finite box size effects or by subtle numerical effects. Our results show notably that, at wavenumbers {\it below} the Nyquist wavenumber, the dispersion increases monotonically in time throughout the simulation, while the same is true above the Nyquist wavenumber once non-linearity sets in. For normalizations typical of cosmological simulations, we find lower bounds on errors at the Nyquist wavenumber of order of a percent, and larger above this scale. The only way this error may be reduced below these levels at these scales, and indeed convergence to the physical limit firmly established, is by extrapolation, at fixed values of the other relevant parameters, to the regime in which the mean comoving interparticle distance becomes less than the force smoothing scale.

26 pages, 15 figures, minor changes, slightly shortened, version to be published in MNRAS

References in corpus (46)

• Simulating the joint evolution of quasars, galaxies and their large-scale distribution
• Stable clustering, the halo model and nonlinear cosmological power spectra
• GADGET: A code for collisionless and gasdynamical cosmological simulations
• The Inner Structure of LambdaCDM Halos I: A Numerical Convergence Study
• Nonlinear evolution of cosmological power spectra
• Local Density Fluctuations, Hyperuniformity, and Order Metrics
• Transients from Initial Conditions in Cosmological Simulations
• Evolution of structure in cold dark matter universes
• Convergence and scatter of cluster density profiles
• Discreteness effects in simulations of Hot/Warm dark matter
• The Halo Mass Function: High-Redshift Evolution and Universality
• Dark matter annihilation in the halo of the Milky Way
• Calibrating the Nonlinear Matter Power Spectrum: Requirements for Future Weak Lensing Surveys
• Initial conditions to cosmological N-body simulations, or how to run an ensemble of simulations
• MLAPM - a C code for cosmological simulations
• The Glass-like Universe: Real-space correlation properties of standard cosmological models
• The Cosmic Code Comparison Project
• Two body relaxation in CDM simulations
• Fundamental Discreteness Limitations of Cosmological N-Body Clustering Simulations
• Simulation techniques for cosmological simulations
• Adhesive Gravitational Clustering
• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• On effects of resolution in dissipationless cosmological simulations
• Discreteness effects in cosmological N-body simulations
• Dependence of the non-linear mass power spectrum on the equation of state of dark energy
• Demonstrating Discreteness and Collision Error in Cosmological N-body Simulations of Dark Matter Gravitational Clustering
• Generating Cosmological Gaussian Random Fields
• Linear perturbative theory of the discrete cosmological N-body problem
• Generation of Primordial Cosmological Perturbations from Statistical Mechanical Models
• Critical Index and Fixed Point in the Transfer of Power in Nonlinear Gravitational Clustering
• Circular velocity profiles of dark matter haloes
• Quantification of discreteness effects in cosmological N-body simulations: I. Initial Conditions
• Effects of the size of cosmological N-Body simulations on physical quantities -- I: Mass Function
• Reliability of the dark matter clustering in cosmological N-body simulation on scales below the mean separation length of particles
• Quantification of discreteness effects in cosmological N-body simulations: II. Evolution up to shell crossing
• Galaxy Formation: Warm Dark Matter, Missing Satellites, and the Angular Momentum Problem
• A Test of the Particle Paradigm in N-Body Simulations
• Gravitational evolution of a perturbed lattice and its fluid limit
• Power-law correlation and discreteness in cosmological N-body simulations
• Gravitational Dynamics of an Infinite Shuffled Lattice of Particles
• Chaos and the continuum limit in the gravitational N-body problem II. Nonintegrable potentials
• Gravitational Dynamics of an Infinite Shuffled Lattice: Particle Coarse-grainings, Non-linear Clustering and the Continuum Limit
• Two Body Relaxation in Simulated Cosmological Haloes
• Gravitational dynamics of an infinite shuffled lattice: early time evolution and universality of non-linear correlations
• Particle linear theory on a self-gravitating perturbed cubic Bravais lattice
• Gravitational collapse in an expanding background and the role of substructure II: Excess power at small scales and its effect of collapse of structures at larger scales

Cited by in corpus (38)

• The Coyote Universe I: Precision Determination of the Nonlinear Matter Power Spectrum
• Disruption of Dark Matter Substructure: Fact or Fiction?
• An accurate physical model for halo concentrations
• Matter power spectrum and the challenge of percent accuracy
• How to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations
• Light propagation in statistically homogeneous and isotropic universes with general matter content
• Improving Initial Conditions for Cosmological $N$-Body Simulations
• Toward an accurate mass function for precision cosmology
• Matter Power Spectrum Covariance Matrix from the DEUS-PUR ΛCDM simulations: Mass Resolution and non-Gaussian Errors
• Cosmic variance limited Baryon Acoustic Oscillations from the DEUS-FUR $\Lambda$CDM simulation
• 1-d gravity in infinite point distributions
• On the starting redshift for cosmological simulations: Focusing on halo properties
• Quantifying resolution in cosmological N-body simulations using self-similarity
• The Fractal Geometry of the Cosmic Web and its Formation
• An alternative method to generate pre-initial conditions for cosmological $N$-body simulations
• Novel Adaptive softening for collisionless N-body simulations: Eliminating spurious halos
• Evolution of isolated overdensities as a control on cosmological N body simulations
• Modelling non-linear evolution using Lagrangian Perturbation Theory (LPT) re-expansions
• Non-linear gravitational clustering of cold matter in an expanding universe: indications from 1D toy models
• The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations
• Exponents of non-linear clustering in scale-free one dimensional cosmological simulations
• Initial conditions for cosmological N-body simulations of the scalar sector of theories of Newtonian, Relativistic and Modified Gravity
• Phase Space dynamics of triaxial collapse: Joint density-velocity evolution
• Dynamics of finite and infinite self-gravitating systems with cold quasi-uniform initial conditions
• Accuracy of power spectra in dissipationless cosmological simulations
• Improving initialization and evolution accuracy of cosmological neutrino simulations
• Stable clustering and the resolution of dissipationless cosmological N-body simulations
• Numerical convergence of pre-initial conditions on dark matter halo properties
• Improved two-point correlation function estimates using glass-like distributions as a reference sample
• A toy model to test the accuracy of cosmological N-body simulations
• Inner structure of dark matter halos at high $z$ in cosmological models with non-power-law primordial spectra
• Smooth halos in the cosmic web
• Numerical Discreteness Errors in Multi-Species Cosmological N-body Simulations
• Impacts of pre-initial conditions on anisotropic separate universe simulations: a boosted tidal response in the epoch of reionization
• Perturbation-theory informed integrators for cosmological simulations
• Glass-Like Random Catalogues for Two-Point Estimates on the Light Cone
• Particle loads for cosmological simulations with equal-mass dark matter and baryonic particles
• A dynamics-based density profile for dark haloes. I. Algorithm and basic results