Juerg Diemand, Ben Moore, Joachim Stadel
Published 2004-02-11, updated 2005-03-23Version 5
We present new results from a series of LCDM simulations of cluster mass halos resolved with high force and mass resolution. These results are compared with recently published simulations from groups using various codes including PKDGRAV, ART, TPM, GRAPE and GADGET. Careful resolution tests show that with 25 million particles within the high resolution region we can resolve to about 0.3% of the virial radius and that convergence in radius is proportional to the mean interparticle separation. The density profiles of 26 high resolution clusters obtained with the different codes and from different initial conditions agree very well. The average logarithmic slope at one percent of the virial radius is $\gamma = 1.26$ with a scatter of $\pm 0.17$. Over the entire resolved regions the density profiles are well fitted by a smooth function that asymptotes to a central cusp $\rho \propto r^{-\gamma}$, where we find $\gamma=1.16\pm 0.14$ from the mean of the fits to our six highest resolution clusters.
Comments: MNRAS accepted version. Due to an error in the initial conditions these simulations have a lower sigma_8 than the published value, 0.7 instead of 0.9. We thank Mike Kuhlen who helped us finding this mistake. See the erratum at http://www-theorie.physik.unizh.ch/~diemand/proerr.pdf . Images and movies available at http://www-theorie.physik.unizh.ch/~diemand/clusters/
Journal: Mon.Not.Roy.Astron.Soc. 353 (2004) 624
The evolution of substructure III: the outskirts of clusters
Internal and External Alignment of the Shapes and Angular Momenta of LCDM Halos
Dark Matter and Baryon Fraction at the Virial Radius in Abell 2256
