No Arabic abstract
In the present work we describe the formalism necessary to derive the properties of dark matter halos beyond two virial radius using the spherical collapse model (without shell crossing), and provide the framework for the theoretical prediction presented in Prada et al. (2005). We show in detail how to obtain within this model the probability distribution for the spherically-averaged enclosed density at any radii P(delta,r). Using this probability distribution, we compute the most probable and mean density profiles, which turns out to differ considerably from each other. We also show how to obtain the typical profile, as well as the probability distribution and mean profile for the spherically averaged radial velocity. Two probability distributions are obtained: a first one is derived using a simple assumption, that is, if Q is the virial radius in Lagrangian coordinates, then the enclosed linear contrast delta_l(q,Q) must satisfy the condition that delta_l(q=Q) = delta_vir, where delta_vir is the linear density contrast within the virial radius Rvir at the moment of virialization. Then we introduce an additional constraint to obtain a more accurate P(delta,r) which reproduces to a higher degree of precision the distribution of the spherically averaged enclosed density found in the simulations. This new constraint is delta_l(q,Q) < delta_vir for all q > Q, which means that there are no radii larger than Rvir where the density contrast is larger than that used to define the virial radius. Finally, we compare in detail our theoretical predictions for the probability distributions with the results found in the simulations.
We present a test to quantify how well some approximate methods, designed to reproduce the mildly non-linear evolution of perturbations, are able to reproduce the clustering of DM halos once the grouping of particles into halos is defined and kept fixed. The following methods have been considered: Lagrangian Perturbation Theory (LPT) up to third order, Truncated LPT, Augmented LPT, MUSCLE and COLA. The test runs as follows: halos are defined by applying a friends-of-friends (FoF) halo finder to the output of an N-body simulation. The approximate methods are then applied to the same initial conditions of the simulation, producing for all particles displacements from their starting position and velocities. The position and velocity of each halo are computed by averaging over the particles that belong to that halo, according to the FoF halo finder. This procedure allows us to perform a well-posed test of how clustering of the matter density and halo density fields are recovered, without asking to the approximate method an accurate reconstruction of halos. We have considered the results at $z=0,0.5,1$, and we have analysed power spectrum in real and redshift space, object-by-object difference in position and velocity, density Probability Distribution Function (PDF) and its moments, phase difference of Fourier modes. We find that higher LPT orders are generally able to better reproduce the clustering of halos, while little or no improvement is found for the matter density field when going to 2LPT and 3LPT. Augmentation provides some improvement when coupled with 2LPT, while its effect is limited when coupled with 3LPT. Little improvement is brought by MUSCLE with respect to Augmentation. The more expensive particle-mesh code COLA outperforms all LPT methods [abridged]
Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10-100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.
We study properties of dark matter halos at high redshifts z=2-10 for a vast range of masses with the emphasis on dwarf halos with masses 10^7-10^9 Msun/h. We find that the density profiles of relaxed dwarf halos are well fitted by the NFW profile and do not have cores. We compute the halo mass function and the halo spin parameter distribution and find that the former is very well reproduced by the Sheth & Tormen model while the latter is well fitted by a lognormal distribution with lambda_0 = 0.042 and sigma_lambda = 0.63. We estimate the distribution of concentrations for halos in mass range that covers six orders of magnitude from 10^7 Msun/h to 10^13} Msun/h, and find that the data are well reproduced by the model of Bullock et al. The extrapolation of our results to z = 0 predicts that present-day isolated dwarf halos should have a very large median concentration of ~ 35. We measure the subhalo circular velocity functions for halos with masses that range from 4.6 x 10^9 Msun/h to 10^13 Msun/h and find that they are similar when normalized to the circular velocity of the parent halo. Dwarf halos studied in this paper are many orders of magnitude smaller than well-studied cluster- and Milky Way-sized halos. Yet, in all respects the dwarfs are just down-scal
We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.
The proper motions of stars in the outskirts of globular clusters are used to estimate cluster velocity dispersion profiles as far as possible within their tidal radii. We use individual color-magnitude diagrams to select high probability cluster stars for 25 metal-poor globular clusters within 20 kpc of the sun, 19 of which have substantial numbers of stars at large radii. Of the 19, 11 clusters have a falling velocity dispersion in the 3-6 half mass radii range, 6 are flat, and 2 plausibly have a rising velocity dispersion. The profiles are all in the range expected from simulated clusters started at high redshift in a zoom-in cosmological simulation. The 11 clusters with falling velocity dispersion profiles are consistent with no dark matter above the Galactic background. The 6 clusters with approximately flat velocity dispersion profiles could have local dark matter, but are ambiguous. The 2 clusters with rising velocity dispersion profiles are consistent with a remnant local dark matter halo, but need membership confirmation and detailed orbital modeling to further test these preliminary results.