No Arabic abstract
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.
We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height mass parameter, $ u(M,z)=delta_{rm crit}(z)/sigma(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.
We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r) profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschingers spherical secondary-infall similarity solution. In this model, Q propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius.
We show that hidden hot dark matter, hidden-sector dark matter with interactions that decouple when it is relativistic, is a viable dark matter candidate provided it has never been in thermal equilibrium with the particles of the standard model. This hidden hot dark matter may reheat to a lower temperature and number density than the visible Universe and thus account, simply with its thermal abundance, for all the dark matter in the Universe while evading the typical constraints on hot dark matter arising from structure formation. We find masses ranging from ~3 keV to ~10 TeV. While never in equilibrium with the standard model, this class of models may have unique observational signatures in the matter power spectrum or via extra-weak interactions with standard model particles.
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.
A self-interacting dark matter halo can experience gravothermal collapse, resulting in a central core with an ultrahigh density. It can further contract and collapse into a black hole, a mechanism proposed to explain the origin of supermassive black holes. We study dynamical instability of the core in general relativity. We use a truncated Maxwell-Boltzmann distribution to model the dark matter distribution and solve the Tolman-Oppenheimer-Volkoff equation. For given model parameters, we obtain a series of equilibrium configurations and examine their dynamical instability based on considerations of total energy, binding energy, fractional binding energy, and adiabatic index. The numerical results from our semi-analytical method are in good agreement with those from fully relativistic N-body simulations. We further show for the instability to occur in the classical regime, the boundary temperature of the core should be at least $10%$ of the mass of dark matter particles; for a $10^9~{rm M_odot}$ seed black hole, the particle mass needs to be larger than a few keV. These results can be used to constrain different collapse models, in particular, those with dissipative dark matter interactions.