1 Simple 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.1 Flat 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:
3 a ˙ 2 a 2 = 8 π G ρ m = 8 π G ρ m 0 ( a 0 a ) 3 ⇔ ( a ˙ ) 2 = 8 π G ρ m 0 a 0 3 / 3 a = H 0 2 Ω m 0 a 0 3 a \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*} 3 a 2 a ˙ 2 = 8 π G ρ m = 8 π G ρ m 0 ( a a 0 ) 3 ⇔ ( a ˙ ) 2 = 8 π G ρ m 0 a 0 3 /3 a = H 0 2 Ω m 0 a a 0 3 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 Ω m 0 = 1 \Omega_m^0=1 Ω m 0 = 1 . Now let’s integrate the differential equation:
a a 0 d a a 0 = H 0 d t ⇒ ( a a 0 ) 3 / 2 = 3 2 H 0 t
\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 0 a a 0 d a = H 0 d t ⇒ ( a 0 a ) 3/2 = 2 3 H 0 t
⇒ a ( t ) a 0 = ( 3 2 H 0 t ) 2 / 3 \Rightarrow \frac{a(t)}{a_0} = \left( \frac{3}{2}H_0 t\right)^{2/3} ⇒ a 0 a ( t ) = ( 2 3 H 0 t ) 2/3 with at the beginning of the Universe t = 0 t=0 t = 0 and a = 0 a=0 a = 0 . We have thus obtained a direct relationship between the scale factor and the age of the Universe.
1.2 Flat 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 ) a 0 = ( 2 H 0 t ) 1 / 2 \frac{a(t)}{a_0} = \left( 2 H_0 t\right)^{1/2} a 0 a ( t ) = ( 2 H 0 t ) 1/2 1.3 Empty 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 :
Ω k 0 = 1 − Ω m 0 − Ω r 0 − Ω Λ 0 ≈ 1 \Omega_k^0 = 1 - \Omega_m^0 - \Omega_r^0 - \Omega_\Lambda^0 \approx 1 Ω k 0 = 1 − Ω m 0 − Ω r 0 − Ω Λ 0 ≈ 1 .
This Universe is therefore hyperbolic [1] . The first Friedmann equation is written :
a ˙ 2 a 2 = H 0 2 Ω k 0 a 0 2 a 2
\frac{\dot{a}^2}{a^2} = H_0^2 \Omega_k^0 \frac{a_0^2}{a^2} a 2 a ˙ 2 = H 0 2 Ω k 0 a 2 a 0 2 then :
a ˙ = a 0 2 H 0 2 Ω k 0 = a 0 H 0
\dot a = \sqrt{a_0^2 H_0^2 \Omega_k^0} = a_0 H_0 a ˙ = a 0 2 H 0 2 Ω k 0 = a 0 H 0 Integration gives a Universe expanding at constant speed:
a ( t ) = a 0 H 0 t a(t) = a_0 H_0 t a ( t ) = a 0 H 0 t 2 Historical 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.1 Einstein’s first model (1917) ¶ For a universe with only non-relativistic matter, P m = 0 P_m=0 P m = 0 . If it is static, then a ¨ = a ˙ = 0 \ddot{a}=\dot{a}=0 a ¨ = a ˙ = 0 and Friedmann’s second equation is written : c 2 k a 2 = c 2 Λ \frac{c^2 k}{a^2}=c^2 \Lambda a 2 c 2 k = c 2 Λ Injecting this equation into the first Friedmann equation gives :
3 c 2 k a 2 = 8 π G ρ m + c 2 Λ ⇔ 2 c 2 Λ = 8 π G ρ m ⇔ Λ = 4 π G ρ m / c 2 > 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} 3 a 2 c 2 k = 8 π G ρ m + c 2 Λ ⇔ 2 c 2 Λ = 8 π G ρ m ⇔ Λ = 4 π G ρ m / c 2 > 0 So, to obtain a static universe, Einstein had to introduce a non-zero cosmological constant. Moreover, we deduce that c 2 k / a 2 > 0 c^2 k / a^2 > 0 c 2 k / a 2 > 0 so k = + 1 k=+1 k = + 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 :
R ≡ a E = Λ − 1 / 2 R\equiv a_E = \Lambda^{-1/2} R ≡ a E = Λ − 1/2 Combining the two Friedmann equations for a universe with cold matter only, we find: 2 a ¨ a = − 8 π G ρ 3 + 2 c 2 Λ 3 \frac{2 \ddot{a}}{a} = - \frac{8\pi G \rho}{3} + \frac{2 c^2 \Lambda}{3} a 2 a ¨ = − 3 8 π Gρ + 3 2 c 2 Λ Let’s consider a perturbation of the radius, a = Λ − 1 / 2 + δ a a=\Lambda^{-1/2} + \delta a a = Λ − 1/2 + δ a , which induces a perturbation of the density of matter on the sphere ρ m = ρ 0 + δ ρ = c 2 Λ / 4 π G + δ ρ \rho_m = \rho_0 + \delta \rho = c^2 \Lambda/4 \pi G + \delta \rho ρ m = ρ 0 + δ ρ = c 2 Λ/4 π G + δ ρ . This last equation becomes
2 δ a ¨ a = − 2 c 2 Λ 3 − 8 π G δ ρ 3 c 2 + 2 c 2 Λ 3 = − 8 π G δ ρ 3 c 2 \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} 2 a δ a ¨ = − 3 2 c 2 Λ − 3 c 2 8 π G δ ρ + 3 2 c 2 Λ = − 3 c 2 8 π G δ ρ Now the quantity of matter is conserved, so we also have :
ρ a 3 = c s t e ⇒ δ ρ ρ = − 3 δ a a \rho a^3 = cste \Rightarrow \frac{\delta \rho}{\rho} = -3 \frac{\delta a}{a} ρ a 3 = cs t e ⇒ ρ δ ρ = − 3 a δ a and :
δ a ¨ = 4 π G ρ 0 δ a ⇔ δ a ¨ − c 2 Λ δ a = 0 \delta \ddot{a} = 4 \pi G \rho_0 \delta a \Leftrightarrow \delta \ddot{a} - c^2 \Lambda \delta a = 0 δ a ¨ = 4 π G ρ 0 δ a ⇔ δ a ¨ − c 2 Λ δ a = 0 i.e. the perturbation δ a \delta a δ a grows exponentially with time (Λ > 0 \Lambda>0 Λ > 0 ). The universe is unstable.
2.2 Universe with only one cosmological constant (de Sitter, 1917) ¶ In a flat universe k = 0 k=0 k = 0 and the first Friedmann equation gives in Sitter’s model : a ˙ 2 a 2 = c 2 Λ 3 ⇔ − 3 c 2 Λ d a d t + 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 a 2 a ˙ 2 = 3 c 2 Λ ⇔ − c 2 Λ 3 d t d a + a = 0 The solution of this differential equation is
a ( t ) = a 0 e s q r t c 2 Λ / 3 t a(t) = a_0 e^{sqrt{c^2 \Lambda/3}t} a ( t ) = a 0 e s q r t c 2 Λ/3 t The de Sitter universe grows exponentially with time.
To show that a de Sitter universe can be considered a static universe, let’s write T 0 = c 2 Λ / 3 T_0 = \sqrt{c^2 \Lambda / 3} T 0 = c 2 Λ/3 and change the coordinates r ′ ( t ) = a 0 e t / T 0 σ r'(t) = a_0 e^{t/T_0}\sigma r ′ ( t ) = a 0 e t / T 0 σ . The FLRW metric then becomes d s 2 = − d t 2 + a 2 ( d σ 2 + σ 2 d θ 2 + σ 2 sin θ d ϕ 2 ) = − d t 2 + [ ( d r ′ − r ′ d t / T 0 ) 2 + r ′ 2 d θ 2 + r ′ 2 sin θ 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] d s 2 = − d t 2 + a 2 ( d σ 2 + σ 2 d θ 2 + σ 2 sin θ d ϕ 2 ) = − d t 2 + [ ( d r ′ − r ′ d t / T 0 ) 2 + r ′2 d θ 2 + r ′2 sin θ d ϕ 2 ] Using the transformation :
t = t ′ + T 0 2 log ( r ′ 2 T 0 2 − 1 ) t = t' + \frac{T_0}{2}\log \left( \frac{r'^2}{T_0^2} - 1 \right) t = t ′ + 2 T 0 log ( T 0 2 r ′2 − 1 ) we find :
d s 2 = ( 1 − r ′ 2 T 0 2 ) d t ′ 2 − ( 1 − r ′ 2 T 0 2 ) − 1 d r ′ 2 − r ′ 2 d θ 2 − r ′ 2 sin θ 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 d s 2 = ( 1 − T 0 2 r ′2 ) d t ′2 − ( 1 − T 0 2 r ′2 ) − 1 d r ′2 − r ′2 d θ 2 − r ′2 sin θ d ϕ 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.3 Eddington-Lemaître model (1927) ¶ Exercise 3 (Lemaître model)
Friedmann and Lemaître were the first cosmologists to independently propose non-static models of universes with arbitrary curvatures. Friedmann’s equations have been extensively studied in this course, but Lemaître was the first to propose the idea that the Universe developed from a primitive atom. His model is based on a universe composed entirely of cold matter, with a cosmological constant and arbitrary spatial curvature (no radiation).
In such a model, show that, just after a big bang at t = 0 t=0 t = 0 , at the beginning of the universe the scale factor increases as : a ( t ) a 0 = ( 3 2 H 0 Ω m 0 t ) 2 / 3 \frac{a(t)}{a_0} =\left( \frac{3}{2}H_0\sqrt{\Omega_m^0}t\right)^{2/3} a 0 a ( t ) = ( 2 3 H 0 Ω m 0 t ) 2/3 As the universe expands, however, the energy density of matter decreases and the cosmological constant eventually dominates. Show that, for large t t t , the scale factor increases as : a ( t ) ∝ e H 0 Ω Λ 0 t a(t) \propto e^{H_0\sqrt{\Omega_\Lambda^0}t} a ( t ) ∝ e H 0 Ω Λ 0 t Calculate a ¨ \ddot{a} a ¨ and show that the expansion of the universe first decelerates, then accelerates. Calculate a ∗ a_* a ∗ the scale factor at the transition. In the single-matter Lemaître model, the first Friedmann equation is written : a ˙ 2 a 2 = H 0 2 [ Ω m 0 ( a 0 a ) 3 + Ω Λ 0 + Ω k 0 ( a 0 a ) 2 ] ⇔ a ˙ 2 = H 0 2 [ Ω m 0 a 0 3 a + Ω Λ 0 a 2 + Ω k 0 a 0 2 ] \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] a 2 a ˙ 2 = H 0 2 [ Ω m 0 ( a a 0 ) 3 + Ω Λ 0 + Ω k 0 ( a a 0 ) 2 ] ⇔ a ˙ 2 = H 0 2 [ Ω m 0 a a 0 3 + Ω Λ 0 a 2 + Ω k 0 a 0 2 ] At t ≈ 0 t\approx 0 t ≈ 0 , the Universe was extremely small so the matter term dominates:
a ˙ 2 ≈ H 0 2 [ Ω m 0 a 0 3 a ] ⇔ a a 0 a ˙ a 0 = H 0 Ω m 0 ⇔ a ( t ) a 0 = ( 3 2 H 0 Ω m 0 t ) 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} a ˙ 2 ≈ H 0 2 [ Ω m 0 a a 0 3 ] ⇔ a 0 a a 0 a ˙ = H 0 Ω m 0 ⇔ a 0 a ( t ) = ( 2 3 H 0 Ω m 0 t ) 2/3 Then, after some time, a a a becomes large and the cosmological constant term dominates: a ˙ 2 ≈ H 0 2 ( Ω Λ 0 a 2 ) ⇔ a ˙ = H 0 Ω m 0 a ( t ) ⇒ a ( t ) ∝ e H 0 Ω Λ 0 t \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} a ˙ 2 ≈ H 0 2 ( Ω Λ 0 a 2 ) ⇔ a ˙ = H 0 Ω m 0 a ( t ) ⇒ a ( t ) ∝ e H 0 Ω Λ 0 t Deriving equation (24) , we find that : 2 a ˙ a ¨ = H 0 2 [ − a ˙ Ω m 0 a 0 3 a 2 + 2 a ˙ a Ω Λ 0 ] ⇔ a ¨ a 0 = H 0 2 2 [ 2 Ω Λ 0 ( a a 0 ) − Ω m 0 ( a 0 a ) 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] 2 a ˙ a ¨ = H 0 2 [ − a ˙ Ω m 0 a 2 a 0 3 + 2 a ˙ a Ω Λ 0 ] ⇔ a 0 a ¨ = 2 H 0 2 [ 2 Ω Λ 0 ( a 0 a ) − Ω m 0 ( a a 0 ) 2 ] When a a a is small, we see that a ¨ \ddot{a} a ¨ is negative and expansion decelerates. However, when a a a is large, a ¨ > 0 \ddot{a}>0 a ¨ > 0 and the expansion of the universe accelerates. The transition occurs at :
a ¨ = 0 ⇔ 0 = H 0 2 2 [ 2 Ω Λ 0 a ∗ a 0 − Ω m 0 a 0 2 a ∗ 2 ] ⇔ a ∗ a 0 = ( Ω m 0 2 Ω Λ 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} a ¨ = 0 ⇔ 0 = 2 H 0 2 [ 2 Ω Λ 0 a 0 a ∗ − a ∗ 2 Ω m 0 a 0 2 ] ⇔ a 0 a ∗ = ( 2 Ω Λ 0 Ω m 0 ) 1/3 For the ΛCDM model with Ω m 0 ≈ 0.3 \Omega_m^0\approx 0.3 Ω m 0 ≈ 0.3 and Ω Λ 0 ≈ 0.7 \Omega_\Lambda^0\approx 0.7 Ω Λ 0 ≈ 0.7 , we have a ∗ / a 0 ≈ 0.6 a_*/a_0 \approx 0.6 a ∗ / a 0 ≈ 0.6 hence a transition redshift to z ≈ 0.67 z\approx 0.67 z ≈ 0.67 .
3 ΛCDM model ¶ The expansion of the Universe is now well described by the flat ΛCDM model (Ω k 0 = 0 \Omega_k^0=0 Ω k 0 = 0 ). The proportions of each of these components are now estimated at Planck2018 :
O m e g a Λ 0 = 0.685 , Ω m 0 = 315 Omega_\Lambda^0 = 0.685,\quad \Omega_m^0=315 O m e g a Λ 0 = 0.685 , Ω m 0 = 315 Cold matter can be separated into two contributions: dark matter Ω c 0 = 0.264 \Omega_{c}^0=0.264 Ω c 0 = 0.264 and baryonic matter [^baryons] Ω b 0 = 0.049 \Omega_b^0=0.049 Ω b 0 = 0.049 .
4 Mechanical analogy ¶ Exercise 4 (Mechanical analogy)
Write the first Friedmann equation as follows: 1 2 Ω k 0 = f ( Ω i 0 , a ) \frac{1}{2}\Omega_k^0 = f(\Omega_i^0,a) 2 1 Ω k 0 = f ( Ω i 0 , a ) Interpret this equation by analogy with the equation for the conservation of mechanical energy of a massive body in one-dimensional motion, and describe the role of each “potential energy” term.
Derive this equation with respect to time and make the analogy with Newton’s law applied to a massive body in one-dimensional motion. Analyze again the role of each “force” term. In what follows, we neglect the radiation component. Plot the potential energy terms as a function of the scale factor a a a and describe the fate of the following universes.
Spherical model (k = + 1 k=+1 k = + 1 ) with Λ = 4 π G ρ m / c 2 \Lambda=4\pi \GN \rho_m / c^2 Λ = 4 π G ρ m / c 2 (show that Einstein’s static model is a special case of this model and is unstable).
Matter-only models with different signs of curvature (the Einstein-de Sitter model corresponds to the case of flat curvature).
ΛCDM models with different signs for the cosmological constant and arbitrary curvatures (calculate the transition scale factor a ∗ a_* a ∗ between expansion deceleration and acceleration).
In terms of Ω i 0 \Omega_i^0 Ω i 0 , the first Friedmann equation is written : H 2 = ( a ˙ a ) 2 = H 0 2 ( Ω m 0 a 3 + Ω r 0 a 4 + Ω Λ 0 + Ω k 0 a 2 ) 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) H 2 = ( a a ˙ ) 2 = H 0 2 ( a 3 Ω m 0 + a 4 Ω r 0 + Ω Λ 0 + a 2 Ω k 0 ) which gives
1 2 Ω k 0 = 1 2 a ˙ 2 H 0 2 − 1 2 Ω m 0 a − 1 2 Ω r 0 a 2 − 1 2 Ω Λ 0 a 2 \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 2 1 Ω k 0 = 2 1 H 0 2 a ˙ 2 − 2 1 a Ω m 0 − 2 1 a 2 Ω r 0 − 2 1 Ω Λ 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:
1 2 Ω k 0 \frac{1}{2}\Omega_k^0 2 1 Ω k 0 is constant with a a a can be identified as the conserved mechanical energy of the massive body.1 2 a ˙ 2 H 0 2 \frac{1}{2}\frac{\dot{a}^2}{H_0^2} 2 1 H 0 2 a ˙ 2 represents the kinetic energy of the solid body.1 2 Ω m 0 a \frac{1}{2}\frac{\Omega_m^0}{a} 2 1 a Ω m 0 resembles a gravitational potential centered around a = 0 a=0 a = 0 .1 2 Ω r 0 a 2 \frac{1}{2}\frac{\Omega_r^0}{a^2} 2 1 a 2 Ω r 0 is another type of attractive potential.1 2 Ω Λ 0 a 2 \frac{1}{2}\Omega_\Lambda^0 a^2 2 1 Ω Λ 0 a 2 is an inverted harmonic potential (repulsive) centered around a = 0 a=0 a = 0 .a ¨ H 0 2 = − 1 2 Ω m 0 a 2 − Ω r 0 a 3 + Ω Λ 0 a \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 H 0 2 a ¨ = − 2 1 a 2 Ω m 0 − a 3 Ω r 0 + Ω Λ 0 a This equation resembles Newton’s law applied to a massive body in one-dimensional motion. Let’s make the analogy:
a ¨ H 0 2 \frac{\ddot{a}}{H_0^2} H 0 2 a ¨ acceleration of the solid body1 2 Ω m 0 a 2 \frac{1}{2}\frac{\Omega_m^0}{a^2 } 2 1 a 2 Ω m 0 gravitational force (attractive)Ω Λ 0 a \Omega_\Lambda^0 a Ω Λ 0 a elastic repulsive forceLet’s define :
V e f f ( a ) = − 1 2 Ω m 0 a − 1 2 Ω Λ 0 a 2 V_{\rm eff}(a) = - \frac{1}{2}\frac{\Omega_m^0}{a} - \frac{1}{2}\Omega_\Lambda^0 a^2 V eff ( a ) = − 2 1 a Ω m 0 − 2 1 Ω Λ 0 a 2 In this universe model, we have Ω Λ 0 = Ω m 0 / 2 \Omega_\Lambda^0 = \Omega_m^0 / 2 Ω Λ 0 = Ω m 0 /2 and : V e f f ( a ) = − 1 2 Ω m 0 a − 1 4 Ω m 0 a 2 V_{\rm eff}(a) = - \frac{1}{2}\frac{\Omega_m^0}{a} - \frac{1}{4}\Omega_m^0 a^2 V eff ( a ) = − 2 1 a Ω m 0 − 4 1 Ω m 0 a 2 d V e f f d a = 0 ⇒ ( 1 a 2 − a ) Ω m 0 = 0 ⇒ a = 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)} d a d V eff = 0 ⇒ ( a 2 1 − a ) Ω m 0 = 0 ⇒ a = 1 (today) At a = 1 a=1 a = 1 or t = 0 t=0 t = 0 , the first Friedmann equation gives :
1 = Ω m 0 + Ω Λ 0 + Ω k 0 = 1 − 3 2 Ω m 0 1 = \Omega_m^0 + \Omega_\Lambda^0 + \Omega_k^0 = 1 - \frac{3}{2}\Omega_m^0 1 = Ω m 0 + Ω Λ 0 + Ω k 0 = 1 − 2 3 Ω m 0 The model is spherical so Ω k 0 = − k c 2 / H 0 2 < 0 \Omega_k^0 = -k c^2 / H_0^2 < 0 Ω k 0 = − k c 2 / H 0 2 < 0 with k = + 1 k=+1 k = + 1 which implies that Ω m 0 > 2 / 3 \Omega_m^0 > 2/3 Ω m 0 > 2/3 . In Einstein’s Universe, Ω m 0 = 1 \Omega_m^0=1 Ω m 0 = 1 .
Figure 1: Potential energies for a spherical universe with Ω m 0 = 1 \Omega_m^0=1 Ω m 0 = 1 .
According to figure Figure 1 , Einstein’s solution at a = a 0 a=a_0 a = a 0 is unstable.
Table 1: Potential energies in the case of matter-only models with different curvatures: (top left), Ω m 0 = 1.5 ⇒ k = + 1 \Omega_m^0=1.5\Rightarrow k=+1 Ω m 0 = 1.5 ⇒ k = + 1 (top right), Ω m 0 = 0.5 ⇒ k = − 1 \Omega_m^0=0.5\Rightarrow k=-1 Ω m 0 = 0.5 ⇒ k = − 1 (bottom)
Ω m 0 = 1 ⇒ k = 0 \Omega_m^0=1\Rightarrow k=0 Ω m 0 = 1 ⇒ k = 0 Ω m 0 = 1.5 ⇒ k = + 1 \Omega_m^0=1.5\Rightarrow k=+1 Ω m 0 = 1.5 ⇒ k = + 1 Ω m 0 = 0.5 ⇒ k = − 1 \Omega_m^0=0.5\Rightarrow k=-1 Ω m 0 = 0.5 ⇒ k = − 1
In these models, curvature is again given by :
1 = Ω m 0 + Ω k 0 ⇒ Ω k 0 = 1 − Ω m 0 ⇒ { k = + 1 if Ω m 0 > 1 k = 0 if Ω m 0 = 1 k = − 1 if Ω m 0 < 1 1 = \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. 1 = Ω m 0 + Ω k 0 ⇒ Ω k 0 = 1 − Ω m 0 ⇒ ⎩ ⎨ ⎧ k = + 1 k = 0 k = − 1 if Ω m 0 > 1 if Ω m 0 = 1 if Ω m 0 < 1 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 t → ∞ t\rightarrow \infty t → ∞ . An expanding hyperbolic universe also expands to infinity.
The transition scale factor is given by : d V e f f d a = 0 → a ∗ = ( Ω m 0 2 Ω Λ 0 ) 1 / 3 \frac{d V_{\rm eff} }{da}= 0 \rightarrow a_* = \left(\frac{\Omega_m^0}{2 \Omega_\Lambda^0}\right)^{1/3} d a d V eff = 0 → a ∗ = ( 2 Ω Λ 0 Ω m 0 ) 1/3 Table 2: Potential energies in the case of single-matter models with different curvatures: (top left), Ω m 0 = 1.5 ⇒ k = + 1 \Omega_m^0=1.5\Rightarrow k=+1 Ω m 0 = 1.5 ⇒ k = + 1 (top right), Ω m 0 = 0.5 ⇒ k = − 1 \Omega_m^0=0.5\Rightarrow k=-1 Ω m 0 = 0.5 ⇒ k = − 1 (bottom)
Ω m 0 = 0.3 , Ω Λ 0 = 0.7 \Omega_m^0=0.3, \Omega_\Lambda^0=0.7 Ω m 0 = 0.3 , Ω Λ 0 = 0.7 Ω m 0 = 0.3 , Ω Λ 0 = 1.5 \Omega_m^0=0.3, \Omega_\Lambda^0=1.5 Ω m 0 = 0.3 , Ω Λ 0 = 1.5 Ω m 0 = 0.3 , Ω Λ 0 = 0.5 \Omega_m^0=0.3, \Omega_\Lambda^0=0.5 Ω m 0 = 0.3 , Ω Λ 0 = 0.5 Ω m 0 = 0.3 , Ω Λ 0 = − 0.7 \Omega_m^0=0.3, \Omega_\Lambda^0=-0.7 Ω m 0 = 0.3 , Ω Λ 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.
5 Age of the Universe ¶ Exercise 5 (Age of the Universe)
Calculate a ( t ) a(t) a ( t ) and the age of the universe for three different models, assuming H 0 ≈ 70 km/s/Mpc H_0\approx 70\,\text{km/s/Mpc} H 0 ≈ 70 km/s/Mpc : Ω m 0 = 1 , Ω Λ 0 = 0 , Ω r 0 ≈ 0 \Omega_m^0=1, \Omega_\Lambda^0=0, \Omega_r^0 \approx 0 Ω m 0 = 1 , Ω Λ 0 = 0 , Ω r 0 ≈ 0 (flat Einstein-de Sitter model)Ω m 0 = Ω Λ 0 ≈ 0 , Ω r 0 ≈ 0 \Omega_m^0=\Omega_\Lambda^0\approx 0, \Omega_r^0 \approx 0 Ω m 0 = Ω Λ 0 ≈ 0 , Ω r 0 ≈ 0 (empty universe)Ω m 0 = Ω r 0 ≈ 0 , Ω Λ 0 = \Omega_m^0=\Omega_r^0\approx 0, \Omega_\Lambda^0 = Ω m 0 = Ω r 0 ≈ 0 , Ω Λ 0 = 1$ (flat model dominated by vacuum)We need 1 1\, 1 Mpc= 3 × 1 0 19 =3\times 10^{19}\, = 3 × 1 0 19 km.
Which universe is the oldest? For a ΛCDM model with 32 % 32\% 32% matter and 68 % 68\% 68% dark energy, cosmologists estimate the age of the universe at 13.8 Gyr 13.8\,\text{Gyr} 13.8 Gyr . This figure can be compared with the age of the oldest stars, which is measured at Hobson et al. , 2006 [p. 410] :
t first stars = 11.5 ± 1.3 Gyr t_{\text{first stars}} = 11.5\pm 1.3\,\text{Gyr} t first stars = 11.5 ± 1.3 Gyr Now assume that H 0 ≈ 500 H_0\approx 500\, H 0 ≈ 500 km/s/Mpc as originally estimated by Edwin Hubble.
Calculate again the age of the universe for the first case, which was the most reasonable model at the time (flat, full of matter...). Conclude. We use the first Friedmann equation : H 2 = ( 1 a d a d t ) 2 = H 0 2 [ Ω m 0 a 3 + Ω Λ 0 + ( 1 − Ω m 0 − Ω Λ 0 ) a 2 ] . 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]. H 2 = ( a 1 d t d a ) 2 = H 0 2 [ a 3 Ω m 0 + Ω Λ 0 + a 2 ( 1 − Ω m 0 − Ω Λ 0 ) ] . because at a ( t 0 ) = 1 a(t_0)=1 a ( t 0 ) = 1 we still have 1 = Ω k 0 + Ω m 0 + Ω Λ 0 1=\Omega_k^0 + \Omega_m^0 + \Omega_\Lambda^0 1 = Ω k 0 + Ω m 0 + Ω Λ 0 .
We have Ω m 0 = 1 , Ω Λ 0 = 0 , Ω r 0 ≈ 0 \Omega_m^0=1, \Omega_\Lambda^0=0, \Omega_r^0 \approx 0 Ω m 0 = 1 , Ω Λ 0 = 0 , Ω r 0 ≈ 0 then : d a d t = a H 0 Ω m 0 a 3 + Ω Λ 0 + ( 1 − Ω m 0 − Ω Λ 0 ) a 2 = a H 0 Ω m 0 a 3 = H 0 a − 1 / 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} d t d a = a H 0 a 3 Ω m 0 + Ω Λ 0 + a 2 ( 1 − Ω m 0 − Ω Λ 0 ) = a H 0 a 3 Ω m 0 = H 0 a − 1/2 With t = 0 t=0 t = 0 at the Big Bang,
a d a = H 0 d t ⇒ ∫ 0 a a ′ d a ′ = H 0 ∫ 0 t d t ′ ⇒ 2 3 a 3 / 2 = H 0 t \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 a d a = H 0 d t ⇒ ∫ 0 a a ′ d a ′ = H 0 ∫ 0 t d t ′ ⇒ 3 2 a 3/2 = H 0 t We thus obtain :
a ( t ) = ( 3 2 H 0 t ) 2 / 3 and t U = 2 3 1 H 0 = 9.1 Gyr a(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} a ( t ) = ( 2 3 H 0 t ) 2/3 and t U = 3 2 H 0 1 = 9.1 Gyr with :
1 H 0 = 1 70 km/s/Mpc = 4.3 × 1 0 17 s = 13.6 Gyr \frac{1}{H_0} = \frac{1}{70\,\text{km/s/Mpc}}= 4.3\times 10^{17}\,\text{s} = 13.6\,\text{Gyr} H 0 1 = 70 km/s/Mpc 1 = 4.3 × 1 0 17 s = 13.6 Gyr Ω m 0 = Ω Λ 0 ≈ 0 , Ω r 0 ≈ 0 \Omega_m^0=\Omega_\Lambda^0\approx 0,\Omega_r^0 \approx 0 Ω m 0 = Ω Λ 0 ≈ 0 , Ω r 0 ≈ 0 then :d a d t = a H 0 ( 1 − Ω m 0 − Ω Λ 0 ) a 2 = a H 0 1 a = H 0 \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 d t d a = a H 0 a 2 ( 1 − Ω m 0 − Ω Λ 0 ) = a H 0 a 1 = H 0 With t = 0 t=0 t = 0 at the Big Bang,
a ( t ) = H 0 t ⇒ t U = 1 H 0 = 13.6 Gyr a(t) = H_0 t \Rightarrow t_U = \frac{1}{H_0} = 13.6\,\text{Gyr} a ( t ) = H 0 t ⇒ t U = H 0 1 = 13.6 Gyr Ω m 0 = Ω r 0 ≈ 0 , Ω Λ 0 = 1 \Omega_m^0=\Omega_r^0\approx 0, \Omega_\Lambda^0 = 1 Ω m 0 = Ω r 0 ≈ 0 , Ω Λ 0 = 1 then:d a d t = a H 0 Ω Λ 0 \frac{da}{dt} = a H_0 \sqrt{\Omega_\Lambda^0} d t d a = a H 0 Ω Λ 0 With t = t 0 t=t_0 t = t 0 now,
d a a = H 0 Ω Λ 0 d t ⇒ ∫ a a 0 d a ′ a ′ = H 0 Ω Λ 0 ∫ t t 0 d t ′ ⇒ a ( t ) = a 0 e H 0 Ω Λ 0 ( t − t 0 ) \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)} a d a = H 0 Ω Λ 0 d t ⇒ ∫ a a 0 a ′ d a ′ = H 0 Ω Λ 0 ∫ t t 0 d t ′ ⇒ a ( t ) = a 0 e H 0 Ω Λ 0 ( t − t 0 ) We have a ( t ) ⇒ 0 a(t)\Rightarrow 0 a ( t ) ⇒ 0 when t → − ∞ t\rightarrow -\infty t → − ∞ so the age of the universe is infinite in Sitter’s model.
Einstein de-Sitter’s value is incompatible with measuring the age of the first stars with such a value of H 0 H_0 H 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.
1 H 0 = 1 500 km/s/Mpc = 4.3 × 1 0 17 s = 1.9 Gyr \frac{1}{H_0} = \frac{1}{500\,\text{km/s/Mpc}}= 4.3\times 10^{17}\,\text{s} = 1.9\,\text{Gyr} H 0 1 = 500 km/s/Mpc 1 = 4.3 × 1 0 17 s = 1.9 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 H 0 H_0 H 0 was corrected.
6 Evolution of cosmological parameters ¶ Figure 2: Evolution of cosmological parameters.
This will embed the output of a notebook cell