Astronomy. Astrophysics. Mass functions and bias of dark matter halos. P. Valageas. 1. Introduction
|
|
- Abigayle Sanders
- 5 years ago
- Views:
Transcription
1 A&A 508, (2009) DOI: / / c ESO 2009 Astronomy & Astrophysics Mass functions and bias of dark matter halos P. Valageas Institut de Physique Théorique, CEA Saclay, Gif-sur-Yvette, France patrick.valageas@cea.fr Received 14 May 2009 / Accepted 8 October 2009 ABSTRACT Aims. We revisit the study of the mass functions and the bias of dark matter halos. Methods. Focusing on the limit of rare massive halos, we point out that exact analytical results can be obtained for the large-mass tail of the halo mass function. This is most easily seen from a steepest-descent approach, that becomes asymptotically exact for rare events. We also revisit the traditional derivation of the bias of massive halos, associated with overdense regions in the primordial density field. Results. We check that the theoretical large-mass cutoff agrees with the mass functions measured in numerical simulations. For halos defined by a nonlinear threshold δ = 200 this corresponds to using a linear threshold δ L 1.59 instead of the traditional value We also provide a fitting formula that matches simulations over all mass scales and obeys the exact large-mass tail. Next, paying attention to the Lagrangian-Eulerian mapping (i.e. corrections associated with the motions of halos), we improve the standard analytical formula for the bias of massive halos. We check that our prediction, which contains no free parameter, agrees reasonably well with numerical simulations. In particular, it recovers the steepening of the dependence on scale of the bias that is observed at higher redshifts, which published fitting formulae did not capture. This behavior mostly arises from nonlinear biasing. Key words. large-scale structure of Universe methods: analytical gravitation 1. Introduction The distribution of nonlinear virialized objects, such as galaxies or clusters of galaxies, is a fundamental test of cosmological models. First, this allows us to check the validity of the standard cosmological scenario for the formation of large-scale structures, where nonlinear objects form thanks to the amplification by gravitational instability of small primordial density fluctuations, built for instance by an early inflationary stage (e.g., Peebles 1993; Peacock 1998). For cold dark matter (CDM) scenarios (Peebles 1982), where the amplitude of the initial perturbations grows at smaller scales, this gives rise to a hierarchical process, where increasingly large and massive objects form as time goes on, as increasingly large scales turn nonlinear. This process has been largely confirmed by observations, which find smaller galaxies at very high redshifts (e.g., Trujillo et al. 2007) while massive clusters of galaxies (which are the largest bound objects in the Universe) appear at low redshifts (e.g., Borgani et al. 2001). Second, on a more quantitative level, statistical properties, such as the mass function and the twopoint correlation of these objects, provide strong constraints on the cosmological parameters (e.g. through the linear growth factor D + (t) of density perturbations) and on the primordial fluctuations (e.g. through the initial density power spectrum P L (k)). For these purposes, the most reliable constraints come from observations of the most massive objects (rare-event tails) at the largest scales. Indeed, in this regime the formation of large-scale structures is dominated by the gravitational dynamics (baryonic physics, which involves intricate processes associated with pressure effects, cooling and heating, mostly occurs at galactic scales and below), which further simplifies as one probes quasi-linear scales or rare events where effects associated with multiple mergings can be neglected. Moreover, in this regime astrophysical objects, such as galaxies or clusters of galaxies, can be directly related to dark matter halos, and their abundance is highly sensitive to cosmological parameters thanks to the steep decline of the high-mass tail of the mass function (e.g., Evrard 1989). Thus, the computation of the halo mass function (and especially its large-mass tail) has been the focus of many works, as it is one of the main properties measured in galaxy and cluster surveys that can be compared with theoretical predictions. Most analytical derivations follow the Press-Schechter approach (Press & Schechter 1974; Blanchard et al. 1992) or its main extension, the excursion set theory of Bond et al. (1991). In this framework, one attempts to estimate the number of virialized objects of mass M from the probability to have a linear density contrast δ L at scale M above some given threshold δ c. Thus, one identifies current nonlinear halos from positive density fluctuations in the initial (linear) density field, on a one-to-one basis. This is rather well justified for rare massive objects, where one can expect such a link to be valid since such halos should have remained well-defined objects until now (as they should have suffered only minor mergers). By contrast, small and typical objects have experienced many mergers and should be sensitive to highly non-local effects (e.g. tidal forces, mergers), so that such a direct link should no longer hold, as can be checked in numerical simulations (Bond et al. 1991). As noticed by Press & Schechter (1974), the simplest procedure only yields half the mass of the Universe in such objects (essentially because only half of the Gaussian initial fluctuations have a positive density contrast, whatever the smoothing scale), which they corrected by an ad-hoc multiplicative factor 2. In this respect the main result of the excursion set theory was to provide an analytical derivation of this missing factor 2, in the simplified case of a top-hat filter in Fourier space (Bond et al. 1991). Then, it arises from Article published by EDP Sciences
2 94 P. Valageas: Mass functions and bias of dark matter halos the fact that objects of mass larger than M are associated with configurations such that the linear density contrast goes above the threshold δ c at some scale M M, which includes cases missed by the Press-Schechter prescription where the linear density contrast decreases below this scale M so that δ L <δ c at scale M ( cloud-in-cloud problem). The characteristic threshold δ c is usually taken as the linear density contrast reached when the spherical collapse dynamics predicts collapse to a zero radius. In an Einstein-de Sitter universe this corresponds to δ c and to a nonlinear density contrast δ 177 (assuming full virialization in half the turn-around radius). This linear threshold only shows a very weak dependence on cosmological parameters. Numerical simulations have shown that the Press-Schechter mass function (PS) is reasonably accurate, especially in view of its simplicity. Thus, it correctly predicts the typical mass scale of virialized halos at any redshift, as could be expected since the large-mass behavior is rather well justified. In addition, it predicts a universal scaling that appears to be satisfied by the mass functions measured in simulations, that is, the dependence on halo mass, redshift and cosmology is fully contained in the ratio ν = δ c /σ(m, z), where σ(m, z) is the rms linear density fluctuation at scale M. However, as compared with numerical results it overestimates the low-mass tail and it underestimates the highmass tail. This has led to many numerical studies which have provided various fitting formulae for the mass function of virialized halos, written in terms of the scaling variable ν (or σ) (Sheth & Tormen 1999; Jenkins et al. 2001; Reed et al. 2003; Warren et al. 2006; Tinkeretal.2008). We may note here that a theoretical model that attempts to improve over the PS mass function is to consider the ellipsoidal collapse dynamics within the excursion-set approach, to take into account the deviations from spherical symmetry for intermediate mass halos (Sheth et al. 2001). Note that, as described in the original paper (and emphasized in Robertson et al. 2009), at large mass this would recover the spherical collapse. However, the halo mass obtained by such methods is generally underestimated for non center-ofmass particles (since analytical computations assume that particles are located at the center of their halo, i.e. they only consider the linear densities within spherical cells centered on the point of interest). To correct for this effect, in practice one treats the threshold δ c as a free parameter, close to 1.6, to build fitting formulae from numerical simulations. For instance, this correction is contained in the parameter a in the exponential cutoff of Sheth et al. (2001). A second property of virialized halos that can be used to constrain cosmological models, beyond their number density, is their two-point correlation. Indeed, observations show that galaxies and clusters do not obey a Poisson distribution but show significant large-scale correlations (e.g., McCracken et al. 2008; Padilla et al. 2004). In particular, their two-point correlation roughly follows the underlying matter correlation, up to a multiplicative factor b 2, called the bias, that grows for more massive and extreme objects. Following the spirit of the Press-Schechter picture, Kaiser (1984) found that this behavior arises in a natural fashion if halos are associated with large overdensities in the Gaussian initial (linear) density field, above the threshold δ c. This was further expanded by Bardeen et al. (1986) and Bond &Myers(1996), who considered the clustering of peaks in the Gaussian linear density field. A simpler derivation, based on a peak-background split argument, and taking care of the mapping from Lagrangian to Eulerian space, was presented in Mo & White (1996). It provides a prediction for the bias b(m) asa function of halo mass, in the limit of large distance, r,that agrees reasonably well with numerical simulations. However, as for the PS mass function, in order to improve the agreement with numerical results various fitting formulae have been proposed (Sheth & Tormen 1999; Hamana et al. 2001; Pillepich et al. 2009). Again, since the ellipsoidal collapse model reduces to the spherical dynamics for rare massive halos, it also requires free parameters to improve its accuracy, but the latter are consistent with those used for the mass function (Sheth et al. 2001). In this article we revisit the derivation of the mass function and the bias of rare massive halos, following the spirit of the Press & Schechter (1974) andkaiser(1984) approaches. That is, we use the fact that large halos can be identified from overdensities in the Gaussian initial (linear) density field. First, we briefly review in Sect. 2 some properties of the growth of linear fluctuations and of the spherical dynamics in ΛCDM cosmologies. Next, we recall in Sect. 3 that in the quasi-linear regime (i.e. at large scales), the probability distribution of the nonlinear density contrast δ r within spherical cells of radius r can be obtained from spherical saddle-points of a specific action S, for moderate values of δ r where shell-crossing does not come into play. We also discuss the properties of these saddle-points as a function of mass, scale and redshift. Then, we point out in Sect. 4 that this provides the exact exponential tail of the halo mass function. This applies to any nonlinear density contrast threshold δ that is used to define halos, provided it is below the upper bound δ + where shell-crossing comes into play. We compare our results with numerical simulations and we give a fitting formula that applies over all mass scales and satisfies the exact large-mass cutoff. Next, we recall in Sect. 5 that these results also provide the density profile of dark matter halos at outer radii (i.e. beyond the virial radius) in the limit of large mass. Finally, we study the bias of massive halos in Sect. 6, paying attention to some details such as the Lagrangian-Eulerian mapping, and we compare our results with numerical simulations. We conclude in Sect Linear perturbations and spherical dynamics We consider in this article a flat CDM cosmology with two components, (i) a non-relativistic component (dark and baryonic matter, which we do not distinguish here) that clusters through gravitational instability; and (ii) an uniform dark energy component that does not cluster at the scales of interest, with an equation-of-state parameter w = p de /ρ de. For the numerical computations we shall focus on a ΛCDM cosmology, where the dark energy is associated with a cosmological constant that is exactly uniform with w = 1. However, our results directly extend to curved universes (i.e. Ω k 0) and to dark energy models with a possibly time-varying w(z), as long as we can neglect the dark energy fluctuations on the scales of interest, which is valid for realistic cases. Focussing on the case of constant w, wefirst recall in this section the equations that describe the dynamics of the background and of linear matter density perturbations, as well as the nonlinear spherical dynamics. The evolution of the scale factor a(t) is determined by the Friedmann equation (Wang & Steinhardt 1998), H 2 (t) H 2 0 =Ω m0 a 3 +Ω de0 a 3 3w, with H(t) = ȧ a, (1) where subscripts 0 denote current values at z = 0, when a = 1, and a dot denotes the derivative with respect to cosmic time t. On the other hand, the density parameters vary with time as Ω m (a) = Ω m0 Ω m0 +Ω de0 a, Ω Ω de0 de(a) = (2) 3w Ω m0 a 3w +Ω de0
3 Next, introducing the matter density contrast, δ(x, t) = (ρ m (x, t) ρ m )/ρ m,whereρ m (t) is the mean matter density, linear density fluctuations grow as δ L + 2 ȧ a δ L 4πGρ m δ L = 0, (3) where the subscript L denotes linear quantities. For numerical purposes it is convenient to use the logarithm of the scale factor as the time variable. Then, using the Friedmann Eq. (1) the linear growth factor D + (t) evolves from Eq. (3)as [ 1 D ] 2 w Ω de D Ω m D + = 0, (4) P. Valageas: Mass functions and bias of dark matter halos 95 where we note with a prime the derivative with respect to ln a, as D + = dd + /dlna. Then, the normalized linear growth factor, g(t), defined as g(t) = D +(t) and g(t = 0) = 1, (5) a(t) obeys [ 5 g ] 2 w Ω de g (1 w)ω de g = 0, (6) with the initial conditions g 1andg 0ata 0. In the following we shall also need the dynamics of spherical density fluctuations. For such spherically symmetric initial conditions, the physical radius r(t), that contains the constant mass M until shell-crossing, evolves as r = 4πG 3 r [ρ m + (1 + 3w)ρ de ] with ρ m = 3M 4πr (7) 3 Note that ρ m is the mean density within the sphere of radius r (and notthe local density at radiusr). Using again the Friedmann Eq. (1)thisreadsas r 3 2 (1 + w Ω de)r + ρ m + (1 + 3w)ρ de r = 0. (8) 2(ρ m + ρ de ) As for the linear growth factor D + (t), it is convenient to introduce the normalized radius y(t)definedas y(t) = r(t) and q(t) a(t), y(t = 0) = 1. (9) q(t) Thus, q(t) is the Lagrangian coordinate of the shell r(t), that is, the physical radius that would enclose the same mass M in a uniform universe with the same cosmology. This also implies that the density ρ m writes as ρ m (t) = ρ m (t) y(t) 3. (10) Then, Eq. (8) leads to [ 1 y ] 2 w Ω de y + Ω ( m y 3 1 ) y = 0. (11) 2 Of course, we can check that in the linear regime, where y L = 1 δ L /3, we recover Eq. (4) for the linear growth of δ L. Then, to obtain the nonlinear density contrast, δ(z) = ρ m /ρ m 1, associated with the linear density contrast, δ L (z), at a given redshift z, wesolveeq.(11) with the initial condition y(z i ) = 1 δ Li /3andy (z i ) = δ Li /3 at some high redshift z i, with δ Li /δ L = D + (z i )/D + (z). For any redshift z this defines a mapping Fig. 1. The function F (δ L ) that describes the spherical dynamics when there is no shell-crossing, at z = 0. The solid line corresponds to the ΛCDM cosmology (with Ω m0 = 0.27) and the dashed line to the Einstein-de Sitter case Ω m0 = 1. The second peak corresponds to a second collapse to the center, but for our purposes we only need F (δ L ) before first collapse. Fig. 2. The linear density contrast δ L = F 1 (δ) associated with the nonlinear density contrast δ, through the spherical dynamics, as a function of the cosmological parameter Ω m (z) at the redshift of interest. We show the three cases δ = 100, 200 and 300. The dashed line is the fit (12) to the case δ = 200. δ L δ = F (δ L ) that fully describes the spherical dynamics before shell-crossing. We compare in Fig. 1 the function F (δ L ) obtained at z = 0 within the ΛCDM cosmology that we consider in this article (solid line) with the Einstein-de Sitter case (dashed line), where it has a well-known parametric form (Peebles 1980). As is well known, we can check that the dependence on Ω m0 is very weak until full collapse to a point, which occurs at slightly lower values of δ L for low Ω m0.weshowinfig.2 the linear density contrasts, δ L = F 1 (δ), associated with three nonlinear density contrasts, δ = 100, 200 and 300, as a function of the cosmological parameter Ω m (z). In agreement with Fig. 1, they show a slight decrease for smaller Ω m.forδ = 200 a simple fit (dashed line) is provided by δ = 200: δ L (Ω m ) 0.24, (12) which agrees with the exact curve to better than for Ω m > 0.2, which covers the range of practical interest.
4 96 P. Valageas: Mass functions and bias of dark matter halos In this article we consider Gaussian initial fluctuations, which are fully defined by the linear density power spectrum δ L (k 1 ) δ L (k 2 ) = δ D (k 1 + k 2 ) P L (k 1 ), (13) where δ D is the Dirac distribution and we normalized the Fourier transform as δ(x) = dk e ik.x δ(k), (14) where x and k are the comoving spatial coordinate and wavenumber. This gives for the linear density two-point correlation C L (x 1, x 2 ) = δ L (x 1 )δ L (x 2 ) = 4π dkk 2 P L (k) sin(k x 2 x 1 ) (15) k x 2 x 1 We also introduce the smoothed density contrast, δ r (x), within the sphere of radius r and volume V around position x, dx δ r (x) = V V δ(x + x ) = dk e ik.x δ(k) W(kr), (16) with a top-hat window that reads in Fourier space as dx sin(kr) kr cos(kr) W(kr) = V V eik.x = 3 (17) (kr) 3 Then, in the linear regime, the cross-correlation of the smoothed linear density contrast at scales r 1 and r 2 and positions x 1 and x 2 = x 1 + x reads as σ 2 r 1,r 2 (x) = δ Lr1 (x 1 )δ Lr2 (x 1 + x) = 4π dkk 2 P L (k) W(kr 1 ) W(kr 2 ) sin(kx) (18) kx In particular, σ r = σ r,r (0) is the usual rms linear density contrast at scale r. 3. Distribution of the density contrast We recall here that in the quasi-linear limit, σ r 0, the distribution P(δ r ) of the nonlinear density contrast δ r at scale r can be derived from a steepest-descent method (Valageas 2002a). In agreement with the alternative approach of Bernardeau (1994a), this shows that rare-event tails are dominated by spherical saddle-points, which we use in Sects. 4 6 to obtain the properties of massive halos. Since the system is statistically homogeneous we can consider the sphere of radius r centered on the origin x = 0. Then, we first introduce the cumulant generating function ϕ(y), e ϕ(y)/σ2 r = e yδ r/σ 2 r = dδ r e yδ r/σ 2 r P(δr ), (19) 1 which determines the distribution P(δ r ) through the inverse Laplace transform +i dy P(δ r ) = i 2πiσ 2 e [yδ r ϕ(y)]/σ 2 r. (20) r In Eq. (19) we rescaled the cumulant generating function by a factor σ 2 r so that it has a finite limit in the quasi-linear limit σ r 0 for the Gaussian initial density fluctuations that we study in this article (Bernardeau et al. 2002). In particular, its expansion at y = 0 reads as ϕ(y) = ( y) p p=2 p! δr p c σr 2(p 1) (21) Then, the average (19) can be written as the path integral e ϕ(y)/σ2 r = (det CL 1 )1/2 Dδ L (q) e S[δ L]/σ 2 r, (22) where CL 1 is the inverse matrix of the two-point correlation (15) and the action S reads as S[δ L ] = yδ r [δ L ] + σ2 r 2 δ L.CL 1.δ L (23) Here δ r [δ L ] is the nonlinear functional that affects to the initial condition defined by the linear density field δ L (q) the nonlinear density contrast δ r within the sphere of radius r, and we note δ L.CL 1.δ L = dq 1 dq 2 δ L (q 1 )CL 1(q 1, q 2 )δ L (q 2 ). Note that the action S does not depend on the normalization of the linear power spectrum since both σ 2 r and C L are proportional to P L. In the quasi-linear limit, σ r 0, as shown in Valageas (2002a) the path integral (22) is dominated by the minimum 1 of the action S[δ L ], σ r 0: ϕ(y) min S[δ L]. (24) δ L (q) Using the spherical symmetry of the top-hat window W, in agreement with the approach of Bernardeau (1994a), one obtains a spherical saddle-point with the radial profile σ 2 q,q δ Lq = δ Lq, (25) σ 2 q where σ q,q = σ q,q (0, 0) and q is the Lagrangian coordinate associated with the Eulerian radius r through the spherical dynamics recalled in Sect. 2, andq is a dummy Lagrangian coordinate along the radial profile (hereafter we denote by the letter q Lagrangian radii, which are associated with the linear density field, whereas r denotes Eulerian radii associated with the nonlinear density contrast δ r ). Thus, the radii q and r are related by q 3 = (1 + δ r ) r 3, with δ r = F (δ Lq ), (26) where the function δ r = F (δ Lq ), that describes the spherical dynamics, was obtained below Eq. (11)andshowninFig.1. We can note that the profile (25) is also the mean conditional profile of the linear density contast δ Lq, under the constraint that is it equal to a given value δ Lq at a given radius q (e.g., Bernardeau 1994a). The reason why the nonlinear dynamics gives back this result is that the nonlinear density contrast δ r only depends on the linear density contrast δ Lq at the associated Lagrangian radius q, through the mapping δ r = F (δ Lq ). Indeed, as long as shell-crossing does not modify the mass enclosed within the shell of Lagrangian coordinate q, its dynamics is independent of the motion of inner and outer shells (thanks to Gauss theorem). Then, in order to obtain the minimum of the action S we could proceed in two steps. First, for arbitrary Lagrangian radius q and linear contrast δ Lq, we minimize S with 1 As described in details in Valageas (2002a), depending on the slope of the linear power spectrum the saddle-point (25) may only be a local minimum or maximum, but it still governs the rare-event limit, in agreement with physical expectations.
5 P. Valageas: Mass functions and bias of dark matter halos 97 Fig. 3. Theradialprofile(25) of the linear density contrast δ Lq of the saddle point of the action S[δ L (q)]. We show the profiles obtained with a ΛCDM cosmology for the masses M = 10 11, 10 12, 10 13, and h 1 M. A larger mass corresponds to a lower ratio δ Lq /δ Lq at large radii q /q > 1. respect to the profile δ Lq over q q. From the previous argument, only the second Gaussian term in (23) variessothat this partial minimization leads to the profile (25) (indeed, for Gaussian integrals the saddle-point method is exact). Second, we minimize over the Lagrangian radius q (or equivalently over δ L or δ r ),whichleadstoeq.(29) below. Here we also use the fact that a spherical saddle-point with respect to spherical fluctuations is automatically a saddle-point with respect to arbitrary non-spherical perturbations and it can be seen that for small y one obtains a minimum (then we assume that at finite y strong deviations from spherical symmetry do not give rise to deeper minima, which seems natural from physical expectations). We refer to Valageas (2002a, 2009b) for more detailed derivations. Note that the shape of the linear profile (25) depends on the shape of the linear power spectrum, whence on the mass scale M of the saddle point for a curved CDM linear power spectrum, but not on redshift. We show in Fig. 3 the profile (25) obtained for several masses M. For a power-law linear power spectrum, of slope n, Eq.(25) leads to δ Lq q (n+3) at large radii, q. Then, since for a CDM cosmology n increases at larger scales, the profile shows a steeper falloff at large radii for larger mass, in agreement with Fig. 3 (in this section we consider a ΛCDM cosmology with (Ω m0, Ω Λ0,σ 8, n s, h) = (0.27, 0.73, 0.79, 0.95, 0.7)). We show in Fig. 4 the Lagrangian map, q r,givenby the spherical dynamics (i.e. the function F (δ L ) where we neglect shell-crossing)for the saddle-point (25), with a nonlinear density contrast δ = 200 at the Eulerian radius r, atredshiftz = 0. Inner shells have already gone once through the center but they have not reached radius r yet. Even though their dynamics is no longer exactlygivenbyeq.(11), an exact computation would give the same property as the increasing mass seen by these particles, as they pass outer shells, should slow them down as compared with the constant-mass dynamics. In agreement with Fig. 3,forlarger masses, which have a larger central linear density contrast, shellcrossing has moved to larger radii (the local maximum of r /r, to the left of r = 0, is higher). From the Lagrangian map, q r, we define the nonlinear density threshold, δ +, as the nonlinear density contrast δ r reached within radius r at the time when inner shells first cross this radius r. Then, up to δ +, the mass within the Lagrangian Fig. 4. The Lagrangian map q r given by the spherical dynamics (neglecting shell-crossing) for the saddle-point (25) with a nonlinear density contrast δ = 200 at the Eulerian radius r. We plot the curves obtained at z = 0 for several masses, as in Fig. 3. Inner shells have already gone once through the center (hence their dynamics is no longer exactly given by Eq. (11)) but they have not crossed the radius r yet. Fig. 5. The nonlinear density contrast δ +, beyond which shell-crossing must be taken into account, as a function of redshift for several mass scales. shell q has remained constant, so that the saddle-point (25) (26) is exact (this slightly underestimates δ + as the expansion of inner shells should be somewhat slowed down by the mass of the outer shells they have overtaken). We show in Fig. 5 the dependence on redshift of this threshold δ +, for several masses. In agreement with Figs. 3, 4, this threshold is smaller for larger masses. It shows a slight decrease at higher redshift as Ω m (z) growsto unity. We can see that for massive clusters at z = 0, which have a mass of order h 1 M, the density contrast δ 200 should separate outer shells with radial accretion from inner shells with a significant transverse velocity dispersion, built by the radialorbit instability that dominates the dynamics after shell-crossing, see Valageas 2002b (in particular the Appendix A). This agrees with numerical simulations, see for instance the lower panel of
6 98 P. Valageas: Mass functions and bias of dark matter halos Fig. 3 in Cuesta et al. (2008), which show that rare massive clusters exhibit a strong radial infall pattern, with a low velocity dispersion, beyond the virial radius (where δ r 200), while inner radii show a large velocity dispersion (even though we can distinguish close to the virial radius the outward velocity associated with the shells that have gone once through the center). Small-mass halos do not show such a clear infall pattern and the velocity dispersion is significant at all radii (e.g., upper panel of Fig. 3 in Cuesta et al. 2008). This expresses the fact that such halos, associated with typical events and moderate density fluctuations, are no longer governed by the spherical saddle-point (25). Indeed, at low mass and small Lagrangian radius q, σ q is no longer very small so that the path integral (22) is no longer dominated by its minimum and one must integrate over all typical initial conditions, including large deviations from spherical symmetry. Next, the amplitude δ Lq and the minimum ϕ are given as functions of y by the implicit equations (Legendre transform) y = τ dτ dg and ϕ = yg + τ2 2, (27) where the function τ(g) is defined from the spherical dynamics through the parametric system (Valageas 2002a; Bernardeau 1994b) σ r G = δ r = F (δ Lq ) and τ = δ Lq (28) σ q The system (27) also reads as ϕ = min G [ yg + τ2 (G) 2 ] (29) Note that the function τ(g), whence the cumulant generating function ϕ(y), depends on the shape of the linear power spectrum through the ratio of linear power σ r /σ q at Eulerian and Lagrangian scales r and q. Finally, from the inverse Laplace transform (20), a second steepest-descent integration over the variable y gives (Valageas 2002a) σ r 0: P(δ r ) e τ2 /(2σ 2 r ) = e δ2 Lq /(2σ2 q). (30) Thus, in the quasi-linear limit, the distribution P(δ r )isgoverned at leading order by the Gaussian weight of the linear fluctuation δ Lq that is associated with the nonlinear density contrast δ r through the spherical dynamics. We can note that Eq. (30) could also be obtained from a Lagrange multiplier method, without introducing the generating function ϕ(y). Indeed, in the rare-event limit, where P(δ) is governed by a single (or a few) initial configuration, we may write rare events: P(δ) max e 1 2 δ L.C 1 L.δ L. (31) {δ L [q] δ r [δ L ]=δ} That is, P(δ) is governed by the maximum of the Gaussian weight e (δ L.C 1 L.δL)/2 subject to the constraint δ r [δ L ] = δ (assuming there are no degenerate maxima). Then, we can obtain this maximum by minimizing the action S[δ L ]/σ 2 r of Eq. (23), where y plays the role of a Lagrange multiplier. This gives the saddle-point (25), and the amplitude δ Lq and the radius q are directly obtained from the constraint δ = F (δ Lq ), as in Eq. (26). Then, we do not need the explicit expression of the Lagrange multiplier y, asthisissufficient to obtain the last expression of the asymptotic tail (30). Nevertheless, it is usefulto introducethe generating function ϕ(y), which makes it clear that the Lagrange multiplier y is also the Laplace conjugate of the nonlinear density contrast δ as in Eq. (19), since it is also of interest by itself, as it yields the density cumulants through the expansion (21). The asymptotic (30) holds in the rare-event limit. This corresponds to both the quasi-linear limit, σ r 0atfixeddensity contrast δ r, and to the low-density limit, δ r 0atfixedσ r,as long as there is no shell-crossing. This latter requirement gives a lower boundary δ, for linear power spectra with a slope n > 1, and an upper boundary δ +, for any linear spectrum (Valageas 2002b). This upper boundary was shown in Fig. 5 for a ΛCDM cosmology, for several masses. Indeed, at large positive density contrast, shell-crossing always occurs, as seen in Figs. 4 and 5. This invalidates the mapping δ L δ = F (δ L ) obtained from Eq. (11), as mass is no longer conserved within the Lagrangian shell q, so that the asymptotic behavior (30) is no longer exact. In fact, as shown in Valageas (2002b), after shell-crossing it is no longer sufficient to follow the spherical dynamics, even if we take into account shell-crossing. Indeed, a strong radial-orbit instability develops so that the sensitivity to initial perturbations actually diverges when particles cross the center of the halo. Then, the functional δ r [δ L (q)] is singular at such spherical states (i.e. it is discontinuous as infinitesimal deviations from spherical symmetry lead to a finite change of δ r ) and the path integral (22) is no longer governed by spherical states that have a zero measure. As noticed above, this also means that, in the limit of rare massive halos, the nonlinear density threshold δ + separates outer shells with a smooth radial flow from inner shells with a significant transverse velocity dispersion. Thus, δ + marks the virialization radius where isotropization of the velocity tensor becomes important, in agreement with numerical simulations (e.g., Cuesta et al. 2008). It is interesting to note that a similar approach can be developed for the adhesion model (Gurbatov et al. 1989), where particles move according to the Zeldovich dynamics (Zeldovich 1970) but do not cross because of an infinitesimal viscosity (i.e. they follow the Burgers dynamics in the inviscid limit). Moreover, in the one-dimensional case, with a linear powerspectrum slope n = 2 orn = 0 (i.e. the linear velocity is a Brownian motion or a white noise), it is possible to derive the exact distribution P(δ r ) by other techniques and to check that it agrees with the asymptotic tail (30), as seen in Valageas (2009a,b). 4. Mass function of collapsed halos The quasi-linear limit (30) ofthedistributionp(δ r ) clearly governs the large-mass tail of the mass function n(m)dm, where we define halos as spherical objects with a fixed density contrast δ r = δ. Indeed, since the Lagrangian radius q is related to the halo mass by M = ρ m 4πq 3 /3, massive halos correspond to large Eulerian and Lagrangian radii, and the limit M corresponds to the quasi-linear limit σ r 0. Then, going from the distribution P(δ r ) to the mass function n(m) can introduce geometrical power-law prefactors since halos are not exactly centeredonthecellsofafixedgrid,asdiscussedinbetancort-rijo & Montero-Dorta (2006), but the exponential cutoff remains the same as in Eq. (30), whence M : ln[n(m)] δ2 L 2σ 2, with δ L = F 1 (δ), (32) where σ = σ q with M = ρ m 4πq 3 /3. A simple approximation for the mass function that satisfies this large-mass falloff can
7 P. Valageas: Mass functions and bias of dark matter halos 99 be obtained from the Press-Schechter method (PS), see Press & Schechter (1974), but using the actual linear threshold δ L = F 1 (δ) associated with the nonlinear threshold δ that defines the halo radius r, rather than the linear threshold δ c associated with the spherical collapse down to a point. This yields for the number of halos in the range [M, M + dm] per unit volume, n(m)dm = ρ m M f (σ) dσ σ, (33) with 2 δ L f (σ) = π σ e δ2 L /(2σ2). (34) The mass function (34) includes the usual prefactor 2 that gives the normalization dσ f (σ) = 1, (35) 0 σ which ensures that all the mass is contained in such halos. However, we must note that the power-law prefactor in (34) has no strong theoretical justification since only the exponential cutoff (32) was exactly derived from (30). In principle, it could be possible to derive subleading terms (whence powerlaw prefactors) in the quasi-linear limit for P(δ r ), by expanding the path integral (22) around the saddle-point (25). Then, one would also need to take more care of the prefactors that would arise as one goes from the distribution P(δ r ) to the mass function n(m). However, this expansion would not allow one to derive the shape of the mass function at low masses, as this regime is far from the quasi-linear limit and involves multiple mergers far from spherical symmetry. Until a new method is devised to handle this regime, one must treat the prefactors as free parameters to be fitted to numerical simulations, in order to describe the mass function from small to large masses. Here we can note that in the Press-Schechter approach (and in most models) the mass function (34) is obtained from the linear density reached within the sphere centered on each mass element. That is, a particle is assumed to be part of a halo of mass larger than M if the sphere of mass M (or a sphere of mass larger than M in the excursion set approach of Bond et al. 1991), centered on this particle, has a linear density contrast above a threshold δ L. As pointed out in Sheth et al. (2001), and discussed in their Sect. 3, since all particles are not located at the center of their parent halo, the correct criteria should rather be that the particle belongs to a sphere, not necessarily centered on this point, of mass greater than M, that has collapsed by the redshift of interest. Within the excursion set approach of Bond et al. (1991), this means that one should consider the first-crossing distributions associated with center-of-mass particles, rather than with randomly chosen particles, as argued in Sheth et al. (2001). Then, the latter suggest that this could modify the numerical factor in the exponential tail of the mass function, that is, lead to a smaller factor δ L in Eq. (32) than the one associated with spherical (or ellipsoidal) collapse dynamics. We point out here that this effect should only give subleading corrections to the tail (32), so that the factor δ L in Eq. (32) remains exactly given by the spherical collapse dynamics. This can be seen from the fact that the same effect would apply to the density probability distribution P(δ r ), as randomly placed cells of radius r are typically not centered on halo profiles. Nevertheless, the results (27) (31) are exact in the quasilinear limit (as can also be checked by the comparison with standard perturbation theory for the cumulant generating function ϕ(y)), and off-center effects would be included in the subleading terms, computed as usual by expanding the path integral (22) around its saddle-point, which would typically give the power-law prefactor to the tail (30). In terms of the halo mass function itself, this point was also studied in Betancort- Rijo & Montero-Dorta (2006), who found as expected that at large mass such a geometrical factor only modifies that powerlaw prefactor. As for the probability distribution P(δ r ), it is interesting to note that a similar high-mass tail can be derived for the adhesion model, where halos are defined as zero-size objects (shocks). Again, in the one-dimensional case, for both n = 2 and n = 0, where the exact mass function can be obtained by other means, one can check that it agrees with the analog of the asymptotic tail (32) (Valageas 2009a,b,c). Thus, we can check that in these two non-trivial examples the leading-order terms for the large-mass decays of P(δ r )andn(m) are exactly set by saddle-point properties and are not modified by the off-center effects discussed above. We can note that the explanation of the large-mass tail (32) by the exact asymptotic result of the steepest-descent method describedinsect.3agrees nicely with numerical simulations. This is most clearly seen in Figs. 3 and 6 of Robertson et al. (2009), who trace back the linear density contrast δ L of the Lagrangian regions that form halos at z = 0. Their results show that the distribution of linear contrasts δ L, measured as a function of mass (or of σ(m)), has a roughly constant lower bound δ L, with δ L 1.6, and an upper bound δ+ Lthat grows with σ(m). We must note however that the difference between 1.59 and (which would be the standard threshold in their ΛCDM cosmology with Ω m0 = 0.27) is too small to discriminate both values from the results shown in their figures, so that these numerical simulations alone do not give the asymptotic value δ L to better than about 0.2. They obtain the same results when they define halos by nonlinear density contrasts δ = 100 or δ = 600, with a lower bound δ L that grows somewhat with δ. Intermsof the approach described above, this behavior expresses the fact that the most probable way to build a massive halo of nonlinear density contrast δ is to start from a Lagrangian region of linear density contrast δ L = F 1 (δ), which obeys the spherical profile (25) and corresponds to the saddle-point of the action (23). As recalled in Sect. 3, the path integral (22) is increasingly sharply peaked around this initial state at large mass scales, which explains why the dispersion of linear density contrasts δ L measured in the simulations decreases with σ(m). At smaller mass, one is sensitive to an increasingly broad region around the saddle-point, which mostly includes non-spherical initial fluctuations. Since these initial conditions are less efficient to concentrate matter in a small region one needs a larger linear density contrast δ L to reach the same nonlinear threshold δ within the Eulerian radius r, which is why the distribution is not symmetric and mostly broadens by increasing its typical upper bound δ + L.In principles, it may be possible to estimate the width of this distribution, at large masses, by expanding the action S[δ L ] around its saddle-point. We compare in Fig. 6 the prediction (34) (solid line labeled as δ L = F 1 (δ)) with results from numerical simulations, for the nonlinear threshold δ = 200 at redshift z = 0. In our case, this corresponds to a linear density contrast δ L 1.59, that is obtained from Fig. 2 or the fit (12). As in Sect. 3, we consider a ΛCDM cosmology with (Ω m0, Ω Λ0,σ 8, n s, h) = (0.27, 0.73, 0.79, 0.95, 0.7). This corresponds to the cosmological parameters of the largest-box numerical simulations of
8 100 P. Valageas: Mass functions and bias of dark matter halos Fig. 6. The mass function at redshift z = 0 of halos defined by the nonlinear density contrast δ = 200. The points are the fits to numerical simulations from Sheth & Tormen (1999), Jenkins et al. (2001), Reed et al. (2003), Warren et al. (2006) and Tinker et al. (2008). The dashed line ( δ c ) is the usual Press-Schechter mass function, while the solid line that goes close to the Press-Schechter prediction at small mass is the mass function (34) with the exact cutoff (32). Here we have δ L The second solid line that agrees with simulations over the whole range is the fit (36), that obeys the same large-mass exponential cutoff. Tinker et al. (2008), which allow the best comparison with the theoretical predictions, as they also define halos by the density threshold δ = 200 with a spherical-overdensity algorithm. The numerical results are the fits to the mass function given in Sheth &Tormen(1999) ( ST ), Jenkins et al. (2001) ( J ), Reed et al. (2003) ( R ), Warren et al. (2006) ( W ) and Tinker et al. (2008) ( T ). Note that these mass functions are defined in slightly different fashions, using either a spherical-overdensity or friendsof-friends algorithm, and density contrast thresholds that vary somewhat about δ = 200. However, they agree rather well, as the dependence on δ is rather weak. This can be understood from Fig. 1, which shows that around δ = 200 the linear density contrast δ L = F 1 (δ) has a very weak dependence on δ. Wealso plot the usual Press-Schechter prediction (dashed line labeled δ c ), that amounts to replace δ L by δ c = in Eq. (34)(since Ω m = 0.27 at z = 0). We can see that using the exact value δ L = F 1 (δ) significantly improves the agreement with numerical simulations at large masses. Note that there are no free parameters in Eq. (34). Of course, at small masses the mass function (34) closely follows the usual Press-Schechter prediction and shows the same level of disagreement with numerical simulations. This is expected since only the exponential cutoff (32) has been exactly derived from Sect. 3. Since at large masses the power-law prefactor 1/σ in Eq. (34) is also unlikely to be correct, as discussed above, we give a simple fit that matches the numerical simulations from small to large masses (solid line that runs through the simulation points) while keeping the exact exponential cutoff: f (σ) = 0.5 [ (0.6 ν) (0.62 ν) 0.5] e ν2 /2, (36) with ν = δ L σ, δ L = F 1 (δ) 1.59 for δ = 200, Ω m = (37) At higher redshift or for other Ω m one simply needs to use the relevant mapping F 1 (δ), shown in Fig. 2 or given by the simple fit (12) forδ = 200. This mass function also satisfies the normalization (35), whatever the value of the threshold δ L. Thus, we suggest that mass functions of virialized halos should be defined with a fixed nonlinear density threshold, such as δ = 200, and fits to numerical simulations should use the exact exponentialcutoff of Eq. (32), with the appropriate linear density contrast δ L = F 1 (δ), rather than treating this as a free parameter. This would automatically ensure that the large mass tail has the right form (up to subleading terms such as power-law prefactors), as emphasized by the reasonable agreement with simulations of Eq. (34) at large masses. Moreover, it is best no to introduce unnecessary free parameters that become partly degenerate. Here we must note that Barkana (2004) had already noticed that, taking the spherical collapse at face value and defining halos by a nonlinear threshold δ (he chose δ = 18π 2 ), one should use the linear threshold δ L = F 1 (δ) for the Press-Schechter mass function. Unfortunately, noticing that the value of δ L givenbyfits to numerical simulations was even lower, he concluded that this was not sufficient to reconcile theoretical predictions with numerical results. As shown by Fig. 6, this is not the case, as the parameter δ L used in fitting formulae is partly degenerate with the exponents of the power-law prefactors, so that it is possible to match numerical simulations while satisfying the large-mass tail (32). Of course, the actual justification of the asymptote (32) is provided by the analysis of Sect. 3, which shows that below the upper bound δ + the spherical collapse is indeed relevant and asymptotically correct at large masses, as it corresponds to the saddle-point of the action (23). Note that the result (32) also implies that the mass function f (σ) is not exactly universal, since the mapping δ δ L = F 1 (δ) shows a (very) weak dependence on cosmological parameters, see Fig. 2. Using the exact tail (32) should also improve the robustness of the mass function with respect to changes of cosmological parameters and redshifts. Next, we compare in Fig. 7 the mass functions obtained at z = 0 for the density thresholds δ = 100 and δ = 300 with the numerical simulations from Tinker et al. (2008), who also considered these density thresholds. We show the usual Press- Schechter mass function (dashed line) and the results obtained from Eqs. (34)and(36) (solid lines), with now δ L = F 1 (100) 1.55 and δ L = F 1 (300) We can see that the large mass tail remains consistent with the simulations, but for the case δ = 100 it seems that the shape of the mass function at low and intermediate masses is modified and cannot be absorbed through the rescaling of δ L. It appears to follow Eq. (34) rather than the fit (36), but this is likely to be a mere coincidence. We should note that Fig. 5 shows that δ + < 300 for massive halos, so that shell-crossing should be taken into account and the tail (32)isno longer exact for δ = 300, although it should still provide a reasonable approximation, as checked in Fig. 7. Thus, even though Tinker et al. (2008) also studied higher density thresholds, we do not consider such cases here as the tail (32) no longer applies. 5. Halo density profile As for the mass function n(m), the analysis of Sect. 3 shows that the density profile of rare massive halos is given by the spherical saddle-point (25), see also Barkana (2004) andpradaetal. (2006). This holds for halos selected by some nonlinear density threshold δ in the limit of rare events, provided shell-crossing has not occurred beyond the associated radius r. In particular, this only applies to the outer part of the halo since in the inner part, at r r, shell-crossing must be taken into account. Then, as discussed in Sect. 3, a strong radial-orbit instability comes
9 P. Valageas: Mass functions and bias of dark matter halos 101 Fig. 7. The mass functions at z = 0 of halos defined by the nonlinear density contrasts δ = 100 (upper panel) andδ = 300 (lower panel). The points show the numerical simulations of Tinker et al. (2008). The dashed line is the usual Press-Schechter mass function, while the solid lines correspond to Eqs. (34) and(36), with now δ L = F 1 (100) 1.55 (upper panel) andδ L = F 1 (300) 1.61 (lower panel). into play and modifies the profile in this inner region, as deviations from spherical symmetry govern the dynamics and the virialization process (Valageas 2002b). We compare in Fig. 8 the nonlinear density profile obtained from Eq. (25) with fits to numerical simulations. We plot the overdensity within radius r,1+ δ r = ρ m (< r )/ρ m, as a function of radius r. This is again obtained from Eq. (25) with the mapping q 3 = (1+δ r )r 3 and δ r = F (δ Lq ). Since this only applies to the limit of rare events, we choose for each mass M aredshiftz such that σ(m, z) = 0.5. Thus, smaller masses are associated with higher redshifts. The results do not significantly depend on the precise value of the criterium used to define rare events, here σ(m, z) = 0.5. The points in Fig. 8 are the results obtained from a Navarro et al. profile (Navarro et al. 1997,NFW), ρ(r ρ s ) = (r /r s )(1 + r /r s ), (38) 2 or an Einasto profile (Einasto 1965), ln[ρ(r )/ρ 2 ] = 2 [ (r /r 2 ) α 1 ]. (39) α For the NFW profile, the characteristic radius r s is obtained from the concentration parameter, c(m, z) = r/r s,whereasin the previous sections r is the Eulerian radius where δ r = 200 (i.e. r = r 200 ). Then, we use in Fig. 8 the two fits obtained in Fig. 8. The density profile of rare massive halos, for several masses M. For each mass the redshift is such that σ(m, z) = 0.5, as the theoretical prediction from Eq. (25) (solid line) only applies to rare events. The second branch that appears at small radii is due to the shell-crossing, hence the theoretical prediction only holds to the right of this branch. The points are fits to numerical simulations, based on an NFW profile (Dolag et al. 2004; Duffy etal.2008) or an Einasto profile (Duffy et al. 2008). Note that we show the mean density within radius r, ρ m (< r )/ρ m, rather than the local density at radius r, so that the NFW and Einasto profiles are integrated once. Dolag et al. (2004)andDuffyetal.(2008) from numerical simulations, of the form c(m, z) = a(m/m 0 ) b (1 + z) c. For the Einasto profile, we use the fits obtained in Duffy etal.(2008) forthe concentration parameter, c(m, z) = r/r 2, and for the exponent α (written as a quadratic polynomial over ν = δ c /σ(m, z)asingao et al. 2008). Then, we integrate over r the profiles (38) (39), to obtain the overdensity profiles ρ m (<r )/ρ m shown in Fig. 8 (as ρ m (<r ) is the mean density within radius r and not the local density at radius r as in Eqs. (38) (39)). We can check in Fig. 8 that our results agree reasonably well with these fits to numerical simulations over the range where both predictions are valid. Note that our prediction has no free parameter, since it is given by the saddle-point profile (25). At large radii we recover the mean density of the Universe, while the numerical profiles (38) (39) go to zero, but this is an artefact of the forms (38) (39), that were designed to give a sharp boundary for the halos and are mostly used for the high-density regions. For a study of numerical simulations at outer radii, see Figs. 8 and 10 of Prada et al. (2006) which recover the mean density of the Universe at large scales. At small radii, the second branch that makes a turn somewhat below r is due to shell-crossing that makes the function ρ m (<r ) bivaluate. Below the maximum shell-crossing radius, for each Eulerian radius r there are two Lagrangian radii q, a large one that corresponds to shells that are still falling in, and a smaller one that corresponds to shells that have already gone once through the center (note that in agreement with Fig. 4, shell-crossing appears at slightly smaller relative radii for smaller mass). Then, the
2. What are the largest objects that could have formed so far? 3. How do the cosmological parameters influence structure formation?
Einführung in die beobachtungsorientierte Kosmologie I / Introduction to observational Cosmology I LMU WS 2009/10 Rene Fassbender, MPE Tel: 30000-3319, rfassben@mpe.mpg.de 1. Cosmological Principles, Newtonian
More informationPhysics 463, Spring 07. Formation and Evolution of Structure: Growth of Inhomogenieties & the Linear Power Spectrum
Physics 463, Spring 07 Lecture 3 Formation and Evolution of Structure: Growth of Inhomogenieties & the Linear Power Spectrum last time: how fluctuations are generated and how the smooth Universe grows
More informationDynamics of gravitational clustering
A&A 382, 412 430 (2002) DOI: 10.1051/0004-6361:20011663 c ESO 2002 Astronomy & Astrophysics Dynamics of gravitational clustering II. Steepest-descent method for the quasi-linear regime P. Valageas Service
More informationAST4320: LECTURE 10 M. DIJKSTRA
AST4320: LECTURE 10 M. DIJKSTRA 1. The Mass Power Spectrum P (k) 1.1. Introduction: the Power Spectrum & Transfer Function. The power spectrum P (k) emerged in several of our previous lectures: It fully
More informationThe mass of a halo. M. White
A&A 367, 27 32 (2001) DOI: 10.1051/0004-6361:20000357 c ESO 2001 Astronomy & Astrophysics The mass of a halo M. White Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA e-mail: mwhite@cfa.harvard.edu
More informationLecture 6: Dark Matter Halos
Lecture 6: Dark Matter Halos Houjun Mo February 23, 2004 Hierarchical clustering, the matter in the universe condenses to form quasiequilibrium objects (dark halos) of increasing masses in the passage
More informationN-body Simulations. Initial conditions: What kind of Dark Matter? How much Dark Matter? Initial density fluctuations P(k) GRAVITY
N-body Simulations N-body Simulations N-body Simulations Initial conditions: What kind of Dark Matter? How much Dark Matter? Initial density fluctuations P(k) GRAVITY Final distribution of dark matter.
More informationThe Los Cabos Lectures
January 2009 The Los Cabos Lectures Dark Matter Halos: 2 Simon White Max Planck Institute for Astrophysics EPS statistics for the standard ΛCDM cosmology Millennium Simulation cosmology: Ωm = 0.25, ΩΛ
More informationarxiv:astro-ph/ v1 27 Nov 2000
A&A manuscript no. (will be inserted by hand later) Your thesaurus codes are: 02 (3.13.18) - methods: N-body simulations ASTRONOMY AND ASTROPHYSICS The mass of a halo Martin White arxiv:astro-ph/0011495v1
More informationAnalytic model for galaxy and dark matter clustering
Analytic model for galaxy and dark matter clustering Uroš Seljak Department of Physics, Jadwin Hall, Princeton University, Princeton, NJ 08544 (January 2000) We investigate an analytic model to compute
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Friday 8 June 2001 1.30 to 4.30 PAPER 41 PHYSICAL COSMOLOGY Answer any THREE questions. The questions carry equal weight. You may not start to read the questions printed on
More informationWeek 3: Sub-horizon perturbations
Week 3: Sub-horizon perturbations February 12, 2017 1 Brief Overview Until now we have considered the evolution of a Universe that is homogeneous. Our Universe is observed to be quite homogeneous on large
More informationPreliminaries. Growth of Structure. Today s measured power spectrum, P(k) Simple 1-D example of today s P(k) Growth in roughness: δρ/ρ. !(r) =!!
Growth of Structure Notes based on Teaching Company lectures, and associated undergraduate text with some additional material added. For a more detailed discussion, see the article by Peacock taken from
More informationarxiv:astro-ph/ v1 26 Jul 2002
Non linear predictions from linear theories in models with Dark Energy R. Mainini, A.V. Macciò & S.A. Bonometto arxiv:astro-ph/0207581v1 26 Jul 2002 Physics Department G. Occhialini, Università degli Studi
More informationPONT 2011 Avignon - 20 Apr 2011 Non-Gaussianities in Halo Clustering properties. Andrea De Simone
PONT 2011 Avignon - 20 Apr 2011 Non-Gaussianities in Halo Clustering properties Andrea De Simone Based on: ADS, M. Maggiore, A. Riotto - 1102.0046 ADS, M. Maggiore, A. Riotto - MNRAS (1007.1903) OUTLINE
More informationHalo models of large scale structure
Halo models of large scale structure Asantha Cooray Theoretical Astrophysics, California Institute of Technology, Pasadena CA 91125. E-mail: asante@caltech.edu arxiv:astro-ph/0206508v1 28 Jun 2002 Ravi
More informationAST Cosmology and extragalactic astronomy Lecture 7
AST4320 - Cosmology and extragalactic astronomy Lecture 7 Press-Schechter Formalism: applications, origin of the `fudge factor 2, modifications Press-Schechter (PS) Formalism: Preface... PS assumes 1.
More informationGalaxies 626. Lecture 3: From the CMBR to the first star
Galaxies 626 Lecture 3: From the CMBR to the first star Galaxies 626 Firstly, some very brief cosmology for background and notation: Summary: Foundations of Cosmology 1. Universe is homogenous and isotropic
More informationNon-linear structure in the Universe Cosmology on the Beach
Non-linear structure in the Universe Cosmology on the Beach Puerto Vallarta January, 2011 Martin White UC Berkeley/LBNL (http://mwhite.berkeley.edu/talks) Strong non-linearity Martin White UCB/LBNL Limited
More informationPrinceton December 2009 The fine-scale structure of dark matter halos
Princeton December 2009 The fine-scale structure of dark matter halos Simon White Max Planck Institute for Astrophysics The dark matter structure of CDM halos A rich galaxy cluster halo Springel et al
More informationpseudo- evolution of halo mass and observable-mass relations Andrey Kravtsov The University of Chicago
pseudo- evolution of halo mass and observable-mass relations Andrey Kravtsov The University of Chicago Baseline model for cluster scaling relations Kaiser 1986, 1991 If mass is defined within a spherical
More informationDark Matter Halos. They will have a large impact on galaxies
Dark Matter Halos Dark matter halos form through gravitational instability. Density perturbations grow linearly until they reach a critical density, after which they turn around from the expansion of the
More informationarxiv:astro-ph/ v1 7 Jan 2000
Biased Estimates of Ω from Comparing Smoothed Predicted Velocity Fields to Unsmoothed Peculiar Velocity Measurements Andreas A. Berlind, Vijay K. Narayanan 1 and David H. Weinberg Department of Astronomy,
More informationChapter 9. Cosmic Structures. 9.1 Quantifying structures Introduction
Chapter 9 Cosmic Structures 9.1 Quantifying structures 9.1.1 Introduction We have seen before that there is a very specific prediction for the power spectrum of density fluctuations in the Universe, characterised
More informationHalo concentration and the dark matter power spectrum
Mon. Not. R. Astron. Soc. 340, 1199 1204 (2003) Halo concentration and the dark matter power spectrum Kevin M. Huffenberger and Uroš Seljak Department of Physics, Jadwin Hall, Princeton University, Princeton,
More informationIntroduction to Cosmology
Introduction to Cosmology João G. Rosa joao.rosa@ua.pt http://gravitation.web.ua.pt/cosmo LECTURE 2 - Newtonian cosmology I As a first approach to the Hot Big Bang model, in this lecture we will consider
More informationConstruction of the one-point PDF of the local aperture mass in weak lensing maps
Astron. Astrophys. 364, 1 16 () ASTRONOMY AND ASTROPHYSICS Construction of the one-point PDF of the local aperture mass in weak lensing maps F. Bernardeau and P. Valageas Service de Physique Théorique,
More informationGalaxy Formation! Lecture Seven: Galaxy Formation! Cosmic History. Big Bang! time! present! ...fluctuations to galaxies!
Galaxy Formation Lecture Seven: Why is the universe populated by galaxies, rather than a uniform sea of stars? Galaxy Formation...fluctuations to galaxies Why are most stars in galaxies with luminosities
More information4. Structure of Dark Matter halos. Hence the halo mass, virial radius, and virial velocity are related by
13-4-12see http://www.strw.leidenuniv.nl/ franx/college/galaxies12 12-c04-1 13-4-12see http://www.strw.leidenuniv.nl/ franx/college/galaxies12 12-c04-2 4. Structure of Dark Matter halos Obviously, we cannot
More informationCluster abundances, clustering and cosmology. Ravi K Sheth (ICTP/UPenn)
Cluster abundances, clustering and cosmology Ravi K Sheth (ICTP/UPenn) Motivation: A biased view of dark matters Observational probes Gravitational Instability The spherical (and tri-axial) collapse models
More informationCosmological Structure Formation
Cosmological Structure Formation Beyond linear theory Lagrangian PT: Zeldovich approximation (Eulerian PT: Spherical collapse) Redshift space distortions Recall linear theory: When radiation dominated
More informationSimulating non-linear structure formation in dark energy cosmologies
Simulating non-linear structure formation in dark energy cosmologies Volker Springel Distribution of WIMPS in the Galaxy Early Dark Energy Models (Margherita Grossi) Coupled Dark Energy (Marco Baldi) Fifth
More information4. Structure of Dark Matter halos. Hence the halo mass, virial radius, and virial velocity are related by
6-4-10see http://www.strw.leidenuniv.nl/ franx/college/galaxies10 10-c04-1 6-4-10see http://www.strw.leidenuniv.nl/ franx/college/galaxies10 10-c04-2 4. Structure of Dark Matter halos Obviously, we cannot
More informationarxiv:astro-ph/ v2 9 Apr 2002
Astronomy & Astrophysics manuscript no. (will be inserted by hand later) The phase-diagram of cosmological baryons P. Valageas 1, R. Schaeffer 1, and Joseph Silk 2 arxiv:astro-ph/0112273v2 9 Apr 2002 1
More informationClusters: Context and Background
Clusters: Context and Background We re about to embark on a subject rather different from what we ve treated before, so it is useful to step back and think again about what we want to accomplish in this
More informationA look into the Dark Universe
Zentrum für Astronomie, Institut für Theoretische Astrophysik Universität Heidelberg Dark Energy dark energy: assumed to drive accelerated cosmic expansion scalar field having negative pressure? early
More informationThe theoretical view of high-z Clusters. Nelson Padilla, PUC, Chile Pucón, November 2009
The theoretical view of high-z Clusters Nelson Padilla, PUC, Chile Pucón, November 2009 The Plan: 1) To see what the observations are telling us using models that agree with the cosmology, and with other
More informationThe fine-scale structure of dark matter halos
COSMO11, Porto, August 2011 The fine-scale structure of dark matter halos Simon White Max-Planck-Institute for Astrophysics COSMO11, Porto, August 2011 Mark Vogelsberger The fine-scale structure of dark
More informationD. f(r) gravity. φ = 1 + f R (R). (48)
5 D. f(r) gravity f(r) gravity is the first modified gravity model proposed as an alternative explanation for the accelerated expansion of the Universe [9]. We write the gravitational action as S = d 4
More informationStructure formation. Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München
Structure formation Yvonne Y. Y. Wong Max-Planck-Institut für Physik, München Structure formation... Random density fluctuations, grow via gravitational instability galaxies, clusters, etc. Initial perturbations
More informationCONSTRAINTS ON THE EFFECTS OF LOCALLY-BIASED GALAXY FORMATION
CONSTRAINTS ON THE EFFECTS OF LOCALLY-BIASED GALAXY FORMATION Robert J. Scherrer NASA/Fermilab Astrophysics Center, Fermi National Accelerator Laboratory, P.O. Box 500, Batavia, IL 60510, and Department
More informationComments on the size of the simulation box in cosmological N-body simulations
Mon. Not. R. Astron. Soc. 358, 1076 1082 (2005) doi:10.1111/j.1365-2966.2005.08858.x Comments on the size of the simulation box in cosmological N-body simulations J. S. Bagla and Suryadeep Ray Harish-Chandra
More informationPhysical Cosmology 12/5/2017
Physical Cosmology 12/5/2017 Alessandro Melchiorri alessandro.melchiorri@roma1.infn.it slides can be found here: oberon.roma1.infn.it/alessandro/cosmo2017 Structure Formation Until now we have assumed
More informationformation of the cosmic large-scale structure
formation of the cosmic large-scale structure Heraeus summer school on cosmology, Heidelberg 2013 Centre for Astronomy Fakultät für Physik und Astronomie, Universität Heidelberg August 23, 2013 outline
More information4 Evolution of density perturbations
Spring term 2014: Dark Matter lecture 3/9 Torsten Bringmann (torsten.bringmann@fys.uio.no) reading: Weinberg, chapters 5-8 4 Evolution of density perturbations 4.1 Statistical description The cosmological
More informationThe bias and mass function of dark matter haloes in non-markovian extension of the excursion set theory
Mon. Not. R. Astron. Soc. 4, 644 65 (0) doi:0./j.365-966.00.787.x The bias and mass function of dark matter haloes in non-markovian extension of the excursion set theory Chung-Pei Ma, Michele Maggiore,
More informationcluster scaling relations and mass definitions
cluster scaling relations and mass definitions Andrey Kravtsov Department of Astronomy & Astrophysics Kavli Institute for Cosmological Physics The University of Chicago Abell 85 SDSS Abell 85 SDSS/ Abell
More informationTheory of galaxy formation
Theory of galaxy formation Bibliography: Galaxy Formation and Evolution (Mo, van den Bosch, White 2011) Lectures given by Frank van den Bosch in Yale http://www.astro.yale.edu/vdbosch/teaching.html Theory
More informationYou may not start to read the questions printed on the subsequent pages until instructed to do so by the Invigilator.
MATHEMATICAL TRIPOS Part III Thursday 3 June, 2004 9 to 12 PAPER 67 PHYSICAL COSMOLOGY Attempt THREE questions. There are four questions in total. The questions carry equal weight. You may not start to
More informationHalo Model. Vinicius Miranda. Cosmology Course - Astro 321, Department of Astronomy and Astrophysics University of Chicago
Introduction program Thanks Department of Astronomy and Astrophysics University of Chicago Cosmology Course - Astro 321, 2011 Introduction program Thanks Outline 1 Introduction 2 Linear Theory Press -
More informationThe imprint of the initial conditions on large-scale structure
Stars, Planets and Galaxies 2018 Harnack House, Berlin The imprint of the initial conditions on large-scale structure Simon White Max Planck Institute for Astrophysics The Planck map of TCMB the initial
More informationA5682: Introduction to Cosmology Course Notes. 11. CMB Anisotropy
Reading: Chapter 8, sections 8.4 and 8.5 11. CMB Anisotropy Gravitational instability and structure formation Today s universe shows structure on scales from individual galaxies to galaxy groups and clusters
More informationThe Current Status of Too Big To Fail problem! based on Warm Dark Matter cosmology
The Current Status of Too Big To Fail problem! based on Warm Dark Matter cosmology 172th Astronomical Seminar Dec.3 2013 Chiba Lab.M2 Yusuke Komuro Key Word s Too Big To Fail TBTF Cold Dark Matter CDM
More informationPhysical Cosmology 18/5/2017
Physical Cosmology 18/5/2017 Alessandro Melchiorri alessandro.melchiorri@roma1.infn.it slides can be found here: oberon.roma1.infn.it/alessandro/cosmo2017 Summary If we consider perturbations in a pressureless
More informationGaia Revue des Exigences préliminaires 1
Gaia Revue des Exigences préliminaires 1 Global top questions 1. Which stars form and have been formed where? - Star formation history of the inner disk - Location and number of spiral arms - Extent of
More informationStructure of Dark Matter Halos
Structure of Dark Matter Halos Dark matter halos profiles: DM only: NFW vs. Einasto Halo concentration: evolution with time Dark matter halos profiles: Effects of baryons Adiabatic contraction Cusps and
More informationNon-perturbative probability distribution function for cosmological counts in cells
CERN-TH-018-4 INR-TH-018-07 Non-perturbative probability distribution function for cosmological counts in cells arxiv:1811.07913v1 [astro-ph.co] 19 Nov 018 Mikhail M. Ivanov 1a,b,c Alexander A. Kaurov
More informationFormation time-distribution of dark matter haloes: theories versus N-body simulations
Mon. Not. R. Astron. Soc. 344, 1327 1333 (2003) Formation time-distribution of dark matter haloes: theories versus N-body simulations W. P. Lin, 1,2 Y. P. Jing 1,2 and Lihwai Lin 3 1 Shanghai Astronomical
More information8.1 Structure Formation: Introduction and the Growth of Density Perturbations
8.1 Structure Formation: Introduction and the Growth of Density Perturbations 1 Structure Formation and Evolution From this (Δρ/ρ ~ 10-6 ) to this (Δρ/ρ ~ 10 +2 ) to this (Δρ/ρ ~ 10 +6 ) 2 Origin of Structure
More informationStructure formation in the concordance cosmology
Structure formation in the Universe, Chamonix, May 2007 Structure formation in the concordance cosmology Simon White Max Planck Institute for Astrophysics WMAP3 team WMAP3 team WMAP3 team WMAP3 team In
More informationR. K. Sheth, H. J. Mo & G. Tormen collapses at time z if the initial overdensity within it exceeds a critical value, sc(z). This critical value depend
Mon. Not. R. Astron. Soc. 000, 000{000 (0000) Printed July 1999 (MN LATEX style le v1.4) Ellipsoidal collapse and an improved model for the number and spatial distribution of dark matter haloes Ravi K.
More informationA A + B. ra + A + 1. We now want to solve the Einstein equations in the following cases:
Lecture 29: Cosmology Cosmology Reading: Weinberg, Ch A metric tensor appropriate to infalling matter In general (see, eg, Weinberg, Ch ) we may write a spherically symmetric, time-dependent metric in
More informationarxiv:astro-ph/ v1 3 Oct 2001
AiBA 2001: High-z Clusters, issing Baryons, and CB Polarization ASP Conference Series, Vol. 999, 2002 L-W Chen, C-P a, K-W Ng and U-L Pen, eds Clustering of Dark atter Halos on the Light-cone arxiv:astro-ph/0110073v1
More informationThe assembly of matter in galaxy clusters
Mon. Not. R. Astron. Soc. 297, 648 656 (1998) The assembly of matter in galaxy clusters Giuseppe Tormen 1;2 * 1 Max Planck Institute für Astrophysik, Karl-Schwarzschild-Strasse 1, 85740 Garching bei München,
More informationAn excursion set model of the cosmic web: The abundance of sheets, filaments and halos
SLAC-PUB-11626 astro-ph/0511365 An excursion set model of the cosmic web: The abundance of sheets, filaments and halos Jiajian Shen Department of Astronomy and Astrophysics, The Pennsylvania State University,
More information5.1 Circular Velocities and Rotation Curves
Chapter 5 otation Curves 5.1 Circular Velocities and otation Curves The circular velocity v circ is the velocity that a star in a galaxy must have to maintain a circular orbit at a specified distance from
More informationCluster Formation in Lagrangian Perturbation Theory. Lavinia Heisenberg
Cluster Formation in Lagrangian Perturbation Theory Zentrum fu r Astronomie (ITA/ZAH) Universita t Heidelberg 18.12.2009, University of Canterbury, Christchurch outline 1 Introduction 2 LPT 3 ELL 4 PINOCCHIO
More informationThe evidence of absence: galaxy voids in the excursion set formalism
Mon. Not. R. Astron. Soc. 366, 467 479 (2006) doi:10.1111/j.1365-2966.2005.09862.x The evidence of absence: galaxy voids in the excursion set formalism Steven R. Furlanetto 1 and Tsvi Piran 1,2 1 Division
More informationOrigin of Structure Formation of Structure. Projected slice of 200,000 galaxies, with thickness of a few degrees.
Origin of Structure Formation of Structure Projected slice of 200,000 galaxies, with thickness of a few degrees. Redshift Surveys Modern survey: Sloan Digital Sky Survey, probes out to nearly 1000 Mpc.
More informationDependence of the inner dark matter profile on the halo mass
Mon. Not. R. Astron. Soc. 344, 1237 1249 (2003) Dependence of the inner dark matter profile on the halo mass Massimo Ricotti Institute of Astronomy, Madingley Road, Cambridge CB3 0HA Accepted 2003 June
More informationarxiv:astro-ph/ v1 20 Sep 2006
Formation of Neutrino Stars from Cosmological Background Neutrinos M. H. Chan, M.-C. Chu Department of Physics, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong, China arxiv:astro-ph/0609564v1
More informationhalo merger histories
Physics 463, Spring 07 Lecture 5 The Growth and Structure of Dark Mater Halos z=7! z=3! z=1! expansion scalefactor z=0! ~milky way M>3e12M! /h c vir ~13 ~virgo cluster M>3e14M! /h, c vir ~6 halo merger
More informationLinear Theory and perturbations Growth
Linear Theory and perturbations Growth The Universe is not homogeneous on small scales. We want to study how seed perturbations (like the ones we see in the Cosmic Microwave Background) evolve in an expanding
More informationFuture evolution of bound supercluster in an accelerating Universe
Future evolution of bound supercluster in an accelerating Universe P.A. Araya-Melo, A. Reisenegger, A. Meza, R. van de Weygaert, R. Dünner and H. Quintana August 27, 2010 Mon. Not. R. Astron. Soc. 399,97-120
More informationProbing Dark Matter Halos with Satellite Kinematics & Weak Lensing
Probing Dark Matter Halos with & Weak Lensing Frank C. van den Bosch (MPIA) Collaborators: Surhud More, Marcello Cacciato UMass, August 2008 Probing Dark Matter Halos - p. 1/35 Galaxy Formation in a Nutshell
More informationClusters: Context and Background
Clusters: Context and Background We reabouttoembarkon asubjectratherdifferentfrom what we vetreatedbefore, soit is useful to step back and think again about what we want to accomplish in this course. We
More informationVoid abundance predic.on and spherical evolu.on model.!!!!! Ixandra Achitouv! Universitäts- Sternwarte München!
Void abundance predic.on and spherical evolu.on model!!!!! Ixandra Achitouv! Universitäts- Sternwarte München! Voids hierarchy - Sheth & Van de Weijgaert (04) Voids Hierarchies: Number of voids dn/dr at
More informationPAPER 71 COSMOLOGY. Attempt THREE questions There are seven questions in total The questions carry equal weight
MATHEMATICAL TRIPOS Part III Friday 31 May 00 9 to 1 PAPER 71 COSMOLOGY Attempt THREE questions There are seven questions in total The questions carry equal weight You may make free use of the information
More informationDwarf Galaxies as Cosmological Probes
Dwarf Galaxies as Cosmological Probes Julio F. Navarro The Ursa Minor dwarf spheroidal First Light First Light The Planck Satellite The Cosmological Paradigm The Clustering of Dark Matter The Millennium
More informationSEMIANALYTICAL DARK MATTER HALOS AND THE JEANS EQUATION
The Astrophysical Journal, 634:756 774, 2005 December 1 # 2005. The American Astronomical Society. All rights reserved. Printed in U.S.A. SEMIANALYTICAL DARK MATTER HALOS AND THE JEANS EQUATION Crystal
More informationGalaxy Formation Seminar 2: Cosmological Structure Formation as Initial Conditions for Galaxy Formation. Prof. Eric Gawiser
Galaxy Formation Seminar 2: Cosmological Structure Formation as Initial Conditions for Galaxy Formation Prof. Eric Gawiser Cosmic Microwave Background anisotropy and Large-scale structure Cosmic Microwave
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationBAO & RSD. Nikhil Padmanabhan Essential Cosmology for the Next Generation VII December 2017
BAO & RSD Nikhil Padmanabhan Essential Cosmology for the Next Generation VII December 2017 Overview Introduction Standard rulers, a spherical collapse picture of BAO, the Kaiser formula, measuring distance
More informationPAPER 73 PHYSICAL COSMOLOGY
MATHEMATICAL TRIPOS Part III Wednesday 4 June 2008 1.30 to 4.30 PAPER 73 PHYSICAL COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationarxiv: v1 [astro-ph.co] 4 Sep 2009
Spherical collapse model with and without curvature Seokcheon Lee 1,2 arxiv:99.826v1 [astro-ph.co] 4 Sep 29 1 Institute of Physics, Academia Sinica, Taipei, Taiwan 11529, R.O.C. 2 Leung Center for Cosmology
More informationThe Iguaçu Lectures. Nonlinear Structure Formation: The growth of galaxies and larger scale structures
April 2006 The Iguaçu Lectures Nonlinear Structure Formation: The growth of galaxies and larger scale structures Simon White Max Planck Institute for Astrophysics z = 0 Dark Matter ROT EVOL Cluster structure
More informationarxiv:astro-ph/ v1 22 Sep 2005
Mass Profiles and Shapes of Cosmological Structures G. Mamon, F. Combes, C. Deffayet, B. Fort (eds) EAS Publications Series, Vol.?, 2005 arxiv:astro-ph/0509665v1 22 Sep 2005 MONDIAN COSMOLOGICAL SIMULATIONS
More informationNewtonian Gravity and Cosmology
Chapter 30 Newtonian Gravity and Cosmology The Universe is mostly empty space, which might suggest that a Newtonian description of gravity (which is valid in the weak gravity limit) is adequate for describing
More informationEnergy and Matter in the Universe
Chapter 17 Energy and Matter in the Universe The history and fate of the Universe ultimately turn on how much matter, energy, and pressure it contains: 1. These components of the stress energy tensor all
More informationThe Caustic Technique An overview
The An overview Ana Laura Serra Torino, July 30, 2010 Why the mass of? Small scales assumption of dynamical equilibrium Mass distribution on intermediate scales (1 10 Mpc/h) Large scales small over densities
More informationTESTING GRAVITY WITH COSMOLOGY
21 IV. TESTING GRAVITY WITH COSMOLOGY We now turn to the different ways with which cosmological observations can constrain modified gravity models. We have already seen that Solar System tests provide
More informationNon-linear dynamics and mass function of cosmic structures
Astron. Astrophys. 25, 49 449 997) ASTRONOMY AND ASTROPHYSICS Non-linear dynamics and mass function of cosmic structures I. Analytical results Edouard Audit, Romain Teyssier 2, and Jean-Michel Alimi Laboratoire
More informationPrimordial Black Holes in Cosmology. Lectures 1 & 2 : What are PBHs? Do they exist? Massimo Ricotti (University of Maryland, USA)
Primordial Black Holes in Cosmology Lectures 1 & 2 : What are PBHs? Do they exist? Massimo Ricotti (University of Maryland, USA) Institute of Cosmos Sciences, University of Barcelona 23/10/2017 What are
More informationIntroduction to Inflation
Introduction to Inflation Miguel Campos MPI für Kernphysik & Heidelberg Universität September 23, 2014 Index (Brief) historic background The Cosmological Principle Big-bang puzzles Flatness Horizons Monopoles
More informationOrbital Motion in Schwarzschild Geometry
Physics 4 Lecture 29 Orbital Motion in Schwarzschild Geometry Lecture 29 Physics 4 Classical Mechanics II November 9th, 2007 We have seen, through the study of the weak field solutions of Einstein s equation
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 4 Stellar orbits and dark matter 1 Using Kepler s laws for stars orbiting the center of a galaxy We will now use Kepler s laws of gravitation on much larger scales. We will study
More informationObservational Cosmology
(C. Porciani / K. Basu) Lecture 7 Cosmology with galaxy clusters (Mass function, clusters surveys) Course website: http://www.astro.uni-bonn.de/~kbasu/astro845.html Outline of the two lecture Galaxy clusters
More informationNEWTONIAN COSMOLOGY. Figure 2.1: All observers see galaxies expanding with the same Hubble law. v A = H 0 r A (2.1)
M. Pettini: Introduction to Cosmology Lecture 2 NEWTONIAN COSMOLOGY The equations that describe the time evolution of an expanding universe which is homogeneous and isotropic can be deduced from Newtonian
More informationLearning from WIMPs. Manuel Drees. Bonn University. Learning from WIMPs p. 1/29
Learning from WIMPs Manuel Drees Bonn University Learning from WIMPs p. 1/29 Contents 1 Introduction Learning from WIMPs p. 2/29 Contents 1 Introduction 2 Learning about the early Universe Learning from
More informationOn an Analytical Framework for Voids: Their abundances, density profiles and local mass functions
Mon. Not. R. Astron. Soc., () Printed February 6 (MN LATEX style file v.) On an Analytical Framework for Voids: Their abundances, density profiles and local mass functions S. G. Patiri, J. Betancort-Rijo,,
More information