## Brazilian Journal of Physics

*On-line version* ISSN 1678-4448

### Braz. J. Phys. vol.31 no.1 São Paulo Mar. 2001

#### http://dx.doi.org/10.1590/S0103-97332001000100007

**Gravitational clustering to all perturbative orders **

Elcio Abdalla^{a} and Roya Mohayaee^{b}

^{a}Instituto de Física, Universidade de São Paulo,

*C.P. 66318, CEP 05315-970, São Paulo, Brazil*

^{b}Dipartimento di Fisica Galileo Galilei (INFN)

*Via Marzolo 8, 35131 Padova - Italy*

**Received on 8 May, 2000 **

We derive the time evolution of the density contrast to all orders of perturbation theory, by solving the Einstein equation for scale-invariant fluctuations. These fluctuations are represented by an infinite series in inverse powers of the radial parameter. In addition to the standard growing modes, we find infinitely many more new growing modes for open and closed universes.

The large scale structure of the universe is believed to have grown, due to gravitational instability, from small primordial density fluctuations. These fluctuations are fully characterized by the density contrast, d(*t*, )=dr(*t*, )/r(*t*). At any given time, the connected correlation functions of the density contrast, determine the power spectrum and the space distribution of the matter in the universe [1, 2, 3, 4, 5]. The time dependence of the density contrast, on the other hand, determines if the inhomogeneities would indeed grow, or would just oscillate or decay, in the first place.

In this article, we are primarily concerned with the time evolution of the density contrast. In almost all standard methods, this is determined by the continuity, Euler and Poisson equations in the Newtonian regime [1, 6, 7, 8, 9, 10] and by the Einstein equation in the relativistic era [1, 6, 7, 8, 9, 10]. In the nonlinear regime, the perturbation is carried to second order and new growing modes containing nonlocal terms are obtained for dust universes [1, 5, 11]. The local nonlinear terms have also been studied to all orders for the flat universe [4, 12]. A single complete treatment covering all possible cases and also extending to an arbitrary perturbative order does not exist.

Our starting point is not the usual equations of fluid hydrodynamics but the Einstein equation which naturally covers the present era too. We expand the Einstein equation around the background Friedmann universe by scale-invariant fluctuations. These fluctuations are expressed as an infinite series in inverse powers of the radial parameter. This method, which assumes spherical symmetry, allows us to obtain the full expression for the growth rate of density inhomogeneities in Newtonian, relativistic, zero and non-zero curvatures, linear and local-nonlinear regimes in a single treatment. Thus in this way, we make major simplifications on current perturbative methods. In a flat universe, the leading growing modes are given by an infinite series in the matter and radiation -dominated eras. For the open and closed universes we find, in addition to the modes presented in the literature [8, 11], infinitely many more new growing modes which are also given by infinite sums. In the nonlinear regime, we also find new growing modes for closed and open universes in the extreme limits of small and large times.

Contrary to the standard analysis, where the density contrast is often taken to be a random Gaussian field, we take it to be a scale-invariant quantity. That the standard growing modes are contained in our results, is a further solid evidence that the time evolution of the density contrast is unaffected by its statistical properties.

Since, the inverse radial parameter now plays the rôle of the perturbation parameter, the order by which the structures at different scales are formed is easily determined. Our perturbation scheme implies a bottom-up or a hierarchical clustering scenario for the formation of structures and thus a universe dominated by cold dark matter [8].

To account for the inhomogeneity of space-time, we use the spherically symmetric metric,

where the scale factor *R _{p}*(

*t*,

*r*) is a function of both time and coordinate. This metric does not assume homogeneity and is contained in the Tolman metric for a pressure-less universe [13], although not restricted to this era. We assume that the inhomogeneities are given by scale-invariant fluctuations and expand the scale factor as

around the background Friedmann-Robertson-Walker scale factor, *R*(*t*). This is a genuine perturbation of the metric and not just a gauge mode. A solution z_{m} of d*g*_{mn}= z_{m},n + z_{n},m where d*g*_{mn} is the metric perturbation, such that d*g*_{0r }= 0 and d*g _{rr}* ~

*t*

^{a}*r*

^{-n}, for a general

*a*and

*n*, cannot be constructed [14].

The time-time and radial-radial components of the Einstein equation are expanded as in (2) and, on using the equation of state *P*=wr, are reduced to the following second-order recursive differential equation [15]:

where summations over the indices *n*, *m*, *l* and *s* are implied and d*R _{n}* represents

*R*+ d

*R*

_{0}when

*n*= 0. We have solved the above equation by a Maple program to very large orders for a flat universe,

*i.e.*

*k*= 0 [18]. The solutions are substituted back in the time-time component of the Einstein equation and a series expression for the density and subsequently the leading-order growing modes for the density contrast are found. The fastest growing modes can be written in the closed forms

in the matter-dominated era, w = 0, and

in the radiation-dominated era. We see that in addition to the usual growing modes, *t*^{2/3} and *t* in the matter and radiation dominated eras respectively, we have infinitely many more growing modes. The higher growing modes are the local nonlinear modes [1, 12]. The above density contrasts also tell us that smaller structures enter nonlinear regime, d > 1, faster than the larger structures. Thus a hierarchical clustering scenario is expected in this scheme.

For the closed and open universes, the recursive equation (3) is highly inhomogeneous and cannot be solved. However, assuming that different perturbative orders evolve independently of each other, we find, in the matter-dominated era, the linearized equation

where y is real for an open and is imaginary, y = -*i*q, for a closed universe. The above equation can be solved exactly. The solutions can be used in the Einstein equation to obtain the density. The growing modes of the density contrast are given by

We see that in addition to the usual growing mode [8] which is given by *n *= 0 terms in the above equation, infinitely many more growing modes exist in the linear regime. The *n* ¹ 0 modes have been overlooked in the previous analyses and can have important consequences for the power spectrum.

In the nonlinear regime, the full inhomogeneous equation (3) for open and closed universes can also be solved analytically at low orders and by a fully algebraic Maple program at higher orders [15] in the limits of small and large times. At small times, the first growing modes d_{+} ~ *t*^{2/3} and d_{+} ~ *t* in the matter and radiation-dominated eras occur at the lowest order in the perturbation series, *i.e.* at the r-independent order [15]. At higher orders the solution contains oscillatory and or decaying modes only. Therefore, open and closed universes behave as a flat universe at small times in the linear regime, as expected. However, this statement is not completely true. In a flat universe the first growing modes appear at the third, *r*^{-3}, rather than the first, order in the perturbation expansion, as expressed by equations (4) and (5). In other words, at small times the inhomogeneities will enter the nonlinear phase faster in the open and closed universes than in the flat universe. Therefore, unlike common expectations, structures can grow faster in an open universe than in a flat universe.

At large times, a new growing mode arises in an open universe. Taking the asymptotic limit of the hypergeometric solutions to the Einstein equation [15] we obtain

This result indicates that, at very large times, inhomogeneities will only grow in an open universe if the radiation pressure is non-vanishing. However, since the baryonic matter cannot grow in the radiation-dominated era due to its strong coupling to the radiation, the above mode is only relevant for the growth of the perturbation in the non-baryonic dark matter component.

We, therefore, conclude that for an open universe, in addition to the growing solution appearing at the zero mode (*i.e.* at *n *= 0 in equation (7)), many more growing modes exist. These modes have not been obtained previously, even though the flat universe modes, (4) and (5), have been well known in the context of the relativistic theory of small perturbations and spherical top hat model (see for example [8]). Our method of finding the evolution of a specific type of perturbation, namely a scale-invariant fluctuation, enables us to obtain the expression (3), which can be solved in general without being restricted to an Einstein-de Sitter dust universe in a linear regime. The extra growing modes presented in equation (7) can have significant consequences for the formation of structure in an open universe. It would also be interesting to obtain the nonlinear-nonlocal growing modes [1]. We shall discuss both these issues in future works.

E.A thanks Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq-Brazil) and Fundação de Amparo a Pesquisa do Estado de São Paulo (FAPESP) for financial support.

**References **

[1] P.J.E. Peebles, *The large scale structure of universe* (Princeton University Press, NJ, 1980). [ Links ]

[2] J.N. Fry, Astrophys. J. **279**, 499 (1984). [ Links ]

[3] M.H. Goroff, B. Grinstein, S.-J. Rey and M.B. Wise, Astrophys. J. **311**, 6 (1986). [ Links ]

[4] R. Scoccimarro, S. Colombi, J.N. Fry, J.A. Frieman, E. Hivon, A. Melott, Astrophys. J. **496**, 586 (1998); [ Links ]R. Scoccimarro and J. Frieman, Astrophys. J. **473**, 620 (1996); [ Links ]R. Scoccimarro and J. Frieman, Astrophys. J. Supp. **105**, 37 (1996). [ Links ]

[5] F.R. Bouchet, S. Colombi, E. Hivon and R. Juszkiewicz, Astronom.& Astrophys. **296**, 575 (1995). [ Links ]

[6] S. Weinberg, *Gravitation and Cosmology* (John Wiley and Sons, NY, 1972). [ Links ]

[7] E.W. Kolb and M.S. Turner, *The early Universe* (Addison Wesley, Redaing, Mass., 1990). [ Links ]

[8] P. Coles and F. Lucchin, *Cosmology: The origin and evolution of cosmic structure* (John Wiley & Sons Ltd, West Sussex, 1995). [ Links ]

[9] V. Sahni and P. Coles, Phys. Rep. **262**, 1 (1995). [ Links ]

[10] F.R. Bouchet, Lecture notes for the Course CXXXII of the Enrico Fermi School on *Dark Matter in the Universe* (Varenna, 1996); Report No. astro-ph/9603013. [ Links ]

[11] H. Martel and W. Freudling, Astrophys. J. **371**, 1 (1991). [ Links ]

[12] P. Fosalba and E. Gaztañaga, Report No. astro-ph/9712095, Report No. astro-ph/9712263. [ Links ]

[13] G. Lemaître, Ann. Soc. Sci. Bruxelles **A53**, 51 (1933); [ Links ]G. Lemaître, Gen. Rel. Grav. **29**, 641 (1997); [ Links ]R.C. Tolman, Proc. Natl. Acad. Sci. **20**, 169 (1934); [ Links ]H. Bondi, Mon. Not. R. Astron. Soc. **107**, 410 (1947). [ Links ]

[14] G.F.R. Ellis and M. Bruni, Phys. Rev. **D 40**, 1804 (1989). [ Links ]

[15] E. Abdalla and R. Mohayaee, Report No. astro-ph/9810146, Report No. astro-ph/9711162. [ Links ]

[16] The complete algebraic Maple programs can be found at the site http://www.fma.if.usp.br/~ fractals/perturbation . [ Links ]