No Arabic abstract
We explore the cosmological halo-to-halo scatter of the distribution of mass within dark matter halos utilizing a well-resolved statistical sample of clusters from the cosmological Millennium simulation. We find that at any radius, the spherically-averaged dark matter density of a halo (corresponding to the smooth-component) and its logarithmic slope are well-described by a Gaussian probability distribution. At small radii (within the scale radius), the density distribution is fully determined by the measured Gaussian distribution in halo concentrations. The variance in the radial distribution of mass in dark matter halos is important for the interpretation of direct and indirect dark matter detection efforts. The scatter in mass profiles imparts approximately a 25 percent cosmological uncertainty in the dark matter density at the Solar neighborhood and a factor of ~3 uncertainty in the expected Galactic dark matter annihilation flux. The aggregate effect of halo-to-halo profile scatter leads to a small (few percent) enhancement in dark matter annihilation background if the Gaussian concentration distribution holds for all halo masses versus a 10 percent enhancement under the assumption of a log-normal concentration distribution. The Gaussian nature of the cluster profile scatter implies that the technique of stacking halos to improve signal to noise should not suffer from bias.
Dark matter-dominated cluster-scale halos act as an important cosmological probe and provide a key testing ground for structure formation theory. Focusing on their mass profiles, we have carried out (gravity-only) simulations of the concordance LCDM cosmology, covering a mass range of 2.10^{12}-2.10^{15} solar mass/h and a redshift range of z=0-2, while satisfying the associated requirements of resolution and statistical control. When fitting to the Navarro-Frenk-White profile, our concentration-mass (c-M) relation differs in normalization and shape in comparison to previous studies that have limited statistics in the upper end of the mass range. We show that the flattening of the c-M relation with redshift is naturally expressed if c is viewed as a function of the peak height parameter, u. Unlike the c-M relation, the slope of the c- u relation is effectively constant over the redshift range z=0-2, while the amplitude varies by ~30% for massive clusters. This relation is, however, not universal: Using a simulation suite covering the allowed wCDM parameter space, we show that the c- u relation varies by about +/- 20% as cosmological parameters are varied. At fixed mass, the c(M) distribution is well-fit by a Gaussian with sigma_c/c = 0.33, independent of the radius at which the concentration is defined, the halo dynamical state, and the underlying cosmology. We compare the LCDM predictions with observations of halo concentrations from strong lensing, weak lensing, galaxy kinematics, and X-ray data, finding good agreement for massive clusters (M > 4.10^{14} solar mass/h), but with some disagreements at lower masses. Because of uncertainty in observational systematics and modeling of baryonic physics, the significance of these discrepancies remains unclear.
The most accurate way to get information on the mass of the MACHOs (Massive Astrophysical Compact Halo Objects) is to use the method of mass moments. For the microlensing events detected so far by the EROS and the MACHO collaborations in the Large Magellanic Cloud the average mass turns out to be 0.08$M_{odot}$. Assuming a spherical standard halo model we find that MACHOs contribute about 20% to the halo dark matter. The eleven events recorded by OGLE, mainly during its first two years of operation, in the galactic bulge lead to an average mass of 0.29$M_{odot}$, whereas forty events detected by MACHO during its first year give 0.16$M_{odot}$, thus suggesting that the lens objects are faint disk stars.
Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $dlogrho/dlog r approx -1$, and some dispersion, which is generally parametrized by the varying index $alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $phi_0$, Einasto $alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.
Several direct detection experiments, including recently CDMS-II, have reported signals consistent with 5 to 10 GeV dark matter (DM) that appear to be in tension with null results from XENON and LUX experiments; these indicate a careful review of the theoretical basis, including the galactic DM velocity distribution function (VDF). We establish a VDF parameter space from DM-only cosmological simulations and illustrate that seemingly contradictory experimental results can be made consistent within this parameter space. Future experimental limits should be reported after they are marginalized over a range of VDF parameters.
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.