Matter power spectrum and the challenge of percent accuracy

Abstract

Abstract. Future galaxy surveys require one percent precision in the theoretical knowledge of the power spectrum over a large range including very nonlinear scales. While this level of accuracy is easily obtained in the linear regime with perturbation theory, it represents a serious challenge for small scales where numerical simulations are required. In this paper we quantify the precision of present-day N -body methods, identifying main potential error sources from the set-up of initial conditions to the measurement of the final power spectrum. We directly compare three widely used N -body codes, Ramses, Pkdgrav3, and Gadget3 which represent three main discretisation techniques: the particle-mesh method, the tree method, and a hybrid combination of the two. For standard run parameters, the codes agree to within one percent at k ≤ 1 h Mpc −1 and to within three percent at k ≤ 10 h Mpc −1 . We also consider the bispectrum and show that the reduced bispectra agree at the sub-percent level for k ≤ 2 h Mpc −1 . In a second step, we quantify potential errors due to initial conditions, box size, and resolution using an extended suite of simulations performed with our fastest code Pkdgrav3. We demonstrate that the simulation box size should not be smaller than L = 0.5 h −1 Gpc to avoid systematic finite-volume effects (while much larger boxes are required to beat down the statistical sample variance). Furthermore, a maximum particle mass of M p = 10 9 h −1 M is required to conservatively obtain one percent precision of the matter power spectrum. As a consequence, numerical simulations covering large survey volumes of upcoming missions such as DES, LSST, and Euclid will need more than a trillion particles to reproduce clustering properties at the targeted accuracy.arXiv:1503.05920v3 [astro-ph.CO]

Inverted initial conditions: Exploring the growth of cosmic structure and voids

“…The corrections reach ∼5% even at small scales, k ¼ 1.0h −1 Mpc, and modest redshifts, z ¼ 1, where one might hope that the box-size effects are minimal. This is consistent with what is found by averaging over hundreds of realizations [43] or comparing to increased box sizes which better sample the large-scale modes [23].…”

The Mira-Titan Universe. II. Matter Power Spectrum Emulation

“…In order to avoid problems due to large neutrino velocities, they started the simulations as late as z in = 24 for the lightest neutrinos, using the Zel'dovich approximation. This leads to systematic inaccuracies in the power spectrum at the few percent level as shown in Heitmann et al (2010) and Schneider et al (2016). (In both papers effects at the 2-3% level were shown with a starting redshift z in = 50 at k ∼ 1hMpc −1 , extrapolating these results would suggest a 5% error due to the late start alone, and even more at higher k).…”

Renormalization-group flow of the effective action of cosmological large-scale structures

“…Technically, this task amounts to solving the collisionless Vlasov-Boltzmann equation for appropriate classes of initial conditions [1]. With recent progress in both computational power and coding techniques, N -body simulations describe by now structure evolution down to small scales where baryonic effects become increasingly important [2][3][4]. Even for pure cold dark matter simulations, perturbative techniques cannot be expected to apply on scales where non-linearities become strong and the process of virialization starts to dominate.…”