No Arabic abstract
We study the central dark matter (DM) cusp evolution in cosmological galactic halos. Models with and without baryons (baryons+DM, hereafter BDM model, and pure DM, PDM model, respectively) are advanced from identical initial conditions. The DM cusp properties are contrasted by a direct comparison of pure DM and baryonic models. We find a divergent evolution between the PDM and BDM models within the inner ~10 kpc region. The PDM model forms a R^{-1} cusp as expected, while the DM in the BDM model forms a larger isothermal cusp R^{-2} instead. The isothermal cusp is stable until z~1 when it gradually levels off. This leveling proceeds from inside out and the final density slope is shallower than -1 within the central 3 kpc (i.e., expected size of the R^{-1} cusp), tending to a flat core within ~2 kpc. This effect cannot be explained by a finite resolution of our code which produces only a 5% difference between the gravitationally softened force and the exact Newtonian force of point masses at 1 kpc from the center. Neither is it related to the energy feedback from stellar evolution or angular momentum transfer from the bar. Instead it can be associated with the action of DM+baryon subhalos heating up the cusp region via dynamical friction and forcing the DM in the cusp to flow out and to `cool down. The process described here is not limited to low z and can be efficient at intermediate and even high z.
Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $dlogrho/dlog r approx -1$, and some dispersion, which is generally parametrized by the varying index $alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $phi_0$, Einasto $alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.
We present an equilibrium statistical mechanical theory of collisionless self-gravitational systems with isotropic velocity distributions. Compared to existing standard theories, we introduce two changes: (1) the number of possible microstates is computed in energy (orbit) space rather than phase space and (2) low occupation numbers are treated more appropriately than using Stirlings approximation. Combined, the two modifications predict that the relaxed parts of collisionless self-gravitating systems, such as dark-matter halos, have a differential energy distribution N(E) ~ [exp(phi_0 - E) - 1], dubbed DARKexp. Such systems have central power-law density cusps rho(r) ~ r^-1, which suggests a statistical mechanical origin of cusps in simulated dark-matter halos.
The multicomponent dark matter model with self-scattering and inter-
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]
We describe an extensive FUSE survey of highly ionized oxygen in the vicinity of the Milky Way that serves as an example of the type of study that would be desirable for other galactic systems. Understanding the origin of hot gas in the vicinity of galaxies and its relationship to the intergalactic medium presents a major observational challenge. Ultraviolet absorption-line spectroscopy is currently the most direct means for comprehensive investigations of the gas in galactic environments, but even with present (and near-term) facilities the number of background objects available to probe nearby galaxy halos and low-redshift cosmological structures is limited. Studying these structures over a range of impact parameters and angular separations would provide fundamental information about the baryonic content of the hot gas, its physical conditions, and its origins. A large space telescope optimized for high resolution spectroscopy in the 900-3200 Angstrom wavelength region at a sensitivity sufficient to observe faint AGNs/QSOs at angular separations of <1 degree would be ideal for such studies.