No Arabic abstract
(Abridged) We apply a very general statistical theorem introduced by Cramer (1936) to study the origin of the deviations of the halo spin PDF from the reference lognormal shape. We find that these deviations originate from correlations between two quantities entering the definition of spin, namely the ratio $J/M^{5/2}$ (which depends only on mass) and the total gravitational binding energy $E$. To reach this conclusion, we have made usage of the results deduced from two high spatial- and mass resolution simulations. Our simulations cover a relatively small volume and produce a sample of more than 16.000 gravitationally bound halos, each traced by at least 300 particles. We verify that our results are stable to different systematics, by comparing our results with those derived by the GIF2 and by a more recent simulation performed by Maccio et al. We find that the spin probability distribution function shows systematic deviations from a lognormal, at all redshifts z <= 1. These deviations depend on mass and redshift: at small masses they change little with redshift, and also the best lognormal fits are more stable. The J-M relationship is well described by a power law of exponent $alpha$ very near to the linear theory prediction (alpha=5/3), but systematically lower than this at z<= 0.3. We argue that the fact that deviations from a lognormal PDF are present only for high-spin halos could point to a role of large-scale tidal fields in the evolution of the spin PDF.
We use a 200 $h^{-1}Mpc$ a side N-body simulation to study the mass accretion history (MAH) of dark matter halos to be accreted by larger halos, which we call infall halos. We define a quantity $a_{rm nf}equiv (1+z_{rm f})/(1+z_{rm peak})$ to characterize the MAH of infall halos, where $z_{rm peak}$ and $z_{rm f}$ are the accretion and formation redshifts, respectively. We find that, at given $z_{rm peak}$, their MAH is bimodal. Infall halos are dominated by a young population at high redshift and by an old population at low redshift. For the young population, the $a_{rm nf}$ distribution is narrow and peaks at about $1.2$, independent of $z_{rm peak}$, while for the old population, the peak position and width of the $a_{rm nf}$ distribution both increases with decreasing $z_{rm peak}$ and are both larger than those of the young population. This bimodal distribution is found to be closely connected to the two phases in the MAHs of halos. While members of the young population are still in the fast accretion phase at $z_{rm peak}$, those of the old population have already entered the slow accretion phase at $z_{rm peak}$. This bimodal distribution is not found for the whole halo population, nor is it seen in halo merger trees generated with the extended Press-Schechter formalism. The infall halo population at $z_{rm peak}$ are, on average, younger than the whole halo population of similar masses identified at the same redshift. We discuss the implications of our findings in connection to the bimodal color distribution of observed galaxies and to the link between central and satellite galaxies.
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 the probability distribution function (PDF) of relative velocity between two different dark matter halos (i.e. pairwise velocity) with a set of high-resolution cosmological $N$-body simulations. We investigate the pairwise velocity PDFs over a wide range of halo masses of $10^{12.5-15}, h^{-1}M_{odot}$ and redshifts of $0<z<1$. At a given set of masses, redshift and the separation length between two halos, our model requires three parameters to set the pairwise velocity PDF, whereas previous non-Gaussian models in the literature assume four or more free parameters. At the length scales of $r=5-40, [h^{-1}, mathrm{Mpc}]$, our model predicts the mean and dispersion of the pairwise velocity for dark matter halos with their masses of $10^{12.5-13.5} , [h^{-1}M_{odot}]$ at $0.3 < z < 1$ with a 5%-level precision, while the model precision reaches a 20% level (mostly a 10% level) for other masses and redshifts explored in the simulations. We demonstrate that our model of the pairwise velocity PDF provides an accurate mapping of the two-point clustering of massive-galaxy-sized halos at the scales of $O(10), h^{-1}mathrm{Mpc}$ between redshift and real space for a given real-space correlation function. For a mass-limited halo sample with their masses greater than $10^{13.5}, h^{-1}M_{odot}$ at $z=0.55$, our model can explain the monopole and quadropole moments of the redshift-space two-point correlations with a precision better than 5% at the scales of $5-40$ and $10-30, h^{-1}mathrm{Mpc}$, respectively. Our model of the pairwise velocity PDF will give a detailed explanation of statistics of massive galaxies at the intermediate scales in redshift surveys, including the non-linear redshift-space distortion effect in two-point correlation functions and the measurements of the kinematic Sunyaev-Zeldovich effect.
Measurements of the total amount of stars locked up in galaxies as a function of host halo mass contain key clues about the efficiency of processes that regulate star formation. We derive the total stellar mass fraction f_star as a function of halo mass M500c from z=0.2 to z=1 using two complementary methods. First, we derive f_star using a statistical Halo Occupation Distribution model jointly constrained by data from lensing, clustering, and the stellar mass function. This method enables us to probe f_star over a much wider halo mass range than with group or cluster catalogs. Second, we derive f_star at group scales using a COSMOS X-ray group catalog and we show that the two methods agree to within 30%. We quantify the systematic uncertainty on f_star using abundance matching methods and we show that the statistical uncertainty on f_star (~10%) is dwarfed by systematic uncertainties associated with stellar mass measurements (~45% excluding IMF uncertainties). Assuming a Chabrier IMF, we find 0.012<f_star<0.025 at M500c=10^13 Msun and 0.0057<f_star<0.015 at M500c=10^14 Msun. These values are significantly lower than previously published estimates. We investigate the cause of this difference and find that previous work has overestimated f_star due to a combination of inaccurate stellar mass estimators and/or because they have assumed that all galaxies in groups are early type galaxies with a constant mass-to-light ratio. Contrary to previous claims, our results suggest that the mean value of f_star is always significantly lower than f_gas for halos above 10^13 Msun. Combining our results with recently published gas mas fractions, we find a shortfall in f_star+f_gas at R500c compared to the cosmic mean. This shortfall varies with halo mass and becomes larger towards lower halos masses.
We revisit the question of what mechanism is responsible for the spins of halos of dark matter. The answer to this question is of high importance for modeling galaxy intrinsic alignment, which can potentially contaminate current and future lensing data. In particular, we show that when the dark matter halos pass nearly by each other in dense environments-- namely halo assemblies-- they swing and spin each other via exerting mutual tidal torques. We show that this has a significant contribution to the spin of dark matter halos comparable to that of calculated by the so-called tidal torque theory (TTT). We use the results of state-of-the-art simulation of Illutris to check the prediction of this theory against the simulation data.