Next Contents Previous

5.1. Linear evolution of perturbations

When the perturbations are small, one can ignore the second term in the right hand side of (65) and replace (1 + delta) by unity in the first term on the right hand side. The resulting equation is valid in the linear regime and hence will be satisfied by each of the Fourier modes deltak(t) obtained by Fourier transforming delta(t, x) with respect to x. Taking delta(t, x) = D(t) f (x), the D(t) satisfies the equation

Equation 66 (66)

The power spectra P(k, t) = < |deltak(t)|2 > at two different redshifts in the linear regime are related by

Equation 67 (67)

where T (called transfer function) depends only on the parameters of the background universe (denoted by `bg') but not on the initial power spectrum and can be computed by solving (66). It is now clear that the only new input which structure formation scenarios require is the specification of the initial perturbation at all relevant scales, which requires one arbitrary function of the wavenumber k = 2pi / lambda.

Let us first consider the transfer function. The rate of growth of small perturbations is essentially decided by two factors: (i) The relative magnitudes of the proper wavelength of perturbation lambdaprop(t) propto a(t) and the Hubble radius dH(t) ident H-1(t) = (adot / a)-1 and (ii) whether the universe is radiation dominated or matter dominated. At sufficiently early epochs, the universe will be radiation dominated and dH(t) propto t will be smaller than the proper wavelength lambdaprop(t) propto t1/2. The density contrast of such modes, which are bigger than the Hubble radius, will grow as a2 until lambdaprop = dH(t). [When this occurs, the perturbation at a given wavelength is said to enter the Hubble radius. One can use (66) with the right hand side replaced by 4pi(1 + w)(1 + 3w) G rho in this case; this leads to D propto t propto a2.] When lambdaprop < dH and the universe is radiation dominated, the perturbation does not grow significantly and increases at best only logarithmically [193]. Later on, when the universe becomes matter dominated for t > teq, the perturbations again begin to grow. It follows from this result that modes with wavelengths greater than deq ident dH(teq) - which enter the Hubble radius only in the matter dominated epoch - continue to grow at all times while modes with wavelengths smaller than deq suffer lack of growth (in comparison with longer wavelength modes) during the period tenter < t < teq. This fact leads to a distortion of the shape of the primordial spectrum by suppressing the growth of small wavelength modes in comparison with longer ones. Very roughly, the shape of T2(k) can be characterized by the behaviour T2(k) propto k-4 for k > keq and T2 approx 1 for k < keq. The wave number keq corresponds to the length scale

Equation 68 (68)

(eg., [44], p.75). The spectrum at wavelengths lambda >> deq is undistorted by the evolution since T2 is essentially unity at these scales. Further evolution can eventually lead to nonlinear structures seen today in the universe.

At late times, we can ignore the effect of radiation in solving (66). The linear perturbation equation (66) has an exact solution (in terms of hyper-geometric functions) for cosmological models with non-relativistic matter and dark energy with a constant w. It is given by

Equation 69 (69)

Figure 15 shows the growth factor for different values of w including the one for cosmological constant (corresponding to w = - 1) and an open model (with w = -1/3.)

Figure 15

Figure 15. The growth factor for different values of w including the one for cosmological constant (corresponding to w = - 1) and an open model (with w = - 1/3).

For small values of a, D approx a which is an exact result for OmegaLambda = 0, OmegaNR = 1 model. The growth rate slows down in the cosmological constant dominated phase (in models with OmegaNR + OmegaLambda = 1 with w = - 1) or in the curvature dominated phase (open models with OmegaNR < 1 corresponding to w = - 1/3). Between the two cases, there is less growth in open models compared to models with cosmological constant.

It is possible to rewrite equation (65) in a different form to find an approximate solution for even variable w(a). Converting the time derivatives into derivatives with respect to a (denoted by a prime) and using the Friedmann equations, we can write (65) as

Equation 70 (70)

In a universe populated by only non relativistic matter and dark energy characterized by an equation of state function w(a), this equation can be recast in a different manner by introducing a time dependent Omega [as in equation (43)] by the relation Q(t) = (8pi G / 3) [rhoNR(t) / H2(t)] so that (dQ/d ln a) = 3wQ(1 - Q). Then equation (65) becomes in terms of the variable f ident (d lndelta / d ln a)

Equation 71 (71)

Unfortunately this equation is not closed in terms of f and Q since it also involves delta = exp[integ(da / a)f]. But in the linear regime, we can ignore the second term on the right hand side and replace (1 + delta) by unity in the first term thereby getting a closed equation:

Equation 72 (72)

This equation has approximate power law solutions [189] of the form f = Qn when |dw / dQ| << 1 / (1 - Q). Substituting this ansatz, we get

Equation 73 (73)

[Note that Q(t) rightarrow 1 at high redshifts, which is anyway the domain of validity of the linear perturbation theory]. This result shows that n is weakly dependent on OmegaNR; further, n approx (4/7) for open Friedmann model with non relativistic matter and n approx (6/11) approx 0.6 in a k = 0 model with cosmological constant.

It is possible to provide simple analytic fitting functions for the transfer function, incorporating all the above effects. For models with a cosmological constant, the transfer function is well fitted by [194]

Equation 74 (74)

where p = k / (Gamma h Mpc-1) and Gamma = OmegaNR h exp[- OmegaB(1 + (2 h)1/2 / OmegaNR)] is called the `shape factor'.

Let us next consider the initial power spectrum P(k, zi). The following points need to be emphasized regarding the initial fluctuation spectrum.

(1) It can be proved that known local physical phenomena, arising from laws tested in the laboratory in a medium with (P/rho) > 0, are incapable producing the initial fluctuations of required magnitude and spectrum (eg., [9], p. 458). The initial fluctuations, therefore, must be treated as arising from physics untested at the moment.

(2) Contrary to claims sometimes made in the literature, inflationary models are not capable of uniquely predicting the initial fluctuations. It is possible to come up with viable inflationary potentials ([196], chapter 3) which are capable of producing any reasonable initial fluctuation.

A prediction of the initial fluctuation spectrum was indeed made by Harrison [197] and Zeldovich [198], who were years ahead of their times. They predicted - based on very general arguments of scale invariance - that the initial fluctuations will have a power spectrum P = Akn with n = 1. Considering the simplicity and importance of this result, we shall briefly recall the arguments leading to the choice of n = 1.

If the power spectrum is P propto kn at some early epoch, then the power per logarithmic band of wave numbers is Delta2 propto k3P(k) propto k(n+3). Further, when the wavelength of the mode is larger than the Hubble radius, dH(t) = (adot / a)-1, during the radiation dominated phase, the perturbation grows as a2 making Delta2 propto a4k(n+3). We need to determine how Delta scales with k when the mode enters the Hubble radius dH(t). The epoch aenter at which this occurs is determined by the relation 2pi aenter / k = dH. Using dH propto t propto a2 in the radiation dominated phase, we get aenter propto k-1 so that

Equation 75 (75)

It follows that the amplitude of fluctuations is independent of scale k at the time of entering the Hubble radius, only if n = 1. This is the essence of Harrison-Zeldovich and which is independent of the inflationary paradigm. It follows that verification of n = 1 by any observation is not a verification of inflation. At best it verifies a far deeper principle of scale invariance. We also note that the power spectrum of gravitational potential Pphi scales as Pphi propto P / k4 propto k(n-4). Hence the fluctuation in the gravitational potential (per decade in k) Delta2phi propto k3 Pphi is proportional to Delta2phi propto k(n-1). This fluctuation in the gravitational potential is also independent of k for n = 1 clearly showing the special nature of this choice.[It is not possible to take n strictly equal to unity without specifying a detailed model; the reason has to do with the fact that scale invariance is always broken at some level and this will lead to a small difference between n and unity]. Given the above description, the basic model of cosmology is based on seven parameters. Of these 5 parameters (H0, OmegaB, OmegaDM, OmegaLambda, OmegaR) determine the background universe and the two parameters (A, n) specify the initial fluctuation spectrum.

The presence of dark energy, with a constant w, will also affect the transfer function and hence the final power spectrum. An approximate fitting formula can be given along the following lines [195]. Let the power spectrum be written in the form

Equation 76 (76)

where AQ is a normalization, TQ is the modified transfer function and gQ = (D / a) is the ratio between linear growth factor in the presence of dark energy compared to that in Omega = 1 model. Writing TQ as the product TQLambda TLambda where TLambda is given by (74), numerical work shows that

Equation 77 (77)

where k is in Mpc-1, and alpha is a scale-independent but time-dependent coefficient well approximated by alpha = (-w)s with

Equation 78 (78)

where the matter density parameter is OmegaNR(a) = OmegaNR / [OmegaNR + (1 - OmegaNR) a-3w]. Similarly, the relative growth factor can be expressed in the form gQLambda ident (gQ / gLambda) = (- w)t with

Equation 79 (79)

Finally the amplitude AQ can be expressed in the form AQ = deltaH2(c / H0)n+3 / (4pi), where

Equation 80 (80)

and alpha0 = alpha(a = 1) of equation (78), and

Equation 81 (81)

This fit is valid for -1 ltapprox w ltapprox - 0.2.

Next Contents Previous