The inner structure of ΛCDM haloes  I. A numerical convergence study
Abstract
We present a comprehensive set of convergence tests which explore the role of various numerical parameters on the equilibrium structure of a simulated dark matter halo. We report on results obtained with two independent, stateoftheart, multistepping, parallel Nbody codes: PKDGRAV and GADGET. We find that convergent mass profiles can be obtained for suitable choices of the gravitational softening, timestep, force accuracy, initial redshift, and particle number. For softenings chosen so that particle discreteness effects are negligible, convergence in the circular velocity is obtained at radii where the following conditions are satisfied: (i) the timestep is much shorter than the local orbital timescale; (ii) accelerations do not exceed a characteristic acceleration imprinted by the gravitational softening; and (iii) enough particles are enclosed so that the collisional relaxation timescale is longer than the age of the Universe. Convergence also requires sufficiently high initial redshift and accurate force computations. Poor spatial, time, or force resolution leads generally to systems with artificially low central density, but may also result in the formation of artificially dense central cusps. We have explored several adaptive timestepping choices and we have obtained the best results when individual timesteps are chosen according to the local acceleration and the gravitational softening (Δt_{i}~ (ɛ/a_{i})^{1/2}), although further experimentation may yield better and more efficient criteria. The most stringent requirement for convergence is typically that imposed on the particle number by the collisional relaxation criterion. This implies that, in order to estimate accurate circular velocities at radii where the density contrast may reach ~10^{6}, the region must enclose of the order of 3000 particles (or more than a few times 10^{6} within the virial radius). Applying these criteria to a galaxysized ΛCDM halo, we find that the spherically averaged density profile becomes progressively shallower from the virial radius inwards, reaching a logarithmic slope shallower than 1.2 at the innermost resolved point, ~0.005 r_{200}, with little evidence for convergence to a powerlaw behaviour in the inner regions.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 January 2003
 DOI:
 10.1046/j.13658711.2003.05925.x
 arXiv:
 arXiv:astroph/0201544
 Bibcode:
 2003MNRAS.338...14P
 Keywords:

 gravitation;
 cosmology: theory;
 dark matter;
 Astrophysics
 EPrint:
 23 pages, 15 figures. Published in MNRAS. Changes : (1) Addition of paragraph at end of Introduction, highlighting more technical sections of paper, (2) Revised summaries at end of major (sub)sections, (3) Modified caption to Figure 15