Volume 552, April 2013
|Number of page(s)||16|
|Section||Cosmology (including clusters of galaxies)|
|Published online||21 March 2013|
A toy model to test the accuracy of cosmological N-body simulations
1 Centro Studi e Ricerche Enrico Fermi, via Panisperna 89 A, Compendio del Viminale, 00184 Rome, Italy
2 Istituto dei Sistemi Complessi CNR, via dei Taurini 19, 00185 Rome, Italy
Received: 3 January 2013
Accepted: 31 January 2013
The evolution of an isolated over-density represents a useful toy model to test the accuracy of a cosmological N-body code in the non-linear regime as it is approximately equivalent to that of a truly isolated cloud of particles, with same density profile and velocity distribution, in a non-expanding background. This is the case as long as the system size is smaller than the simulation box side, so that its interaction with the infinite copies can be neglected. In such a situation, the over-density rapidly undergoes to a global collapse forming a quasi stationary state in virial equilibrium. However, by evolving the system with a cosmological code (GADGET) for a sufficiently long time, a clear deviation from such quasi-equilibrium configuration is observed. This occurs in a time tLI that depends on the values of the simulation numerical parameters such as the softening length and the time-stepping accuracy, i.e. it is a numerical artifact related to the limited spatial and temporal resolutions. The analysis of the Layzer-Irvine cosmic energy equation confirms that this deviation corresponds to an unphysical dynamical regime. By varying the numerical parameters of the simulation and the physical parameters of the system we show that the unphysical behaviour originates from badly integrated close scatterings of high-velocity particles. We find that, while the structure may remain virialized in the unphysical regime, its density and velocity profiles are modified with respect to the quasi-equilibrium configurations, converging, however, to well defined shapes, the former characterised by a Navarro Frenk White-type behaviour.
Key words: methods: numerical / galaxies: halos / galaxies: formation / large-scale structure of Universe / dark matter
© ESO, 2013
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.