No Arabic abstract
In the present paper, we improve the Extended Secondary Infall Model (ESIM) of Williams et al. (2004) to obtain further insights on the cusp/core problem. The model takes into account the effect of ordered and random angular momentum, dynamical friction and baryon adiabatic contraction in order to obtain a secondary infall model more close to the collapse reality. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on cluster of galaxies scales. The results obtained suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The NFW profile can be reobtained, in our model only if the system is constituted just by dark matter and the magnitude of angular momentum and dynamical friction are reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from de Blok & Bosma (2002) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. On scales smaller than $simeq 10^{11} h^{-1} M_{odot}$ the slope $alpha simeq 0$ and on cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $alpha simeq 0.6$ for a cluster of $simeq 10^{14} h^{-1} M_{odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.
In this paper, in the framework of the secondary infall model, the correlation between the central surface density and the halo core radius of galaxy, and cluster of galaxies, dark matter haloes was analyzed, this having recently been studied on a wide range of scales. We used Del Popolo (2009) secondary infall model taking into account ordered and random angular momentum, dynamical friction, and dark matter (DM) adiabatic contraction to calculate the density profile of haloes, and then these profiles are used to determine the surface density of DM haloes. The main result is that $r_ast$ (the halo characteristic radius) is not an universal quantity as claimed by Donato et al. (2009) and Gentile et al. (2009). On the contrary, we find a correlation with the halo mass $M_{200}$ in agreement with Cardone & Tortora (2010), Boyarsky at al. (2009) and Napolitano et al. (2010), but with a significantly smaller scatter, namely $0.16 pm 0.05$. We also consider the baryon column density finding this latter being indeed a constant for low mass systems such as dwarfs, but correlating with mass with a slope $alpha= 0.18 pm 0.05$. In the case of the surface density of dark matter for a system composed only of dark matter, as in dissipationless simulations, we get $alpha=0.20 pm 0.05$. These results leave little room for the recently claimed universality of (dark and stellar) column density.
(Abridged) We study the outer density profiles of dark matter haloes predicted by a generalized secondary infall model and observed in a N-body cosmological simulation of a Lambda CDM model. We find substantial systematic variations in shapes and concentrations of the halo profiles as well as a strong correlation of the profiles with the environment. In the N-body simulation, the average outer slope of the density profiles, beta (rhopropto r^{-beta}), of isolated haloes is approx 2.9; 68% of these haloes have values of beta between 2.5 and 3.8. Haloes in dense environments of clusters are more concentrated and exhibit a broad distribution of beta with values larger than for isolated haloes . Contrary to what one may expect, the haloes contained within groups and galaxy systems are less concentrated and have flatter outer density profiles than the isolated haloes. The concentration decreases with M_h, but its scatter for a given mass is substantial. The mass and circular velocity of the haloes are strongly correlated: M_h propto V_m^{alpha} with alpha ~ 3.3 (isolated) and ~3.5 (haloes in clusters). For M_h=10^12M_sun the rms deviations from these relations are Delta logM_h=0.12 and 0.18, respectively. Approximately 30% of the haloes are contained within larger haloes or have massive companions (larger than ~0.3 the mass of the current halo) within 3 virial radii. The remaining 70% of the haloes are isolated objects. The distribution of beta as well as the concentration-mass and M_h-V_m relations for the isolated haloes agree very well with the predictions of our seminumerical approach which is based on a generalization of the secondary infall model and on the extended Press-Schechter formalism.
We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html
We use N-body simulations of dark matter haloes in cold dark matter (CDM) and a large set of different warm dark matter (WDM) cosmologies to demonstrate that the spherically averaged density profile of dark matter haloes has a shape that depends on the power spectrum of matter perturbations. Density profiles are steeper in WDM but become shallower at scales less than one percent of the virial radius. Virialization isotropizes the velocity dispersion in the inner regions of the halo but does not erase the memory of the initial conditions in phase space. The location of the observed deviations from CDM in the density profile and in phase space can be directly related to the ratio between the halo mass and the filtering mass and are most evident in small mass haloes, even for a 34 keV thermal relic WDM. The rearrangement of mass within the haloes supports analytic models of halo structure that include angular momentum. We also find evidence of a dependence of the slope of the inner density profile in CDM cosmologies on the halo mass with more massive haloes exhibiting steeper profiles, in agreement with the model predictions and with previous simulation results. Our work complements recent studies of microhaloes near the filtering scale in CDM and strongly argue against a universal shape for the density profile.
We present N-body simulations of a new class of self-interacting dark matter models, which do not violate any astrophysical constraints due to a non-power-law velocity dependence of the transfer cross section which is motivated by a Yukawa-like new gauge boson interaction. Specifically, we focus on the formation of a Milky Way-like dark matter halo taken from the Aquarius project and re-simulate it for a couple of representative cases in the allowed parameter space of this new model. We find that for these cases, the main halo only develops a small core (~1 kpc) followed by a density profile identical to that of the standard cold dark matter scenario outside of that radius. Neither the subhalo mass function nor the radial number density of subhaloes are altered in these models but there is a significant change in the inner density structure of subhaloes resulting in the formation of a large density core. As a consequence, the inner circular velocity profiles of the most massive subhaloes differ significantly from the cold dark matter predictions and we demonstrate that they are compatible with the observational data of the brightest Milky Way dSphs in such a velocity-dependent self-interacting dark matter scenario. Specifically, and contrary to the cold dark matter case, there are no subhaloes that are more concentrated than what is inferred from the kinematics of the Milky Way dSphs. We conclude that these models offer an interesting alternative to the cold dark matter model that can reduce the recently reported tension between the brightest Milky Way satellites and the dense subhaloes found in cold dark matter simulations.