Evolution of the density profiles of dark matter haloes

We use numerical simulations in a LambdaCDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to seven million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10(11) -10(14) M-., down toapproximately 0.5 per cent of the virial radius, and from redshift of 5 to the present, covering a larger range in parameter space than previous works. We confirm that the cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. Motivated by the diversity and evolution of halo profile shapes, we fit our haloes to the two-parameter profile, rho proportional to 1 {(c(gamma)r /r(vir))(gamma)[1 + (c(gamma)r /r(vir))](3-gamma)}, where the steepness of the cusp is given by the asymptotic inner slope parameter,, and its radial extent is described by the concentration parameter, c (with c defined as the virial radius divided by the concentration radius). In our simulations, we find gamma similar or equal to 1.4 - 0.08 log(10)(M/M-*) for haloes of 0.01 M-*-1000 M-*, with a large scatter of Deltagamma similar to +/-0.3, where M-* is the redshift dependent characteristic mass of collapsing haloes; and c(gamma) similar or equal to 8.(M/M)(-0.15), with a large M/M-* dependent scatter roughly equal to +/- c(gamma). Our redshift zero haloes have inner slope parameters ranging approximately from r(-1) (Navarro, Frank and White) to r(-1.5) (Moore et al.), with a median of roughly r(-1.3). This two-parameter profile fit works well for all types of haloes in our simulations, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power-law cosmologies of P proportional to k(n), with n = (0, -1, -2, -2.7). Here we find that the concentration radius and the inner cusp slope are both a function of n, with larger concentration radii and shallower cusps for steeper power spectra. We have completed a thorough resolution study and find that the minimum resolved radius is well described by the mean interparticle separation over a range of masses and redshifts. The trend of steeper and more concentrated cusps for smaller M/M-* haloes clearly shows that dwarf-sized LambdaCDM haloes have, on average, significantly steeper density profiles within the inner few per cent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac. uk/ similar toreed/profile.html.