Transients from Initial Conditions in Cosmological Simulations

Abstract

We study the impact of setting initial conditions in numerical simulations using the standard procedure based on the Zel'dovich approximation (ZA). As it is well known from perturbation theory, ZA initial conditions have incorrect second and higher-order growth and therefore excite long-lived transients in the evolution of the statistical properties of density and velocity fields. We also study the improvement brought by using more accurate initial conditions based on second-order Lagrangian perturbation theory (2LPT). We show that 2LPT initial conditions reduce transients significantly and thus are much more appropriate for numerical simulations devoted to precision cosmology. Using controlled numerical experiments with ZA and 2LPT initial conditions we show that simulations started at redshift z_i=49 using the ZA underestimate the power spectrum in the nonlinear regime by about 2,4,8 % at z=0,1,3 respectively, whereas the mass function of dark matter halos is underestimated by 5% at m=10^15 M_sun/h (z=0) and 10% at m=2x10^14M_sun/h (z=1). The clustering of halos is also affected to the few percent level at z=0. These systematics effects are typically larger than statistical uncertainties in recent mass function and power spectrum fitting formulae extracted from numerical simulations. At large scales, the measured transients in higher-order correlations can be understood from first principle calculations based on perturbation theory.

14 pages, 14 figures, code to generate 2LPT initial conditions available at http://cosmo.nyu.edu/roman/2LPT . Typos corrected, Fig.13 symbols consistent with Fig.11,12

References in corpus (3)

• Cosmological parameters from combining the Lyman-alpha forest with CMB, galaxy clustering and SN constraints
• Cosmology and the Bispectrum
• The Lyman-alpha forest and WMAP year three

Cited by in corpus (25)

• Nonlinear Evolution of Baryon Acoustic Oscillations
• The morphology of HII regions during reionization
• The Scale Dependence of Halo and Galaxy Bias: Effects in Real Space
• The Halo Mass Function: High-Redshift Evolution and Universality
• Large-Scale Anisotropic Correlation Function of SDSS Luminous Red Galaxies
• The detectability of baryonic acoustic oscillations in future galaxy surveys
• Multi-Point Propagators in Cosmological Gravitational Instability
• `Eppur Si Muove': On The Motion of the Acoustic Peak in the Correlation Function
• The Effect of Thermal Neutrino Motion on the Non-linear Cosmological Matter Power Spectrum
• Probing Reionization with the 21 cm-Galaxy Cross Power Spectrum
• The simulated 21 cm signal during the epoch of reionization : full modeling of the Ly-alpha pumping
• Accurate estimators of power spectra in N-body simulations
• How Very Massive Metal Free Stars Start Cosmological Reionization
• Perturbation Theory Reloaded II: Non-linear Bias, Baryon Acoustic Oscillations and Millennium Simulation In Real Space
• Evaluating backreaction with the peak model of structure formation
• An analytic model for the bispectrum of galaxies in redshift space
• Higher Order Contributions to the 21 cm Power Spectrum
• Angular Signatures of Annihilating Dark Matter in the Cosmic Gamma-Ray Background
• Impact of Scale Dependent Bias and Nonlinear Structure Growth on the ISW Effect: Angular Power Spectra
• Towards quantitative control on discreteness error in the non-linear regime of cosmological N body simulations
• Simulations of Baryon Acoustic Oscillations I: Growth of Large-Scale Density Fluctuations
• Using the Zeldovich dynamics to test expansion schemes
• Subhaloes in Scale-Free Cosmologies
• Transients from initial conditions based on Lagrangian perturbation theory in N-body simulations
• Baryon oscillations in galaxy and matter power-spectrum covariance matrices