2.II.15D

Cosmology | Part II, 2005

(a) Consider a homogeneous and isotropic universe with scale factor a(t)a(t) and filled with mass density ρ(t)\rho(t). Show how the conservation of kinetic energy plus gravitational potential energy for a test particle on the edge of a spherical region in this universe can be used to derive the Friedmann equation

(a˙a)2+kc2a2=8πG3ρ(*)\tag{*} \left(\frac{\dot{a}}{a}\right)^{2}+\frac{k c^{2}}{a^{2}}=\frac{8 \pi G}{3} \rho

where kk is a constant. State clearly any assumptions you have made.

(b) Now suppose that the universe was filled throughout its history with radiation with equation of state P=ρc2/3P=\rho c^{2} / 3. Using the fluid conservation equation and the definition of the relative density Ω\Omega, show that the density of this radiation can be expressed as

ρ=3H028πGΩ0a4,\rho=\frac{3 H_{0}^{2}}{8 \pi G} \frac{\Omega_{0}}{a^{4}},

where H0H_{0} is the Hubble parameter today and Ω0\Omega_{0} is the relative density today (t=t0)\left(t=t_{0}\right) and a0a(t0)=1a_{0} \equiv a\left(t_{0}\right)=1 is assumed. Show also that kc2=H02(Ω01)k c^{2}=H_{0}^{2}\left(\Omega_{0}-1\right) and hence rewrite the Friedmann equation ()(*) as

(a˙a)2=H02Ω0(1a4βa2)(†)\tag{†} \left(\frac{\dot{a}}{a}\right)^{2}=H_{0}^{2} \Omega_{0}\left(\frac{1}{a^{4}}-\frac{\beta}{a^{2}}\right)

where β(Ω01)/Ω0\beta \equiv\left(\Omega_{0}-1\right) / \Omega_{0}.

(c) Now consider a closed model with k>0k>0 (or Ω>1)\Omega>1). Rewrite ( \dagger ) using the new time variable τ\tau defined by

dtdτ=a\frac{d t}{d \tau}=a

Hence, or otherwise, solve ()(†) to find the parametric solution

a(τ)=1β(sinατ),t(τ)=1αβ(1cosατ),a(\tau)=\frac{1}{\sqrt{\beta}}(\sin \alpha \tau), \quad t(\tau)=\frac{1}{\alpha \sqrt{\beta}}(1-\cos \alpha \tau),

where αH0(Ω01).\alpha \equiv H_{0} \sqrt{\left(\Omega_{0}-1\right)} . \quad Recall that dx/1x2=sin1x.]\left.\int d x / \sqrt{1-x^{2}}=\sin ^{-1} x .\right]

Using the solution for a(τ)a(\tau), find the value of the new time variable τ=τ0\tau=\tau_{0} today and hence deduce that the age of the universe in this model is

t0=H01Ω01Ω01t_{0}=H_{0}^{-1} \frac{\sqrt{\Omega_{0}}-1}{\Omega_{0}-1}

Typos? Please submit corrections to this page on GitHub.