EDP Sciences
Free Access
Issue
A&A
Volume 502, Number 3, August II 2009
Page(s) 733 - 747
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/200811404
Published online 22 June 2009

Density profiles of dark matter haloes on galactic and cluster scales

A. Del Popolo1,2,3 - P. Kroupa2

1 - Dipartimento di Fisica e Astronomia, Universitá di Catania, Viale Andrea Doria 6, 95125 Catania, Italy
2 - Argelander-Institut für Astronomie, Auf dem Hügel 71, 53121 Bonn, Germany
3 - Istanbul Technical University, Ayazaga Campus, Faculty of Science and Letters, 34469 Maslak/ISTANBUL, Turkey

Received 24 November 2008 / Accepted 7 May 2009

Abstract
Aims. In the present paper, we improve the ``extended secondary infall model'' (ESIM) of Williams and collaborators to obtain further insights into the cusp/core problem.
Methods. A secondary infall model close to the collapse reality is obtained by simultaneously taking into account effects that till now have been studied separately, namely ordered and random angular momentum, dynamical friction, and baryon adiabatic contraction. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on galaxy cluster scales.
Results. Our results imply that angular momentum and dynamical friction are able, on galactic scales, of overcoming the competing effect of adiabatic contraction and eliminating the Cusp. The NFW profile is not the standard outcome of the model, and it can be recovered in our model only if the system consists entirely of dark matter and the magnitude of angular momentum and dynamical friction are lower than the values predicted by the model itself. Comparison of the rotation curves of LSB galaxies with the results of our model are in good agreement. On scales smaller than $\simeq$ $ 10^{11}~ h^{-1} ~M_{\odot}$, the slope is $\alpha \simeq 0$, and on cluster scales we observe a similar evolution in the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.7$ for a cluster of $\simeq$ $10^{14}~ h^{-1}~ M_{\odot}$. The total mass profile differs from that of dark matter showing a central cusp that is reproduced by a NFW model.

Key words: cosmology: large-scale structure of Universe - cosmology: theory - galaxies: formation

1 Introduction

The structure of dark matter haloes is of fundamental importance on the study of the formation and evolution of galaxies and clusters of galaxies. At the simplest level, dark matter haloes form when, in the early Universe, the matter within (and surrounding) an overdense region experiences deceleration because of the gravitational force, decouples from the Hubble flow, collapses, and eventually, virializes. From the theoretical point of view, the structure of dark matter haloes can be studied both analytically and numerically. A significant part of the analytical work completed so far was based on the secondary infall model (SIM) introduced by Gunn & Gott (1972), Gott (1975), and Gunn (1977). Calculations based on this model predict that the density profile of the virialized halo should scale as $\rho \propto r^{-9/4}$. Self-similar solutions were found by Fillmore & Goldreich (1984) and Bertschinger (1985), while Hoffman & Shaham (1985, hereafter HS) studied density profiles around density peaks. More recently, modifications of the self-similar collapse model to include more realistic dynamics of the growth process have been proposed (e.g., Avila-Reese et al. 1998; Nusser & Sheth 1999; Henriksen & Widrow 1999; Subramanian et al. 2000; Del Popolo et al. 2000, hereafter DP2000). Ryden & Gunn (1987, hereafter RG87) were the first to relax the assumption of purely radial self-similar collapse by including non-radial motions arising from secondary perturbations in the halo. Numerous authors have emphasized the effect of an isotropic velocity dispersion (thus of non-radial motion) in the core of collisionless haloes. One common result of the previous studies (RG87; White & Zaritsky 1992; Avila-Reese et al. 1998; Hiotelis 2002; Nusser 2001; Le Delliou & Enriksen 2003; Ascasibar et al. 2004; Williams et al. 2004) is that a larger amount of angular momentum leads to shallower final density profiles in the inner region of the halo. Moreover, baryons have been invoked both to shallow (El-Zant et al. 2001, hereafter EZ01, 2003; Romano-Diaz et al. 2008) and to steepen (Blumenthal et al. 1986) the dark matter profile.

As previously reported, the structure of dark matter haloes can also be studied by numerical simulations. Quinn et al. (1986) pioneered the use of N-body simulations to study halo formation and confirmed HS results. More recent studies (Dubinski & Carlberg 1991; Lemson 1995; Cole & Lacey 1996; Navarro et al. 1996, 1997 (NFW); Moore et al. 1998; Jing & Suto 2000; Klypin et al. 2001; Bullock et al. 2001; Power et al. 2003; and Navarro et al. 2004) found that although the spherically-averaged density profiles of the N-body dark matter haloes are similar, regardless of the mass of the halo or the cosmological model, their profiles differ significantly from the single power laws predicted by the theoretical studies. The N-body profiles are characterized by an r-3 decline at large radii and a cuspy profile of the form $\rho(r)\propto r^{-\alpha}$, where $\alpha < 2$ close to the center. The true value of the inner density slope $\alpha$ is a matter of controversy, with NFW suggesting $\alpha=1$, but with Moore et al. (1998), Ghigna et al. (2000), and Fukushige & Makino (2001) arguing for $\alpha=1.5$, while Jing & Suto (2000) and Klypin et al. (2001) claimed that the true value of $\alpha$may depend on halo mass, merger history, and substructure. Power et al. (2003) pointed out that the logarithmic slope becomes increasingly shallow inwards, with little sign of approaching an asymptotic value at the resolved radii. In that case, the precise value of $\alpha$, at a given cut-off scale, would not be particularly meaningful. This result was later confirmed by Hayashi et al. (2003) and Fukushige et al. (2004), and it is predicted by several analytical models (e.g., Taylor & Navarro 2001; Hoeft et al. 2003). Finally, Navarro et al. (2004) proposed a new fitting formula with a logarithmic slope that decreases inward more gradually than the NFW profile.

While numerical simulations universally produce a cuspy density profile, observed rotation curves of dwarf spiral and LSB galaxies seem to indicate that the shape of the density profile on small scales is significantly shallower than found in numerical simulations (Flores & Primak 1994; Moore 1994; Burkert 1995; Kravtsov et al. 1998; Salucci & Burkert 2000; Borriello & Salucci 2001; de Blok et al. 2001; de Blok & Bosma 2002; Marchesini et al. 2002; de Blok 2003; de Blok et al. 2003). It seems that the data generally favor logarithmic density slopes close to 0.2 (de Blok 2003; de Blok et al. 2003; Spekkens et al. 2005).

On cluster scales, X-ray analyses have led to wide ranging results, from $\alpha=0.6$ (Ettori et al. 2002) to $\alpha=1.2$ (Lewis et al. 2003) or even $\alpha=1.9$ (Arabadjis et al. 2002). Ricotti's (2003) N-body simulations suggest that the density profile of DM haloes is not universal (in agreement with Jing & Suto 2000; Subramanian et al. 2000; Simon et al. 2003b; Cen et al. 2004; Ricotti & Wilkinson 2004; Ricotti et al. 2007), but that there are instead shallower cores in dwarf galaxies and steeper cores in clusters.

The discrepancy between simulations and observations has often been accused of representing a genuine crisis of the CDM scenario and has become known as the ``cusp/core'' problem. Since LSB galaxies are thought to be ideal for comparing with theory, because their dynamics are dominated by dark matter with little contribution from baryons (Bothun et al. 1997), the discrepancy with simulations is particularly troublesome. The significance of this disagreement, however, remains controversial and different solutions have been proposed. A number of authors attribute the problem to a real failure of the CDM model, or to that of simulations (de Blok et al. 2001a, 2001b; Borriello & Salucci 2001; de Blok et al. 2003). This has led to suggestions that dark matter properties may deviate from standard CDM and several alternatives have been suggested, such as warm (Colin et al. 2000; Sommer-Larsen & Dolgov 2001), repulsive (Goodman 2000), fluid (Peebles 2000), fuzzy (Hu et al. 2000), decaying (Cen 2001), annihilating (Kaplinghat et al. 2000), or self-interacting (Spergel & Steinhardt 2000; Yoshida et al. 2000; Davé et al. 2001) dark matter. Others argue that the inconsistency may reflect the finite resolution of the observations that has not been properly accounted for in the analysis of the HI rotation curves (van den Bosch et al.  2000; van den Bosch & Swaters 2001; Rhee et al. 2004). Alternatively, it has been suggested that stellar feedback from the first generation of stars formed in galaxies was so efficient that the remaining gas was expelled on a timescale comparable to, or less than, the local dynamical timescale. The dark matter subsequently adjusted to form an approximately constant density core (e.g., Gelato & Sommer-Larsen 1999). This is however unlikely to affect cluster cusps.

Simulations, observations, and semi-analytic models agree about the outer parts of the haloes' structure, and the disagreement in the inner regions could be connected to limits in numerical simulations (de Blok 2003; Taylor et al. 2004) or dissipationless N-body simulations being unable to account for the effects of baryons on dark matter evolution. Interestingly, the amount of central substructure seen in the semi-analytic haloes is consistent with the amount of substructure inferred from strong lensing experiments (Taylor et al. 2004). Thus, the semi-analytic haloes might provide a more accurate picture of the spatial distribution of substructure around galaxies even if an analytical method, no matter how sophisticated, will never be able to capture the full extent of complexity of a non-linear process.

Nevertheless the large amount of work carried out by many researchers to date using N-body simulations has met with limited success in elucidating the physics of halo formation. The reason is that the point of force of numerical simulations (namely to capture the full extent of the complexity of a non-linear process) is also their weakness: numerical simulations yield little physical insight beyond empirical findings precisely because they are so rich in dynamical processes, which are hard to disentangle and interpret in terms of the underlying physics. Analytical and semi-analytical models are more flexible than N-body simulations (see Williams et al. 2004). So even if analytical models such as SIM consider collapse and virialization of halos that are spherically symmetric, that have suffered no major mergers, and that have suffered quiescent accretion, their investigation can provide important results. One of the most widely used type of semi-analytical model is SIM, whose most often questioned assumptions are spherical symmetry and the absence of peculiar velocities (non-radial motions): in the ``real'' collapse, accretion does not happen in spherical shells but by aggregation of subclumps of matter that have already collapsed and a large fraction of observed clusters of galaxies exhibit significant substructure (Kriessler et al. 1995). Motions are not purely radial, especially when the perturbation detaches from the general expansion. Nevertheless, the SIM provides good results in describing the formation of dark matter haloes, because in energy space the collapse is ordered and gentle, differently from the chaotic collapse that is seen in N-body simulations (Zaroubi et al. 1996). This is confirmed by other studies (Tóth & Ostriker 1992; Huss et al. 1999a,b; Moore et al. 1999). We should also add that analytical and semi-analytical models have some advantages over N-body simulations: a) they are flexible because one can study the effects of physical processes one at a time; b) one can incorporate many physical effects at least in a schematic manner; c) they are computationally efficient (it takes about 10 s to compute the density profile of a given object at a given epoch on a desktop PC, Ascasibar et al. 2007).

In this paper, we present an analytical model for haloes formation based on SIM. The model is an extension of the Williams et al. (2004) model that take account of ordered angular momentum, dynamical friction, and adiabatic contraction. This model derives the initial shape of proto-haloes from the fluctuation spectrum at high redshifts, and haloes are endowed with secondary perturbations, which impart random motions to halo particles, as in Williams et al. (2004), and ordered angular momentum. The statistical properties of the secondary perturbations are derived from the same fluctuation power spectrum, and therefore their effects on particle random velocities are treated self-consistently. The ordered angular momentum acquired by the proto-structures is obtained from the tidal interaction of the proto-structure with the neighboring ones and using the theory of random fields. Dynamical friction is calculated according to Kandrup's (1980) approach. Adiabatic contraction is taken into account using the Blumenthal et al. (1986) model as modified by Gnedin et al. (2004).

The plan of the paper is the following: in Sect. 2, we introduce the model that shall be used to calculate the density profiles. In Sect. 3, we show the results of the model, while Sect. 4 is devoted to conclusions.

2 Model

As shown by the spherical collapse model[*] of Gunn & Gott (1972), a bound mass shell of initial comoving radius xi will expand to a maximum radius $x_{\rm m}$ (named apapsis or turnaround radius $x_{\rm ta}$). As successive shells expand to their maximum radius, they acquire angular momentum and then contract on orbits determined by the angular momentum. Dissipative physics and the process of violent relaxation will eventually intervene and convert the kinetic energy of collapse into random motions (virialization). Following RG87 and R88a,b, we restrict ourselves to haloes that form primarily via nearly-smooth accretion of matter, treating collapse and virialization of haloes that are spherically symmetric, that have suffered no major mergers, and that have experienced only quiescent accretion of somewhat lumpy material. The most important feature of the original RG87 method is that the dynamical evolution is carried out while conserving angular and radial momenta of individual halo shells. The initial shape of the proto-haloes is derived from the fluctuation spectrum at high redshift, and halos are endowed with secondary perturbations, which impart random motions to halo particles. The statistical properties of the secondary perturbations are derived from the same fluctuation power spectrum of the primary perturbations, and therefore their effects on particle random velocities are treated self-consistently. Ordered angular momentum is also treated self-consistently as shown in R88a.

It is possible to subdivide the halo into four regions, starting from the outside. In region (1), furthest from the center of the initial density peak, dark matter particles begin to feel the gravitational tug of the central peak and start to fall behind the Hubble flow. Further, in region (2), the particles start to decouple from the Hubble flow and begin to collapse. In region (3), the central density peak dominates the motion of particles; this is the region of infall and shell crossing. Finally, in region (4), the central part of the density peak, virialization occurs, or has already been reached. In region (4), virialization is well underway or is complete. Many analytical solutions in regions (1)-(3) ignore the non-radial motions of the particles; all orbits are assumed to be purely radial. The calculations carried out in RG87 demonstrate that non-radial motions are very important in determining the outcome in this region. Non-radial motions provide angular momentum to particles, which prevents them from penetrating to the very center of the halo, as confirmed by e.g., Avila-Reese et al. (1998), Subramanian et al. (2000), and Hozumi et al. (2000).

In most popular cosmological scenarios, the density field soon after recombination can be represented by a Gaussian random field. High density contrast peaks in the field will eventually achieve overdensities of order 1 and enter a non-linear stage of evolution. These peaks will then collapse to form bound structures. We start with one of these peaks, and, for simplicity assume that it is spherically symmetric. The peak is divided into a very small central core and many spherically symmetric concentric mass shells, each labeled by its initial comoving distance from the center, x.

The first step in obtaining the density profile, is to calculate the initial density profile produced by a primordial fluctuation, $\overline \delta_i (x_i)$. A well known result is the expression for the radial density profile of a fluctuation centered on a primordial peak of arbitrary height $\nu$

$\displaystyle \langle \delta (x) \rangle = \delta_0(x)=\frac{\nu \xi (x)}{\xi (...
...mma ^2\xi (x)+\frac{%
R_{\ast }^2}3\nabla ^2\xi(x) \right] \cdot \xi (0)^{-1/2}$     (1)

(Bardeen et al. 1986; RG87), where x is the comoving separation, $\nu=\delta(0)/\sigma $[*] is the height of a density peak, $\xi (r)$ is the two-point correlation function

\begin{displaymath}\xi(x)= \frac{1}{2 \pi^2 x} \int_0^{\infty} P(k) k \sin(k x) {\rm d} k ,
\end{displaymath} (2)

$\gamma $ and $R_{\ast}$ are two spectral parameters (see BBKS), and $\vartheta (\nu \gamma ,\gamma )$ is defined in BBKS.

The CDM spectrum used in this paper is that of BBKS (their equation (G3)), which has the transfer function

$\displaystyle T(k) = \frac{\ln \left( 1+2.34 q\right)}{2.34 q}
\left[1+3.89q+
(16.1 q)^2+(5.46 q)^3+(6.71)^4\right]^{-1/4}$     (3)

where $q=\frac{k\theta^{1/2}}{\Omega_{\rm X} h^2 {\rm Mpc^{-1}}}$. Here $\theta=\rho_{\rm er}/(1.68 \rho_{\rm\gamma})$represents the ratio of the energy density in relativistic particles to that in photons ($\theta=1$ corresponds to photons and three flavors of relativistic neutrinos). The spectrum is connected to the transfer function by the equation

\begin{displaymath}P(k)=P_{\rm CDM} {\rm e}^{-1/2 k^2 R_f^2},
\end{displaymath} (4)

where Rf is the smoothing (filtering) scale and $P_{\rm CDM}$ is given by

\begin{displaymath}P_{\rm CDM}= A k T^2(k)
\end{displaymath} (5)

being A is the normalization constant. We normalized the spectrum by assuming that the mass variance of the density field

\begin{displaymath}\sigma^2(M)=\frac{1}{2 \pi^2} \int_0^\infty {\rm d}k k^2 P(k) W^2(kR)
\end{displaymath} (6)

convolved with the top hat window

\begin{displaymath}W(kR)=\frac{3}{(kR)^3} (\sin kR-kR \cos kR)
\end{displaymath} (7)

of radius 8 h-1 Mpc-1 is $\sigma _{8}=0.76$ (Romano-Diaz et al. 2008). Throughout the paper, we adopt a $\Lambda$CDM cosmology with WMAP3 parameters, $\Omega_{\rm m}=1-\Omega_{\Lambda}=0.24$, $\Omega_{\Lambda}=0.76$, $\Omega_b=0.043$, and h=0.73, where h is the Hubble constant in units of 100 km s-1 Mpc-1. The density profile of an initial, pre-collapse halo, at early times, $\delta_i(x)$, is related to $\delta_0(x)$ by means of the linear growth factor D(z) (Peebles 1980) given by:

\begin{displaymath}\delta_i(x)=\delta_0(x)/D(z_i).
\end{displaymath} (8)

Let us divide the main smooth spherically symmetric halo, $\delta(x)$, into many concentric mass shells. Each shell is uniquely labeled by x, its initial comoving radius. The halo density differs from the average background density, i.e., density interior to any shell is greater than critical at all times. Each shell evolution is divided into two stages. Initially a shell expands with the Hubble flow, but with a slight deceleration arising from the central mass concentration. Eventually the shells outward radial velocity decreases to zero, after which the shell collapses, by some distance, back towards the center of the halo. The dividing moment is called the turnaround, and at any given time corresponds to the line dividing regions (1) and (2). In a halo with a declining density profile, the turnaround happens at progressively later cosmic times for shells at greater distances from the center.

In such cases, the time evolution of a shell until reaching turnaround is given by a set of parametric equations (Gunn & Gott 1972; Peebles 1980)

\begin{displaymath}r(\theta)={1\over 2}x{\bar\delta_0}^{-1}(1-\cos\theta),
\end{displaymath} (9)


\begin{displaymath}t(\theta)={3\over 4}t_0{\bar\delta_0}^{-3/2}(\theta-\sin\theta),
\end{displaymath} (10)

where r is particle's proper radius, t is cosmic time, and $\bar\delta_0(x)$is the initial average fractional density excess inside the shell (Eq. (13) of RG87),

\begin{displaymath}\bar\delta_0(x)={3\over x^2}\; \int_0^x \delta_0(y)\;y^2\;{\rm d}y,
\end{displaymath} (11)

and t0 is the present time. The mass of a shell, and the mass within a shell are constant; these relations, together with Eqs. (9) and (10) can be used to compute the fractional overdensity for any shell, parameterized by its $\delta_0/{\bar\delta_0}$, at a cosmic time corresponding to $\theta$, (Eq. (23) of RG87)
$\displaystyle \delta(\theta)+1 = {{9(\theta-\sin\theta)^2}\over{2(1-\cos\theta)...
...\sin\theta(\theta-\sin\theta)}\over{2(1-\cos\theta)^2}}\Bigr]
\Biggr)^{-1}\cdot$     (12)

This can be used to construct the density run of a halo at a fixed cosmic time t/t0.

The mean density distribution about a peak is spherically symmetric. However, the initial density peak will in general have a triaxial shape, leading to non-spherical collapse. As shown by R88a, since the quadrupole moment of the protogalaxy is largely due to the outermost shells, where triaxiality is insignificant, it is justifiable to ignore the triaxility in computations. In other terms, the dynamics of the halo collapse are dictated by the potential, which, being a double integral over all space, is much rounder than the mass distribution. Therefore, the effects of intrinsic triaxiality of initial density peaks are smaller than those due to the secondary perturbations, and so can be ignored. Furthermore, initial triaxiality is less severe in larger, 2-4$\sigma$ peaks (BBKS), which are the subject of the present study. With this caveat the smooth part of our density peak is still described by Eq. (1).

Moreover, in reality, the initial density peak will not be smooth, but contain many smaller scale positive and negative perturbations that originate in the same Gaussian random field producing the main peak. These secondary perturbations will perturb the motion of the dark matter particles from their otherwise purely radial orbits.

So, in order to investigate effects such as tidal torques and non-radial collapse, it is necessary to consider the non-spherical portion of the density distribution. In addition to the smooth halo, RG87 and Williams et al. (2004) considered contributions from the secondary perturbations produced by the same Gaussian random field that created the main halo. The overall initial density profile, evolved linearly to the present day, can be written as

\begin{displaymath}\rho(\vec{ x})=\rho_0[1+\delta_0(x)][1+\epsilon_0(\vec{ x})],
\end{displaymath} (13)

where $\rho_0$ is the present-day background density, the density excess due to the main halo is $\delta_0(x)$, and is assumed to be spherically symmetric, and $\epsilon_0(\vec{ x})$ is the density excess caused by the random secondary perturbations. In a statistical sense, the growth rate depends on $\delta_0/{\bar\delta_0}$, or, equivalently, x, either of which can be used to parameterize the strength of the tidal field in a spherically symmetric halo. The final expression for the growth rate is given by (Eq. (28) of RG87 and Williams et al. 2004)

\begin{displaymath}\epsilon(x,\theta)={\epsilon_0(x)\over{\bar\delta_0}}
{{f_2(\...
...}\over{f_1(\theta)-[\delta_0(x)/\bar\delta_0(x)]f_2(\theta)}},
\end{displaymath} (14)

where $\epsilon_0(x)$ is the amplitude of the initial perturbation, given by

\begin{displaymath}\epsilon_0(\vec{ x})={1\over{(2\pi)^3}}
\int {\rm d}^3k\; \e...
...uad \quad
\langle \vert\epsilon_{\vec{ }}\vert^2 \rangle =P(k)
\end{displaymath} (15)

and f1 and f2 are functions of the ``time'' parameter $\theta$: $f_1(\theta)=16{-}16\cos\theta+\sin^2\theta-9\theta\sin\theta$, and $f_2(\theta)=12{-}12\cos\theta+3\sin^2\theta-9\theta\sin\theta$.

At any given time and proper position, the random acceleration due to perturbation field $\epsilon(\vec{ r})$ is given by (RG87 Eq. (40))

                       $\displaystyle \vec{ g}(\vec{ r},t)$ = $\displaystyle \vec{ g}_{\rm tot}(\vec{ r},t)-\vec{ g}_b(\vec{ r},t)$  
  $\textstyle \approx$ $\displaystyle G\int {\rm d}^3r^\prime
{{\rho_b (\vec{ r}^\prime,t)\epsilon(\vec{ r},t)}\over
{\vert\vec{ r}^\prime-\vec{ r}\vert^3}}(\vec{ r}^\prime-\vec{ r}),$ (16)

which is an integral over all space, and $\rho_b(\vec{ r},t)=\rho_0 (t) [1+\delta(\vec{ r},t)]$ is the background density due to the main halo, which is related to the total density given by Eq. (13) by, $\rho(\vec{ r},t)=\rho_b(\vec{ r},t)[1+\epsilon(\vec{ r},t)]$, and $\rho_0 (t)$ is the average density of the Universe at epoch t. Since the density distribution in the main peak and perturbation field and growth rate of perturbations are functions of radial position only, we use x instead of $\vec{ x}$, and r instead of $\vec{ r}$. Because we aim to determine the rms value of acceleration, we will replace g with vector g . A major simplification is accomplished by decoupling the time dependence of acceleration, i.e., the rate of growth of acceleration, from its spatial variation. Let the initial acceleration field due to secondary perturbations be denoted by g0(x). Then the dimensionless rate of growth of acceleration is given by,

Fg(x,t)=g(x,t)/g0(x)=g(r,t)/g0(r). (17)

With these, the proper displacement of a particle can be evaluated as
                            dp(x,t) = $\displaystyle \int_0^t {\rm d}t_1 \int_0^{t_1} {\rm d}t_2\; g(x,t_2)$  
  = $\displaystyle g_0(x) \int_0^t\! {\rm d}t_1
\int_0^{t_1}\! {\rm d}t_2\; F_g(x,t_2)= \Delta g_0(x) t_0^2 F_r(x,t).$ (18)

We now describe the two functions, Fg(x,t) and g0(x), separately. The acceleration growth rate becomes (RG87 Eq. (44)[*])
                          $\displaystyle F_g(x,\theta)$ = $\displaystyle {g(x,\theta)\over g_0(x)}=
{{\rho_0[1+\delta(x,\theta)]\epsilon(x,\theta)r(x,\theta)}
\over{\rho_0\epsilon_0 x}}$  
  = $\displaystyle 8\bar\delta_0{{f_2(\theta)}\over
{[f_1(\theta)-{{\delta_0(x)}\over{\bar\delta_0(x)}}\;f_2(\theta)]^2}}\cdot$ (19)

The spatial dependence of acceleration is given by

\begin{displaymath}\Delta g_0(d,x)=
4 G \rho_0
\Bigl[
\int P_{\rm s}(k){\rm e}^{-kd}\left(1-{{\sin kx}\over{kx}}\right){\rm d}k\Bigr]^{1/2}.
\end{displaymath} (20)

The comoving displacement is given by

\begin{displaymath}d(x,t)=d_{\rm p}(x,t)/a(t)=\Delta g_0(x) t_0^2 F_r(x,t)/a(t).
\end{displaymath} (21)

In the early part of the halo evolution most of the shells are still expanding. During this time, secondary perturbations grow, and so does the acceleration, the velocity, and the displacement caused by these perturbations to the particles in the shell. Because the secondary peaks are randomly distributed within the halo, they displace the dark matter particles in random directions from their original positions. This can be visualized as a shell with an internal velocity dispersion, resulting in a ``swollen'' shell. If we concentrate on a single particle, its orbit, viewed from the rest frame of the parent shell, will oscillate between an inner and an outer radius of that shell (i.e., pericenter, $r_{\rm p}$ and apocenter, $r_{\rm a}$). In general, this orbit will not be closed and will resemble a rosette. In a real situation, the ``swelling'' of any given shell will gradually increase as the shell expands away from the center of the halo, and the influence of the secondary peaks grows, but for simplicity, as long as a given shell is expanding its dark matter particle positions and velocities are not corrected for the effects of secondary perturbations. So, the contributions are evaluated analytically but are not imparted to the shell until it reaches turnaround.

The magnitude of the extra velocity imparted to a typical dark matter particle at the time of turnaround (Eq. (48) of RG87),

\begin{displaymath}\vert\Delta \vec{v}_{\rm rms}(x,t_c/2)\vert=F_{\rm v}(x,t_c/2)\Delta g_0[d(x,t_c/2),x]t_0,
\end{displaymath} (22)

where the velocity growth factor $F_{\rm v}(x,t)$ is

\begin{displaymath}F_{\rm v}(x,t)= \int_0^\theta F_g(x,\theta_1) \frac{{\rm d}t}{{\rm d} \theta_1} {\rm d} \theta_1
\end{displaymath} (23)

and the tangential, $v_{\rm tan}$, and radial, $v_{\rm rad}$, components of velocity are

\begin{displaymath}(\Delta v_{\rm tan})^2 = {2\over 3}\vert\Delta \vec{v}_{\rm r...
...\rm rad})^2={1\over 3}\vert\Delta \vec{v}_{\rm rms}(x)\vert^2.
\end{displaymath} (24)

The mean radial velocity at t=tc/2 is zero, while the mean tangential velocity is given by

\begin{displaymath}j_\theta(x)=\left (\frac{2}{3} \right)^{1/2} r_{\rm m} \Delta v(x,t_{\rm c}/2).
\end{displaymath} (25)

It is important to note that $j_\theta$ is inversely correlated with the height of the peak in $\delta$, since both $\Delta v$ and $r_{\rm m}$ (the maximum expansion radius) decrease with an increase in the height $\nu$ of the peak. The random angular momentum is important to the particle (shell) orbits since the larger it is, the larger is the orbital ellipticity, and the shell then penetrates less to the center, resulting in a flattening of the inner density profile. In several of the previously quoted papers (even Williams et al. 2004), only the random angular momentum was taken into consideration, while we also take account of ordered angular momentum, dynamical friction, and adiabatic contraction.

The collapse starts from the innermost shell, the one adjacent to the core. When the first shell reaches its $r_{\rm m}$, it collapses and finds its $r_{\rm a}$ and $r_{\rm p}$ within the overall halo potential. It is assumed that the potential changes slowly compared to the dynamical timescales of the shells, so that every shell conserves its adiabatic invariants, the radial and tangential momenta


 
    $\displaystyle j_\theta(x)=\Delta v_{\rm tan}\; r_{\rm m}$ (26)
    $\displaystyle j_r(x)=\int_{r_{\rm p}}^{r_{\rm a}} v_{\rm rad}~ {\rm d}r$ (27)

throughout the collapse. This is an important assumption in the RG87 formalism, it is crucial to the computation of dynamics of shell crossing. Until the moment when a given shell reaches its maximum expansion radius $r_{\rm m}(x)$ at a time t=tc(x)/2 corresponding to $\theta=\pi$, the shell is assumed to be thin, its radial extent determined by the initial shell separation. At $r_{\rm m}$, the average dark matter particle in the shell acquires its additional random velocity, given by Eqs. (22) and (24). In either case, the apocenter and pericenter can be calculated from the particle's energy integral and $\Delta v_{\rm tan}$ and $\Delta v_{\rm rad}$. The energy integral,

\begin{displaymath}E= \psi(r_{\rm m})-\frac{1}{2} \left[\left(\frac{{\rm d}r}{{\rm d}t}\right)^2+\left(\frac{j_\theta}{r}\right)^2 \right]
\end{displaymath} (28)

is conserved if $\psi$ does not change too rapidly changing with time, an approximation which is good outside the immediate region of the core.

The radial velocity of a particle is

\begin{displaymath}v_{\rm rad}=\left[2(E-\psi(r))-(j_\theta/r)^2\right]^{1/2},
\end{displaymath} (29)

(RG87). We note that RG87 define their energy integral and potential as the negatives of the conventional definitions of these quantities. We use the conventional definitions, i.e., the potential is a negative quantity inside the halo, and the energy integral of a bound particle is negative.

The addition of a random radial velocity ensures that the apocenter distance is greater than $r_{\rm m}$, the radius of maximum expansion in the absence of velocity perturbations. A measure of the mean radial momentum of the particle is the radial action, (RG87)

\begin{displaymath}j_r= \int_{r_{\rm p}}^{r_{\rm a}} v_r {\rm d}r= \overline{v}_r {\rm d}r(r_{\rm a}-r_{\rm p}).
\end{displaymath} (30)

The radial distribution of mass within a shell between apocenter and pericenter is not uniform. The density in the radial range dr around r is proportional to the amount of time the particle spends there, (Eq. (53) of RG87)

\begin{displaymath}P(r)\;{\rm d}r={{v_{\rm rad}^{-1}\;{\rm d}r}\over{\int_{r_{\rm p}}^{r_{\rm a}}v_{\rm rad}^{-1}\;{\rm d}r}}\cdot
\end{displaymath} (31)

If $M_{\rm shell}$ is the mass of the added shell, then the total mass distribution of the core and shell together is:

\begin{displaymath}M_1(r)=M(r)+M_{\rm shell} \int_{r_{\rm p}}^r P(r') {\rm d} r'.
\end{displaymath} (32)

From this mass function, the new potential $\psi_1(r)$ is calculated. The next mass shell is added and its probability distribution is calculated in the potential $\psi_1(r)$. However, the mass distribution of a newly added shell overlaps with shells that have collapsed earlier; that is, the pericenter of the shell is at a smaller radius than the apocenters of some fraction of the previously added shells. Thus, after adding each new shell, we must recompute the orbits for each shell with which it overlaps. Since the potential does not change violently, the adiabatic invariants of the orbits are conserved. In this case, the adiabatic invariants are the angular momentum  $j_\theta$ and the radial action jr. By repeatedly adding shells in this manner, while adjusting the orbits so that the angular momentum and radial action of the orbits are conserved, a self-consistent mass distribution is built up[*].

The collapse of a perturbation taking into account angular momentum, and dynamical friction may be calculated by solving the equation for the radial acceleration (Kashlinsky 1986, 1987; Antonuccio-Delogu & Colafrancesco 1994; Peebles 1993)

\begin{displaymath}\frac{{\rm d}v_r}{{\rm d}t}=\frac{h^2(r,\nu )+j^2(r, \nu)}{r^3}-G(r) -\mu \frac{{\rm d}r}{{\rm d}t}
\end{displaymath} (33)

where $h(r,\nu )$ is the ordered specific angular momentum generated by tidal torques, $j(r, \nu)$ the random angular momentum (see RG87), G(r) the acceleration, and $\mu $ the coefficient of dynamical friction.

In the peculiar case of $\mu=0$, Eq. (33) can be integrated to obtain the square of the velocity

\begin{displaymath}v(r)^2=2 \left[\epsilon -G \int_0^r \frac{m_T(y)}{y^2} {\rm d} y +\int_0^r \frac{h^2}{y^3} {\rm d}y
\right],
\end{displaymath} (34)

where $\epsilon$ is the specific binding energy of the shell that can be obtained from the previous equation at turnaround when, d $r/{\rm d}t=0$.

If $\mu\neq 0$, the previous equation must be substituted with

\begin{displaymath}\frac{{\rm d} v^2}{{\rm d} t}+2 \mu v^2=2 \left[
\frac{h^2+j^2}{r^3} -G\frac{m_T}{r^2}
\right] v
\end{displaymath} (35)

which can be solved numerically for v. In the previous equation, we have two unknown quantities, the specific ordered angular momentum, h, and the coefficient of dynamical friction, $\mu $.

The ordered angular momentum is connected to the tidal interaction of the proto-structure with the neighboring ones. The explanation of the galaxy spin gain by tidal torques was pioneered by Hoyle (1949). Peebles (1969) performed the first detailed calculation of the acquisition of angular momentum in the early stages of protogalactic evolution. More recent analytic computations (White 1984; Hoffman 1986, R88; Eisenstein & Loeb 1995; Catelan & Theuns 1996) and numerical simulations (Barnes & Efstathiou 1987) re-investigated the role of tidal torques in generating the angular momenta of galaxies.

In the present paper, we take into account both types of angular momentum: random j, and ordered, h. As described in RG87, to calculate the ordered angular momentum, one has first to obtain the rms torque, $\tau (r)$, on a mass shell and then calculate the total specific angular momentum, $h(r,\nu )$, acquired during expansion by integrating the torque over time (Ryden 1988a, hereafter R88, Eq. (35))

$\displaystyle h(r,\nu )=\frac 13\left( \frac 34\right) ^{2/3}
\frac{\tau _{\rm ...
...\vartheta )\frac{\delta _{\rm o}}{\overline{\delta _{\rm o}}}}
{\rm d}\vartheta$     (36)

where $M_{\rm sh}= 4 \pi \rho_{\rm b}\left[1+\delta(x)\right] x^2 \delta x$ is the mass in a thin spherical shell of internal radius x, $\tau_{\rm o}$ is the tidal torque at time t0, while the mean overdensity inside the shell, $\overline{%
\delta }(r)$, is given by Eq. (11).

As remarked by Del Popolo et al. (2001), the angular momentum obtained from Eq. (36) is evaluated at the time of maximum expansion $t_{\rm m}$. With the BBKS power spectrum (filtered on a galactic scale), for a $\nu=3$ peak of mass $\simeq$ $2 \times 10^{11} ~M_{\odot}$, the model gives a value of $h=2.5 \times 10^{74}$  $\frac{\rm g ~cm^2}{\rm s}$.

It is interesting to compare this result with a different method such as that of Catelan & Theuns (1996), who calculated the angular momentum at maximum expansion time (see their Eqs. (31)-(32)) and compared it with previous theoretical and observational estimates. Assuming the same value of mass and $\nu$ used to obtain our previously quoted result and the same distribution of angular momentum as adopted by Catelan & Theuns (1996), we obtain a value for the angular momentum of $h=2.0 \times 10^{74}$  $\frac{\rm g ~cm^2}{\rm s}$ in agreement with our result.

The total angular momentum of a system is often expressed in terms of the dimensionless spin parameter

\begin{displaymath}\lambda=\frac{L \vert E\vert^{1/2}}{GM^{5/2}},
\end{displaymath} (37)

where L is the angular momentum, summed over shells, and E is the binding energy of the halo. Expressing the values of angular momentum calculated as previously reported in terms of the spin parameter, we get a value of $\lambda \simeq 0.04$ for a 3$\sigma$ peak of $\simeq$ $ 10^{12}~M_{\odot}$[*]. This is comparable with values found in the literature (Barnes & Efstathiou 1987; Vivtska et al. 2002). Using the parameters in Vivitska et al. (2002), the maximum of the distribution of $\lambda$ (well approximated by a log-normal distribution) is $\lambda=0.035$, while there is a 90% probability that $\lambda$ is in the range 0.02-0.1.

Since our paper deals also with LSB galaxies, we recall that LSB galaxies are more angular-momentum-dominated than normal galaxies of the same luminosity (McGaugh & de Blok 1998). These galaxies are characterized by high values of $\lambda$, since, as shown in Boissier et al. (2002), 35% of all galaxies (in number) with $0.06<\lambda<0.21$ are LSBs. According to several analytical papers (Nusser 2001; Hiotelis 2002; Le Delliou & Henriksen 2003; Ascasibar et al. 2004; Williams et al. 2004), one therefore expects that these objects should typically have shallower density cusps. Moreover the cores of these galaxies are also much less dense than the simulations indicate. In the case of LSB galaxies, we assumed a value of $\lambda=0.06$ as a conservative estimate. In agreement with what has been reported, larger values of $\lambda$ should correspond to additional flattening of the rotation curves in Fig. 4. A more complete description of this calculation can be found in Del Popolo (2006).

We note that after turnaround, not only the random angular momentum, j, will contribute to the non-radial motions in the protostructure but also the ordered angular momentum, h. In fact, as successive shells expand to their maximum radius, they acquire angular momentum and then collapse on orbits determined by their angular momentum. Actual peaks are not spherical, thus the infall of matter will not be purely radial. Random substructure in the region surrounding the peak will divert infalling matter onto non-radial orbits. The role of these random motions is of fundamental importance to structure formation.

Dynamical friction was taken into account by introducing the dynamical friction force to the equation of motions. In ordered to calculate $\mu $, we recall that in hierarchical universes, matter is concentrated in lumps, and the lumps into groups and so on, which act as gravitational field generators. One can calculate the stochastic force generated by these field generators and then following Kandrup (1980) the dynamical friction force per unit mass

                              $\displaystyle \mathcal{F}$ = $\displaystyle -\mu v=-\frac{4.44[Gm_{\rm a}n_{\rm ac}]^{1/2} \log \left( 1.12N^{2/3}\right)}N \frac v{a^{3/2}}$  
  = $\displaystyle -\beta_{\rm o} \frac v{a^{3/2}}$ (38)

where $N=\frac{4\pi }3R_{\rm sys}^3n_{\rm a}$, is the total number of field particles, $R_{\rm sys}$ the system radius, $m_{\rm a}$ and $n_{\rm a}$ are, respectively, the average mass and the number density of the field particles, $n_{\rm ac}=n_{\rm a}\times a^3$ is the comoving number of field particles, and a is the expansion parameter, connected to the proper radius of a shell by

r(ri,t)=ria(ri,t). (39)

The number and mass of the field generators is then calculated using the theory of Gaussian random fields (BBKS). The acceleration obtained from Eq. (38) is used in the equation of motion. The reader is again referred to Del Popolo (2006) for a more complete description.

The shape of the central density profile is influenced by baryonic collapse: baryons drag dark matter in, causing the so-called adiabatic contraction (AC) and a steepening of the dark matter density slope. By considering a spherically symmetric protostructure, before dissipation, that consists of a baryonic fraction F with $F \ll 1$ of dissipational baryons and a fraction 1-F of dissipationless dark matter particles constituting the halo, Blumenthal et al. (1986) described an approximate analytical model to calculate the effects of AC. One assumes that the dissipational baryons and the halo particles are well mixed initially (i.e., the ratio of their densities is F throughout the protostructure). This is because the original angular momentum of the dark matter halo comes from gravitational (tidal) interactions with its environment. Thus, the dark matter and gas experience the same torque in the process of halo assembly and should initially have (almost) the same specific angular momentum (Klypin et al. 2002). A usual assumption is, therefore, that initially baryons had the same density profile as the dark matter (see the previous discussion and Mo et al. 1998; Cardone & Sereno 2005; Treu & Koopmans 2002; Keeton 2001).

In real systems, baryons assume a distribution in the central regions of the halo determined by the competition between dissipation and star formation. The dissipative infall of baryons ceases either when most of the baryon gas is converted into stars, ending dissipation, or when a rotationally supported gas disk is formed. As the baryons dissipate energy, they fall toward the center of the halo. If the baryons conserve their angular momentum as they fall inward, then the baryons will dissipate energy until they form a rotationally supported disk, whose shape is determined by the angular momentum distribution h(M) of the baryons. Summarizing, the two methods for halting dissipative infall produce a spheroidal distribution (if star formation ends the infall) or in a disk of stars (if angular momentum ends the infall) (R88a). Following usual practice, the final baryon distribution is assumed to be a disk for spiral galaxies (Blumenthal et al. 1986; Flores et al. 1993; Mo et al. 1998; Klypin et al. 2002; Cardone & Sereno 2005). In our calculations, we shall assume the Klypin et al. (2002) model for the baryon distribution (see their Sect. 2.1), when dealing with mass scales typical of spiral galaxies. In the case of elliptical galaxies and clusters, a typical assumption (that we use in this paper) is that baryons collapse to a Hernquist configuration (Rix et al. 1997; Keeton 2001; Treu & Koopmans 2002). One objection that one could advance is that clearly if the final baryon distribution is assumed to be a disk, the inner gravity potential is clearly non-spherical. The usual assumption that is made at this step (see RG87; R88a,b; Mo et al. 1998; Klypin et al. 2001; Cardone et al. 2006) is that the final configuration (in this case the disk) is assembled slowly, so that we can assume that the halo responds adiabatically to the modification of the gravitational potential and remains spherical while contracting. The angular momentum of the dark matter particles is then conserved and a particle that is initially at radius r, ends up at a radius r', as described in the following.

The scale of the Hernquist profile is fixed by the competition between dissipation and star formation, when most of the baryon gas is converted into stars (see RG87; Gnedin et al. 2004).

As achieved by RG87, we take into account adiabatic compression improving its description as outlined by Gnedin et al. (2004, see the following). The Blumenthal et al. (1986), RG87, and R88ab models of adiabatic contraction can be described as follows. Since the baryon fraction F is much less than one, the adiabatic invariants of the dark matter orbits are conserved (Blumenthal et al. 1986; RG87). If the matter is on circular orbits, the invariant is r M(r); if the dark matter is on radial orbits in a self-similar distribution, the invariant is $r_{\rm m} M(r_{\rm m})$ (Blumenthal et al. 1986). Since the actual orbits are neither circular nor radial, the two adiabatic invariants that are preserved are the angular momentum $j_\theta$ and jr, of the orbits (RG87). In the absence of dissipation, the mass distribution M0(r) has a value of F, which is constant throughout the structure. A self-consistent final solution, in which the baryons form a disk (or Hernquist configuration), conserving h, and the dark matter is compressed, conserving $j_\theta$ and jr, is found by an iterative procedure. In the zeroth-order approximation, the baryons have their predissipation distribution $M_{\rm B0}=F M_0(r)$, and the dark matter has the distribution $M_{\rm D0}=(1-F) M_0(r)$. While the dark matter distribution is held constant, the baryons fall inward, preserving their angular momentum until they are on circular orbits. Baryons originally at radius r and with specific angular momentum h will end up, in the stiff halo, at a radius,

\begin{displaymath}r'=\frac{h^2}{G[M_{\rm B0}(r)+M_{\rm D0}(r')]}\cdot
\end{displaymath} (40)

The baryons now have a new, more compact distribution $M_{\rm B1}(r)$, and the central concentration of the baryons will draw the dark matter inward. The new potential of the baryons plus the dark matter is given by

\begin{displaymath}\Phi_1(r)=-G \int_r^\infty \frac{M_{\rm D0}(s)+M_{\rm B1}(s)}{s^2} {\rm d}s.
\end{displaymath} (41)

After fixing the value of $j_\theta$ for each mass shell of the halo, the value of the apocenter is adjusted until the orbits of the dark matter in the potential $\Phi_1(r)$ have the same value of jr, which they had in the predissipation distribution. The distribution of mass $M_{\rm D1}(r)$ of the dark matter in the potential $\Phi_1$ is built orbit by orbit, ensuring that $j_\theta$ and jr are conserved. In the second iteration one has:

\begin{displaymath}r''=\frac{h^2}{G[M_{\rm B1}(r)+M_{\rm D1}(r'')]},
\end{displaymath} (42)

and the new potential $\Phi_2$ is calculated and the dark matter orbits are adjusted to preserve $j_\theta$ and jr. The iteration continues until convergence. As shown by Gnedin et al. (2004) and Klypin et al. (2002), the previous model can be improved in two ways. It is worth stressing that there is some debate of the validity of the adiabatic compression formalism, with authors like Jesseit et al. (2002) finding substantial agreement between the final dark matter distribution in numerically simulated haloes and that predicted by the adiabatic compression approach. On the other hand, this result has been contradicted on the basis of a set of higher resolution numerical simulations by Gnedin et al. (2004). According to these authors, the standard adiabatic compression formalism systematically overpredicts the dark matter density profile in the inner 5% of the virial radius, and the adiabatic compression formalism overpredicts the dark matter density less than $\simeq$$ 10\%$ at $r/r_{\rm v} \simeq 0.1$, while the error quickly decreases for larger values of  $r/r_{\rm v}$. Gustafsson et al. (2006) confirmed that the Blumenthal et al. (1986) model overestimates the central dark matter density. They also showed that the modified model proposed by Gnedin et al. (2004), even if it is a considerable improvement, it is not perfect. It is also found that the contraction parameters in their model depend not only on the orbital structure of the dark-matter-only halos, but also on the stellar feedback prescription that is most relevant to the baryonic distribution. As a caveat, we recall that according to Romano-Diaz et al. (2008), the study of Gustafsson et al. (2006), even if focused on the AC, had insufficient resolution to achieve conclusive results. In the following, we shall take into account the Gnedin et al. (2004) simple modification, which describes numerical results more accurately than Blumenthal et al. (2006). They proposed a modified adiabatic contraction model based on conservation of the product of the current radius and the mass enclosed within the orbit-averaged radius,

\begin{displaymath}M(\bar{r})r= {\rm const.},
\end{displaymath} (43)

where the orbit-averaged radius is

\begin{displaymath}\bar{r} = {2 \over T_r} \int_{r_p}^{r_{\rm a}} r ~ {{\rm d}r \over v_r},
\end{displaymath} (44)

where $T_{\rm r}$ is the radial period, $r_{\rm a}$ is the apocenter radius, and $r_{\rm p}$ the pericenter radius. The previous, classical AC model assumes no angular momentum exchange between different components (e.g., baryons and dark matter). However, transfer of angular momentum from the baryons to the dark matter can be produced by dynamical friction, and the change of the AC model due to this effect can be taken into account by means of a simple model by Klypin et al. (2002).

 \begin{figure}
\subfigure[]{\includegraphics[width=8.5cm,clip]{11404fig1a.eps} }
\subfigure[]{\includegraphics[width=8.5cm,clip]{11404fig1b.eps} }\end{figure} Figure 1:

Panel a) Dark matter haloes generated with the model of Sect. 2. In this case we do not take into account baryon collapse. The upper dashed-line represents the results of the Aquarius Experiment for a halo of $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$, the dot-dashed line represents the NFW profile for c=10, calculated by means of Eqs. (46)-(48) expressing the scaling radius $r_{\rm s}$ in terms of the virial radius. The short-dashed line, dashed line and solid line represent, respectively, the density profile obtained by means of the model of the present paper for $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$, $M \simeq 10^{10}~ h^{-1} ~M_{\odot}$, and $M \simeq 10^{8} ~h^{-1}~ M_{\odot}$ without inclusion of baryonic infall. Panel  b) same of panel  a) but taking account of baryonic infall.

Open with DEXTER

The universal baryon fraction has been inferred by observations involving different physical processes (Turner 2002). The power spectrum of matter inhomogeneities in large-scale structure is sensitive to the quoted F; the Two Degree Field Galaxy Redshift Survey reported a value of $0.15 \pm 0.07$ (Percival et al. 2001). Measurements of the angular power spectrum of the CMBR also provided a very significant estimate. The combined analysis in Jaffe et al. (2001) of several data sets gives F= 0.186+0.010-0.008. By using WMAP collaboration data (Spergel et al. 2003), one obtains $F \simeq 0.16$ ( $\Omega_{\rm b}=0.047 \pm 0.006$; $\Omega_{\rm m}=0.29 \pm 0.07$). After baryons cool and form stars, the baryon-to-total mass at the virial radius does not need to equal the previous quoted universal baryon fraction and may deviate from it depending on feedback effects, hierarchical formation details, and heating by the extragalactic UVB flux. In the AC calculation, one has to take into account that for large galaxies, only about half of the baryons, in principle available within the virial radius, goes into the central galaxy. For dwarf galaxies, this amount may be much smaller because of feedback effects, hierarchical formation details, and heating by the extragalactic UVB flux. In the following, in the case of large galaxies, we use one half of the universal baryonic fraction, and for dwarfs, values obtainable from Hoeft et al. (2007) (their Fig. 1). For example, in the case of a dwarf galaxy of $10^{10}~M_{\odot}$ it is about 0.06.

3 Results and discussion

After fixing the initial conditions and describing how to calculate angular momentum, dynamical friction, and adiabatic contraction, we can use the model in Sect. 2 to obtain the density profile of haloes. In Fig. 1, we plot the profiles obtained with our model and those predicted by numerical simulations. In Fig. 1, the short-dashed line, the dashed line, and the solid line represents, respectively, the density profile for haloes of $10^{12}~ h^{-1}~ M_{\odot}$, $10^{10} ~h^{-1} ~M_{\odot}$, and $10^{8}~ h^{-1}~ M_{\odot}$ calculated by means of the model of this paper. In order to compare the results of the model with those of N-body simulations, we plotted the NFW profile (dot-dashed line) for haloes with masses equal to $10^{12}~ h^{-1}~ M_{\odot}$, and c=10, and the results for the same mass obtained in the Aquarius Project by Navarro et al. (2008) (upper dashed-lines). Navarro et al. (2008) showed that density profiles deviate slightly but systematically from the NFW model, and are approximated more accurately by a fitting formula where the logarithmic slope is a power-law of radius, the Einasto profile

\begin{displaymath}\ln(\rho(r)/\rho_{-2})=-2/\alpha [(r/r_{-2})^{\alpha}-1],
\end{displaymath} (45)

where r-2 and $\rho_{-2}$ are connected to the scaling radius and density of the NFW profile (see the following) by $r_{-2}=r_{\rm s}$ and $\rho_{-2}=\rho_{\rm s}/4$.

Density profiles become monotonically shallower inwards, down to the innermost resolved point, with no indication that they approach power-law behavior. The innermost slope they measure is slightly shallower than -1, a result supported by estimates of the maximum possible asymptotic inner slope. Shallower cusps, such as the r-0.75 behavior predicted by the model of Taylor & Navarro (2001), cannot yet be excluded. Similarly, Graham et al. (2006) showed that Einasto and Prugniel-Simien models describe simulated dark matter halos better than a NFW model. A comparison of the quoted models with data from real galaxies gives an extrapolated inner (0.01-1 kpc) logarithmic profile slope ranging from approximately -0.2 to approximately -1.5, with a typical value at 0.1 kpc around -0.7.

The NFW profile for the given mass was calculated by means of the relationships connecting the concentration parameter, c, and the virial mass, $M_{\rm v}$, to the shape of NFW profile. We used the following equation for c

\begin{displaymath}c \simeq 13.6 \left( \frac{M_{\rm v}}{10^{11}~ M_{\odot}} \right)^{-0.13},
\end{displaymath} (46)

(Gentile et al. 2007), and the usual one for the NFW profile,

\begin{displaymath}\rho(r)=\frac{\rho_{\rm s}}{r/r_{\rm s}(1+r/r_{\rm s})^2}=\frac{\rho_{\rm b} \delta_{\rm v}}{r/r_{\rm s}(1+r/r_{\rm s})^2}
\end{displaymath} (47)

where

\begin{displaymath}\delta_{\rm v}=\frac{\Delta_{\rm v}}{3}
\frac{c^3}{\log(1+c)-c/(1+c)},
\end{displaymath} (48)

and $\Delta_{\rm v}$ is the virial overdensity (see Bryan & Norman 1998). The scaling radius, $r_{\rm s}$, of the NFW profile is connected to the virial radius, concentration parameter and virial mass by $c=r_{\rm v}/r_{\rm s}$, where,

\begin{displaymath}r_{\rm s} \simeq 8.8 \left( \frac{M_{\rm v}}{10^{11}~ M_{\odot}} \right)^{0.46} ~ {\rm kpc }
\end{displaymath} (49)

(Gentile et al. 2007). The residuals from the most robust Einasto fits are extremely small, and show no sign of systematic deviations down to the innermost resolved point. NFW residuals are then less than 30% over the full radial range and below 20% within r-2.

Figure 1 shows that halos generated using our method are different in character from the profiles predicted by numerical simulations, such as those of NFW or the Aquarius project. This result is unsurprising since the types of evolution that numerical N-body and our halos undergo are rather different: the former are produced as a cumulative result of many minor and major mergers of smaller sub-halos, while the latter are the product of quiescent accretion of lumpy matter onto the primary halo. The NFW profiles change slope dramatically from $\alpha=1$ to $\alpha=3$ at the characteristic radius $r_{\rm s}$. For example, in the case of the $10^{12}~ h^{-1}~ M_{\odot}$, the NFW profile has a characteristic scale length, equal to $0.1 ~r_{\rm v}$, beyond which the density profile steepens, so that much of the mass accumulates within 10% of the virial radius. We note that Navarro et al. (2004), Stadel et al. (2008), and the Aquarius Experiment start to show that the spherically-averaged density profile is not described, as in the NFW model, by just an outer power-law $\rho \propto r^{-3}$ and inner power-law $\rho \propto r^{-1}$, but becomes progressively shallower inwards, and at the innermost resolved radius, the logarithmic slope is $\gamma \equiv -{\rm d} \ln \rho /{\rm d} \ln r \preceq 1$, ruling out claims of a steep $\rho \propto r^{-1.2}$ central cusp, and not excluding shallower cusps as given by the r-0.75 behavior predicted by the model of Taylor & Navarro (2001). In other words, simulations of higher resolution are merging to the conclusion that the inner slope is not so steep as it appeared in older simulations and so in some cases even marginally consistent the DM halo of NGC 4605 ( $\rho \propto r^{-0.65}$) (Simon et al. 2003). In the following, we discuss another reason why density profiles obtained in N-body simulations differ from those obtained in our model.

NFW profiles change slope rapidly from $\alpha=1$ to $\alpha=3$ at the characteristic radius $r_{\rm s}$. For example, in the case of the $10^{12}~ h^{-1}~ M_{\odot}$ model, the NFW profile has a characteristic scale length, equal to $0.1 ~r_{\rm v}$, beyond which the density profile steepens, so that much of the mass accumulates within 10% of the virial radius. The haloes obtained using the model in Sect. 2 differ in character from the profiles predicted by numerical simulations, such as those of NFW. Within the virial radius the log-log density slope changes gradually and the slopes of the inner part of haloes flatten with decreasing mass. Our results show a steepening of the density profile with increasing mass of slopes $\alpha \simeq 0$ for $M \simeq 10^{8}{-} 10^{10} ~h^{-1}~ M_{\odot}$, and $\alpha \simeq 0.8$ for $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$. The profiles with $M \simeq 10^{8}{-} 10^{10} ~h^{-1}~ M_{\odot}$ are fitted well by a Burkert profile, which has been considered to be a good fit to the dark matter rotation curves inferred from observations (e.g., Salucci & Burkert 2000). The functional form of this profile is characterized by

\begin{displaymath}\rho(r)= \frac{\rho_{\rm o}}{(1+r/r_{\rm o})[1+(r/r_{\rm o})^2]},
\end{displaymath} (50)

where $\rho_{\rm o} \simeq \rho_{\rm s}$ and $r_{\rm o} \simeq r_{\rm s}$ (EZ01). The dark matter within the core is given by $M_{\rm o}=1.6 \rho_{\rm o} r_{\rm o}^3$. Although the dark matter parameters $r_{\rm o}$, $\rho_{\rm o}$ and $M_{\rm o}$ are in principle independent, the observations reveal a clear connection (Burkert 1995) given by

\begin{displaymath}M_{\rm o}= 4.3 \times 10^7 \left (r_{\rm o}/{\rm kpc}\right)^{7/3} ~M_{\odot}.
\end{displaymath} (51)

In panel a of Fig. 1, we plot results for the model of this paper in the case baryonic infall is not taken into account, and in panel b when baryonic infall is taken into account. Comparing panels a and b, we observe that the presence of baryons leads to a steeper profile in agreement with previous studies (Blumenthal et al. 1986; Flores et al. 1993; Klypin et al. 2002; Gnedin et al. 2004; Gustafsson et al. 2006). It is important to note that the effect of baryonic infall is more evident at early redshifts, as shown in Fig. 4a. As shown in Fig. 4a, from z=3 to z=2 the profile becomes even steeper than the NFW cusp. Later, the effect of angular momentum and dynamical friction reduces the slope and erases the cusp. However, as shown in Fig. 1 the initial steepening of the profile remains visible in the profiles calculated by the model in the present paper ( $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$, $M \simeq 10^{10}~ h^{-1} ~M_{\odot}$, and $M \simeq 10^{8} ~h^{-1}~ M_{\odot}$), since the profile that takes into account baryonic infall (panel b) are steeper than those that do not take it into account (panel a).

The model of the present paper gives just one profile for a given halo mass. In reality, N-body simulations and analytical models that differ from those of Ryden-Williams and that of this paper (e.g., the extended Press-Schechter formalism) predict an ensemble of profiles for a given halo mass, depending on the different formation histories. For example, using the extended Press-Schechter (EPS) conditional probabilities for halo progenitor masses one can construct detailed histories of the mass assembly of dark matter haloes. One is interested in computing the mass accretion histories (hereafter MAHs) of dark matter haloes, defined as $\Psi(M_0,z)=M(z)/M_0$, where M(z) is defined as the mass of the main progenitor halo, and M0 is the halo mass at z = 0. One then defines the average mass accretion history of a halo of mass M0 as

\begin{displaymath}\Psi(M_0,z)= 1/N \sum_{i=1}^N \Psi_i (M_0,z),
\end{displaymath} (52)

where the summation is over an ensemble of N random realizations (see van den Bosch & Frank 2002). One can then convert the mass growth curves M(z) in density profiles as for example shown by Nusser & Sheth (1999). The single profile that we obtain in the present paper must be considered to be the average profile of all simulations. Two important things must be noted: a) less massive haloes are less concentrated; b) the halo's inner slope is shallower for smaller mass.

The first point can be explained as follows: higher peaks (larger $\nu$), which are progenitors of more massive haloes, have greater density contrast at their center, and so shells do not expand far before beginning to collapse. This reduces j and h and allows haloes to become more concentrated. An alternative explanation is connected to the quoted angular momentum-density anticorrelation demonstrated by Hoffman (1986): $j \propto \nu ^{-3/2}$ (and similarly for h). So, density peaks with low (high) values of $\nu$ acquire a larger (smaller) amount of angular momentum than high $\nu$ peaks and consequently the halo becomes less (more) concentrated. We note that the quoted trend of increased central concentration as a function of mass applies only to halos that started out as peaks in the density field smoothed by a fixed Rf scale. Our conclusions do not mean that, for example, clusters of galaxies should be far more centrally concentrated than galaxies, since different smoothing scales would apply in the two cases. Point (b) can be explained in a similar way to (a), as described previously. Less massive objects are generated by peaks of smaller $\nu$, which acquire more angular momentum (h and j). The angular momentum sets the shape of the density profile at the inner regions. For pure radial orbits, the core is dominated by particles from the outer shells. As the angular momentum increases, these particles remains closer to the maximum radius, resulting in a shallower density profile. Particles with smaller angular momentum will be able to enter the core but with a reduced radial velocity compared with the purely radial SIM. For some particles, the angular momentum is so large that they will never fall into the core (their rotational kinetic energy causes them to become unbound). Summarizing, particles of larger angular momenta are prevented from coming close to the halo's center and so contributing to the central density, which has the effect of flattening the density profile. This result agrees with the previrialization conjecture (Peebles & Groth 1976; Davis & Peebles 1977; Peebles 1990), according to which initial asphericities and tidal interactions between neighboring density fluctuations induce significant non-radial motions that oppose the collapse. In order to reproduce the NFW profile, we performed an experiment similar to that performed by Williams et al. (2004), namely we reduced the magnitude of the h and j angular momentum, dynamical friction, $\mu $, and we considered the system to consist only of dark matter (F=0). The experiment was performed on a halo of mass $10^{12}~ h^{-1}~ M_{\odot}$, and to reproduce NFW profile with c=10 and mass $\simeq$ $10^{12}~ h^{-1}~ M_{\odot}$, we had to reduce the magnitude of h by a factor of 2, and both j and $\mu $ by a factor of 2.5. The result of this experiment is represented by the dashed line in Fig. 2, which closely reproduces the NFW profile (dot-dashed line), and the dotted line is the profile for the $10^{12}~ h^{-1}~ M_{\odot}$ halo obtained by means of our model[*].

 \begin{figure}
\par\hspace*{1mm}
\psfig{file=11404fig2.eps,width=8.5cm}
\end{figure} Figure 2:

Dark matter haloes generated with the model of Sect. 2. The dotted-dashed line is the NFW profile, the short-dashed line is the profile for the $10^{12}~ h^{-1}~ M_{\odot}$ halo obtained by means of our model. The dashed line represents the density profile obtained from the $10^{12}~ h^{-1}~ M_{\odot}$ halo by reducing the magnitude of h, j and $\mu $ as described in the text.

Open with DEXTER

Williams et al. (2004) similarly had to reduce the random velocities, which is equivalent to reducing the angular momentum, to obtain a NFW profile. After each reduction in the random velocities, the profiles became steeper at the center. This effect can be interpreted, as the central density being accumulated by shells whose pericenters are very close to the center of the halo. The correlation between increasing angular momentum and the reduction in the inner slopes of haloes has also been noted by several other authors (Avila-Reese et al. 1998, 2001; Subramanian et al. 1999; Nusser 2001; Hiotelis 2002; Le Delliou & Henriksen 2003; Ascasibar et al. 2004). As previously mentioned, there could be, according to Williams et al. (2004), another reason for the difference in density profiles obtained by N-body simulations and by models such as those in this paper (or William's). As shown in their Fig. 1 (lower-left panel), the specific angular momentum distribution (SAM) in NFW-like halo is more centrally concentrated than the SAM distribution of the reference halo obtained with their model, model which is similar to that of this paper, and is closer to those of typical halos emerging from numerical simulations. This may suggest, in agreement with Williams et al. (2004), that haloes in N-body simulations lose a considerable amount of angular momentum between 0.1 and 1 rv. Since virialization proceeds from inside out, this means that the angular momentum loss takes place during the later stages of the haloes' evolution, rather than very early on. This appears to be confirmed by the so-called angular momentum catastrophe, namely that dark matter halos generated by gas-dynamical simulations are too small and have too little angular momentum compared to the halos of real disk galaxies, possibly because it was lost during repeated collisions through dynamical friction or other mechanisms (van den Bosch et al. 2002; Navarro & Steinmetz 2000). The problem can be solved by introducing stellar feedback processes (Weil et al. 1998), but part of the angular momentum problem seems to be caused by numerical effects, most likely related to the shock-capturing, artificial viscosity used in smoothed particle hydrodynamics (SPH) simulations (Sommer-Larsen & Dolgov 2001).

We discussed the effect of changing the magnitude of angular momentum but did not consider the effect of changing the magnitude of $\mu $ (dynamical friction), which is similar to changing the magnitude of angular momentum: an increase in the term $\mu $ produces shallower profiles in much the same way as larger values of angular momentum. This is expected from Del Popolo (2006) (Fig. 1), showing that dynamical friction influences the dynamics of collapse in a similar way to that of angular momentum, slowing down the collapse of outer shells and compelling the particles to remain closer to the maximum radius.

 \begin{figure}
\par\psfig{file=11404fig3.eps,width=9.0cm}
\end{figure} Figure 3:

Comparison of the rotation curves obtained with the model in Sect. 2 (solid lines) with the rotation curves of four LSB galaxies studied by de Blok & Bosma (2002). The dashed line represents the fit with NFW model (see Sect. 3 for details).

Open with DEXTER

At this point, we emphasize that the flattening of the density profile is caused by two main effects: angular momentum and dynamical friction. The other effect that we have considered, namely baryonic infall, counteracts the effect of the density profile flattening produced by angular momentum and dynamical friction. The role of baryonic infall is predominant at early times and produces a steepening of the profile, while at later times its effect is overwhelmed by the effect of angular momentum and dynamical friction. The result is similar to that described by Williams et al. (2004), apart from that here the effect of ordered angular momentum and dynamical friction adds to that of the random angular momentum studied by Williams et al. (2004), producing a more significant flattening of the density profiles. In Fig. 3, we plot the rotation curves obtained by our model and we compare them to four LSB galaxies taken from the high-resolution data of de Blok & Bosma (2002). Of the 26 LSBs that are presented in that paper we picked the ones with high inclination angles, smooth rotation curves, and good agreement between HI, Ha, and optical data. In all four cases, the data are compared with the rotation curve obtained using our model (solid line) and rotation curves obtained from the NFW profile (dotted lines), given by

\begin{displaymath}V(r)=V_{\rm v}
\left \{
\frac{\ln (1+cx)-cx/(1+cx)}
{x [\ln(1+c)-c/(1+c)]}
\right \}^{1/2},
\end{displaymath} (53)

where $x=r/r_{\rm v}$ and $V_{\rm v}$ is the virial velocity[*]. The concentration parameter c was chosen to be consistent with the NFW predictions, calculating it from the mass of the galaxy through Eq. (46). For all the four galaxies the rotation curves obtained with the model of the present paper (solid lines) are a considerably better fit than NFW (dotted lines). Figure 3 shows that the typical rotation velocities of NFW haloes are higher than those of the rotation curves obtained using our model, in which more massive haloes tend to be more centrally concentrated and have flatter rotation curves. Less massive haloes are less concentrated, and have more slowly rising rotation curves. In contrast, NFW rotation curves rise very steeply and, as a consequence, NFW fits to dwarf galaxy rotation curves have too low concentration parameters (van den Bosch & Swaters 2001) compared to N-body predictions. The NFW profiles fail to reproduce the velocities and the shape of the observed rotation curves, predicting too high velocities in the central part of haloes, and even leaving c as a free parameter, instead of using Eq. (46), there is no appreciable improvement in the fit. Using Eq. (46), one obtains very low values of c. The result is similar to that described by Gentile et al. (2004): data are far more accurately described by core-like profiles, such as the Burkert profile, which predicts flatter rotation curves than NFW profiles given by
$\displaystyle V(r)= \left
(\frac{2 \pi G \rho_0 r_{\rm o}^3}{r} \right)^{1/2}
\...
...rm o}) \sqrt(1+(r/r_{\rm o})^2)
\right]-
\arctan (r/r_{\rm o})
\right \}^{1/2},$     (54)

Our rotation curves are very similar to those generated by the Burkert profile and the residuals and discrepant points for our rotation curves are close to those given in Gentile et al. (2004) for the Burkert's fit to their data. Our haloes appear to be a closer match to the haloes of spiral and dwarf galaxies, than N-body haloes. This may indicate that the haloes of real late-type disk galaxies undergo a formation scenario similar to the one depicted by our method, i.e., collapse proceeds through a quiescent accretion of lumpy material and minor mergers, rather than through a merger-driven formation process characteristic of hierarchical models. Before going on, we notice that of the four rotation curves, NGC 100 seems to be almost equally well fitted by the two type of models. In some other cases, the difference between the NFW rotation curves and those of the present paper is even smaller than for NGC100 and so it would be worthwhile to discuss halo-to-halo scatter. During the hierarchical assembly of dark matter haloes, the inner regions of early virialized objects often survive accretion onto a larger system, thus giving rise to a population of subhaloes. Depending on their orbits and masses, these subhaloes therefore either merge, are disrupted, or survive to the present day. To fully describe, in a statistical sense, the non-linear distribution of mass in the Universe, it is essential that halo substructure is taken into account. These subhaloes appear ubiquitous in high resolution cosmological simulations and provide the source of fluctuations. One can study the evolution of halos under the influence of a generation of subhaloes, while real halos grow continuously by accretion and mergers. Fluctuations caused by subhaloes in parent halos are important for understanding the time evolution of dark matter density profiles and the halo-to-halo scatter of the inner cusp seen in recent ultrahigh resolution cosmological simulations. This scatter may be explained by subhalo accretion histories: when we allow for a population of subhaloes of varying concentration and mass, the total inner profile of dark matter can either steepen or flatten.

As noted in the introduction, while on galactic scales many studies predict central cores and it seems that the cusp/core problem is a real problem that cannot be attributed to systematic errors in the data (de Blok et al. 2003), on cluster scales the situation is less clear with slopes ranging from $\alpha=0.5$ (e.g., Sand et al. 2002, 2004) to $\alpha=1.9$ (Arabadjis et al. 2002). In order to study the problem on cluster scales, we calculated the density profile evolution of dark matter and that of the total matter distribution for a halo of $\simeq$ $10^{14}~ h^{-1}~ M_{\odot}$. We used the method of Sect. 2 to calculate the density profile and then we repeated the experiment shown in Sect. 3 to reproduce the NFW profile for the quoted mass that we considered as the initial condition[*]. We then evolved this profile using the method of Sect. 2 and taking the redshift dependence in the model into account using the technique described in Del Popolo (2001) (the reader is referred to the quoted paper for more details). The goal was to study the evolution of a NFW profile similar to El-Zant et al. (2001, 2004), Tonini et al. (2006), and Romano-Diaz et al. (2008), in a way.

Before showing the results, we recall that in real haloes angular momentum is lost in the final stages of collapse by the transfer of angular momentum from the subclumps to external material by dynamical friction (Quinn & Zurek 1988; R88b; Klypin et al. 2002). Once the baryons condense to form stars and galaxies, they experience a dynamical friction force from the less massive dark matter particles as they move through the halo. Energy and angular momentum is transferred to dark matter, increasing its random motion. Moreover, angular momentum acquired in the expansion phase gives rise to non-radial motions in the collapse phase. The spherical approximation that we use, which ignores the possibility of large substructures forming, neglects the transfer of angular momentum by means of the interactions of subclumps[*]. In the present paper, the effect of dynamical friction was taken into account as an additive term in the equation of motion, Eq. (33), without changing the spherical symmetry of the problem as done by Peebles (1993) and Del Popolo (2006). The effect of the dynamical friction term in the equation of motion is very similar to that of the term of angular momentum: an increase in the term $\mu $ produces shallower profiles in a similar way to larger values of angular momentum. This is expected from Fig. 1 of Del Popolo (2006), showing that dynamical friction influences the dynamics of collapse in a similar way to that of angular momentum slowing down the collapse of outer shells and so compelling the particles to remain closer to the maximum radius.

 \begin{figure}
{\hbox{(a)
\psfig{figure=11404fig4a.ps,width=9.0cm} }}
{\hbox{(b)
\psfig{figure=11404fig4b.ps,width=9.0cm} }}\end{figure} Figure 4:

Panel a) Density profile evolution of a $10^{14}~ h^{-1}~ M_{\odot}$ halo. The solid line represents the NFW initial profile at z=3. Notice how the baryonic infall at redshift 2<z<3 steepens the NFW cusp. The profile at z=2, z=1.5, z=1 and z=0 is represented by the dot-short-dashed line, dotted line, short-dashed line, and long-dashed line respectively. Panel  b). Density profile evolution of a $10^{14}~ h^{-1}~ M_{\odot}$ halo for the total mass. The solid line represents the NFW initial profile at z=3. The profile at z=0 is represented by the dotted line.

Open with DEXTER

Figure 4a plots the evolution in a density profile of a mass $10^{14}~ h^{-1}~ M_{\odot}$. The solid line represents the initial density profile at z=3, which subsequently steepens due to baryon settling in virialized dark matter haloes (AC) (dot-short-dashed line) at z=2. Moreover, angular momentum acquired in the expansion phase gives rise to non-radial motions in the collapse phase. The effect of angular momenta and dynamical friction overcomes that of the AC and the profile starts to flatten (z=1.5 dotted line; z=1 short-dashed line; z=0 long-dashed line). The final dark matter profile (long-dashed line) is characterized by a log-log slope of $\alpha \simeq 0.7$ at $0.01 r_{\rm s}$. The situation is therefore similar to that of haloes on galactic scales but the slope is steeper than for dwarf galaxies. In Fig. 4b, we plot the evolution in the total density profile at initial redshift, z=3, (solid line) and at the final redshift z=0 (dotted line). The plot shows that the total density profile slightly changes with time and the cusp is not ``erased'' as in the case of the dark matter profile. This result also implies that the baryonic component becomes steeper than the original NFW profile. The behavior of total mass agrees with that inferred from X-ray observations by Chandra and XMM (Buote 2003, 2004; Lewis et al. 2003), weak lensing (Dahle et al. 2003), and strong lensing (Bartelmann 2003), which are consistent with a cusp of value $\alpha=1$ or larger. The behavior of the dark matter halo agrees with that inferred by the analysis of Sand et al. (2002, 2004), who fitted the baryonic and dark matter profiles only in the very inner part of the cluster MS 2137-23 within $\simeq$ 50  h-1 kpc with a generalized NFW profile

\begin{displaymath}\rho(r)= \frac{\rho_{\rm b} \delta_{\rm v}}{(r/r_{\rm s})^{\alpha}(1+r/r_{\rm s})^{3-\alpha}}
\end{displaymath} (55)

obtaining a nearly flat core $\alpha=0.35$. The steepening of the baryonic component is consistent with the results of Brunzendorf & Meusinger (1999), who found that the projected galaxy distribution in the Perseus cluster diverges as r-1.

The results previously reported have several implications on the effort to test predictions of the CDM model observationally. The test that has, as previously mentioned a number of times, is the density distribution in the inner regions of galaxies and clusters.

On galactic scales, as our results show, the infall of baryons at early times steepen the cusp because of AC (in agreement with previous studies: Blumenthal et al. 1986; Gnedin et al. 2004; Gustafsson et al. 2006), but later the cusp is erased by dynamical friction and non-radial motions effects. This results agree with other analyses that study separately the effects of dynamical friction (e.g., EZ01) and those of angular momentum (non-radial motions) (e.g., Nusser 2001; Hiotelis 2002; Ascasibar et al. 2004; Williams et al. 2004), and with the recent SPH simulations of Romano-Diaz et al. (2008).

On larger scales the situation is different. The analysis of the density distribution for bright galaxies is complicated by the uncertain contribution of stars to the total mass profile (Treu & Koopmans 2002; Mamon & Lokas 2004). Some analyses tend to favor inner slopes shallower than predicted by CDM (e.g., Gentile et al. 2004), but others measure inner profile slopes that are at least marginally consistent with predictions (Treu & Koopmans 2002, 2004; Koopmans & Treu 2003; Jimenez et al. 2003). As previously reported, our results exhibit a steepening in the density profile with increasing mass for a density profile of haloes of ${\rm mass } > 10^{12} ~h^{-1}~ M_{\odot}$ with slopes >0.8. This agrees with recent N-body simulations having a logarithmic slope that decreases inward more gradually than the NFW profile (Hayashi et al. 2003; Navarro et al. 2004; Stadel et al. 2008). In the case of Stadel et al. (2008) the logarithmic slope is 0.8 at $0.05 \%$ of $r_{\rm v}$.

The density distribution of galaxy clusters can, in principle, provide a cleaner test of the models because the effects of the baryons and gas on the dark matter distribution are expected to be smaller and simpler. However, in this case observations also predict slopes ranging from $\alpha \leq 0.5$ ( $\alpha=0.35$, Sand et al. 2002, 2004) to values larger than one (Arabadjis et al. 2002) with in some cases different results being found for the same object. This is the case for the cluster MS2137-23 studied by Sand et al. (2004), who found a shallow density slope ($\leq $0.5), while Dalal & Keeton (2003), Bartelman & Meneghetti (2004) and Gavazzi (2003) contested Sand's results, which according to them is neglecting lens ellipticity, and found consistency between the inner slope and NFW profile. Our result on cluster scales differs for the dark matter and total mass distribution: the first tends to be less cuspy than the NFW profile in agreement with some observations (e.g., Sand et al. 2002, 2004), while the second is a bit more cuspy than the NFW profile (in agreement with Brunzendorf & Meusinger 1999) .

In the present paper, we did not take into account the baryon-DM interaction. As previously reported, once the baryons condense to form stars and galaxies, energy and angular momentum is transferred to dark matter, increasing its random motion.

In future work, we will develop a semianalytic model to follow the evolution of the baryonic component, and its interaction with DM, and we will compare these results with those of full hydrodynamical simulations. This should provide fresh insights into galaxy formation, and directly address possible small-scale challenges to the CDM theory. The role of the interaction between DM and baryons could bring further changes in the final density profile, which we shortly discuss. We note that discrepancies between observations and CDM only become apparent on length scales on which the baryons begin to play a role not only for the cusp/core problem but also for the missing dwarfs problem. We also note that if, indeed, most star-forming galaxies in the early universe lost their DM cusps because of stellar feedback, the missing dwarfs (satellites) problem could also be solved. Dwarf galaxies without a central cusp have a lower average core density than cuspy ones, and are hence much easier to disrupt tidally during the hierarchical assembly of larger galaxies (Mashchenko & Sills 2005). As a consequence, the removal of galactic cusps by stellar feedback in the early universe would result in fewer satellites today (Mashchenko et al. 2006). This again indicates that baryon physics could be one of the missing pieces of the puzzle, and will very likely make a major contribution toward a solution. If this is true, it would be unwise to ignore the conclusions to which data are leading us, namely that small scales tells us more about galaxy formation than it does about CDM[*]. In other words, the centers of galaxies are special places, the only places where we can study dark matter under peculiar conditions.

4 Conclusions

In this paper, we have studied the cusp/core problem by improving the ESIM model of Williams et al. (2004). We have taken into account, simultaneously, the effects of ordered and random angular momentum, dynamical friction and adiabatic contraction. The improved SIM of the present paper, taking account the previous effects predicts haloes characterized by a log-log density slope that changes gradually within the virial radius and slopes of the inner haloes that flatten with decreasing mass. As in previous papers, AC steepens the density profiles. The density profiles of structure having masses smaller than $ 10^{11}~ h^{-1} ~M_{\odot}$ are well fitted by Burkert profiles. The comparison of some of the rotational curves given by de Blok & Bosma (2002) with the rotational curves obtained by means of our model, provides a good agreement. In the case of galaxy clusters, the density profile evolution is similar to that observed on galactic scales with the difference that the final slope is observed to be steeper than in the dwarf galaxies case. However, the total mass profile is still cuspy and evolves slightly with time. The behavior of the dark matter halo is in agreement with the analysis of Sand et al. (2002, 2004), who fitted the baryonic and dark matter profiles only to the very inner part of the cluster MS 2137-23 within $\simeq$ 50  h-1 kpc with a generalized NFW profile. The behavior of the total mass is in agreement with X-ray observations by Chandra and XMM (Buote 2003, 2004; Lewis et al. 2003) weak lensing (Dahle et al. 2003), and strong lensing (Bartelmann 2003), which are consistent with a cusp of $\alpha=1$ or of higher value. As also reported by Williams et al. (2004), numerically generated halos must have lost a considerable amount of their angular momentum in the outer parts, roughly between 0.1 and $1 ~R_{\rm v}$, possibly by dynamical friction or other mechanisms. A future model that takes into account the interaction between DM and baryons would help us to better understand the problem.

Acknowledgements
We would like to thank Massimo Ricotti, Nicos Hiotelis, and Antonaldo Diaferio for their very helpful suggestions and comments. Antonio Del Popolo acknowledges the financial support from the German Research Fondation (DFG) under grant NO KR 1635/16-1.

References

Footnotes

... model[*]
A slightly overdense sphere, embedded in the Universe, is a useful non-linear model, as it behaves exactly as a closed sub-universe because of Birkhoff's theorem. The sphere is divided into spherical ``shells''. A spherical ``shell'' may be defined as the set of particles at a given radius that are all at the same phase in their orbits (see Le Delliou & Henriksen 2003).
...$\nu=\delta(0)/\sigma $[*]
$\sigma$ is the mass variance filtered on a scale Rf.
... Eq. (44)[*]
Note a typo in their paper: $f_1(\theta)$ in the numerator of Eq. (44) should be $f_2(\theta)$.
... up[*]
In other words, the mass particles, once position and velocities are assigned, are allowed to follow the appropriate orbit in the gravitational potential of the previously collapsed matter. As each mass shell is added, those previously added shells with which it overlaps have their orbits adjusted so that the angular momentum  $j_\theta$ and radial moments jr integrated from pericenter to apocenter, are conserved.
...$ 10^{12}~M_{\odot}$[*]
Notice that there is a very mild dependence of $\lambda$ on peak height, or mass.
... model[*]
Note that the density profile of the Aquarius Experiment can be reproduced in the same manner with a small change in the h and j parameters used for the NFW profile.
... velocity[*]
The value of the characteristic velocity $V_{\rm v}$ of the halo is defined in the same way as the virial radius $r_{\rm v}$.
... condition[*]
We could have started directly from a NFW profile for the given mass calculated by means of the NFW fitting formula.
... subclumps[*]
An approximate analytical approach to the problem of exchange of angular momentum between the dark matter particles and infalling baryons was developed by Klypin et al. (2002).
... CDM[*]
Other possibilities are that the observed dark matter cores are telling us that dark matter has pressure at small scales or something unexpected.

All Figures

  \begin{figure}
\subfigure[]{\includegraphics[width=8.5cm,clip]{11404fig1a.eps} }
\subfigure[]{\includegraphics[width=8.5cm,clip]{11404fig1b.eps} }\end{figure} Figure 1:

Panel a) Dark matter haloes generated with the model of Sect. 2. In this case we do not take into account baryon collapse. The upper dashed-line represents the results of the Aquarius Experiment for a halo of $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$, the dot-dashed line represents the NFW profile for c=10, calculated by means of Eqs. (46)-(48) expressing the scaling radius $r_{\rm s}$ in terms of the virial radius. The short-dashed line, dashed line and solid line represent, respectively, the density profile obtained by means of the model of the present paper for $M \simeq 10^{12}~ h^{-1} ~M_{\odot}$, $M \simeq 10^{10}~ h^{-1} ~M_{\odot}$, and $M \simeq 10^{8} ~h^{-1}~ M_{\odot}$ without inclusion of baryonic infall. Panel  b) same of panel  a) but taking account of baryonic infall.

Open with DEXTER
In the text

  \begin{figure}
\par\hspace*{1mm}
\psfig{file=11404fig2.eps,width=8.5cm}
\end{figure} Figure 2:

Dark matter haloes generated with the model of Sect. 2. The dotted-dashed line is the NFW profile, the short-dashed line is the profile for the $10^{12}~ h^{-1}~ M_{\odot}$ halo obtained by means of our model. The dashed line represents the density profile obtained from the $10^{12}~ h^{-1}~ M_{\odot}$ halo by reducing the magnitude of h, j and $\mu $ as described in the text.

Open with DEXTER
In the text

  \begin{figure}
\par\psfig{file=11404fig3.eps,width=9.0cm}
\end{figure} Figure 3:

Comparison of the rotation curves obtained with the model in Sect. 2 (solid lines) with the rotation curves of four LSB galaxies studied by de Blok & Bosma (2002). The dashed line represents the fit with NFW model (see Sect. 3 for details).

Open with DEXTER
In the text

  \begin{figure}
{\hbox{(a)
\psfig{figure=11404fig4a.ps,width=9.0cm} }}
{\hbox{(b)
\psfig{figure=11404fig4b.ps,width=9.0cm} }}\end{figure} Figure 4:

Panel a) Density profile evolution of a $10^{14}~ h^{-1}~ M_{\odot}$ halo. The solid line represents the NFW initial profile at z=3. Notice how the baryonic infall at redshift 2<z<3 steepens the NFW cusp. The profile at z=2, z=1.5, z=1 and z=0 is represented by the dot-short-dashed line, dotted line, short-dashed line, and long-dashed line respectively. Panel  b). Density profile evolution of a $10^{14}~ h^{-1}~ M_{\odot}$ halo for the total mass. The solid line represents the NFW initial profile at z=3. The profile at z=0 is represented by the dotted line.

Open with DEXTER
In the text


Copyright ESO 2009

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.