Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTesting approximate predictions of displacements of cosmological dark matter halos
104
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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]
rate research
Read More
We explore the cosmological constraints on the parameter w_dm of the dark matter barotropic equation of state (EoS) to investigate the warmness of the dark matter fluid. The model is composed by the dark matter and dark energy fluids in addition to the radiation and baryon components. We constrain the values of w_dm using the latest cosmological observations that measure the expansion history of the Universe. When w_dm is estimated together with the parameter w_de of the barotropic EoS of dark energy we found that the cosmological data favor a value of w_dm = 0.006 +- 0.001, suggesting a -warm- dark matter, and w_de= -1.11 +- 0.03$ that corresponds to a phantom dark energy, instead of favoring a cold dark matter and a cosmological constant (w_dm = 0, w_de = -1). When w_dm is estimated alone but assuming w_de = -1, -1.1, -0.9, we found w_dm = 0.009 +- 0.002, 0.006 +- 0.002, 0.012 +- 0.002 respectively, where the errors are at 3 sigma (99.73%), i.e., w_dm > 0 with at least 99.73% of confidence level. When (w_dm, Omega_dm0) are constrained together, the best fit to data corresponds to (w_dm=0.005 +- 0.001, Omega_dm0 = 0.223 +- 0.008) and with the assumption of w_de = -1.1 instead of a cosmological constant (i.e., w_de = -1). With these results we found evidence of w_dm > 0 suggesting a -warm- dark matter, independent of the assumed value for w_{rm de}, but where values w_de < -1 are preferred by the observations instead of the cosmological constant. These constraints on w_dm are consistent with perturbative analyses done in previous works.
We compare the predicted conditional mass function (CMF) of dark matter halos from two theoretical prescriptions against numerical N-body simulations, both in overdense and underdense regions and at different Eulerian scales ranging from $5$ to $30,h^{-1},$Mpc. In particular, we consider in detail a locally-implemented rescaling of the unconditional mass function (UMF) already discussed in the literature, and also a generalization of the standard rescaling method described in the extended Press-Schechter formalism. First, we test the consistency of these two rescalings by verifying the normalization of the CMF at different scales, and showing that none of the proposed cases provides a normalized CMF. In order to satisfy the normalization condition, we include a modification in the rescaling procedure. After this modification, the resulting CMF generally provides a better description of numerical results. We finally present an analytical fit to the ratio between the CMF and the UMF (also known as the matter-to-halo bias function) in underdense regions, which could be of special interest to speed-up the computation of the halo abundance when studying void statistics. In this case, the CMF prescription based on the locally-implemented rescaling provides a slightly better description of the numerical results when compared to the standard rescaling.
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.
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.
The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the splashback shell formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Spheroidal Halos) and demonstrate that the code identifies splashback shells correctly and measures their properties with an accuracy of $<5%$ for halos with more than 50,000 particles and mass accretion rates of $Gamma_textrm{DK14}>0.5$. Using SHELLFISH, we present the first estimates for several basic properties of individual splashback shells, such as radius, $R_textrm{sp}$, mass, and overdensity, and provide fits to the distribution of these quantities as functions of $Gamma_textrm{DK14}$, $ u_textrm{200m}$, and $z.$ We confirm previous findings that $R_textrm{sp}$ decreases with increasing $Gamma_textrm{DK14}$, but we show that independent of accretion rate, it also decreases with increasing $ u_textrm{200m}$. We also study the 3D structures of these shells and find that they generally have non-ellipsoidal oval shapes. We find that splashback radii estimated by SHELLFISH are $20%-30%$ larger than those estimated in previous studies from stacked density profiles at high accretion rates. We demonstrate that the latter are biased low due to the contribution of high-mass subhalos to these profiles and show that using the median instead of mean density in each radial bin mitigates the effect of substructure on density profiles and removes the bias.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
