Dependence of the inner dark matter profile on the halo mass

Abstract

I compare the density profile of dark matter (DM) haloes in cold dark matter (CDM) N-body simulations with 1, 32, 256 and 1024 h - 1 Mpc box sizes. I compare the profiles when the most massive haloes are composed of about 10 5 DM particles. The DM density profiles of the haloes in the 1 h - 1 Mpc box at redshift z ∼ 10 show systematically shallower cores with respect to the corresponding haloes in the 32 h - 1 Mpc simulation at z ∼ 3 that have masses, M d m , typical of the Milky Way and are fitted by a Navarro-Frenk-White (NFW) profile. The DM density profiles of the haloes in the 256 h - 1 Mpc box at z ∼ 0 are consistent with having steeper cores than the corresponding haloes in the 32 h - 1 Mpc simulation, but higher mass resolution simulations are needed to strengthen this result. Combined, these results suggest that the density profile of DM haloes is not universal, presenting shallower cores in dwarf galaxies and steeper cores in clusters. More work is needed to validate this finding at z = 0. Physically, the result sustains the hypothesis that the mass function of the accreting satellites determines the inner slope of the DM profile. However, the result can also be interpreted as a trend with the dynamical state in the assembly process of haloes of different mass. In comoving coordinates, r, the profile ρ d m with X = c Δ r/r Δ and a = (9 + 3n)/(5 + n) 1.3 + (M 1 / 6 dm,14 - 1)/(M 1 / 6 dm,14 + 1), provides a good fit to all the DM haloes from dwarf galaxies to clusters at any redshift with the same concentration parameter c Δ ∼ 7. Here, r Δ (M d m ) is the virial radius, n is the effective spectral index of the initial power spectrum of density perturbations and M d m , 1 4 = M d m /(3 × 10 1 4 M O .). The slope, y, 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 γ 3 as in the NEW profile, but γ 4 at a > 1 when Ω Λ ∼ 1 and the universe is inflating. The shape of the DM profiles presents a significant scatter around the mean. It is therefore important to analyse a significant statistical sample of haloes in order to determine the mean profile. I compare the DM profiles in the 1 h - 1 Mpc box with the same simulation including stars, baryons and radiative transfer presented by Ricotti, Gnedin and Shull. Radiative feedback effects produce a larger scatter in the density profile shapes but, on average, do not affect the shape of the DM profiles significantly.