
Cosmology. The Origin and Evolution of Cosmic Structure - Coles P., Lucchin F
..pdf
28 First Principles
For a ‘dust’ universe (p = 0), which is a good approximation to our Universe at the present time, Equations (1.12.9) and (1.12.6) give
|
K |
|
|
Kc2 |
|
|
Λ = |
|
, |
ρ = |
|
. |
(1.12.10) |
a2 |
4πGa2 |
Since ρ > 0, we must have K = 1 and therefore Λ > 0. The value of Λ which makes the universe static is just
ΛE = |
4πGρ |
|
c2 . |
(1.12.11) |
The model we have just described is called the Einstein universe. This universe is static (but unfortunately unstable, as one can show), has positive curvature and a curvature radius
aE = ΛE−1/2 = |
c |
|
(4πGρ)1/2 . |
(1.12.12) |
After the discovery of the expansion of the Universe in the late 1920s there was no longer any reason to seek static solutions to the field equations. The motivation which had led Einstein to introduce his cosmological constant term therefore subsided. Einstein subsequently regarded the Λ-term as the biggest mistake he had made in his life. Since then, however, Λ has not died but has been the subject of much interest and serious study on both conceptual and observational grounds. The situation here is reminiscent of Aladdin and the genie: after he released the genie from the lamp, it took on a life of its own. For more than 60 years the genie lingered, providing neither compelling observational evidence of its existence nor strong theoretical reasons for it to be taken seriously. However, observations do now suggest that it may have been there all along. We shall return to this resurgence of Λ in the next chapter and also in Chapter 7, but in the meantime we shall restrict ourselves to brief comments on two particularly important models involving the cosmological constant, because we shall encounter them again when we discuss inflation.
The de Sitter universe (de Sitter 1917) is a cosmological model in which the universe is empty (p = 0; ρ = 0) and flat (K = 0). From Equations (1.12.6) we get
p˜ = −ρc˜ 2 |
|
Λc4 |
|
|
= − |
|
, |
(1.12.13) |
|
8πG |
||||
which, on substitution in (1.12.8), gives |
|
|
|
|
a˙2 = 31 |
Λc2a2; |
(1.12.14) |
this equation implies that Λ is positive. Equation (1.12.14) has a solution of the form
a = A exp[(31 Λ)1/2ct], |
(1.12.15) |
corresponding to a Hubble parameter H = a/a˙ = c(Λ/3)1/2, which is actually constant in time. In the de Sitter vacuum universe, test particles move away from
Friedmann Models |
29 |
each other because of the repulsive gravitational e ect of the positive cosmological constant.
The de Sitter model was only of marginal historical interest until the last 20 years or so. In recent years, however, it has been a major component of inflationary universe models in which, for a certain interval of time, the expansion assumes an exponential character of the type (1.12.15). In such a universe the equation of state of the fluid is of the form p −ρc2 due to quantum e ects which we discuss in Chapter 7.
In the Lemaître model (1927) the universe has positive spatial curvature (K = 1). One can demonstrate that the expansion parameter in this case is always increasing, but there is a period in which it remains practically constant. This model was invoked around 1970 to explain the apparent concentration of quasars at a redshift of z 2. Subsequent data have, however, shown that this is not the explanation for the redshift evolution of quasars, so this model is again of only marginal historical interest.
1.13 Friedmann Models
Having dealt with a few special cases, we now introduce the standard cosmological models described by the solutions (1.10.3) and (1.10.5). Their name derives from A. Friedmann, who derived their properties in 1922. His work was not well known at that time partly because his models were not static, and the discovery of the Hubble expansion was still some way in the future. His work was in any case not widely circulated in the western scientific literature. Independently, and somewhat later, the Belgian priest George Lemaître obtained essentially the same results and his work achieved more immediate attention, especially in England where he was championed by Eddington. When the work of Lemaître (1927) was published, Hubble’s observations were just becoming known, so in the West Lemaître is often credited with being the father of the Big Bang cosmology, although that honour should probably be conferred on Friedmann.
The Friedmann models are so important that we shall devote the next chapter to their behaviour. Here we shall just whet the readers appetite with some basic properties. First, we assume a perfect fluid with some density ρ and pressure p. The form of equation of state giving p as a function of ρ does not concern us for now; we discuss it in Section 2.1. For the moment we also ignore the cosmological constant.
The equations we need to solve are (1.10.3) and (1.10.5), which we rewrite here for completeness:
a¨ = − |
4π |
G ρ |
+ 3 |
p |
a |
(1.13.1) |
||||
3 |
c2 |
|||||||||
and |
|
|
|
|
|
|
|
|
|
|
a˙2 + Kc2 |
= |
8πG |
ρa2, |
(1.13.2) |
||||||
|
|
3 |
30 |
First Principles |
|
||
as well as the Equation (1.10.6) |
|
|||
|
|
p |
|
|
|
d(ρa3) = −3 |
|
a2 da. |
(1.13.3) |
|
c2 |
The Equations (1.13.1)–(1.13.3) allow one, at least in principle, to calculate the time evolution of a(t) as well as ρ(t) and p(t) if we know the equation of state.
Let us focus for now on Equation (1.13.3), which can be rewritten in a convenient
form for a = a0: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
a˙ |
2 |
8π |
Gρ |
a |
2 |
|
|
ρ0 |
Kc2 |
|
|||||
|
|
− |
|
|
|
= H02 |
1 − |
|
|
= H02(1 |
− Ω0) = − |
|
, |
(1.13.4) |
||
a0 |
3 |
a0 |
ρ0c |
a02 |
||||||||||||
where H0 = a˙0/a0, Ω0 is the (present) density parameter and |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
3H02 |
|
|
|
|
||
|
|
|
|
|
|
|
|
ρ0c = |
|
. |
|
|
|
(1.13.5) |
||
|
|
|
|
|
|
|
|
8πG |
|
|
|
The su x ‘0’ refers here to a generic reference time t0 which is also used in the particular case where t is the present time. Equation (1.13.5) is a reminder of the importance of ρ0c: if ρ0 < ρ0c, then K = −1, while if ρ0 > ρ0c, K = 1; K = 0 corresponds to the ‘critical’ case when ρ0 = ρ0c.
Let us now include the cosmological constant term Λ. In Section 1.12 we showed how one can treat the cosmological constant as a form of fluid with a strange equation of state, as well as a modification of the law of gravity. In that sense, Λ can be thought of as belonging either on the left-hand or right-hand side of the Einstein equations. Either way, the upshot is that Equations (1.13.1) and (1.13.2)
become |
|
|
|
|
|
|
|
|
|
|
|
|
|
a¨ = − |
4π |
G ρ + 3 |
p |
a |
+ |
Λc2a |
(1.13.6) |
||||||
3 |
c2 |
|
3 |
|
|
||||||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
a˙2 + Kc2 |
|
8πG |
ρa2 + |
Λc2a2 |
|
||||||||
= |
|
|
|
|
|
, |
|
(1.13.7) |
|||||
3 |
|
|
3 |
|
|
respectively. If we ignore the original terms in p and ρ we can see that Equation (1.13.7) can be written in a form similar to Equation (1.13.4):
|
a˙ |
2 |
Λc2 |
|
Λ |
|
|
|
Kc2 |
|
|||
|
− |
|
= H02 |
1 − |
|
= H02(1 |
− Ω0Λ) = − |
|
. |
(1.13.8) |
|||
a0 |
3 |
Λc |
a02 |
||||||||||
In this case the ‘critical’ value is |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Λc = |
3H02 |
|
|
|
|
||
|
|
|
|
|
|
|
, |
|
|
|
(1.13.9) |
||
|
|
|
|
|
|
c2 |
|
|
|
so that Ω0Λ = Λc2/3H02.
If we now reinstate the ‘ordinary’ matter we began with, we can see that the curvature is zero as long as Ω0Λ + Ω0 = 1. The cosmological constant therefore breaks the relationship between the matter density and curvature. Even if Ω0 < 1, a suitably chosen value of Ω0Λ = 1 − Ω0 can be invoked to ensure flat space sections.

Friedmann Models |
31 |
Bibliographic Notes on Chapter 1
The classic papers of Einstein (1917), de Sitter (1917), Friedmann (1922) and Lemaître (1927) are all well worth reading for historical insights. A particularly erudite overview of the role of observation in expanding world models is given by Sandage (1988). More detailed discussions of the basic background, including the role of general relativity in cosmology, can be found in Berry (1989), Harrison (1981), Kenyon (1990), Landau and Lifshitz (1975), Milne (1935), Misner et al. (1972), Narlikar (1993), Peebles (1993), Peacock (1999), Raychaudhuri (1979), Roos (1994), Wald (1984), Weinberg (1972) and Zel’dovich and Novikov (1983).
Problems
1.Suppose that it is discovered that Newton’s law of gravitation is incorrect, and that the force F on a test particle of mass m due to a body of mass M has an additional term that does not depend on M and is proportional to the separation r:
F = −GMmr2 + Amr3 .
Assuming that Newton’s sphere theorem continues to hold, derive the appropriate form of the Friedmann equation in this case and comment on your result.
2. The most general form of a space–time four-metric in the synchronous gauge is
ds2 = c2 dt2 − gαβ dxα dxβ = c2 dt2 − dl2,
where gαβ is the three-metric of the spatial hypersurfaces. By writing the equation of the three-space as that of a constrained surface in four dimensions, show that the most general form of the three-metric compatible with homogeneity and isotropy is given by the Robertson–Walker form.
3. Show that the special-relativistic formula for the Doppler shift,
1 + z = |
|
|
|
|
1 |
+ v/c , |
|||
|
|
1 |
|
|
− v/c
reduces to z v/c in the limit of small velocities. Invert the formula to give v/c in terms of z. Calculate the recession velocity of a galaxy at z = 5 using the specialrelativistic formula.
4. A model is constructed with Ω0 < 1, Λ ≠ 0 and k = 0. Show in this case that
q0 = 32 Ω0 − 1.
5. An object has luminosity distance dL and angular-diameter distance dA. Show that
dA |
= |
1 |
, |
dL |
(1 + z)2 |
independent of cosmology.
2
The Friedmann
Models
2.1 Perfect Fluid Models
In this chapter we shall consider a set of homogeneous and isotropic model universes that contain a relatively simple form of matter. In Section 1.13 we explained how a perfect fluid, described by an energy–momentum tensor of the type (1.10.2), forms the basis of the so-called Friedmann models. The ideal perfect fluid is, in fact, quite a realistic approximation in many situations. For example, if the mean free path between particle collisions is much less than the scales of physical interest, then the fluid may be treated as perfect. It should also be noted that the form (1.10.2) is also required for compatibility with the Cosmological Principle: anisotropic pressure is not permitted. To say more about the cosmological solutions, however, we need to say more about the relationship between p and ρ. In other words we need to specify an equation of state.
As we mentioned in the last section of the previous chapter, we need to specify an equation of state for our fluid in the form p = p(ρ). In many cases of physical interest, the appropriate equation of state can be cast, either exactly or approximately, in the form
p = wρc2 = (Γ − 1)ρc2, |
(2.1.1) |
where the parameter w is a constant which lies in the range
0 w 1. |
(2.1.2) |
We do not use the parameter Γ = 1 + w further in this book, but we have defined it here as it is used by other authors. The allowed range of w given in (2.1.2) is

34 The Friedmann Models
often called the Zel’dovich interval. We shall restrict ourselves for the rest of this chapter to cosmological models containing a perfect fluid with equation of state satisfying this condition.
The case with w = 0 represents dust (pressureless material). This is also a good approximation to the behaviour of any form of non-relativistic fluid or gas. Of course, gas of particles at some temperature T does exert pressure but the typical thermal energy of a particle is approximately kBT (kB is the Boltzmann constant), whereas its rest mass is mpc2, usually very much larger. The relativistic e ect of pressure is usually therefore negligible. In more detail, an ideal gas of nonrelativistic particles of mass mp, temperature T, density ρm and adiabatic index γ exerts a pressure
|
kBT |
|
|
kBT |
|
ρc2 |
|
p = nkBT = |
|
ρmc2 |
= |
|
|
|
= w(T)ρc2, (2.1.3) |
mpc2 |
mpc2 |
1 + (kBT/((γ − 1)mpc2)) |
where ρc2 is the energy density; a non-relativistic gas has w(T) 1 according to Equation (2.1.3) and will therefore be well approximated by a fluid of dust.
At the other extreme, a fluid of non-degenerate, ultrarelativistic particles in thermal equilibrium has an equation of state of the type
p = 31 ρc2. |
(2.1.4) |
For instance, this is the case for a gas of photons. A fluid with an equation of state of the type (2.1.4) is usually called a radiative fluid, though it may comprise relativistic particles of any form.
It is interesting to note that the parameter w is also related to the adiabatic sound speed of the fluid
vs = |
∂p |
1/2 |
|
|
|
S |
, |
(2.1.5) |
|
∂ρ |
where S denotes the entropy. In a dust fluid v = 0 and a radiative fluid has
√ s
vs = c/ 3. Note that the case w > 1 is impossible, because it would imply that vs > c. If w < 0, then it is no longer related to the sound speed, which would have to be imaginary. These two cases form the limits in (2.1.2). There are, however, physically important situations in which matter behaves like a fluid with w < 0, as we shall see later.
We shall restrict ourselves to the case where w is constant in time. We shall also assume that normal matter, described by an equation of state of the form (2.1.3), can be taken to have w(T) 0. From Equations (2.1.1) and (1.13.3) we can easily
obtain the relation |
|
ρa3(1+w) = const. = ρ0wa03(1+w). |
(2.1.6) |
In this equation and hereafter we use the su x ‘0’ to denote a reference time, usually the present. In particular we have, for a dust universe (w = 0) or a matter universe described by (2.1.3),
ρa3 ≡ ρma3 = const. = ρ0ma03 |
(2.1.7) |
Perfect Fluid Models |
35 |
(which simply represents the conservation of mass), and for a radiative universe (w = 13 )
ρa4 ≡ ρra4 = const. = ρ0ra04. |
(2.1.8) |
If one replaces the expansion parameter a with the redshift z, one finds, for dust and non-relativistic matter,
ρm = ρ0m(1 + z)3, |
(2.1.9) |
and, for radiation and relativistic matter,
ρr = ρ0r(1 + z)4. |
(2.1.10) |
The di erence between (2.1.9) and (2.1.10) can be understood quite straightforwardly if one considers a comoving box containing, say, N particles. Let us assume that, as the box expands, particles are neither created nor destroyed. If the particles are non-relativistic (i.e. if the box contains ‘dust’), then the density simply decreases as the cube of the scale factor, equivalent to (2.1.9). On the other hand, if the particles are relativistic, then they behave like photons: not only is their number-density diluted by a factor a3, but also the wavelength of each particle is increased by a factor a resulting in a redshift z. Since the energy of the particles is inversely proportional to their wavelength the total energy must decrease as the fourth power of the scale factor.
Notice the peculiar case in which w = −1 in (2.1.6), which we demonstrated to be the perfect fluid equivalent of a cosmological constant. The energy density does not vary as the universe expands for this kind of fluid.
Models of the Universe made from fluids with −13 < w < 1 have the property that they possess a point in time where a vanishes and the density diverges. This instant is called the Big Bang singularity. To see how this singularity arises, let us rewrite Equation (1.13.4) of the previous chapter using (2.1.6). Introducing the density parameter
|
|
|
|
|
|
|
Ω0w = |
|
|
ρ0w |
|
|
(2.1.11) |
|||
|
|
|
|
|
|
|
|
|
ρ0c |
|
||||||
allows us to obtain the equation |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
a˙ |
2 |
|
|
|
|
|
a0 |
|
|
1+3w |
+ (1 − Ω0w) |
|
|||
|
|
|
= H02 Ω0w |
|
|
|
|
|
|
(2.1.12) |
||||||
a0 |
a |
|
|
|||||||||||||
or, alternatively, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a0 |
2 |
|
|
|
a0 |
1+3w |
|
||||
H2(t) = H02 |
|
|
|
Ω0w |
|
|
+ (1 − Ω0w) , |
(2.1.13) |
||||||||
a |
a |
where H(t) = a/a˙ is the Hubble parameter at a generic time t. Suppose at some generic time, t (for example the present time, t0), the universe is expanding, so that a(t)˙ > 0. From Equation (1.13.1), we can see that a¨ < 0 for all t, provided

36 The Friedmann Models
Figure 2.1 The concavity of a(t) ensures that, if a(t)˙ > 0 for some time t, then there must be a singularity a finite time in the past, i.e. a point when a = 0. It also ensures that the age of the Universe, t0, is less than the Hubble time, 1/H0.
(ρ +3p/c2) > 0 or, in other words, (1+3w) > 0 since ρ > 0. This establishes that the graph of a(t) is necessarily concave. One can see therefore that a(t) must be equal to zero at some finite time in the past, and we can label this time t = 0 (see Figure 2.1). Since a(0) = 0 at this point, the density ρ diverges, as does the Hubble expansion parameter. One can see also that, because a(t) is a concave function, the time between the singularity and the epoch t must always be less than the characteristic expansion time of the Universe, τH = 1/H = a/a˙.
The Big Bang singularity is unavoidable in all homogeneous and isotropic models containing fluids with equation-of-state parameter w > −13 , which includes the Zel’dovich interval (2.1.2). It can be avoided, for example, in models with a non-zero cosmological constant, or if the universe is dominated by ‘matter’ with an e ective equation-of-state parameter w < −13 . One might suspect that the singularity may simply be a consequence of the special symmetry of the Friedmann models, and that inhomogeneous and/or anisotropic models would not display such a feature. However, this is not the case, as was shown by the classic work of Hawking an Penrose. We shall return to the unavoidability of the Big Bang singularity later, in Chapter 6.
Note that the expansion of the universe described in the Big Bang model is not due in any way to the e ect of pressure, which always acts to decelerate the expansion, but is a result of the initial conditions describing a homogeneous and isotropic universe. Another type of initial condition compatible with the Cosmological Principle are those which lead to an isotropic collapse of the universe towards a singularity like a time-reversed Big Bang, often called a Big Crunch.
2.2 Flat Models
In this section we shall find the solution to Equation (2.1.12) appropriate to a flat universe, i.e. with Ωw = 1. When w = 0 this solution is known as the Einstein– de Sitter universe; we shall also give this name to solutions with other values of

Flat Models |
37 |
w ≠ 0. For Ωw = 1, Equation (2.1.12) becomes
|
a˙ |
2 |
|
|
a0 |
1+3w |
|
|
|||
|
|
= H02 |
|
|
|
|
|
|
= H02(1 + z)1+3w, |
(2.2.1) |
|
a0 |
a |
|
|
|
|||||||
which one can immediately integrate to obtain |
|
||||||||||
|
|
|
a(t) = a0 |
|
t |
|
2/3(1+w). |
(2.2.2) |
|||
|
|
|
t0 |
This equation shows that the expansion of an Einstein–de Sitter universe lasts an indefinite time into the future; Equation (2.2.2) is equivalent to the relation
t = t0(1 + z)−3(1+w)/2, |
(2.2.3) |
which relates cosmic time t to redshift z. From Equations (2.2.2), (2.2.3) and (2.1.6), we can derive
H ≡ |
a˙ |
= |
|
|
|
2 |
|
|
= H0 |
t0 |
|
= H0(1 + z)3(1+w)/2, |
(2.2.4 a) |
||||||||||||
|
3 |
( |
1 |
|
|
|
t |
||||||||||||||||||
|
a |
|
|
|
+ w)t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
aa¨ |
|
|
|
|
1 + 3w |
|
|
|
const |
|
|
|
|
|
|
(2.2.4 |
|
) |
|||||
q ≡ − a˙2 |
= |
|
|
|
|
= |
. = q0, |
|
|
|
b |
||||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
t0w ≡ t0 = |
|
|
|
|
|
2 |
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
(2.2.4 c) |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
3(1 + w)H0 |
|
|
|
|
|
|
|
|
|
||||||||||||||||
ρ = ρ0w |
t |
− |
2 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
(2.2.4 d) |
||||||||
= |
|
|
|
|
|
|
|
|
; |
|
|
|
|||||||||||||
|
|
|
|
|
|
||||||||||||||||||||
t0 |
6(1 + w)2πGt2 |
|
|
|
|||||||||||||||||||||
in the last expression we have made use of the relation |
|
|
|
||||||||||||||||||||||
2 |
|
|
2 |
|
|
|
3H02 |
|
|
|
2 |
|
|
|
2 |
1 |
|
|
|
|
|||||
ρ0wt0 ≡ ρ0ct0w |
= |
|
|
|
|
= |
|
. |
(2.2.5) |
||||||||||||||||
8πG |
3(1 + w)H0 |
|
6(1 + w)2πG |
Useful special cases of the above relationship are dust, or matter-dominated universes (w = 0),
a(t) = a0 |
|
t |
2/3 |
, |
(2.2.6 a) |
|||||
t0 |
||||||||||
t = t0(1 + z)−3/2, |
(2.2.6 b) |
|||||||||
H = |
2 |
= H0(1 + z)3/2, |
(2.2.6 c) |
|||||||
|
||||||||||
3t |
||||||||||
q0 |
= |
21 , |
|
|
|
|
|
|
|
(2.2.6 d) |
t0m |
≡ t0 = |
|
2 |
, |
|
(2.2.6 e) |
||||
|
|
|
||||||||
3H0 |
|
|||||||||
ρm = |
|
1 |
|
|
; |
|
(2.2.6 f) |
|||
|
|
|||||||||
6πGt2 |
|