No Arabic abstract
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.
In this work, we compare large scale structure observables for stellar mass selected samples at $z=0$, as predicted by two galaxy models, the hydrodynamical simulation IllustrisTNG and the Santa-Cruz semi-analytic model (SC-SAM). Although both models have been independently calibrated to match observations, rather than each other, we find good agreement between the two models for two-point clustering and galaxy assembly bias signatures. The models also show a qualitatively similar response of occupancy and clustering to secondary halo paramaters other than mass, such as formation history and concentration, although with some quantitative differences. Thus, our results demonstrate that the galaxy-halo relationships in SC-SAM and TNG are quite similar to first order. However, we also find areas in which the models differ. For example, we note a strong correlation between halo gas content and environment in TNG, which is lacking in the SC-SAM, as well as differences in the occupancy predictions for low-mass haloes. Moreover, we show that higher-order statistics, such as cumulants of the density field, help to accurately describe the galaxy distribution and discriminate between models that show degenerate behavior for two-point statistics. Our results suggest that SAMs are a promising cost-effective and intuitive method for generating mock catalogues for next generation cosmological surveys.
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.
We present a new empirical model for the mass assembly of dark matter halos. We approximate the growth of individual halos as a simple power-law function of time, where the power-law index smoothly decreases as the halo transitions from the fast-accretion regime at early times, to the slow-accretion regime at late times. Using large samples of halo merger trees taken from high-resolution cosmological simulations, we demonstrate that our 3-parameter model can approximate halo growth with a typical accuracy of 0.1 dex for t > 1 Gyr for all halos of present-day mass greater than 10^11Msun, including subhalos and host halos in gravity-only simulations, as well as in the TNG hydrodynamical simulation. We additionally present a new model for the assembly of halo populations, which not only reproduces average mass growth across time, but also faithfully captures the diversity with which halos assemble their mass. Our python implementation is based on the autodiff library JAX, and so our model self-consistently captures the mean and variance of halo mass accretion rate across cosmic time. We show that the connection between halo assembly and the large-scale density field, known as halo assembly bias, is accurately captured by our model, and that residual errors in our approximations to halo assembly history exhibit a negligible residual correlation with the density field. Our publicly available source code can be used to generate Monte Carlo realizations of cosmologically representative halo histories; our differentiable implementation facilitates the incorporation of our model into existing analytical halo model frameworks.
The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2times10^{12}M_{odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_sle 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5times10^{12}M_{odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to 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.