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 and Particle Astrophysics, National Taiwan University, Taipei, Taiwan 1617, R.O.C. Abstract We investigate a spherical collapse model with and without the spatial curvature. We obtain the exact solutions of dynamical quantities such as the ratio of the scale factor to its value at the turnaround epoch (x a a ta ) and the ratio of the overdensity radius to its value at the turnaround time (y R R ta ) with general cosmological parameters. The exact solutions of the overdensity at the turnaround epoch = ρ cluster ρ m for the different models are also obtained. Thus, we are able z=zta to obtain the nonlinear overdensity 1 + δ NL = ( x y ) at any epoch for the given model. The ) nonlinear overdensity at the virial epoch vir = ( xvir y vir is obtained by using the virial theorem and the energy conservation. We obtain that the nonlinear overdensity of the Einstein de Sitter (EdS) Universe EdS vir is 18π 2( 2 1 2π 4) + 147 instead of the well known value 18π 2 178. In the open Universe, perturbations are virialized earlier than in flat Universe and thus clusters are denser at the virial epoch. Also the critical density threshold from the linear theory δlin EdS at the virialized epoch is obtained as 2 (9π + 6)2 1.58 instead of 2 (12π)2 1.69. This value is same for the close and the open Universes. We also show that there is the minimum z ta value with the given cosmological parameters for each model. We find that the observed quantities at high redshifts are less sensitive between different models. Thus, the high redshift cluster (z vir.5) is not a good object to probe the curvature of the Universe. The low redshift cluster (z vir.7 i.e. z ta.6) shows the stronger model dependence feature. However, it might be still too small to be distinguished. Because we obtain the analytic forms of dynamical quantities, we are able to estimate the abundances of both virialized and non-virialized clusters at any epoch. Also the temperature and luminosity functions are able to be computed at any epoch. Thus, these analytic forms of x, y, and provide the accurate tools for probing the curvature of the Universe. Unfortunately, the current concordance model prefers the almost flat Universe and the above results might be restricted by the academic interests only. However, mathematical structure of the physical evolution equations of the curved space is identical with that of the dark energy with the equation of state ω de = 1. Thus, we are able to extend these analytic solutions to the general dark energy model and they will provide the useful tools for probing the properties of dark energy.
Background evolution equations of physical quantities in a FRW Universe with the matter are given by H 2 = (ȧ a ) 2 = 8πG ρ m k a 2 = 8πG ρ cr, (.1) ä a = 4πG ρ m, (.2) (ȧ ) ρ m + ρ m =, (.) a where a is the scale factor, ρ m is the energy densities of the matter, ρ cr is the critical energy density, and k is chosen to be +1,, or 1 for spaces of constant positive, zero, or negative spatial curvature. In terms of the ratio of the matter density to the critical density Ω m, the above Friedmann equation (.1) becomes which is valid for all times. k H 2 a 2 Ω k = Ω m 1, (.4) We consider a spherical perturbation in the matter density. ρ cluster is the matter density within the spherical overdensity radius R. The flatness condition is not held because of the perturbation in the matter (The curvature is not a constant inside the overdensity patch). Thus, we have the another equations governing the dynamics of the spherical perturbation [1] R R = 4πG ρ cluster, (.5) ( Ṙ ) ρ cluster + ρ cluster =, (.6) R where ρ cluster is the energy density of the clustering matter. The radius of the overdensity R evolves slower than the scale factor a and reaches its maximum size R ta at the turnaround epoch z ta and then the system begins to collapse. Cosmological parameters and the curvature of the Universe can be constrained from the growth of large scale structure and the abundance of rich clusters of galaxies. There have been numerous works related to this [2,, 4, 5, 6, 7, 8, 9, 1, 11, 12]. Most of them reach to the similar conclusions based on the conventional approximate solutions of the background scale factor and of the overdensity radius. It is natural to expect that the correct values for the virial radius and the nonlinear overdensity obtained from the exact solutions might be different from those obtained from the conventional approximate solutions. In order to investigate this, we need to obtain the exact solutions of physical quantities. Now we adopt the notations in Ref. [11] to investigate the evolutions of a and R x = a a ta, (.7) y = R R ta, (.8) 1
H tta.5..25.2.15.8.9 1. 1.1 1.2 m Τta.68.67.66.65.64.8.9 1. 1.1 1.2 Figure 1: t ta and τ ta for the different values of z ta. a) H t ta verses Ω m for the different values of z ta =.6,1.2, and 2. (from top to bottom). b) τ ta verses Ω m for the same values of z ta as in the left panel. m where a ta is the scale factor of the background evolution at z ta. Then the equations (.1) and (.5) are rewritten as dx = x dτ 1 Q 1 ta, (.9) d 2 y dτ 2 = 1 2 y 2, (.1) where dτ = H(x ta ) Ω m (x ta )dt H ta Ωmta dt, = Ωm Ω k zta Ωmta Ω kta = Ωmta Ω = Ω m mta 1 Ω 1(1 + z ta), m = ρ cluster ρ m zta, and x ta x(z ta ) = 1 from the equation (.7). Ω m and Ω k represent the present value of energy density contrasts of the matter and the curvature term, respectively. Equations (.9) and (.1) can be solved analytically and we will obtain them. The analytic solution of Eq. (.9) is given by x dx = x 1 Q 1 ta τ dτ 2 [ 1 2F x 2, 2, 5 2, x ] = τ, (.11) where F is the hypergeometric function and we use the boundary condition x = when τ = (see Appendix for details). From this equation, the exact turnaround time τ ta is given by τ ta = 2 [ 1 F 2, 2, 5 2, Ω m 1 Ω (1 + z ta ) 1] = H ta Ωmta t ta = H Ω m(1 + z ta ) 2 tta, (.12) m where we use the fact that x ta = 1, the relations = Ω m Ω m 1(1 + z ta) and τ = H ta Ωmta t. This exact analytic form of the turnaround time will be used to investigate the other quantities. As expected, τ ta (t ta ) depends on Ω m (i.e. Ω k ) and z ta as given in the Eq. (.12). We show these properties of τ ta (t ta ) in Fig. 1. In the left panel of Fig. 1, we show t ta dependence 2
6. 5.8 Ζ 5.6 5.4 5.2.8.9 1. 1.1 1.2 Figure 2: verses Ω m when we choose z ta =.6, 1.2, and 2.. m on the values of Ω m for the different choice of z ta models. The solid, dashed, and dotdashed lines (from top to bottom) correspond to z ta =.6,1.2, and 2., respectively. Eq. (.11) is the evolution of the background scale factor a and we can interpret it as the age of the Universe is a decreasing function of Ω m. Larger Ω m implies faster deceleration, which corresponds a more rapidly expanding universe early on. Also larger z ta means the earlier formation of the structure and thus gives the smaller t ta. We also show the Ω m dependence of τ ta for the values of z ta in the right panel of Fig. 1. Because τ ta = H ta Ωmta t ta, τ ta becomes larger for the larger values of Ω m. The evolution of y given in Eq. (.1) described as the exact analytic solution (see Appendix) ArcSin[ y] y(1 y) = τ,when τ τ ta, (.1) y(1 y) ArcSin[ y] + π 2 = (τ τ ta ),when τ τ ta, (.14) where τ and τ ta are given in Eqs. (.11) and (.12). can be obtained from this analytic solutions Eq. (.1) (or equally from Eq. (.14)) by using the fact that y ta = 1 = ( π 2τ ta ) 2 = ( ) 2 ( π F 4 [ 1 2, 2, 5 2, (Ω m 1) Ω m (1 + z ta ) 1]) 2, (.15) [ ] where we use Eq. (.12). When Ω m = 1, F 1 2, 2, 5 2, = 1 and thus = ( π 4 )2. This factor ( π 4 )2 is the well known value of for the Einstein de Sitter (EdS) Universe (Ω m = 1) [1, 1]. The general value of for open or closed Universe is given by Eq. (.15). We show the behavior of in Fig. 2. As shown in the equation (.15), is inversely proportional to τ 2 ta. Thus, decreases as Ω m increases. Because is the ratio of ρ cluster to ρ m at z = z ta, it means that the smaller overdensity
2. 1.8 t t ta,eds 1.6 1.4 1.2..2.4.6.8 1. Figure : t normalized to the turnaround time for the EdS Universe verses y for close, EdS, and open Universes (from top to bottom). y takes longer time to turnaround and collapse. In Fig. 2, the solid, dashed, and dotdashed lines correspond to z ta =.6, 1.2, and 2., respectively. After we obtain the value of, we are able to obtain the values of x and y at any τ without any ambiguity because the analytic solutions of x, y, and given in Eqs. (.11), (.1), (.14), and (.15) are exact. In order to better understand the above results, it is useful to investigate the virialized times for the different models. From the equation (.14), we are able to obtain the collapsing time of each model normalized to the turnaround time for the EdS Universe t = 2 ( y(1 y) ArcSin[ ) 1 y] + π F[ t ta,eds π 2, 2, 5 2,Qta 1], (.16) where t ta,eds is the turnaround time for the EdS Universe. We show this in Fig.. Perturbations is virialized earlier in the EdS Universe than in the close Universe, and even earlier in the open Universe. The dashed, solid, and dotdashed lines represent t normalized to the turnaround time of the EdS Universe, t ta,eds of close, EdS, and open Universes, respectively. The vertical dotted line depicts the ratio of the virialized radius to the turnaround one y vir = 1 2. If only the matter virializes, then from the virial theorem and the energy conservation we have ( U cluster (z ta ) = U cluster + R ) U cluster = 1 2 R zvir 2 U cluster(z vir ), (.17) where U cluster = GM 2 5 R is the potential energy associated with the spherical mass overdensity [1, 4, 5, 14]. From the above equation we are able to obtain y vir = R vir R ta for any model y vir = 1 2. (.18) 4
After replacing y vir = 1 2 into Eq. (.14), we obtain ( τ vir = 2 + 1 ) τ ta = (1 + 2 ) [ 1 F π π 2, 2, 5 2, 1 ]. (.19) x vir is obtained from Eq. (.11) with Eq. (.19) x [ 2 1 vir F 2, 2, 5 2, x ] ( vir = 2 π) + 1 [ 1 F 2, 2, 5 2, 1 ]. (.2) Even though we obtain the analytic expression of x vir in Eq. (.2), generally this equation can be solved non algebraic way. Thus, it might be useful to have the approximate analytic form of x vir for the wide ranges of cosmological parameters (Ω m and z ta) if one analyze the form by interpolating x vir values. We use the exact form in Eq. (.2) though. Before we move further for the nonlinear overdensity, we investigate the common mistake for vir in the EdS Universe. Eqs. (.11), (.14), and (.15) can be simplified in EdS Universe because Q 1 ta = when Ω m = 1. Then the above equations become 2 x 2 = τ, (.21) y(1 y) ArcSin[ y] + π 2 = (τ 2 ) = π ( ) x 2 1, (.22) 2 = ( π 4 )2, (.2) where we use F[ 1 2ω de, 1 2,1 1 2ω de,] = 1 and τ ta = 2. The commonly used assumption to obtain the nonlinear overdensity c in the EdS Universe is that τ c = 2τ ta = 4 which is the collapsing time for y c = even though one uses y vir = 1 2 to obtain c. By using this assumption, x c = 2 2 from the equation (.21). Thus, one obtains EdS c = ( ) xc = y vir ( ) 2 ( 2 ) π 2 4 2 1 = 18π 2 178. (.24) However, we know the exact relation between x and y and thus we do not need to use the above assumption. If we insert y vir = 1 2 in the equation (.22), then we obtain the correct τ vir = 1 + 2 π < τ c. With this correct value of τ vir, we obtain the correct x vir = ( 2 + 1 π )2 by using the Eq. (.21). The correct value of the nonlinear overdensity vir for the EdS universe is EdS vir = ( ) xvir = y vir ( ) 2 ( ( π 4 2 + 1 π 2 1 ) 2 ) ( = 18π 2 4 + 1 ) 2 147. (.25) 2π Thus, the minimum overdensity for the flat Universe is about 147 instead of 178. We show this in tables 1. As z ta increases c approaches to 147 instead of 178. Also for the closed universe vir can be smaller than this value. From the equation (.7) we are able to obtain the virial epoch from the given z ta z vir = 1 + z ta x vir 1 (.26) 5
z ta =.6 z ta = 1.2 z ta = 2. Ω m x vir vir z vir x vir vir z vir x vir vir z vir.8 6.7 1.511 168.6 5.9 1.55 162.46 5.8 1.51 158 1..95 5.66 1.494 151.7 5.6 1.49 15.47 5.61 1.492 149 1.1 1. 5.55 1.49 147.7 5.55 1.49 147.48 5.55 1.49 147 1.1 1.5 5.45 1.485 14.8 5.48 1.487 144.48 5.5 1.487 145 1.2 1.2 5.2 1.474 1.9 5. 1.478 17.49 5.7 1.481 14 1. Table 1:, x vir, vir, and z vir with the given values of cosmological parameters for the open, flat, and closed Universes. Perturbations reach turnaround and virialization earlier in the flat Universe than in the close Universe, and even earlier in the open Universe. Thus, clusters are denser at virial epoch in flat Universe than in the close Universe and even denser in the open Universe. The ratio of cluster to background density at the virialized epoch z vir for open and closed Universes becomes vir = ρ cluster ρ m ( xvir = zvir y vir ) ( ) 2 = 18π 2 4 + 1 2π F[ 1 2, 2, 2, x vir ] 2 147 F[ 1 2, 2, 2, x vir ] 2, (.27) where we use Eqs. (.), (.6), and (.15). We also use the fact that y vir = 1 2 independent of. From the equation (.27), we are able to investigate the linear perturbation at early epoch τ 1 + δ lin (see Appendix for details). τ 1 + δ NL 1 + δ lin = 1 + ( ) 2 5 2 τ = 1 + )2 t Ωmta 5(. (.28) t ta It is equal to the famous result for EdS Universe ( ) = π 4 and Ω mta = 1), δ lin = 5( 2 π t 4 t ta [1]. Again, there is a common mistake for the value of the critical density threshold δ lin (z vir ). In EdS Universe, τ c = 4 and one obtains δ lin(z c ) from Eq. (.28) δ lin (z c ) = ( π ) 2 ( 5 4 2 4 )2 = 2 (12π)2 1.69. (.29) ( However, we show that the correct virialized epoch gives τ vir = 2 π) + 1 τ ta from Eq. (.19). Thus, the correct δ lin (z vir ) value is δ lin (z vir ) = 5 ( ) 2 ( ( 4 + 9π 8 )2 ) 1 = 2 (6 + 9π)2 1.58, (.) where we use = π 2τ ta. The above result given in Eq. (.) is true for with and without curvature and thus it is valid for the close, flat, and open Universes. After we obtain the critical 6
density threshold δ lin (z vir ) we are able to obtain δ lin at any epoch by using the relation δ lin (z) = D g(z) D g (z vir ) δ lin(z vir ) (.1) where D g is the linear growth factor. There is the exact analytic form of D g for the dark energy with the equation of state ω de = 1 [15]. Mathematically, this form of D g is identical with that of curved space and thus we can adopt this form of D g in these models. From the analytic forms of dynamical quantities x, y, and, it is straightforward to estimate the abundances of both virialized and non-virialized clusters at any epoch. Also the temperature and luminosity functions are able to be computed at any epoch [16]. Thus, these analytic forms provide an accurate tool for probing the effect of the curvature on the clustering. As we mentioned, the mathematical structure of the physical evolution equations of the curved space is identical with that of the dark energy with the equation of state ω de = 1. Thus, these analytic solutions give the guideline for the extension of them to those of the general dark energy model and they will provide the useful tools for probing the properties of dark energy [17]. A Appendix First, we derive the exact solution of τ as a function of x given in Eq. (.11). After we replace the variables Z = x and T = x x, the integral in Eq. (.11) becomes [18] x dx = x 2 x 1 Q 1 ta 1 T 2(1 1 ZT) 1 2dT = 2 [ 1 2F x 2, 2, 5 2, x ], (A.1) where F is the hypergeometric function and we use the gamma function relation Γ[1 + b] = bγ[b]. We also derive the exact analytic solution of τ as a function of y given in Eq. (.1). This equation is solved as [19] dy 2 dτ 2 = 1 2 y 2 f(y) y dy c 1 + 2 y f(y )dy = c 2 ± τ. (A.2) We use 2 y f(y )dy = y 1 to obtain y dy c1 + y 1 = c 2 ± τ. (A.) After replacing Z = c 1 y and T = y y, the LHS of Eq. (A.) becomes y y 2 (1 y ) dy = y 1 T 2(1 1 c 1 2dT yt) 1 = 2 1 y 2F [ 1 2, 2, 5 2, c 1 y ]. (A.4) 7
We use the relation c 1 yf[ 1 2, 2, 5 2, c 1 y] = 2 F[1 2, 1 2, 2, c 1 y] 2 (1 c 1 y)f[ 1 2, 2, 2, c 1 y], F[ 1 2, 1 2, 2,y2 ] = y 1 ArcSin[y], and F[a,b,b,y] = (1 y) a to obtain ( y c 1 y c 1 [ c 1 y ArcSin c1 y ( [ c 1 y ArcSin c1 y] ] 1 c 1y 1 c 1y ) ) = c 2l + τ when τ τ ta, (A.5) = c 2g τ when τ τ ta, (A.6) where ArcSin represents arcsine function. We choose the positive (negative) sign in τ to be consistent with the fact that R (i.e. y) increases (decreases) as τ increases before (after) it reaches to the maximum radius R ta ( y ta = 1) at τ ta. Integral constants c 1, c 2l, and c 2g are obtained from the boundary conditions dy dτ ta =, y(τ = ) =, and y ta = 1 dy dτ = ± ta c 1 + y ta = c 1 =, (A.7) = c 2l (A.8) π 2 = c 2g τ ta, (A.9) where we use the fact that ArcSin[1] = π 2. τ ta is obtained from Eq. (.12) τ ta = 2 F [ 2, 1 2, 5 2, 1 ]. (A.1) From Eqs. (A.8), (A.9), and (A.1), we obtain the exact analytic solution of y(τ) ArcSin[ y] y(1 y) = τ,when τ τ ta. (A.11) y(1 y) ArcSin[ y] + π 2 = (τ τ ta ),when τ τ ta. (A.12) [ y ] If we use the relationships ArcSin y(1 y) + π 4 [ ] [ ] + π 4 = 1 2 ArcTan 1 2y 1 2 y(1 y) = 2 y 2F 1 2, 2, 5 2,y, then we are able to rewrite the above equations (A.11) and (A.12) as y(1 y) 1 [ 1 2y 2 ArcTan 2 = 2 [ y 1 2F 2, 2, 5 ] 2,y y(1 y) + 1 2 ArcTan [ 1 2y 2 1 ] + π y(1 y) 4 = τ,when τ τ ta, (A.1) 1 ] + π y(1 y) 4 = π 2 2 y 2 F [ 1 2, 2, 5 2,y ] = (τ τ ta ),when τ τ ta. (A.14) Thus, all of these three representations are equal to each other. We also check the perturbation is linear at early epoch τ δ lin. We use equations (A.1) and (A.11) because τ τ ta in this case. Thus when we choose the limit τ, we obtain 2 x 2 τ 8
and 2 y 2 + 1 5 y 5 2 τ. Thus we obtain ( ( x ) τ 1 ( 1+δ NL = 1+ ) ) 2 y 1 y 1+ 5 y δ NL δ lin = ( 5 ( ) 2 ) 2 2 τ (A.15) If we use τ = H ta Ωmta t and H ta = 2 t 1 ta in the matter dominated epoch, then we obtain ( δ lin = )2 5 t Ωmta t ta. In EdS Universe ( = π 4 and Ω mta = 1), this gives the famous result δ lin = 5( )2 5( π t 4 t ta. For the close and open Universe, δ lin (z vir ) = )2 2 τ vir = 2 (6 + 9π)2 which is equal to that of the EdS one. We thanks K. Umetsu for useful comment. References [1] P. J. E. Peebles, Large-Scale Structure of the Universe (Princeton University Press, 198). [2] J. E. Gunn and J. R. Gott III, Astrophys. J. 176, 1 (197). [] P. J. E. Peebles, Astrophys. J. 284, 49 (1984). [4] O. Lahav, P. B. Lilje, J. R. Primack, and M. J. Rees, Mon. Not. R. Astron. Soc. 251, 128 (1991). [5] P. B. Lilje, Astrophys. J. 86, L (1992). [6] C. Lacey and S. Cole, Mon. Not. R. Astron. Soc. 262, 627 (199). [7] J. D. Barrow and P. Saich, Mon. Not. R. Astron. Soc. 262, 717 (199). [8] P. T. P. Viana and A. R. Liddle, Mon. Not. R. Astron. Soc. 281, 2 (1996) [arxiv:astroph/95117]. [9] V. R. Eke, S. Cole, and C. S. Frenk, Mon. Not. R. Astron. Soc. 282, 26 (1996) [arxiv:astroph/96188]. [1] T. Kitayama and Y. Suto, Astrophys. J. 469, 48 (1996) [arxiv:astro-ph/964141]. [11] L. Wang and P. J. Steinhardt, Astrophys. J. 58, 48 (1998) [arxiv:astro-ph/98415]. [12] C. Horellou and J. Berge, Mon. Not. R. Astron. Soc. 6, 19 (25) [arxiv:astroph/54465]. [1] T. Kihara, Publ. Astron. Soc. Japan. 2, 22 (1968). [14] I. Maor and O. Lahav, JCAP 57, (25) [arxiv:astro-ph/558]. 9
[15] S. Lee and K.-W. Ng, [arxiv:95.1522]; S. Lee, [arxiv:95.474]; S. Lee and K.-W. Ng, [arxiv:96.164]; S. Lee and K.-W. Ng, [arxiv:97.218]. [16] S. Lee and K.-W. Ng, in preparation. [17] S. Lee and K.-W. Ng, in preparation. [18] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications, New York (197). [19] A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations, 2nd Edition, (Chapman and Hall/CRC, Boca Raton, 2). 1