No Arabic abstract
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.
Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $sim10^6$M$_{odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.
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.
We derive a model for Sunyaev--Zeldovich data from a galaxy cluster which uses an Einasto profile to model the clusters dark matter component. This model is similar to the physical models for clusters previously used by the Arcminute Microkelvin Imager (AMI) consortium, which model the dark matter using a Navarro-Frenk-White (NFW) profile, but the Einasto profile provides an extra degree of freedom. We thus present a comparison between two physical models which differ only in the way they model dark matter: one which uses an NFW profile (PM I) and one that uses an Einasto profile (PM II). We illustrate the differences between the models by plotting physical properties of clusters as a function of cluster radius. We generate AMI simulations of clusters which are textit{created} and textit{analysed} with both models. From this we find that for 14 of the 16 simulations, the Bayesian evidence gives no preference to either of the models according to the Jeffreys scale, and for the other two simulations, weak preference in favour of the correct model. However, for the mass estimates obtained from the analyses, the values were within $1sigma$ of the input values for 14 out of 16 of the clusters when using the correct model, but only in 6 out of 16 cases when the incorrect model was used to analyse the data. Finally we apply the models to real data from cluster A611 obtained with AMI, and find the mass estimates to be consistent with one another except in the case of when PM II is applied using an extreme value for the Einasto shape parameter.
We consider the dynamics in and near galaxy clusters. Gas, dark matter and galaxies are presently falling into the clusters between approximately 1 and 5 virial radii. At very large distances, beyond 10 virial radii, all matter is following the Hubble flow, and inside the virial radius the matter particles have on average zero radial velocity. The cosmological parameters are imprinted on the infall profile of the gas, however, no method exists, which allows a measurement of it. We consider the results of two cosmological simulations (using the numerical codes RAMSES and Gadget) and find that the gas and dark matter radial velocities are very similar. We derive the relevant dynamical equations, in particular the generalized hydrostatic equilibrium equation, including both the expansion of the Universe and the cosmological background. This generalized gas equation is the main new contribution of this paper. We combine these generalized equations with the results of the numerical simulations to estimate the contribution to the measured cluster masses from the radial velocity: inside the virial radius it is negligible, and inside two virial radii the effect is below 40%, in agreement the earlier analyses for DM. We point out how the infall velocity in principle may be observable, by measuring the gas properties to distance of about two virial radii, however, this is practically not possible today.
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.