Skip to article frontmatterSkip to article content

1Simple models

Now that we have a model to describe the dynamics of the Universe, let’s calculate its evolution in a few simple cases for practice.

1.1Flat Universe, matter only

Let’s start with the case of a flat Universe, containing only non-relativistic matter. This is the so-called Einstein-Sitter model. It’s the simplest we could think of in 1930. The first Friedmann equation is written:

3a˙2a2=8πGρm=8πGρm0(a0a)3(a˙)2=8πGρm0a03/3a=H02Ωm0a03a \begin{align*} & 3 \frac{\dot{a}^2}{a^2} = 8 \pi \GN \rho_m = 8 \pi \GN \rho_m^0 \left(\frac{a_0}{a}\right)^{3} & \Leftrightarrow (\dot a)^2 = 8 \pi \GN \rho_m^0 a_0^3 / 3 a = H_0^2 \Omega_m^0 \frac{a_0^3}{a} \end{align*}

Before integrating this differential equation, let’s remember that the energy density parameters are linked by a closure relation (70). Therefore, in a flat Universe with only matter, we have Ωm0=1\Omega_m^0=1. Now let’s integrate the differential equation:

aa0daa0=H0dt(aa0)3/2=32H0t \sqrt{\frac{a}{a_0}} \frac{\dd a}{a_0} = H_0 \dd t \Rightarrow \left(\frac{a}{a_0}\right)^{3/2} = \frac{3}{2}H_0 t

a(t)a0=(32H0t)2/3\Rightarrow \frac{a(t)}{a_0} = \left( \frac{3}{2}H_0 t\right)^{2/3}

with at the beginning of the Universe t=0t=0 and a=0a=0. We have thus obtained a direct relationship between the scale factor and the age of the Universe.

1.2Flat Universe, radiation only

Using similar reasoning, we show that for a flat Universe dominated by radiation, we have a different evolution of the scale factor:

a(t)a0=(2H0t)1/2\frac{a(t)}{a_0} = \left( 2 H_0 t\right)^{1/2}

1.3Empty Universe (Milne)

Let’s assume that the Universe is empty, or at least with a total energy density very low compared to the critical density. Then the Universe must be curved, since in this case :

Ωk0=1Ωm0Ωr0ΩΛ01\Omega_k^0 = 1 - \Omega_m^0 - \Omega_r^0 - \Omega_\Lambda^0 \approx 1

. This Universe is therefore hyperbolic [1]. The first Friedmann equation is written :

a˙2a2=H02Ωk0a02a2 \frac{\dot{a}^2}{a^2} = H_0^2 \Omega_k^0 \frac{a_0^2}{a^2}

then :

a˙=a02H02Ωk0=a0H0 \dot a = \sqrt{a_0^2 H_0^2 \Omega_k^0} = a_0 H_0

Integration gives a Universe expanding at constant speed:

a(t)=a0H0ta(t) = a_0 H_0 t

2Historical models

Modern cosmology was born of General Relativity. Since the writing of these equations, scientists have begun to describe the universe mathematically as a physical system. Numerous models have been proposed to describe the different histories of the universe. In this chapter, we will review some of them.

2.1Einstein’s first model (1917)

Solution to Exercise 1
  1. For a universe with only non-relativistic matter, Pm=0P_m=0. If it is static, then a¨=a˙=0\ddot{a}=\dot{a}=0 and Friedmann’s second equation is written :
c2ka2=c2Λ\frac{c^2 k}{a^2}=c^2 \Lambda

Injecting this equation into the first Friedmann equation gives :

3c2ka2=8πGρm+c2Λ2c2Λ=8πGρm Λ=4πGρm/c2>0\begin{align} 3 \frac{c^2 k}{a^2} = 8 \pi G \rho_m + c^2 \Lambda & \Leftrightarrow 2 c^2 \Lambda = 8 \pi G \rho_m \ & \Leftrightarrow \Lambda = 4 \pi G \rho_m / c^2 > 0 \end{align}

So, to obtain a static universe, Einstein had to introduce a non-zero cosmological constant. Moreover, we deduce that c2k/a2>0c^2 k / a^2 > 0 so k=+1k=+1: the static universe must be spherical. Einstein originally argued that the universe should be bounded using Mach’s principle[2], so the solution of a spherical universe seemed satisfactory from this point of view.

The scale factor can be associated with the radius of the spherical universe. Using the second Friedmann equation, we find that the value of the radius is :

RaE=Λ1/2R\equiv a_E = \Lambda^{-1/2}
  1. Combining the two Friedmann equations for a universe with cold matter only, we find:
2a¨a=8πGρ3+2c2Λ3\frac{2 \ddot{a}}{a} = - \frac{8\pi G \rho}{3} + \frac{2 c^2 \Lambda}{3}

Let’s consider a perturbation of the radius, a=Λ1/2+δaa=\Lambda^{-1/2} + \delta a, which induces a perturbation of the density of matter on the sphere ρm=ρ0+δρ=c2Λ/4πG+δρ\rho_m = \rho_0 + \delta \rho = c^2 \Lambda/4 \pi G + \delta \rho. This last equation becomes

2δa¨a=2c2Λ38πGδρ3c2+2c2Λ3=8πGδρ3c2\begin{align} 2 \frac{\delta \ddot{a}}{a} = - \frac{2 c^2 \Lambda}{3} - \frac{8 \pi G \delta \rho}{3 c^2}+ \frac{2 c^2 \Lambda}{3} = - \frac{8 \pi G \delta \rho}{3 c^2} \end{align}

Now the quantity of matter is conserved, so we also have :

ρa3=csteδρρ=3δaa\rho a^3 = cste \Rightarrow \frac{\delta \rho}{\rho} = -3 \frac{\delta a}{a}

and :

δa¨=4πGρ0δaδa¨c2Λδa=0\delta \ddot{a} = 4 \pi G \rho_0 \delta a \Leftrightarrow \delta \ddot{a} - c^2 \Lambda \delta a = 0

i.e. the perturbation δa\delta a grows exponentially with time (Λ>0\Lambda>0). The universe is unstable.

2.2Universe with only one cosmological constant (de Sitter, 1917)

Solution to Exercise 2
  1. In a flat universe k=0k=0 and the first Friedmann equation gives in Sitter’s model :
a˙2a2=c2Λ33c2Λdadt+a=0\frac{\dot{a}^2}{a^2} = \frac{c^2 \Lambda}{3} \Leftrightarrow -\sqrt{\frac{3}{c^2 \Lambda}}\frac{\dd a }{\dd t}+a=0

The solution of this differential equation is

a(t)=a0esqrtc2Λ/3ta(t) = a_0 e^{sqrt{c^2 \Lambda/3}t}

The de Sitter universe grows exponentially with time.

  1. To show that a de Sitter universe can be considered a static universe, let’s write T0=c2Λ/3T_0 = \sqrt{c^2 \Lambda / 3} and change the coordinates r(t)=a0et/T0σr'(t) = a_0 e^{t/T_0}\sigma. The FLRW metric then becomes
ds2=dt2+a2(dσ2+σ2dθ2+σ2sinθdϕ2)=dt2+[(drrdt/T0)2+r2dθ2+r2sinθdϕ2]\dd s^2 =- \dd t^2 + a^2 (\dd \sigma^2 + \sigma^2 \dd \theta^2 + \sigma^2 \sin\theta \dd \phi^2) =- \dd t^2 + \left[\left(\dd r' - r' \dd t / T_0\right)^2 + r'^2 \dd \theta^2 + r'^2 \sin\theta \dd \phi^2\right]

Using the transformation :

t=t+T02log(r2T021)t = t' + \frac{T_0}{2}\log \left( \frac{r'^2}{T_0^2} - 1 \right)

we find :

ds2=(1r2T02)dt2(1r2T02)1dr2r2dθ2r2sinθdϕ2\dd s^2 = \left(1-\frac{r'^2}{T_0^2}\right) \dd t'^2 - \left(1-\frac{r'^2}{T_0^2}\right)^{-1} \dd r'^2 - r'^2 \dd \theta^2 -r'^2 \sin\theta \dd \phi^2

Historically, Sitter’s model was discovered as a static universe with this coordinate system, whose shape closely resembles Schwarzschild’s solution. As soon as the idea of an expanding universe was accepted by the scientific community, Sitter’s universe was considered in its dynamical form, as an exponentially expanding universe dominated by the cosmological constant.

2.3Eddington-Lemaître model (1927)

Solution to Exercise 3
  1. In the single-matter Lemaître model, the first Friedmann equation is written :
a˙2a2=H02[Ωm0(a0a)3+ΩΛ0+Ωk0(a0a)2]a˙2=H02[Ωm0a03a+ΩΛ0a2+Ωk0a02]\frac{\dot{a}^2}{a^2} = H_0^2\left[\Omega_m^0 \left(\frac{a_0}{a}\right)^{3} + \Omega_\Lambda^0 + \Omega_k^0 \left(\frac{a_0}{a}\right)^{2}\right] \Leftrightarrow \dot{a}^2 = H_0^2\left[\Omega_m^0 \frac{a_0^3}{a} + \Omega_\Lambda^0 a^2 + \Omega_k^0 a_0^2 \right]

At t0t\approx 0, the Universe was extremely small so the matter term dominates:

a˙2H02[Ωm0a03a]aa0a˙a0=H0Ωm0a(t)a0=(32H0Ωm0t)2/3\dot a^2 \approx H_0^2\left[\Omega_m^0 \frac{a_0^3}{a}\right] \Leftrightarrow \sqrt{\frac{a}{a_0}}\frac{\dot{a}}{a_0}= H_0 \sqrt{\Omega_m^0} \Leftrightarrow \frac{a(t)}{a_0} = \left(\frac{3}{2}H_0\sqrt{\Omega_m^0}t\right)^{2/3}
  1. Then, after some time, aa becomes large and the cosmological constant term dominates:
a˙2H02(ΩΛ0a2)a˙=H0Ωm0a(t)a(t)eH0ΩΛ0t\dot{a}^2 \approx H_0^2\left(\Omega_\Lambda^0 a^2\right) \Leftrightarrow \dot{a}= H_0 \sqrt{\Omega_m^0} a(t) \Rightarrow a(t) \propto e^{H_0\sqrt{\Omega_\Lambda^0}t}
  1. Deriving equation (24), we find that :
2a˙a¨=H02[a˙Ωm0a03a2+2a˙aΩΛ0]a¨a0=H022[2ΩΛ0(aa0)Ωm0(a0a)2]2\dot{a}\ddot{a} = H_0^2\left[-\dot{a}\Omega_m^0 \frac{a_0^3 }{a^{2}} + 2 \dot{a} a \Omega_\Lambda^0 \right] \Leftrightarrow \frac{\ddot{a}}{a_0} = \frac{H_0^2}{2}\left[2 \Omega_\Lambda^0 \left(\frac{a}{a_0}\right) - \Omega_m^0\left(\frac{a_0}{a}\right)^2\right]

When aa is small, we see that a¨\ddot{a} is negative and expansion decelerates. However, when aa is large, a¨>0\ddot{a}>0 and the expansion of the universe accelerates. The transition occurs at :

a¨=00=H022[2ΩΛ0aa0Ωm0a02a2]aa0=(Ωm02ΩΛ0)1/3\ddot{a}=0 \Leftrightarrow 0=\frac{H_0^2}{2}\left[2 \Omega_\Lambda^0 \frac{a_*}{a_0} - \frac{\Omega_m^0a_0^2}{a_*^2}\right] \Leftrightarrow \frac{a_*}{a_0} = \left( \frac{\Omega_m^0}{2\Omega_\Lambda^0}\right)^{1/3}

For the ΛCDM model with Ωm00.3\Omega_m^0\approx 0.3 and ΩΛ00.7\Omega_\Lambda^0\approx 0.7, we have a/a00.6a_*/a_0 \approx 0.6 hence a transition redshift to z0.67z\approx 0.67.

3ΛCDM model

The expansion of the Universe is now well described by the flat ΛCDM model (Ωk0=0\Omega_k^0=0). The proportions of each of these components are now estimated at Planck2018:

OmegaΛ0=0.685,Ωm0=315Omega_\Lambda^0 = 0.685,\quad \Omega_m^0=315

Cold matter can be separated into two contributions: dark matter Ωc0=0.264\Omega_{c}^0=0.264 and baryonic matter [^baryons] Ωb0=0.049\Omega_b^0=0.049.

4Mechanical analogy

Solution to Exercise 4
  1. In terms of Ωi0\Omega_i^0, the first Friedmann equation is written :
H2=(a˙a)2=H02(Ωm0a3+Ωr0a4+ΩΛ0+Ωk0a2)H^2 = \left(\frac{\dot{a}}{a}\right)^2 = H_0^2 \left( \frac{\Omega_m^0}{a^3} + \frac{\Omega_r^0}{a^4} + \Omega_\Lambda^0 + \frac{\Omega_k^0}{a^2} \right)

which gives

12Ωk0=12a˙2H0212Ωm0a12Ωr0a212ΩΛ0a2\frac{1}{2}\Omega_k^0 = \frac{1}{2}\frac{\dot{a}^2}{H_0^2} - \frac{1}{2}\frac{\Omega_m^0}{a} - \frac{1}{2}\frac{\Omega_r^0}{a^2} - \frac{1}{2}\Omega_\Lambda^0 a^2

This last equation resembles the equation for the conservation of mechanical energy for a massive body in one-dimensional motion. Let’s make the analogy:

  • 12Ωk0\frac{1}{2}\Omega_k^0 is constant with aa can be identified as the conserved mechanical energy of the massive body.
  • 12a˙2H02\frac{1}{2}\frac{\dot{a}^2}{H_0^2} represents the kinetic energy of the solid body.
  • 12Ωm0a\frac{1}{2}\frac{\Omega_m^0}{a} resembles a gravitational potential centered around a=0a=0.
  • 12Ωr0a2\frac{1}{2}\frac{\Omega_r^0}{a^2} is another type of attractive potential.
  • 12ΩΛ0a2\frac{1}{2}\Omega_\Lambda^0 a^2 is an inverted harmonic potential (repulsive) centered around a=0a=0.
  1. a¨H02=12Ωm0a2Ωr0a3+ΩΛ0a\frac{\ddot{a}}{H_0^2} = - \frac{1}{2}\frac{\Omega_m^0}{a^2 } -\frac{\Omega_r^0}{a^3 } + \Omega_\Lambda^0 a

This equation resembles Newton’s law applied to a massive body in one-dimensional motion. Let’s make the analogy:

  • a¨H02\frac{\ddot{a}}{H_0^2} acceleration of the solid body
  • 12Ωm0a2\frac{1}{2}\frac{\Omega_m^0}{a^2 } gravitational force (attractive)
  • ΩΛ0a\Omega_\Lambda^0 a elastic repulsive force

Let’s define :

Veff(a)=12Ωm0a12ΩΛ0a2V_{\rm eff}(a) = - \frac{1}{2}\frac{\Omega_m^0}{a} - \frac{1}{2}\Omega_\Lambda^0 a^2
  1. In this universe model, we have ΩΛ0=Ωm0/2\Omega_\Lambda^0 = \Omega_m^0 / 2 and :
Veff(a)=12Ωm0a14Ωm0a2V_{\rm eff}(a) = - \frac{1}{2}\frac{\Omega_m^0}{a} - \frac{1}{4}\Omega_m^0 a^2
dVeffda=0(1a2a)Ωm0=0a=1 (today)\frac{\dd V_{\rm eff} }{\dd a}= 0 \Rightarrow \left(\frac{1}{a^2}-a\right)\Omega_m^0 = 0 \Rightarrow a=1\text{ (today)}

At a=1a=1 or t=0t=0, the first Friedmann equation gives :

1=Ωm0+ΩΛ0+Ωk0=132Ωm01 = \Omega_m^0 + \Omega_\Lambda^0 + \Omega_k^0 = 1 - \frac{3}{2}\Omega_m^0

The model is spherical so Ωk0=kc2/H02<0\Omega_k^0 = -k c^2 / H_0^2 < 0 with k=+1k=+1 which implies that Ωm0>2/3\Omega_m^0 > 2/3. In Einstein’s Universe, Ωm0=1\Omega_m^0=1.

<Figure size 640x480 with 1 Axes>

Figure 1:Potential energies for a spherical universe with Ωm0=1\Omega_m^0=1.

According to figure Figure 1, Einstein’s solution at a=a0a=a_0 is unstable.

Table 1:Potential energies in the case of matter-only models with different curvatures: (top left), Ωm0=1.5k=+1\Omega_m^0=1.5\Rightarrow k=+1 (top right), Ωm0=0.5k=1\Omega_m^0=0.5\Rightarrow k=-1 (bottom)

Ωm0=1k=0\Omega_m^0=1\Rightarrow k=0Ωm0=1.5k=+1\Omega_m^0=1.5\Rightarrow k=+1Ωm0=0.5k=1\Omega_m^0=0.5\Rightarrow k=-1

In these models, curvature is again given by :

1=Ωm0+Ωk0Ωk0=1Ωm0{k=+1 if Ωm0>1k=0 if Ωm0=1k=1 if Ωm0<11 = \Omega_m^0 + \Omega_k^0 \Rightarrow \Omega_k^0 = 1 -\Omega_m^0 \Rightarrow \left\lbrace\begin{array}{ll} k=+1 & \text{ if } \Omega_m^0 > 1\\\\ k=0 & \text{ if } \Omega_m^0 = 1 \\\\ k=-1 & \text{ if } \Omega_m^0 < 1 \end{array}\right.

By analyzing the three plots in figure Table 1, we can say that a spherical universe composed solely of matter will necessarily collapse at some point, whatever its initial conditions (Einstein’s need to add the cosmological constant). An expanding flat universe expands indefinitely and asymptotically stops expanding at tt\rightarrow \infty. An expanding hyperbolic universe also expands to infinity.

  1. The transition scale factor is given by :
dVeffda=0a=(Ωm02ΩΛ0)1/3\frac{d V_{\rm eff} }{da}= 0 \rightarrow a_* = \left(\frac{\Omega_m^0}{2 \Omega_\Lambda^0}\right)^{1/3}

Table 2:Potential energies in the case of single-matter models with different curvatures: (top left), Ωm0=1.5k=+1\Omega_m^0=1.5\Rightarrow k=+1 (top right), Ωm0=0.5k=1\Omega_m^0=0.5\Rightarrow k=-1 (bottom)

Ωm0=0.3,ΩΛ0=0.7\Omega_m^0=0.3, \Omega_\Lambda^0=0.7Ωm0=0.3,ΩΛ0=1.5\Omega_m^0=0.3, \Omega_\Lambda^0=1.5
Ωm0=0.3,ΩΛ0=0.5\Omega_m^0=0.3, \Omega_\Lambda^0=0.5Ωm0=0.3,ΩΛ0=0.7\Omega_m^0=0.3, \Omega_\Lambda^0=-0.7

Depending on the parameter values, the transition scale occurs in the future or in the past. If the cosmological constant is positive, expanding universes have a decelerated expansion and, after the transition scale, an accelerated expansion. If the cosmological constant is negative, the universe must collapse after a certain time.

So why is the Universe expanding today? That depends entirely on the initial conditions, and in particular on whether the universe was born of a Big Bang. And why was there a Big Bang? We can let our imaginations run wild: brane collisions, God, pan-dimensional mice... but the answer is not (yet) given by the physical sciences.

5Age of the Universe

Solution to Exercise 5
  1. We use the first Friedmann equation :
H2=(1adadt)2=H02[Ωm0a3+ΩΛ0+(1Ωm0ΩΛ0)a2].H^2 = \left(\frac{1}{a}\frac{da}{dt}\right)^2 = H_0^2 \left[ \frac{\Omega_m^0}{a^3}+ \Omega_\Lambda^0+ \frac{\left(1-\Omega_m^0-\Omega_\Lambda^0\right)}{a^2}\right].

because at a(t0)=1a(t_0)=1 we still have 1=Ωk0+Ωm0+ΩΛ01=\Omega_k^0 + \Omega_m^0 + \Omega_\Lambda^0.

  • We have Ωm0=1,ΩΛ0=0,Ωr00\Omega_m^0=1, \Omega_\Lambda^0=0, \Omega_r^0 \approx 0 then :
dadt=aH0Ωm0a3+ΩΛ0+(1Ωm0ΩΛ0)a2=aH0Ωm0a3=H0a1/2\frac{\dd a}{\dd t} = a H_0 \sqrt{\frac{\Omega_m^0}{a^3}+\Omega_\Lambda^0+ \frac{\left(1-\Omega_m^0-\Omega_\Lambda^0\right)}{a^2}} = a H_0 \sqrt{\frac{\Omega_m^0}{a^3}} = H_0 a^{-1/2}

With t=0t=0 at the Big Bang,

ada=H0dt0aada=H00tdt23a3/2=H0t\sqrt{a}\dd a = H_0 \dd t \Rightarrow \int_0^a \sqrt{a'}\dd a' = H_0 \int_0^t \dd t' \Rightarrow \frac{2}{3}a^{3/2} = H_0 t

We thus obtain :

a(t)=(32H0t)2/3andtU=231H0=9.1Gyra(t) = \left( \frac{3}{2}H_0 t\right)^{2/3}\quad\text{and}\quad t_U = \frac{2}{3}\frac{1}{H_0} = 9.1\,\text{Gyr}

with :

1H0=170km/s/Mpc=4.3×1017s=13.6Gyr\frac{1}{H_0} = \frac{1}{70\,\text{km/s/Mpc}}= 4.3\times 10^{17}\,\text{s} = 13.6\,\text{Gyr}
  • Ωm0=ΩΛ00,Ωr00\Omega_m^0=\Omega_\Lambda^0\approx 0,\Omega_r^0 \approx 0 then :
dadt=aH0(1Ωm0ΩΛ0)a2=aH01a=H0\frac{\dd a}{\dd t} = a H_0 \sqrt{\frac{\left(1-\Omega_m^0-\Omega_\Lambda^0\right)}{a^2}} = a H_0 \frac{\sqrt{1}}{a} = H_0

With t=0t=0 at the Big Bang,

a(t)=H0ttU=1H0=13.6Gyra(t) = H_0 t \Rightarrow t_U = \frac{1}{H_0} = 13.6\,\text{Gyr}
  • Ωm0=Ωr00,ΩΛ0=1\Omega_m^0=\Omega_r^0\approx 0, \Omega_\Lambda^0 = 1 then:
dadt=aH0ΩΛ0\frac{da}{dt} = a H_0 \sqrt{\Omega_\Lambda^0}

With t=t0t=t_0 now,

daa=H0ΩΛ0dtaa0daa=H0ΩΛ0tt0dta(t)=a0eH0ΩΛ0(tt0)\frac{\dd a}{a} = H_0 \sqrt{\Omega_\Lambda^0} \dd t \Rightarrow \int_a^{a_0} \frac{\dd a'}{a'} = H_0 \sqrt{\Omega_\Lambda^0} \int_t^{t_0} \dd t' \Rightarrow a(t) = a_0 e^{H_0 \sqrt{\Omega_\Lambda^0} (t- t_0)}

We have a(t)0a(t)\Rightarrow 0 when tt\rightarrow -\infty so the age of the universe is infinite in Sitter’s model.

  1. Einstein de-Sitter’s value is incompatible with measuring the age of the first stars with such a value of H0H_0. However, the ΛCDM model agrees with the “empty universe” model. However, in the “empty universe”, the assumption that matter plays no role is very strong.

1H0=1500km/s/Mpc=4.3×1017s=1.9Gyr\frac{1}{H_0} = \frac{1}{500\,\text{km/s/Mpc}}= 4.3\times 10^{17}\,\text{s} = 1.9\,\text{Gyr}

The first measurement of the Hubble constant was incompatible with measurements of the Earth’s age. Early cosmological models were therefore affected by this fact until the measurement of the constant H0H_0 was corrected.

6Evolution of cosmological parameters

Loading...

Figure 2:Evolution of cosmological parameters.

This will embed the output of a notebook cell

Footnotes
  1. The sign of kk is the inverse of the sign of Ωk0\Omega_k^0.

  2. In theoretical physics, particularly in discussions of gravitational theories, the Mach principle is the name given by Einstein to an imprecise hypothesis often attributed to physicist and philosopher Ernst Mach. The idea is that local inertial frames are determined by the large-scale distribution of matter.

References
  1. Hobson, M. P., Efstathiou, G. P., & Lasenby, A. N. (2006). General Relativity: An Introduction for Physicists. Cambridge University Press. https://books.google.fr/books?id=5dryXCWR7EIC