No Arabic abstract
High-resolution cosmological N-body simulations were performed in order to study the substructure of Milky Way-like galactic halos and the density profiles of halos in a warm dark matter scenario. The results favor this scenario with respect to the cold dark matter one.
Although high-resolution N-body simulations make robust empirical predictions for the density distribution within cold dark matter halos, these studies have yielded little physical insight into the origins of the distribution. We investigate the problem using analytic and semi-analytic approaches. Simple analytic considerations suggest that the inner slope of dark matter halos cannot be steeper than alpha=2 (rho ~ r^-alpha), with alpha=1.5-1.7 being a more realistic upper limit. Our analysis suggests that any number of effects, both real (angular momentum from tidal torques, secondary perturbations) and artificial (two-body interactions, the accuracy of the numerical integrator, round-off errors), will result in shallower slopes. We also find that the halos should exhibit a well-defined relation between r_peri/r_apo and j_theta/j_r. We derive this relation analytically and speculate that it may be universal. Using a semi-analytic scheme based on Ryden & Gunn (1987), we further explore the relationship between the specific angular momentum distribution in a halo and its density profile. For now we restrict ourselves to halos that form primarily via nearly-smooth accretion of matter, and only consider the specific angular momentum generated by secondary perturbations associated with the cold dark matter spectrum of density fluctuations. Compared to those formed in N-body simulations, our ``semi-analytic halos are more extended, have flatter rotation curves and have higher specific angular momentum, even though we have not yet taken into account the effects of tidal torques. Whether the density profiles of numerical halos is indeed the result of loss in angular momentum outside the central region, and whether this loss is a feature of hierarchical merging and major mergers in particular, is under investigation.
The shapes of individual self-gravitating structures of an ensemble of identical, collisionless particles have remained elusive for decades. In particular, a reason why mass density profiles like the Navarro-Frenk-White or the Einasto profile are good fits to simulation- and observation-based dark matter halos has not been found. Given the class of three dimensional, spherically symmetric power-law probability density distributions to locate individual particles in the ensemble mentioned above, we derive the constraining equation for the power-law index for the most and least likely joint ensemble configurations. We find that any dark matter halo can be partitioned into three regions: a core, an intermediate part, and an outskirts part up to boundary radius $r_mathrm{max}$. The power-law index of the core is determined by the mean radius of the particle distribution within the core. The intermediate region becomes isothermal in the limit of infinitely many particles. The slope of the mass density profile far from the centre is determined by the choice of $r_mathrm{max}$ with respect to the outmost halo particle, such that two typical limiting cases arise that explain the $r^{-3}$-slope for galaxy-cluster outskirts and the $r^{-4}$-slope for galactic outskirts. Hence, we succeed in deriving the mass density profiles of common fitting functions from a general viewpoint. These results also allow to find a simple explanation for the cusp-core-problem and to separate the halo description from its dynamics.
We present a study of unprecedented statistical power regarding the halo-to-halo variance of dark matter substructure. Using a combination of N-body simulations and a semi-analytical model, we investigate the variance in subhalo mass fractions and subhalo occupation numbers, with an emphasis on how these statistics scale with halo formation time. We demonstrate that the subhalo mass fraction, f_sub, is mainly a function of halo formation time, with earlier forming haloes having less substructure. At fixed formation redshift, the average f_sub is virtually independent of halo mass, and the mass dependence of f_sub is therefore mainly a manifestation of more massive haloes assembling later. We compare observational constraints on f_sub from gravitational lensing to our model predictions and simulation results. Although the inferred f_sub are substantially higher than the median LCDM predictions, they fall within the 95th percentile due to halo-to-halo variance. We show that while the halo occupation distribution of subhaloes, P(N|M), is super-Poissonian for large <N>, a well established result, it becomes sub-Poissonian for <N> < 2. Ignoring the non-Poissonity results in systematic errors of the clustering of galaxies of a few percent, and with a complicated scale- and luminosity-dependence. Earlier-formed haloes have P(N|M) closer to a Poisson distribution, suggesting that the dynamical evolution of subhaloes drives the statistics towards Poissonian. Contrary to a recent claim, the non-Poissonity of subhalo occupation statistics does not vanish by selecting haloes with fixed mass and fixed formation redshift. Finally, we use subhalo occupation statistics to put loose constraints on the mass and formation redshift of the Milky Way halo. Using observational constraints on the V_max of the most massive satellites, we infer that 0.25<M_vir/10^12M_sun/h<1.4 and 0.1<z_f<1.4 at 90% confidence.
We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the coldest model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $leq 10^9,h^{-1},M_odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{mathrm{max}}^{z=0} leq50~mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.
We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.