1 Can inhomogeneities accelerate the cosmic volume expansion? 1 Tomohiro Kai, 1 Hiroshi Kozaki, 1 Ken-ichi Nakao, 2 Yasusada Nambu and 1 Chul-Moon Yoo arxiv:gr-qc/0605120v1 22 May 2006 1 Department of Physics, Graduate School of Science, Osaka City University Osaka 558-8585, Japan and 2 Department of Physics, Graduate School of Science, Nagoya University, Chikusa, Nagoya 464-8602, Japan If the expanding and contracting regions coexist in the universe, the speed of cosmic expansion can be accelerated. We construct simple inhomogeneous dust-filled universe models of which the speed of cosmic volume expansion is accelerated for finite periods. These models are made by removing spherical domains from the Einstein-de Sitter universe and filling each domain by the Lemaître-Tolman-Bondi dust sphere with the same gravitational mass that contained in the removed region. This is an exact solution of Einstein equations. We find that the the acceleration of the cosmic volume expansion is realized when the size of the contracting region is smaller than the horizon radius of the Einstein-de Sitter universe. 1. Introduction The accelerated expansion of the universe indicated by the observational data of the luminosity distance of Type Ia supernovae and the cosmic microwave background radiation is a great mystery in modern cosmology. The acceleration of a homogeneous and isotropic Friedmann-Robertson-Walker (FRW) universe implies that the existence of the dark energy component of the matter fields which violate energy conditions and we do not know the true character of them. Conversely, if the universe is not homogeneous and isotropic, the acceleration of the cosmic volume expansion does not necessarily mean the existence of the dark energy. Thus to explain the cosmic acceleration without introducing exotic matter fields, the argument about the effects of the inhomogeneities on the cosmic expansion law has been arisen. 1) 5) This approach aims to realize the apparent acceleration of our universe by using backreaction effect due to inhomogeneities in the background FRW universe. However, in our universe, it might be possible to treat the small scale non-linear structures by the Newtonian approximation whereas the amplitude of the metric fluctuations for the general relativistic large scale structure is so small that the linear perturbation approximation might be workable, and therefore the second order effects on the global expansion law of the universe do not seem to be significant. 5) Apart from the magnitude of inhomogeneities, the perturbative analysis of the backreaction effect on the FRW universe shows that spatial inhomogeneities behave as a positive spatial curvature and it reduces the expansion rate of a FRW universe. 6) 8) Further, recently Kasai, Asada and Futamase conjectured that the nonlinear backreaction cannot accelerate the cosmic volume expansion even if the inhomogeneities are large. 9) typeset using PTPTEX.cls Ver.0.9
2 Tomohiro Kai et al. On the other hand, as pointed out by one of the present authors, Nambu, and his collaborator, Tanimoto, 10) non-perturbative feature of inhomogeneity is crucial to obtain a workable model that explains the accelerated expansion of the universe by the backreaction effect. They proposed a model of the universe containing both expanding and contracting regions and the spatially averaged scale factor defined through the three volume of a spatial domain can show accelerated expansion if the sizes of each region are properly chosen. As they have not presented a correct example of the accelerated universe, it is necessary to check whether their mechanism really works. In this paper, we study the comoving volume of the highly inhomogeneous dust filled universe and show that, in principle, the acceleration of the cosmic volume expansion can be realized due to the inhomogeneities. This paper is organized as follows. In Sec. II, we give a qualitative scenario in which inhomogeneities lead to the acceleration of the cosmic volume expansion. In Sec. III, in order to get quantitative information about the acceleration of the cosmic volume expansion, we study the stuffed-up Swiss-cheese model. Then Sec. IV is devoted to summary and discussion. In this paper, we adopt c = G = 1 unit and basically follow the convention of the Riemann and metric tensors and the notation in the textbook by Hawking and Ellis. 11) 2. A scenario: qualitative argument We consider an inhomogeneous dust filled universe which is initially expanding everywhere. We assume that the dust begins contracting in some domains and continues expanding elsewhere (see FIG. 1). Fig. 1. An inhomogeneous universe with expanding and contracting regions. Such a situation agrees with our conventional picture of the real universe; the structures in our universe (stars, galaxies, etc) are formed within the contracting domains in the expanding universe. Further we assume that the universe is almost periodic and thus homogeneous in average sense. We adopt the dust-comoving Gaus-
Can inhomogeneities accelerate the cosmic volume expansion? 3 sian normal coordinate system in which the line element is given by ds 2 = dt 2 + γ ij (t,x k )dx i dx j, (2.1) where i,j,k = 1,2,3 represent the spatial components, and γ ij (t,x k ) is the induced metric on the spacelike hypersurface t =constant which is orthogonal to the trajectories of dust particles. The stress-energy tensor of dust is given by T µν = ρu µ u ν, (2.2) where ρ is the rest mass density of dust and u µ = δ µ 0 is the 4-velocity of dust particles. We consider a compact domain D in the t =constant hypersurface, which is dustcomoving and thus contains fixed rest mass of dust. We assume that the domain D is so large that the spatial periodicity is recognizable in this domain. The volume V D of D is defined by V D = γ d 3 x, (2.3) D where γ is the determinant of the spatial metric γ ij. Following Ref., 9),12) we define the spatially averaged scale factor a D (t) of the domain by 3ȧ D /a D = V D /V D, where the dot means the derivative with respect to t. This definition is equivalent to a D (t) V 1/3 D (t). (2.4) Since we assume that the expanding and contracting domains coexist within this domain D, we rewrite the volume V D in the form V D = V e + V c, (2.5) where V e and V c are the comoving volume of the union of the expanding domains and that of the contracting domains, respectively. The time derivative of the volume is written as V D = V e + V c = V e V c. (2.6) We see from the above equation that even though V D is initially positive, once V c becomes dominant, VD can almost vanish or be negative. In other words, due to the appearances of the contracting domains in D, the cosmic volume expansion of D can strongly slow down or the volume V D of D can decrease. However, the contracting domains will collapse and then their contributions to the volume V D necessarily become negligible; in the realistic situations, almost stationary structures like as stars, black holes and galaxies, etc are formed. This means that V e will be dominant in V D and the volume V D begins increasing again. Here it should be noted that the acceleration of volume expansion is realized around the end of the slowdown or contraction stage. In this period, the second order time derivative of the effective scale factor a D (t) becomes positive and thus the acceleration of cosmic expansion will be realized. This is the mechanism of the acceleration of the cosmic volume expansion pointed out by Nambu and Tanimoto. 10)
4 Tomohiro Kai et al. 3. Stuffed-Up Swiss-Cheese Model In order to quantitatively study the effects of inhomogeneities in the universe, we construct a specific model. As mentioned, the coexistence of the contracting and expanding regions is the origin of the acceleration of the cosmic volume expansion. The model given by Nambu and Tanimoto (NT) realizes such a situation. 10) NT model is a spherically symmetric but inhomogeneous dust universe in which shell crossing singularities necessarily occur at the boundaries between the contracting and expanding regions. After the shell crossing singularity forms, the calculation of the volume becomes impossible in this model. In this paper, to overcome this shortcoming in NT model, we construct a model in which the distribution of the contracting and expanding regions is different from NT model. First we consider the Einstein-de Sitter universe (EdS) and remove disconnected spherical domains from it so as to guarantee the spatial periodicity. Next each removed region is stuffed up with the spherically symmetric but inhomogeneous dust whose Misner-Sharp mass 13) is the same as that of the removed EdS regions. The dynamics of the inhomogeneous dust ball is described by the so called Lemaître-Tolman-Bondi (LTB) solution. We call this model the stuffed-up Swiss-Cheese (SSC) model. We assume that the LTB regions are initially expanding but eventually begin to contract and form singularities and black holes. The line element of the LTB region is given by ds 2 = dt 2 + Y 2 (t,χ) 1 k(χ)χ 2dχ2 + Y 2 (t,χ)(dθ 2 + sin 2 θdϕ 2 ), (3.1) where the prime denotes the differentiation with respect to the radial coordinate χ which is assumed to be non-negative. This coordinate system is dust-comoving Gaussian normal coordinate as in the previous section. Einstein equations lead to the equations for the areal radius Y (t,χ) and the rest mass density ρ(t,χ) of the dust; Ẏ 2 = k(χ)χ 2 + 2M(χ), Y (3.2) ρ = M (χ) 4πY Y 2, (3.3) where k(χ) and M(χ) are arbitrary functions of the radial coordinate. We assume the non-negativity of ρ and monotonicity of Y with respect to χ, i.e., Y > 0 in the regular region. This assumption leads M 0 and thus we can set M(χ) = 4πρ 0 3 χ3, (3.4) where ρ 0 is a non-negative arbitrary constant. The above choice of M(χ) does not loose any generality. Equations (3.1)-(3.3) are invariant under the rescaling of the radial coordinate χ, χ χ = χ(χ). (3.5)
Can inhomogeneities accelerate the cosmic volume expansion? 5 By this property, we can choose the above form of M(χ) generally. The solution of Eq. (3.2) is given by Y = 4πρ ( 0 3k(χ) ( t t i (χ) = 4πρ 0 3k(χ) 1 cos η ( k(χ) η )) χ, (3.6) 1 k(χ) sin ( k(χ) η ) ), (3.7) where t i (χ) is an arbitrary function. We can use this expression of the solution for any sign of k(χ). Note that t i (χ) is the time when a shell focusing singularity appears, where shell focusing singularity means Y = 0 for χ > 0 and Y = 0 at χ = 0. In this paper, we consider a region of t > t i, and hence the time t = t i corresponds to the Big Bang singularity. Hereafter, we assume a simultaneous Big Bang t i = 0. By using the terminology of the cosmological perturbation theory, there are only growing modes near Big Bang singularity in this setting. We denote the boundary between the LTB region and the EdS region by χ = χ b. Then dividing the LTB region into four regions [0,χ 1 ), [χ 1,χ 2 ), [χ 2,χ 3 ) and [χ 3,χ b ). We use the following spatial profile of the curvature function: k 0 for 0 χ < χ 1 { k 0 (χ 2 χ 2 } 2 )2 2χ 2 χ 2 1 + χ 2 χ2 1 + χ2 2 for χ 1 χ < χ 2 2 k(χ) = k 0 ( χ 2 2χ 2 1 + χ 2 ), (3.8) 2 for χ 2 χ < χ 3 k 0 ( χ 2 2χ 2 1 + χ 2 ) { ( χ 2 χ 2 ) 2 } 2 3 2 χ 2 b 1 for χ 3 χ < χ b χ2 3 where k 0 is constant. In order to guarantee 1 kχ 2 > 0, the following inequality should hold, κ := k 0 2 (χ2 1 + χ 2 2) < 1. (3.9) Since we are interested in the case that the LTB regions will contract, we assume k 0 > 0. The singularity formed at the origin χ = 0 can be null naked singularity, whereas that formed in the domain χ > 0 is spacelike. 14) 16) In the case of the present model, since the innermost region [0,χ 1 ) is a Friedman universe with a positive spatial curvature, the singularity formed at χ = 0 is spacelike and thus not naked. We see from Eq. (3.7) that the singularity formation time is given by t = 8π 2 ρ 0 /(3k 3/2 ). Since k is monotonically decreasing with respect to χ, the singularity formation time is monotonically increasing with respect to χ and thus the shell crossing singularity does not occur in the present model. Therefore the singularities formed by the gravitational collapse in the LTB regions are necessarily spacelike in the present SSC model. There is no causal influence of the spacelike singularities on the regular
6 Tomohiro Kai et al. domains. This property allows us to ignore the contribution of the singularity in the calculation of the volume V D of the domain D and thus we can see the evolution of the cosmic volume expansion even after the singularity formation. Further it should be noted that the present choice of k(χ) guarantees that the LTB region is smoothly matched to the EdS universe. The line element of the EdS region is written in the form ds 2 = dt 2 + a 2 e (t)[ dχ 2 + χ 2 ( dθ 2 + sin 2 θdϕ 2)], (3.10) where a e (t) = ( 6πρ 0 t 2) 1/3. (3.11) We assume that the universe has the periodicity of the comoving interval 2l in the original Einstein-de Sitter universe. By this assumption, a cubic region of the comoving edge length 2l corresponds to the compact domain D considered in the previous section. We define the scale factor a(t) of this inhomogeneous universe as 3.1. One scale model a D (t) V 1/3 D. (3.12) 2l We first consider the model with one LTB region of χ b = l in the domain D (see FIG. 2). Fig. 2. The compact domain D of the one-scale SSC model is depicted. Thus in this model, the identical LTB regions are ranged in close order. The volume V D of the domain D is given by ( V D = 8 4π ) ae(t)l a 3 3 e(t)l 3 Y 2 dy + 4π. (3.13) 1 k(χ)χ 2 The integral in the right hand side of the above equation corresponds to the comoving volume of the LTB region. It should be noted that χ in this integral is the function 0
Can inhomogeneities accelerate the cosmic volume expansion? 7 of the areal radius Y at time t through Eqs. (3.6) and (3.7). Thus this integration covers the only region which has not yet collapsed to the singularity, i.e., the domain of Y > 0 at time t. In FIG. 3, we show several examples of the evolution of the scale factor a D (t). We set H 0 l = 6/2,H 0 χ 1 = 6/8,H 0 χ 2 = 6/4 and H 0 χ 3 = 3 6/8 with H 0 = 8πρ0 /3. We determined these values so that the areal radius of the boundary of the LTB region is equal to the horizon radius of the EdS universe at t = l. Each curves represent the scale factor as functions of t for various values of κ defined by Eq. (3.9). As mentioned in the preceding section, slowdown of the volume expansion or the contraction of the domain D really occurs due to the contraction of the LTB region. After the LTB region collapses to a singularity covered by a horizon and forms a black hole, the volume expansion of the domain D, which has been braked by the contracting LTB region, is accelerated until the speed of the volume expansion becomes the same as that of the EdS universe. Thus at the time around the end of the slowdown or contraction stage, the acceleration period appears. The necessary condition for the appearance of the acceleration period is that the volume of the LTB region becomes dominant in V D. If the parameter κ is almost equal to unity, the contribution of χ 2 χ < χ 3 is dominant in the volume of the LTB region and further makes the volume of the LTB region dominant in V D. Thus it is expected that κ almost equal to unity realizes the acceleration of the cosmic volume expansion, and we can see from FIG. 3 that this expectation is true. Fig. 3. The left panel shows the evolution of the averaged scale factor a D of one scale SSC model for various values of κ. The right panel shows the acceleration of the averaged scale factor a D. The size of the LTB region relative to the cosmological horizon is very important. If there are many contracting LTB regions within the cosmological horizon of the EdS universe in the acceleration period, the SSC model might mimic our real universe. By contrast, if the size of the LTB region is comparable to the cosmological horizon in the acceleration period, this SSC model might be very different from our universe. Thus we investigate which of two cases corresponds to this SSC model. It should be again noted that when the volume of the domain χ < χ 3 is dominant in the volume V D and is contracting, the volume expansion of the domain D strongly slows down
8 Tomohiro Kai et al. or the domain D contracts. The acceleration of the volume expansion is realized when the volume of the domain χ < χ 3 becomes so small that its contribution to V D is negligible. In other words, the acceleration period appears around the time at which the dust in the domain χ < χ 3 collapses to the singularity and this time is given by t ac := 8π2 ρ 0 χ 3 3 3 κ 3/2. (3.14) At t = t ac, the ratio of the areal radius of the boundary of a LTB region a e (t)l to the Hubble horizon radius of EdS region H 1 (t) := a e /ȧ e = 3t/2 becomes a e (t ac )l H 1 (t ac ) = ( 2 3π ) 1/3 κ 1/2 l χ 3 0.60κ 1/2 ( l χ 3 ). (3.15) Here it should be noted that κ has to be almost equal to unity so that the acceleration of the cosmic volume expansion is realized. Therefore the size of the LTB region is smaller than the horizon radius of EdS region in the acceleration period. This means that the sub-horizon scale inhomogeneities can mimic our universe with the acceleration period. 3.2. Two scale model As the second model, we consider the situation that there are two kinds of LTB regions; one has the radius χ b = l and the other has the radius χ b = ( 3 1)l. We assume that the large LTB region is put in the cubic domain D in the same manner as the first model, whereas an eighth of the small LTB region is put on each vertex of D (see FIG. 4). Thus in this model, the two kinds of the LTB regions are ranged in closed order in the universe. Fig. 4. The compact domain D of the two scale SSC model is depicted. In this model, the volume of compact domain D is given by [ V D = 8 4π { ( ) }] 3 1 + 3 1 a e (t)l 3 3
Can inhomogeneities accelerate the cosmic volume expansion? 9 ae(t)l + 4π 0 Y 2 dy 1 k(χ L )χ 2 L ( 3 1)ae(t)l + 4π 0 Y 2 dy 1 k(χ S )χ 2 S, (3.16) where χ L = χ L (Y ) and χ S = χ S (Y ) are the radial coordinates as functions of the areal radius Y at time t in the large and small LTB regions through Eqs. (3.6) and (3.7), respectively. We show one example of the evolution of the scale factor a D (t) in FIG. 5. As expected, there are two acceleration periods. The first acceleration period is realized when the small LTB regions collapse to the singularities, whereas the second acceleration period appears when the large LTB region collapses. Thus the coexistence of various scales of LTB regions can lead more than one acceleration periods. Fig. 5. The left panel shows the evolution of the averaged scale factor of two scale SSC model. We set κ = 0.999688 in the large LTB region and κ = 0.999995 in the small LTB region. Except for these values of κ, we used the same values of parameters as the one scale SSC model. 4. Summary and Discussion By using the stuffed-up Swiss-Cheese model, we investigated how the cosmic volume expansion is accelerated by inhomogeneities in the dust-filled universe. We found that the existence of contracting region with the sub-horizon scale in the EdS region will be necessary condition to realize the acceleration of the cosmic volume expansion. About the time at which the central contracting region collapse to form a singularity, the acceleration of the averaged scale factor becomes positive and the acceleration phase of the spatially averaged universe appears. Such a situation might be the same as our real universe; the structures (e.g. galaxies, clusters of galaxies) will form within the contracting regions in our expanding universe. For the model with various scales of inhomogeneities, the temporal variation of the cosmic volume expansion is nontrivial and the complete analysis is left as a future work. The result obtained in this paper shows that the averaged scale factor for the inhomogeneous dust filled universe can minic the cosmic acceleration. However, the analysis of the backreaction effect for the dust filled universe based on the second order cosmological perturbation shows no signs of acceleration of the volume expan-
10 Tomohiro Kai et al. sion. 6),7),9) Thus, the appearance of the acceleration phase in our model is due to highly non-linear effect expected to appear in higher order of the perturbation more than the third order. Here we should note that the issue of the cosmic acceleration has arisen from the observational data of the distance-redshift relation and that we have not directly observed the acceleration of the volume expansion. If the universe is homogeneous and isotropic, then we can conclude from the observational data that the cosmic volume expansion accelerates. By contrast, in the inhomogeneous universe model considered in this paper, the cosmic expansion of the comoving volume is accelerated but it is non-trivial whether the distance-redshift relation is similar to that of the homogeneous and isotropic universe with accelerated volume expansion. There are several studies about the inhomogeneous but isotropic dust filled universe whose distance-redshift relation agrees with the observational data. 17) 20) In these models, there is no contracting region and therefore the cosmic volume expansion will not be accelerated in the sense of this paper. Thus there is a possibility that the acceleration of the cosmic volume expansion due to inhomogeneities might not lead the distanceredshift relation similar to that obtained by the observation of type Ia SNe. This is also a future work. Acknowledgments We are grateful to colleagues in the astrophysics and gravity group of Osaka City University for helpful discussion and criticism. References 1) E. W. Kolb, S. Matarrese, A. Notari and A. Riotto, hep-th/0503117 (2005). 2) E. W. Kolb, S. Matarrese, A. Riotto, astro-ph/0506534 (2005). 3) É. É. Flanagan, Phys. Rev. D 71 (2005), 103521. 4) C. M. Hirata and U. Seljak, Phys. Rev. D 72 (2005), 083501. 5) A. Ishibashi and R. M. Wald, Class. Quantum Grav. 23 (2006), 235. 6) H. Russ, M. H. Soffel, M. Kasai and G. Börner, Phys. Rev. D 56 (1997), 2044. 7) Y. Nambu, Phys. Rev. D 62 (2000), 104010. 8) H. Kozaki and K. Nakao, Phys. Rev. D 66 (2002), 104008. 9) M. Kasai, H. Asada and T. Futamase, Prog. Theor. Phys. 115 (2006), 827. 10) Y. Nambu and M. Tanimoto, gr-qc/0507057 (2005). 11) S. W. Hawking and G. F. R. Ellis, The large scale structure of space-time, Cambridge University Press (1973). 12) M. Kasai, Phys. Rev. Lett. 69 (1992), 2330. 13) C. W. Misner and D. H. Sharp, Phys. Rev. 136 (1964), B571. 14) D. Christodoulou, Commun. Math. Phys. 93 (1984), 171. 15) R. P. A. C. Newman, Class. Quantum Grav. 3 (1986), 527. 16) P. S. Joshi and I. H. Dwivedi, Phys. Rev. D 47 (1993), 5357. 17) M. Celerier, Astron. Astrophys. 353 (2000), 63. 18) K. Tomita, Mon. Not. R. Astron. Soc. 326 (2001), 287. 19) H. Iguchi, T. Nakamura and K. Nakao, Prog. Theor. Phys. 108 (2002), 809. 20) R. A. Vanderveld, É. É. Flanagan and I. Wasserman, astro-ph/0602476 (2006).