1308.0057 (Gerasimos Rigopoulos et al.)

On the accuracy of N-body simulations at very large scales    [PDF]

Gerasimos Rigopoulos, Wessel Valkenburg

We examine the deviation of Cold Dark Matter particle trajectories from the Newtonian result as the size of the region under study becomes comparable to or exceeds the particle horizon. To first order in the gravitational potential, the general relativistic result coincides with the Zel'dovich approximation and hence the Newtonian prediction on all scales. At second order, General Relativity predicts corrections which overtake the corresponding second order Newtonian terms above a certain scale of the order of the Hubble radius. However, since second order corrections are very much suppressed on such scales, we conclude that simulations which exceed the particle horizon but use Newtonian equations to evolve the particles, reproduce the correct trajectories very well. At z=49, the dominant relativistic corrections to the power spectrum on scales close to the horizon are at most of the order of $10^{-3}$, while the differences in the positions of real space features are affected at a level below $10^{-6}$. Similar conclusions hold at z=0. In the process, we clarify the relation of N-body results to relativistic considerations.

View original: http://arxiv.org/abs/1308.0057