EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 518 (July-August 2010) A&A, 518 (2010) A38 Full HTML
Free Access
Issue
A&A
Volume 518, July-August 2010
Herschel: the first science highlights
Article Number A38
Number of page(s) 8
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/201014147
Published online 31 August 2010
A&A 518, A38 (2010)

Triaxial collapse and virialisation of dark-matter haloes

C. Angrick - M. Bartelmann

Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany

Received 27 January 2010 / Accepted 13 May 2010

Abstract
We reconsider the ellipsoidal-collapse model and extend it in two ways: we modify the treatment of the external gravitational shear field, introducing a hybrid model in between linear and non-linear evolution, and we introduce a virialisation criterion derived from the tensor virial theorem to replace the ad hoc criterion employed so far. We compute the collapse parameters $\delta _{\rm c}$ and $\Delta _{\rm v}$ and find that they increase with ellipticity e and decrease with prolaticity p. We marginalise them over the appropriate distribution of e and p and show the marginalised results as functions of halo mass and virialisation redshift. While the hybrid model for the external shear gives results very similar to those obtained from the non-linear model, ellipsoidal collapse changes the collapse parameters typically by $(20\ldots50)\%$, in a way increasing with decreasing halo mass and decreasing virialisation redshift. We qualitatively confirm the dependence on mass and virialisation redshift of a fitting formula for $\delta _{\rm c}$, but find noticeable quantitative differences in particular at low mass and high redshift. The derived mass function is in good agreement with mass functions recently proposed in the literature.

Key words: cosmology: theory - methods: analytical - dark matter - cosmological parameters - galaxies: clusters: general

1 Introduction

The spherical-collapse model (e.g. Mota & van de Bruck 2004; Schäfer & Koyama 2008; Engineer et al. 2000; Bartelmann et al. 2006; Wang & Steinhardt 1998) is a fundamental ingredient in the theory of cosmic structure formation. Following the collapse of a slightly overdense, homogeneous sphere, it allows the derivation of two essential parameters; the overdensity $\Delta _{\rm v}$ of a virialised halo compared to the mean or the critical cosmic density, and the critical linear density contrast $\delta _{\rm c}$. The former is important because it allows relating sizes to masses of virialised structures, and the latter because it establishes a link between linear structure formation and the population statistics of collapsed haloes. Despite fundamental doubts as to the validity of such a simplified model for accurate cosmological predictions, the parameters derived from the spherical-collapse model or variants thereof allow surprisingly far-reaching predictions such as the halo mass function or the correlation properties of haloes, which are confirmed by numerical simulations.

The statistics of a Gaussian random field implies that spherical collapse should not occur. In fact, the probability distribution for the eigenvalues $\lambda_i$ of the Zel'dovich deformation tensor (Doroshkevich 1970) shows that spherical collapse has a vanishing probability,

                $\displaystyle p(\lambda_1,\lambda_2,\lambda_3)$ = $\displaystyle \frac{15^3}{8\pi\sqrt{5}\sigma^6}\exp\left(-\frac{3\delta^2}{\sigma^2}+\frac{15I}{2\sigma^2}\right)$  
    $\displaystyle ~\times(\lambda_1-\lambda_2)(\lambda_2-\lambda_3)(\lambda_1-\lambda_3)$ (1)

because it requires $\lambda_1=\lambda_2=\lambda_3$. In this equation, $\delta\equiv\lambda_1+\lambda_2+\lambda_3$, $I\equiv\lambda_1\lambda_2+\lambda_2\lambda_3+\lambda_1\lambda_3$, and $\sigma^2$ denotes the variance of the matter power spectrum smoothed on a scale that corresponds to a halo mass M. Modifications of the original Press-Schechter mass function (Press & Schechter 1974; Bond et al. 1991) motivated by ellipsoidal collapse (Sheth & Tormen 2002) substantially improve the agreement of analytic predictions on the halo population with numerical simulations (e.g. the Millennium simulation by Springel et al. 2005).

Ellipsoidal collapse was analysed many times before (see Bond & Myers 1996; Eisenstein & Loeb 1995; Bartelmann et al. 1993, for examples). Several authors have worked with the model by Bond & Myers (1996), generalising it for different cosmologies and introducing the scale factor a as a time variable (Monaco 1998; Sheth et al. 2001; Ohta et al. 2004; Monaco 1997,1995; Sheth & Tormen 2002). We are reconsidering it here for two reasons.

First, we want to analyse how different assumptions on the treatment of the environment of a halo impact on the parameters $\Delta _{\rm v}$ and $\delta _{\rm c}$. Previous assumptions have been that the principal-axis system of the homogeneous ellipsoid is either identical with that of the external gravitational shear field, or that the two eigensystems do not coincide, introducing rotation and a deviation from the homogeneous mass profile. We introduce another assumption here, letting the eigensystems of the collapsing ellipsoid and its surrounding shear field follow each other until turn-around of the major principal axis and then decoupling both.

Second, we want to stop the collapse along any of the principal axes according to a physically motivated virialisation condition. Virialisation must be invoked to prevent the axes from collapsing to zero, and thus to be able to follow the entire collapse of an ellipsoid, i.e. the collapse of its three principal axes. Conventionally, the collapse of each axis is stopped when ai(a)=0.177a, where ai is the scale factor of the ith axis and a is the background scale factor (Sheth et al. 2001; Bond & Myers 1996). In this way, $\Delta=a^3/(a_1a_2a_3)=178$ for spherical collapse in an Einstein-de Sitter universe at the time when the third axis virialises. However, the value 0.177 has no fundamental physical motivation, and there is no guarantee for it not to be different for ellipsoidal rather than spherical collapse, or whencosmologies other than EdS are to be considered. We introduce a general virialisation condition based on the tensor virial theorem that avoids introducting such an uncalibrated factor. We find substantial changes on both $\Delta _{\rm v}$ and $\delta _{\rm c}$ from both modifications and point out several discrepancies of our results with earlier studies.

We introduce the ellipsoidal-collapse model including our modifications in Sect. 2, present our results in Sect. 3 and conclude in Sect. 4.

2 The model

in this section, we shall briefly review the ellipsoidal-collapse model of Bond & Myers (1996) for cosmologies with a cosmological constant introducing the scale factor a as time variable. Furthermore, we shall present a physically motivated virialisation condition to stop the collapse of each axis, and show how to find the proper initial ellipticity and prolaticity as a function of mass and virialisation redshift.

2.1 The evolution equations

Let $a_i=R_i/R_{\rm pk}$ be the dimension-less principal axes of the ellipsoid, where Ri with $1\leq i\leq3$ are its dimensional semi-major axes, and $R_{\rm pk}$ the size of a spherical top-hat corresponding to a mass $M=(4\pi/3) \rho_{\rm b} R_{\rm pk}^3$ with the cosmological background density $\rho_{\rm b}$. The evolution of the three principal axes ai with time t in a cosmology with a cosmological constant $ \Lambda=(8\pi G/c^2)\rho_\Lambda$ is given by

\begin{displaymath} \frac{{\rm d}^2 a_i}{{\rm d}t^2}=\frac{8}{3}\pi G \rho_\Lamb... ...\delta}{3}+\frac{b_i}{2}\delta+\lambda_{{{\rm ext},~i}}\right) \end{displaymath} (2)

(Bond & Myers 1996), where G is the gravitational constant and c the speed of light. The density contrast of the ellipsoid with respect to the background density is $\delta=(\rho-\rho_{\rm b})/\rho_{\rm b}=a^3/(a_1a_2a_3)-1$. The parameters bi and $\lambda_{{{\rm ext},~i}}$denote the internal and external contributions to the gravitational tidal shear which occur because of the deviation from sphericity. Generally, the total tidal field is described by the tidal field tensor $\tens{T}$ with the elements $T_{ij}=\partial^2\Phi_{\rm P}/(\partial x_i\partial x_j)=T_{{\rm int,}ij}+T_{{\rm ext,}ij}$, where $\Phi_{\rm P}$ denotes the peculiar gravitational potential, and $T_{{\rm int,}ij}$ and $T_{{\rm ext,}ij}$ are the internal and external contributions to the shear, respectively. After a transformation into the ellipsoid's eigensystem, which is the same as the eigensystem of $\tens{T}$ in this model, the internal shear can be evaluated as

\begin{displaymath} b_i(t)= a_1(t)a_2(t)a_3(t)\int_0^\infty\frac{{\rm d}\tau}{[a_i^2(t)+1]\prod_{k=1}^3[a_k^2(t)+1]^{1/2}}-\frac{2}{3}, \end{displaymath} (3)

while the external shear can be approximated by

\begin{displaymath} \lambda_{{{\rm ext},~i}}(t)\equiv\left\{\begin{array}{rl}\fr... ...)& {\rm (non{\mbox -}linear ~approx.)},\\ \end{array}\right. \end{displaymath} (4)

where D+ is the linear growth factor, and the $\lambda_i$ are the eigenvalues of the Zel'dovich deformation tensor. See the first appendix of Bond & Myers (1996) for details of the calculation. Here and in the following the index ``0'' refers to initial values.

In linear approximation, the environment into which the ellipsoid is embedded evolves completely independently of it, whereas in the non-linear approximation it is tightly coupled to the ellipsoid. In Sect. 3.1 we shall introduce the hybrid model as a third approximation for the external shear.

We can rewrite Eq. (2) by using the scale factor a as time variable using Friedmann's equation $\dot a^2=a^2H_0^2E^2(a)$, where H0 is Hubble's constant. The expansion function of the universe, E(a), introduces the dimension-less density parameters of matter and dark energy today, $\Omega_{\rm m}$ and $\Omega_\Lambda$, respectively. This gives

\begin{displaymath} \frac{{\rm d}^2a_i}{{\rm d}a^2}+\left[\frac{1}{a}+\frac{E'(a... ...5 E^2(a)}C_i(a)-\frac{\Omega_\Lambda}{a^2 E^2(a)}\right]a_i=0, \end{displaymath} (5)

where a prime denotes differentiation with respect to a, and $C_i\equiv(1+\delta)/3+b_i/2+\lambda_{{{\rm ext},~i}}$. Equation (5) defines a set of three coupled second-order differential equations, for which we need six independent initial conditions compatible with the Zel'dovich approximation for early times. These are provided by
                              ai(a0) = $\displaystyle a_0[1-\lambda_i(a_0)],$ (6)
$\displaystyle \left.\frac{{\rm d}a_i}{{\rm d}a}\right\vert _{a_0}$ = $\displaystyle 1-\lambda_i(a_0)-\left.\frac{{\rm d}\!\ln D_+}{{\rm d}\!\ln a}\right\vert _{a_0}\lambda_i(a_0)\approx 1-2\lambda_i(a_0),$ (7)

since $D_+(a)\approx a$ for $a\ll1$. Choosing $a_0=2\times10^{-5}$, this is comfortably fulfilled. In Appendix A, we compare Eqs. (5)-(7) with the results presented by Monaco (1997) and see that they differ.

In the following, we shall assume that the eigenvalues $\lambda_i$ of the Zel'dovich tensor are ordered as $\lambda_1\geq\lambda_2\geq\lambda_3$, which implies that the ellipsoid first collapses along the direction 1, $a_1\rightarrow0$ first. At that time $\delta\rightarrow\infty$, and the collapse of the remaining two axes can no longer be followed so that we have to add a virialisation condition for each axis preventing their collapse to zero.

2.2 The virialisation condition

Conventionally, the collapse of each axis is stopped when ai(a)=0.177a, where ai is the scale factor of the ith axis and a the background scale factor (e.g. Sheth et al. 2001; Bond & Myers 1996). In this way, $\Delta=a^3/(a_1a_2a_3)=178$for spherical collapse in an Einstein-de Sitter (EdS) universe at the time when the third axis is assumed to virialise. However, the value 0.177 has no fundamental motivation in the physics of the collapse, and there is no guarantee that it remain unchanged in the case of ellipsoidal instead of spherical collapse, or in other cosmologies than Einstein-de Sitter.

We shall thus follow a different approach and present a physically well-motivated virialisation condition in the following to stop the collapse of each axis individually. We start from the tensor virial theorem (see Binney & Tremaine 1987, p. 213, 280),

\begin{displaymath} \frac{1}{2}\frac{{\rm d}^2I_{ij}}{{\rm d}t^2}=2K_{ij}+\Pi_{ij}+W_{ij}+V_{ij}, \end{displaymath} (8)

where $\tens{I}$ is the moment of inertia tensor, $\tens{K}$ and $\tens{\Pi}$ are the contributions to the kinetic energy tensor coming from ordered and random motions, respectively, $\tens{W}$ is the potential energy tensor, and $\tens{V}$ is the external potential energy tensor. Their elements are generally defined as
                               Iij $\textstyle \equiv$ $\displaystyle \int_\mathcal{V}{\rm d}^3x~\rho x_i x_j, \qquad W_{ij}\equiv -\int_\mathcal{V}{\rm d}^3x~\rho x_i\frac{\partial\Phi}{\partial x_j},$  
Kij $\textstyle \equiv$ $\displaystyle \frac{1}{2}\int_\mathcal{V}{\rm d}^3x~\rho\overline{v}_i\overline{v}_j, \qquad \Pi_{ij}\equiv\int_\mathcal{V}{\rm d}^3x~\rho\sigma_{ij}^2,$ (9)

and

\begin{displaymath} V_{ij}\equiv-\frac{1}{2}\int_\mathcal{V}{\rm d}^3 x~\rho\lef... ...al x_j}+x_j\frac{\partial\Phi_{\rm ext}}{\partial x_i}\right), \end{displaymath} (10)

where $\rho$ is the density of the fluid, $\Phi$ and $\Phi_{\rm ext}$ are the gravitational potentials of the ellipsoid itself and its surroundings, respectively, $\mathcal{V}$ is the volume which is integrated over, $\sigma_{ij}^2\equiv\overline{v_i v_j}-\overline{v}_i\overline{v}_j$ are the velocity dispersions, and the bar indicates averaging over $\mathcal{V}$.

We now specialise to the case of a homogeneous ellipsoid. For a stable mass configuration, the left-hand side of Eq. (8) has to vanish for each component of the inertial tensor individually. Since the ellipsoid is assumed to be at rest and the ellipsoid's eigensystem is chosen as a reference frame, $\overline{v}_i=0$ so that Kij=0 and $\Pi_{ij}=\int_\mathcal{V}{\rm d}^3 x~\rho \langle v_i^2\rangle \delta_{ij}$, with the Kronecker symbol $\delta_{ij}$. Note again that in this framework the eigensystems of the overdense ellipsoid and the gravitational tidal field are identical. For a homogeneous ellipsoid $v_i(x_i)=(\dot{a}_i/a_i)x_i$, thus

\begin{displaymath} \Pi_{ij}=\frac{1}{5}\dot{a}_i^2 M \delta_{ij}, \end{displaymath} (11)

with the mass M of the ellipsoid. The sum Wij+Vij can be evaluated using $-\nabla(\Phi+\Phi_{\rm ext})=\ddot{\vec{x}}$ and Eq. (2) to be

\begin{displaymath} W_{ij}+V_{ij}=\!\int_\mathcal{V}\!{\rm d}^3 x~ \rho \frac{{\... ...G}{3}\rho_\Lambda\!-\!4\pi G\rho_{\rm b}C_i\right)\delta_{ij}. \end{displaymath} (12)

Requiring that Eq. (8) is fulfilled for each axis separately and introducing the scale factor a as time variable yields the virialisation conditions for the three axes ai,

\begin{displaymath} \left(\frac{a_i'}{a_i}\right)^2=\frac{1}{a^2 E^2(a)}\left(\frac{3\Omega_{\rm m}}{2a^3}C_i-\Omega_\Lambda\right). \end{displaymath} (13)

When this condition is fulfilled for an axis together with $\dot{a}_i<0$, its collapse is stopped and its size is frozen in.

We emphasise that the former equation is consistent with the virialisation condition for spherical collapse in the EdS universe, where $R_{\rm v}/R_{\rm ta}=0.5$. The subscripts ``v'' and ``ta'' denote the time of virialisation and turn-around, respectively. For the special case of EdS, Eq. (13) becomes

\begin{displaymath} R=\frac{G M}{\dot{R}^2}, \end{displaymath} (14)

with $a_1=a_2=a_3\equiv R$. Using the parametric solution

\begin{displaymath} R=\frac{R_0}{2\delta_0}(1-\cos\theta), \qquad t=\frac{3 t_0}{4\delta_0^{3/2}}(\theta-\sin\theta) \end{displaymath} (15)

(Engineer et al. 2000) together with the relations $R_0/\delta_0=R_{\rm ta}$ and R03=9 GMt02/2 indeed gives the expected result that Eq. (14) is satisfied when $R=R_{\rm ta}/2$.

For bound objects in the $\Lambda $CDM model, the virialisation condition is

\begin{displaymath} \langle E_{\rm kin}\rangle=-\frac{1}{2}\langle E_{\rm pot}\rangle+\langle E_\Lambda\rangle, \end{displaymath} (16)

where $\langle E_{\rm kin}\rangle$ is the average kinetic energy of the halo, $\langle E_{\rm pot}\rangle$ its average potential energy, and $\langle E_\Lambda\rangle$ the average effective potential energy contributed by the cosmological constant. Using $\langle E_{\rm pot}\rangle=-3GM^2/(5R)$ and $\langle E_\Lambda\rangle=-\Lambda MR^2/10$ for the homogeneous sphere as well as energy conservation between turn-around and virialisation, one arrives at a cubic equation for $R_{\rm v}/R_{\rm ta}$,

\begin{displaymath} 2\eta\left(\frac{R_{\rm v}}{R_{\rm ta}}\right)^3-(2+\eta)\frac{R_{\rm v}}{R_{\rm ta}}+1=0, \end{displaymath} (17)

with $\eta=\Lambda R_{\rm v}^3/(3GM)$ (Lahav et al. 1991). The relevant solution of this equation agrees precisely with the condition derived from the tensor virial theorem, Eq. (13), as Fig. 1 shows. The small deviation occurs because we choose the time of virialisation rather than collapse as a reference ( $z=z_{\rm v}$) when using the tensor virial theorem. If we use the collapse time ( $z=z_{\rm col}$) instead, both conditions yield identical results.

\begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig1.eps}\vspace*{-2mm} \end{figure} Figure 1:

Comparison of the virialisation conditions from the scalar and tensor virial theorems. The curves differ slightly because they refer to different times. While the tensor virial theorem was evaluated at virialisation, the scaler virial theorem was evaluated at collapse. If we choose the collapse of the dark-matter halo as reference for z, both yield exactly the same result.

Open with DEXTER

2.3 The initial ellipticity and prolaticity

Equations (6), (7) imply that one has to choose initial values for the Zel'dovich deformation tensor to define the initial deviation of the principal axes and their time derivatives from the background. We shall show how they are chosen appropriately so that they comply with the assumed Gaussian nature of the Universe's initial conditions and represent a statistical average of haloes with the same mass M and virialisation redshift $z_{\rm v}$ but different shapes.

\begin{figure} \par\mbox{\includegraphics[width=6cm,clip]{figures/14147fig2a}\vs... ...ncludegraphics[width=6cm,clip]{figures/14147fig2c} } \vspace*{-2mm} \end{figure} Figure 2:

Evolution of the principal axes and the external shear for a $10^{14}\ M_\odot /h$ dark-matter halo with $e=\langle e\rangle$ and p=0 in the reference $\Lambda $CDM model. Left panel: linear approximation. Central panel: non-linear approximation. Right panel: hybrid approximation.

Open with DEXTER

Starting from the probability distribution for the eigenvalues of the Zel'dovich deformation tensor, Eq. (1), the conditional probability distribution for the ellipticity $e\geq0$ and the prolaticity $-e\leq p\leq e$, defined as $e\equiv(\lambda_1-\lambda_3)/(2\delta)$ and $p\equiv(\lambda_1-2\lambda_2+\lambda_3)/(2\delta)$, respectively, was derived by Sheth et al. (2001) to be

\begin{displaymath} g(e,p\vert\delta)=\frac{1125}{\sqrt{10\pi}}e(e^2-p^2)\left(\... ...{5}{2}\left(\frac{\delta}{\sigma}\right)^2(3e^2+p^2)\right]\!. \end{displaymath} (18)

To compute a statistical average of any quantity $\xi(e,p)$ for a given halo that depends on e and p, one should marginalise over the distribution g(e,p),

\begin{displaymath} \langle \xi\rangle=\int_0^\infty{\rm d}e\int_{-e}^e{\rm d}p~ \xi(e,p) g(e,p). \end{displaymath} (19)

However, a Taylor expansion of $\langle\xi\rangle$ up to second order in e and p around their mean values $\langle e\rangle$ and $\langle p \rangle$ under the distribution g(e,p) gives

\begin{displaymath} \langle \xi\rangle=\xi\left(\langle e\rangle,\langle p\rangl... ...}\right\vert _{\langle e \rangle,\langle p\rangle} \sigma_p^2, \end{displaymath} (20)

where $\sigma_e^2$ and $\sigma_p^2$ are the variances for e and p according to the distribution g. Up to first order, $\langle\xi\rangle=\xi\left(\langle e\rangle,\langle p \rangle\right)$. Deviations occur only at second order. The expectation values for e and p as well as their variances are given by
                           $\displaystyle \langle e\rangle$ = $\displaystyle \frac{3\sigma}{\sqrt{10\pi}\delta}, \hspace*{1.5cm} \langle p \rangle=0,$ (21)
$\displaystyle \sigma_e^2$ = $\displaystyle \frac{(19\pi-54)\sigma^2}{60\pi\delta^2}, \qquad \sigma_p^2=\frac{\sigma^2}{20\delta^2}\cdot$ (22)

Since the variances are $\propto(\sigma/\delta)^2$, and $\sigma/\delta<1$, using the approximation $\langle\xi\rangle\approx\xi\left(\langle e\rangle,\langle p\rangle\right)$ introduces only a small error, which is $\sim$1% for $\delta _{\rm c}$ and $\sim$3% for $\Delta _{\rm v}$ for $\Lambda $CDM and OCDM. For EdS, the error is larger and $\sim$3% for $\delta _{\rm c}$ and $\sim$10% for $\Delta _{\rm v}$. But since the latter cosmology is scientifically only of low relevance and usually serves as a reference model only, the usage of the former approximation is well justified. Instead of sampling $\xi(e,p)$ at several points for e and p, one only has to evaluate it once for $\langle e\rangle$ and $\langle p\rangle=0$.

Generally, the eigenvalues $\lambda_i$ are related to e and p by (see e.g. Bardeen et al. 1986; Bond & Myers 1996)

                       $\displaystyle \lambda_1$ = $\displaystyle \frac{\delta}{3}(1+3e+p)=\frac{\delta}{3}+\frac{\sigma}{\sqrt{10\pi}},$ (23)
$\displaystyle \lambda_2$ = $\displaystyle \frac{\delta}{3}(1-2p)=\frac{\delta}{3},$ (24)
$\displaystyle \lambda_3$ = $\displaystyle \frac{\delta}{3}(1-3e+p)=\frac{\delta}{3}-\frac{\sigma}{\sqrt{10\pi}},$ (25)

where we have set $e=\langle e\rangle$ and p=0 in the last step.

3 Results

In this Section, we show the results of the ellipsoidal-collapse model for the parameters $\delta _{\rm c}$ and $\Delta _{\rm v}$ for three different cosmologies and discuss how they are affected by the choice of the external-shear model.

3.1 Influence of the external shear

Figure 2 shows the evolution of the three principal axes and the eigenvalues of the external shear for a reference $\Lambda $CDM model with $\Omega_{\rm m}=0.3$, $\Omega_\Lambda=0.7$, and $\sigma_8=0.8$ for three different models of the external shear and $e=\langle e\rangle$, p=0. The shear in the linear approximation evolves completely smoothly over the entire collapse time. The evolution of each eigenvalue is given by D+(a) so that according to Eq. (4), $\lambda_{\rm ext,1}>0$, $\lambda_{\rm ext,2}=0$, and $\lambda_{\rm ext,3}=-\lambda_{\rm ext,1}<0$ at all times.

This is different in the non-linear approximation: At early times, the evolution of the $\lambda_{{{\rm ext},~i}}$ is the same, but soon thereafter they start evolving non-linearly and $\lambda_{\rm ext,2}$ becomes slightly negative. Noticeably there are steps in the evolution of the external shear whenever an axis virialises because $\lambda_{{\rm ext},~i}\propto a_1a_2a_3$, and the evolution of this volume factor changes after virialisation of each axis. In the right panel, we introduce the hybrid approximation: Initially, the evolution of $\lambda_{{{\rm ext},~i}}$ is described by the non-linear model. When one of the axes turns around, however, the corresponding eigenvalue of the external shear continues evolving linearly, i.e. its value at turn-around is then scaled by $D_+(a)/D_+(a_{\rm ta})$.

We believe that the hybrid model is the preferred model for the evolution of the external shear since it takes into account that the evolution of the ellipsoid itself and its vicinity should be tightly coupled in the beginning. At turn-around, however, they are definitely decoupled so that choosing this moment to switch from non-linear to linear evolution seems appropriate. Hence, we will use the hybrid model for the evolution of the external shear in the following.

\begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig3a}\hspace... ...[width=4.3cm,clip]{figures/14147fig3d}\vspace*{-1mm} \vspace*{-4mm} \end{figure} Figure 3:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) for three different models of the external shear. Both quantities were calculated for a $10^{14}\ M_\odot /h$ dark-matter halo in the reference $\Lambda $CDM model. Top panels: $e=\langle e\rangle$ as a function of virialisation redshift $z_{\rm v}$. Bottom panels: e=0.2 for all $z_{\rm v}$. In both cases, p=0.

Open with DEXTER

The influence of the external-shear model on the parameters $\delta _{\rm c}$ and $\Delta _{\rm v}$ is shown in Fig. 3. For a given mass and a given virialisation redshift $z_{\rm v}$, the initial overdensity is chosen such that the third axis of the ellipsoid finally virialises at $z_{\rm v}$. Both parameters can then be calculated by

\begin{displaymath}% \delta_{\rm c}=\frac{D_+(a_v)}{D_+(a_0)}\sum_{i=1}^3\lambda... ...{\rm m}(\vec{a}_v) \vec{a}_v^3}{a_1(a_v)a_2(a_v)a_3(a_v)}\cdot \end{displaymath} (26)

In this case, $\Delta _{\rm v}$ is defined with respect to the critical density as it can be primarily found in the literature. For the definition with respect to the background density one simply has to omit $\Omega_{\rm m}(a_{\rm v})$.

Using $\langle e\rangle$ as a function of virialisation redshift and p=0 (top panels of Fig. 3), the dependence of $\delta _{\rm c}$ and $\Delta _{\rm v}$ on $z_{\rm v}$is almost the same for the non-linear and the hybrid models. However, both differ from the linear model, showing that the external shear is most important at the beginning of the ellipsoid's evolution. At that time, the non-linear and the hybrid models agree. While $\delta_{\rm c}(z_{\rm v})$ is always smaller in the linear compared to the other two models, the curves for $\Delta_{\rm v}(z_{\rm v})$ cross. This reflects the circumstance that the initial overdensity in the linear model is different from that in the two other models, leading to a different initial ellipticity and therefore to a completely different evolution history. This can be seen in the bottom panels of Fig. 3, for which we have chosen e=0.2 independently of $z_{\rm v}$. They also clearly show that a varying initial ellipticity drives primarily the evolution of $\delta _{\rm c}$, whereas $\Delta _{\rm v}$ also strongly varies for fixed e. In this case, both $\delta _{\rm c}$ and $\Delta _{\rm v}$ are smaller in the linear-shear model compared to the non-linear and the hybrid models.

3.2 Parameters as function of mass and redshift

Before we present general results for the parameters $\delta _{\rm c}$ and $\Delta _{\rm v}$, we should comment on a subtle but very important issue: Whenever we want to compare our results with the ordinary spherical-collapse model, we have to keep in mind that we calculate all quantities at the time when the third axis virialises. Thus, we also have to compare these quantities those from the spherical-collapse model that are also calculated at the time of virialisation and not of collapse, i.e. when $R=R_{\rm ta}/2$ and not R=0 for EdS. This leads to slightly lower reference values of $\delta _{\rm c}$ and $\Delta _{\rm v}$ since $z_{\rm col}<z_{\rm v}$. Here the subscript ``col'' denotes collapse. Using the parametric solutions of Ohta et al. (2004) for the linear and the non-linear overdensity,

\begin{displaymath} \delta_{\rm l}=\frac{3}{5}\left[\frac{3}{4}(\theta-\sin\thet... ...nl}=\frac{9}{2}\frac{(\theta-\sin\theta)^2}{(1-\cos\theta)^3}, \end{displaymath} (27)

respectively, and $\theta=3\pi/2$ at virialisation, we find that $\delta _{\rm c}=1.583$ and $\Delta _{\rm v}=147$ for the EdS universe at $R=R_{\rm v}$ independent of $z_{\rm v}$. Recently, Lee & Ng (2009) arrived at the same values when accounting for the time of virialisation instead of collapse.

The top panels of Fig. 4 show $\delta _{\rm c}$ and $\Delta _{\rm v}$ for three different cosmologies for e=p=0, i.e. spherical systems. The OCDM cosmology is the same as our reference $\Lambda $CDM model except that $\Omega_\Lambda=0$. For the EdS model we set $\Omega_{\rm m}=1$ and $\Omega_\Lambda=0$. Indeed, for the EdS universe the constant values derived analytically are also reproduced by solving Eq. (5) numerically. This demonstrates again that Eqs. (5), (13) are fully consistent with the well-known spherical-collapse model. Note that the qualitative behaviour of $\Delta_{\rm v}(z_{\rm v})$ is the same as $\Delta_{\rm v}(z_{\rm col})$ (compare e.g. with Bartelmann et al. 2006). However, there is a difference for the critical linear overdensity whose shape as a function of $z_{\rm v}$ differs substantially from the shape as a function of $z_{\rm col}$. This should illustrate that the time chosen in the model when virialisation actually occurs ($z_{\rm v}$ or $z_{\rm col}$) can already have substantial impact on the qualitative behaviour of relevant quantities as a function of redshift, hence it is not necessarily a consequence of ellipsoidal collapse alone.

The bottom panels of Fig. 4 show $\delta _{\rm c}$ and $\Delta _{\rm v}$ for a triaxial halo with e=0.2 and p=-0.1. Also in this case both parameters are independent of $z_{\rm v}$ for the EdS universe, although $\delta _{\rm c}$ changes from 1.583 to 2.058, while $\Delta _{\rm v}$stays almost the same, 148 instead of 147. Interestingly, the most drastic changes in the shapes of both parameters occur for the OCDM model, for which the total density is only approximately a third of the critical density, while for the $\Lambda $CDM model their changes are small.

\begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig4a}\hspace... ...pace*{1.5mm} \includegraphics[width=4.3cm,clip]{figures/14147fig4d} \end{figure} Figure 4:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) for three different cosmological models and fixed initial e and p. Top panels: spherical case (e=p=0). Bottom panels: triaxial case with e=0.2 and p=-0.1.

Open with DEXTER

\begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig5a}\hspace... ...ace*{1.5mm} \includegraphics[width=4.3cm,clip]{figures/14147fig5d} \end{figure} Figure 5:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) as a function of halo mass and virialisation redshift for three different cosmologies marginalised over e and p. Top panels: dependence on virialisation redshift $z_{\rm v}$ for a halo with $M=10^{14}~M_\odot/h$. Bottom panels: dependence on mass M for a halo with $z_{\rm v}=0$. Thin dashed lines show the reference values from spherical collapse in the EdS universe, $\delta _{\rm c}=1.583$ and $\Delta _{\rm v}=147$.

Open with DEXTER

The influence of the halo mass M and the virialisation redshift $z_{\rm v}$ on both $\delta _{\rm c}$ and $\Delta _{\rm v}$ are illustrated in Fig. 5. For both decreasing mass and decreasing redshift, $\delta _{\rm c}$ is a monotonically decreasing function, approaching the reference values from the spherical-collapse model in the EdS universe for large masses and high virialisation redshifts, the situation for $\Delta _{\rm v}$ is much more complicated. For all three models, it has a minimum at redshifts 4-5 and at a mass of $\sim$ $10^{16}\ M_\odot/h$. For smaller values of mass and redshift, it is a monotonically decreasing function of both M and $z_{\rm v}$ for the OCDM and EdS models. It is also monotonically decreasing as a function of M in the $\Lambda $CDM model, but reaches a maximum at $z_{\rm v}\sim 1$. This is a direct result of the definition of $\Delta _{\rm v}$ with respect to the critical density. If it was defined with respect to the background density, the factor $\Omega_{\rm m}(a_{\rm v})$ would not appear in Eq. (26), and all three curves would increase with decreasing $z_{\rm v}$.

In Fig. 5 one can clearly see that the intervals that are covered for both $\delta _{\rm c}$ and $\Delta _{\rm v}$ are the largest for the EdS model, indicating that there is a stronger dependence of the ellipsoid's evolution on the total amount of matter in the Universe compared to the size of the cosmological constant if varying initial ellipticities are taken into account.

For either $M\rightarrow\infty$ or $z_{\rm v}\rightarrow\infty$, both $\delta _{\rm c}$ and $\Delta _{\rm v}$ must reach the reference values for spherical collapse in the EdS universe since the initial ellipticity $\langle e\rangle\propto\sigma/\delta$ and $\sigma$ both decrease with increasing mass, and $\delta$ has to be higher the earlier the structure is required to collapse. This expected behaviour can be clearly seen for $\delta _{\rm c}$, whereas for $\Delta _{\rm v}$ this happens finally for very large M and $z_{\rm v}$. We should stress again in this context that a crucial portion of the dependence on mass and virialisation redshift is driven by the change in the initial ellipticity, comparing Figs. 4 and 5.

\begin{figure} \par\includegraphics[width=4.4cm,clip]{figures/14147fig6a}\hspace*{1.5mm} \includegraphics[width=4.2cm,clip]{figures/14147fig6b} \end{figure} Figure 6:

Comparison of the fitting formula of Sheth et al. (2001) for $\delta _{\rm c}(M,z)$ with the results of this work for the reference $\Lambda $CDM model. Left panel: $\delta _{\rm c}$ as a function of mass and z=0. Right panel: $\delta _{\rm c}$ as a function of redshift and $M=10^{14}~M_\odot/h$.

Open with DEXTER

We compare in Fig. 6 the results of this work for $\delta _{\rm c}$ as a function of mass and redshift with the fitting formula of Sheth et al. (2001) given by

\begin{displaymath} \delta_{\rm c}(M,z)=\delta_{\rm c,sph}(z)\left\{1+0.47\left[... ...sigma^2(M,z)}{\delta_{\rm c,sph}^2(z)}\right]^{0.615}\right\}, \end{displaymath} (28)

where $\delta_{\rm c,sph}(z)$ is the redshift-dependent linear overdensity of the usual spherical-collapse model. The left panel shows that the dependence on mass is similar for both the fitting formula and the result of this work. However, the fit by Sheth et al. (2001) is less steep as a function of mass resulting from the differences of their underlying work to ours: first, they use the linear-shear model instead of the hybrid model, leading to a larger $\delta _{\rm c}$ for all redshifts as shown in Fig. 3. Second, the collapse of each axis is stopped using the artificial condition ai=0.177a, leading to $\delta_{\rm c}=1.686$ as a reference value at high mass. Thus, the fitting formula provides a larger value for large masses, while our approach leads to a $\delta _{\rm c}=1.583$ for $M\rightarrow\infty$ due to the virialisation condition that we apply. Third, for a given initial overdensity $\delta$, Sheth et al. (2001) use the most probable value $e_{\rm mp}=(\sigma/\delta)/\!$$\sqrt{5}$ instead of the expectation value $\langle e\rangle$ given by Eq. (21), which leads to initial ellipticities that are slightly too low so that the asymptotic limit for high redshifts is reached earlier than in our case.

The right panel of Fig. 6 shows a similar behaviour for $\delta _{\rm c}$as a function of redshift. The differences between the fitting formula and the result of our work again occur due to the different virialisation time and condition, and the difference between the most probable and the expectation value of e.

For applications, fitting formulae for both $\delta _{\rm c}$ and $\Delta _{\rm v}$ may be useful. We provide expressions here which are inspired by Eq. (28) of Sheth et al. (2001). For $\delta _{\rm c}(M,z)$ we suggest

\begin{displaymath} \delta_{\rm c}(M,z)=\delta_{\rm c,sph}(z_{\rm v})\left\{1+b\... ...z_{\rm v})}{\delta_{\rm c,sph}^2(z_{\rm v})}\right]^c\right\}, \end{displaymath} (29)

where b=0.6536, c=0.6387, and both $\delta_{\rm c,sph}$ and $\sigma^2$ are cosmology-dependent quantities. Note again that the redshift of virialisation, $z_{\rm v}$, has to be chosen as reference for z for both $\delta_{\rm c,sph}$ and $\sigma^2$. For a spatially-flat $\Lambda $CDM model with matter density parameter in the range $\Omega_{\rm m}\in[0.2,0.4]$, $M\in[10^{11},10^{15}]\ M_\odot\ h^{-1}$, and $z\in[0,10]$, the maximal error is $\sim$1.8% with a mean error of $\sim$0.4%.

A similar functional dependence can be found for $\Delta_{\rm v}(M,z)$. Only a small correction term has to be added to arrive at a satisfactory accuracy. We find

$\displaystyle \Delta_{\rm v}(M,z)=\Delta_{\rm v,sph}(z_{\rm v})\left[a\!+\!b~\sigma^{2c}(M,z_{\rm v})\!+\!d(1+z_{\rm v})^{2/5}\log^{9/4}(M)\right]    

with a=0.3819, b=0.5379, c=0.7589, and $d=3.456\times10^{-4}$. In the same range as above, the maximal error is $\sim$5% with a mean error of $\sim$1%.

\begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig7} \end{figure} Figure 7:

Influence of the initial ellipticity e and prolaticity p on the parameters $\delta _{\rm c}$ (top panels) and $\Delta _{\rm v}$ (bottom panels) for three different cosmological models. Left panels: $\Lambda $CDM. Central panels: OCDM. Right panels: EdS. The large circle in each figure indicates the value of the respective parameter for $e=\langle e\rangle$ and p=0. The small circles indicate values for a combination of $e\in[\langle e\rangle \pm \sigma_e/2,\ \langle e\rangle \pm \sigma_e]$ and $p\in [\pm e/4,\ \pm e/2]$ (see Eqs. (21), (22)). For a fixed e both $\delta _{\rm c}$ and $\Delta _{\rm v}$ grow for decreasing p.

Open with DEXTER

3.3 Influence of initial ellipticity and prolaticity

In Fig. 7 we plot both $\delta _{\rm c}$ and $\Delta _{\rm v}$ as a function of the initial ellipticity e and prolaticity p centered around their expectation values given by Eq. (21) for three different cosmologies. For increasing e and decreasing p, both parameters grow qualitatively in the same way as already reported by Sheth et al. (2001, cf. their Fig. 1). Quantitative deviations arise from the differences in the applied algorithm as discussed in Sect. 3.2. For a given mass and virialisation redshift, the initial overdensity for the EdS universe is larger compared to both $\Lambda $CDM and OCDM due to a shorter physical time interval that corresponds to the same redshift interval, resulting in a larger $\langle e\rangle$ and $\sigma_e$, but also in larger curvatures of $\delta _{\rm c}$ and $\Delta _{\rm v}$ with respect to e and p. These are the sources of the larger error in the approximation $\langle\xi\rangle\approx\xi(\langle e \rangle,\langle p\rangle)$ discussed in Sect. 2.3. Since the redshift-time relation is not very different between $\Lambda $CDM and OCDM, the dependences of $\delta _{\rm c}$ and $\Delta _{\rm v}$ on e and p are comparable.

3.4 Mass function

Using Eq. (29), we are able to construct the mass function of dark-matter haloes using the extended Press-Schechter formalism developed by Bond et al. (1991) and Lacey & Cole (1993), which is based on the first-upcrossing distribution of the density contrast $\delta$ as a function of the ``time variable'' $S\equiv\sigma^2(M)$. We shall proceed similarly as Sheth & Tormen (1999,2002) and define the scaled variable $\nu\equiv\delta_{\rm c,sph}^2/S$ to derive the mass function for our standard $\Lambda $CDM cosmology.

As Sheth & Tormen (2002) pointed out, expressing the first-upcrossing distribution f as a function of $\nu$ has the advantage that it is only necessary to calculate $f(\nu)$ for a barrier of height $B(\nu,z)$ at one arbitrary redshift to infer the mass function n(M) at any other redshift by a simple rescaling. For a given first-upcrossing distribution $f(\nu)$, the differential mass function can be calculated using the relation

\begin{displaymath} n(M)=\frac{\rho_{\rm b}}{M^2}\frac{{\rm d}\!\ln\nu}{{\rm d}\!\ln M}\nu f(\nu), \end{displaymath} (30)

where $\rho_{\rm b}$ is the background density of the Universe.

First, we want to find an accurate fit to the first-upcrossing distribution of a moving barrier which is given by the mass-dependent linear overdensity parameter of the ellipsoidal collapse,

\begin{displaymath} B(\nu)=\delta_{\rm c,~sph}\left(1+0.6536~\nu^{-0.6387}\right) \end{displaymath} (31)

(see Eq. (29)). The parameter $\delta_{\rm c,sph}$ is evaluated at $z_{\rm v}=0$. We ran one million random walks and recorded the first-upcrossing values for $\nu\in[0.01,20]$ in 100 equidistant bins in logarithmic space. The resulting distribution $\nu f(\nu)$ is nicely expressed by the function

\begin{displaymath} \nu f(\nu)=A\left[1+(a\nu)^{-p}\right]\sqrt{\frac{a\nu}{2\pi... ...left(-\frac{a\nu}{2} \frac{B(\nu)}{\delta_{\rm c,sph}}\right), \end{displaymath} (32)

thus, our suggested fitting formula is a mixture of the functional forms proposed by Sheth & Tormen (1999) and Sheth & Tormen (2002). The remaining best-fit parameters are A=0.357, p=0.212 and a=1.171. The result is shown in Fig. 8.

\begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig8} \end{figure} Figure 8:

Comparison between the first-upcrossing distribution of the moving barrier, Eq. (32), inferred from an ensemble of one million random walks, and the fitting formula, Eq. (33).

Open with DEXTER

Second, to find a viable mass function from the first-upcrossing distribution, we proceed as Sheth & Tormen (1999) and Sheth et al. (2001), normalise $f(\nu)$ to unity and rescale the variable a such that we are in agreement with the standard Sheth & Tormen mass function and a mass function based on N-body simulations proposed by Courtin et al. (2010). The latter is based on a first-upcrossing distribution that has the same functional form as that proposed by Sheth & Tormen (1999), but slightly different best-fit parameters,

\begin{displaymath} \nu f_{\rm Courtin}(\nu)=\tilde{A}\left[1+(\tilde{a}\nu)^{-\... ...{\tilde{a}\nu}{2\pi}}\exp\left(-\frac{\tilde{a}\nu}{2}\right), \end{displaymath} (33)

with $\tilde{A}=0.348$, $\tilde{a}=0.695$, and $\tilde{p}=0.1$. Note that in their definition of $\nu$, the linear density contrast $\delta_{\rm c,sph}$ has to be taken at collapse. Normalising the first-upcrossing distribution based on the moving barrier of our ellipsoidal-collapse model to unity yields a rescaled parameter $A\rightarrow A'=1.364~A$. We compare the resulting mass function with those by Sheth & Tormen (1999) and Courtin et al. (2010) for three different redshifts in Fig. 9. The parameter a was rescaled by $a\rightarrow a'=0.625~a$. Deviations from the Sheth & Tormen mass function at high masses occur at large redshifts which is compatible with the Courtin et al. mass function. Overall, our proposed mass function lies in between these two, suggesting that differences between N-body simulations and the Sheth & Tormen mass function might be due to an imprecise treatment of the ellipsoidal-collapse dynamics.

\begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig9} \end{figure} Figure 9:

Comparison of the Sheth & Tormen and the Courtin et al. mass functions with the mass function derived from the rescaled upcrossing distribution, Eq. (33), based on our treatment of the ellipsoidal-collapse dynamics.

Open with DEXTER

4 Conclusions

We have reconsidered the collapse of a homogeneous, triaxial ellipsoid in an expanding background universe and extended the treatment of Bond & Myers (1996) in two ways:

  • We have introduced a physically motivated criterion for the onset of virialisation along each principal axis of a collapsing ellipsoid. We derive this criterion from the tensor virial theorem, demanding that the inertial tensor stabilises and its second time derivative vanishes. This approach is a generalisation of the usual virialisation condition for the spherical-collapse scenario and thus fully consistent with it. It replaces the conventional requirement that virialisation is assumed along an axis i when the dimension-less semi-major axis ai of the ellipsoid along that axis reaches ai(a)=0.177a.

  • We have introduced a hybrid model for the influence of the external gravitational shear field acting on the ellipsoid, supplementing the two different models of Bond & Myers (1996). The principal axes of the shear field evolve with those of the halo before it turns around and decouples from the background expansion, and then continues to evolve linearly while the halo evolves non-linearly and collapses. We have shown that the differences between the hybrid and the non-linear model are relatively small.

For a given initial ellipticity e and prolaticity p, and for a specified virialisation redshift $z_{\rm v}$, the ellipsoidal-collapse model then gives a unique answer for the linear density contrast $\delta _{\rm c}$ at virialisation, as well as for the overdensity $\Delta _{\rm v}$ at that time. The probability distribution for e and p, conditional on the density contrast $\delta$, is determined by the probability distribution of the eigenvalues of the Zel'dovich deformation tensor, as shown by Sheth et al. (2001). It is characterised by variance of the matter-density fluctuations on a scale fixed by the halo mass required. We have shown that the marginalisation over e and p can be simplified by evaluating the ellipsoidal collapse at their mean values $\langle e\rangle$ and $\langle p\rangle=0$.

Our main results are as follows:

  • The collapse parameters $\delta _{\rm c}$ and $\Delta _{\rm v}$ depend only weakly on the model for the external gravitational shear. The hybrid model and the non-linear model by Bond & Myers (1996) give approximately the same results.

  • When supplied with our virialisation condition derived from the tensor virial theorem, the ellipsoidal-collapse model returns values $\delta _{\rm c}$ and $\Delta _{\rm v}$ that differ substantially from those obtained with the spherical-collapse model. Depending on halo mass and redshift, deviations of order $(20\ldots50)\%$ are common. After marginalisation over e and p, $\delta _{\rm c}$ and $\Delta _{\rm v}$ increase with decreasing halo mass and with decreasing virialisation redshift.

  • Both parameters increase with increasing initial ellipticity e and decrease with increasing prolaticity p, as already suggested by Sheth et al. (2001).

  • Our results for $\delta _{\rm c}$ qualitatively confirm the dependence on halo mass and virialisation redshift given by a fitting formula by Sheth et al. (2001). Deviations in particular at low mass and high redshift occur due to the differences in the virialisation condition, the model for the external shear, and the marginalisation over e.

  • The mass function based on our refined treatment of the ellipsoidal-collapse model is in good agreement with those proposed by Sheth & Tormen (1999) for low redshifts and Courtin et al. (2010) for high redshifts. This suggests that differences between the Sheth & Tormen mass function and results from N-body simulations at large redshifts may occur due to an imprecise treatment of the ellipsoidal-collapse dynamics.

We now plan to introduce the ellipsoidal-collapse model into the approach presented by Angrick & Bartelmann (2009) to improve the prediction of the X-ray temperature function for galaxy clusters based on the gravitational potential rather than on the density contrast.

Acknowledgements
C.A. wants to thank Francesco Pace for fruitful discussions, and the Deutsche Forschungsgemeinschaft for financial support under grant number BA 1369/12-1, the Heidelberg Graduate School of Fundamental Physics, and the IMPRS for Astronomy & Cosmic Physics at the University of Heidelberg.

Appendix A: Comparison to a previous study

Here we compare our results for the evolution equations of the collapsing ellipsoid with those presented by Monaco (1997) for a flat $\Lambda $CDM and an OCDM model.

Starting from Eq. (5), we can replace E(a) and E'(a) for a flat $\Lambda $CDM model using

\begin{displaymath} E(a)=\sqrt{\Omega_{\rm m}a^{-3}+(1-\Omega_{\rm m})},\quad E'(a)=-\frac{3\Omega_{\rm m}a^{-4}}{2E(a)} \end{displaymath} (A.1)

since $\Omega_\Lambda=1-\Omega_{\rm m}$ and the curvature parameter $\Omega_{\rm k}=0$. This gives

\begin{displaymath} \frac{{\rm d}^2 a_i}{{\rm d}a^2}\mathbf{-}\frac{1\!-\!2(\Ome... ...^2\left[1\!+\!(\Omega_{\rm m}^{-1}\!-\!1)a^3\right]}a_i\!=\!0. \end{displaymath} (A.2)

This equation differs from Eq. (B11) of Monaco (1997) in the second and third term: a factor a2 in the denominator of both terms was omitted. Additionally, the vacuum term $\propto$ $ \rho_\Lambda$ as well as a factor 3 were not included in the nominator of the third term.

For an OCDM model, we have

$\displaystyle E(a)=\sqrt{\Omega_{\rm m}a^{-3}+(1-\Omega_{\rm m})a^{-2}},$     (A.3)
$\displaystyle E'(a)=-\frac{3\Omega_{\rm m}a^{-4}+2(1-\Omega_{\rm m}) a^{-3}}{2E(a)}$     (A.4)

since $\Omega_\Lambda=0$ and $\Omega_{\rm k}=(1-\Omega_{\rm m})a^{-2}$. Inserting this again into Eq. (5) yields
$\displaystyle \frac{{\rm d}^2 a_i}{{\rm d}a^2}-\left\{2a\left[1+(\Omega_{\rm m}... ...}a} +3\left\{2a^2\left[1+(\Omega_{\rm m}^{-1}-1)a\right]\right\}^{-1}C_i a_i=0.$     (A.5)

Equation (B12) of Monaco (1997) is again slightly different: the factor 3 in the last term was omitted.

There is one last difference concerning the initial conditions: Comparing Eq. (B17) of Monaco (1997) with Eq. (6) of this paper, one can find an additional factor a0 in front of $\lambda_i(a_0)$ which should be dropped.

References

All Figures

  \begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig1.eps}\vspace*{-2mm} \end{figure} Figure 1:

Comparison of the virialisation conditions from the scalar and tensor virial theorems. The curves differ slightly because they refer to different times. While the tensor virial theorem was evaluated at virialisation, the scaler virial theorem was evaluated at collapse. If we choose the collapse of the dark-matter halo as reference for z, both yield exactly the same result.

Open with DEXTER
In the text

  \begin{figure} \par\mbox{\includegraphics[width=6cm,clip]{figures/14147fig2a}\vs... ...ncludegraphics[width=6cm,clip]{figures/14147fig2c} } \vspace*{-2mm} \end{figure} Figure 2:

Evolution of the principal axes and the external shear for a $10^{14}\ M_\odot /h$ dark-matter halo with $e=\langle e\rangle$ and p=0 in the reference $\Lambda $CDM model. Left panel: linear approximation. Central panel: non-linear approximation. Right panel: hybrid approximation.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig3a}\hspace... ...[width=4.3cm,clip]{figures/14147fig3d}\vspace*{-1mm} \vspace*{-4mm} \end{figure} Figure 3:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) for three different models of the external shear. Both quantities were calculated for a $10^{14}\ M_\odot /h$ dark-matter halo in the reference $\Lambda $CDM model. Top panels: $e=\langle e\rangle$ as a function of virialisation redshift $z_{\rm v}$. Bottom panels: e=0.2 for all $z_{\rm v}$. In both cases, p=0.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig4a}\hspace... ...pace*{1.5mm} \includegraphics[width=4.3cm,clip]{figures/14147fig4d} \end{figure} Figure 4:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) for three different cosmological models and fixed initial e and p. Top panels: spherical case (e=p=0). Bottom panels: triaxial case with e=0.2 and p=-0.1.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=4.3cm,clip]{figures/14147fig5a}\hspace... ...ace*{1.5mm} \includegraphics[width=4.3cm,clip]{figures/14147fig5d} \end{figure} Figure 5:

$\delta _{\rm c}$ (left panels) and $\Delta _{\rm v}$ (right panels) as a function of halo mass and virialisation redshift for three different cosmologies marginalised over e and p. Top panels: dependence on virialisation redshift $z_{\rm v}$ for a halo with $M=10^{14}~M_\odot/h$. Bottom panels: dependence on mass M for a halo with $z_{\rm v}=0$. Thin dashed lines show the reference values from spherical collapse in the EdS universe, $\delta _{\rm c}=1.583$ and $\Delta _{\rm v}=147$.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=4.4cm,clip]{figures/14147fig6a}\hspace*{1.5mm} \includegraphics[width=4.2cm,clip]{figures/14147fig6b} \end{figure} Figure 6:

Comparison of the fitting formula of Sheth et al. (2001) for $\delta _{\rm c}(M,z)$ with the results of this work for the reference $\Lambda $CDM model. Left panel: $\delta _{\rm c}$ as a function of mass and z=0. Right panel: $\delta _{\rm c}$ as a function of redshift and $M=10^{14}~M_\odot/h$.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig7} \end{figure} Figure 7:

Influence of the initial ellipticity e and prolaticity p on the parameters $\delta _{\rm c}$ (top panels) and $\Delta _{\rm v}$ (bottom panels) for three different cosmological models. Left panels: $\Lambda $CDM. Central panels: OCDM. Right panels: EdS. The large circle in each figure indicates the value of the respective parameter for $e=\langle e\rangle$ and p=0. The small circles indicate values for a combination of $e\in[\langle e\rangle \pm \sigma_e/2,\ \langle e\rangle \pm \sigma_e]$ and $p\in [\pm e/4,\ \pm e/2]$ (see Eqs. (21), (22)). For a fixed e both $\delta _{\rm c}$ and $\Delta _{\rm v}$ grow for decreasing p.

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig8} \end{figure} Figure 8:

Comparison between the first-upcrossing distribution of the moving barrier, Eq. (32), inferred from an ensemble of one million random walks, and the fitting formula, Eq. (33).

Open with DEXTER
In the text

  \begin{figure} \par\includegraphics[width=9cm,clip]{figures/14147fig9} \end{figure} Figure 9:

Comparison of the Sheth & Tormen and the Courtin et al. mass functions with the mass function derived from the rescaled upcrossing distribution, Eq. (33), based on our treatment of the ellipsoidal-collapse dynamics.

Open with DEXTER
In the text

Copyright ESO 2010

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.