Dependence of the Inner DM Profile on the Halo Mass

I compare the density profile of dark matter (DM) halos in cold dark matter (CDM) N-body simulations with 1 Mpc, 32 Mpc, 256 Mpc and 1024 Mpc box sizes. In dimensionless units the simulations differ only for the initial power spectrum of density perturbations. I compare the profiles when the most massive halos are composed of about 10^5 DM particles. The DM density profiles of the halos in the 1 Mpc box show systematically shallower cores with respect to the corresponding halos in the 32 Mpc simulation that have masses, M_{dm}, typical of the Milky Way and are fitted by a NFW profile. The DM density profiles of the halos in the 256 Mpc box are consistent with having steeper cores than the corresponding halos in the 32 Mpc simulation, but higher mass resolution simulations are needed to strengthen this result. Combined, these results indicate that the density profile of DM halos is not universal, presenting shallower cores in dwarf galaxies and steeper cores in clusters. Physically the result sustains the hypothesis that the mass function of the accreting satellites determines the inner slope of the DM profile. In comoving coordinates, r, the profile \rho_{dm} \propto 1/(X^\alpha(1+X)^{3-\alpha}), with X=c_\Delta r/r_\Delta, r_\Delta is the virial radius and \alpha =\alpha(M_{dm}), provides a good fit to all the DM halos from dwarf galaxies to clusters at any redshift with the same concentration parameter c_\Delta ~ 7. The slope, \gamma, of the outer parts of the halo appears to depend on the acceleration of the universe: when the scale parameter is a=(1+z)^{-1} < 1, the slope is \gamma ~ 3 as in the NFW profile, but \gamma ~ 4 at a > 1 when \Omega_\Lambda ~ 1 and the universe is inflating.[abridged]