Skip to article frontmatterSkip to article content

1Cosmological principle

Figure 1:Galaxy distribution compiled by the eBOSS survey. Every dot in this “pie” diagram is a galaxy, color coded by type: green for nearby galaxies, magenta and red for old red galaxies, blue for young blue galaxies, yellow and white for quasars. Credit: A. Raichoor (EPFL) / A. Ross (Ohio State Univ.) / SDSS Collaboration

Temperature map of the Cosmic Microwave Background (CMB) observed by the Planck satellite. The relative difference observed between the  temperature of hot (red) and cold (blue) spots relative to the mean the mean is of the order of \delta \theta / \theta \approx 10^{-5}.

Figure 2:Temperature map of the Cosmic Microwave Background (CMB) observed by the Planck satellite. The relative difference observed between the temperature of hot (red) and cold (blue) spots relative to the mean the mean is of the order of δθ/θ105\delta \theta / \theta \approx 10^{-5}.

To be able to build a model of the Universe, i.e. a theoretical construct capable of describing the contents of the Universe and its evolution, we need to be able to Einstein’s equation of General Relativity. But how much detail is needed to describe the Universe sufficiently well on large scales? Parametrizing Einstein’s equation to include the scale of the solar system is both illusory and unnecessary. What is the the structure of the Universe on larger scales? Here, nature has given us a wonderful gift, which will considerably simplify the writing of a cosmological model based on the equations of General Relativity. of General Relativity.

  1. the Universe is homogeneous: the metric therefore does not depend on an observer’s position in space, so no position is particular in the Universe. This assertion, based on the Copernican principle, is only statistically true, as we can observe that matter has formed lumps (planets, stars, galaxies, etc.) in the middle of large voids. However, observation of the Universe on large scales shows that the Universe is indeed globally homogeneous on scales larger than 100100\,Mpc[1] (see Figure 1 and for example Scrimgeour et al. (2012) for a measure of the Universe’s homogeneity by counting galaxies).

  2. the Universe is isotropic: no direction is privileged. This means that observations made in two different directions across the sky are equivalent. This is well verified by observations of the cosmic microwave background (CMB), whose temperature is measured to be identical at 2.725pm0.0022.725pm0.002\,K in all directions of space Mather et al., 1999. Only temperature fluctuations of the order of 10-5,K are detected in this image of the young Universe (see Figure 2 and for example The Planck Collaboration (2013) for a verification of the isotropy principle using the Sunyaev-Zeldovich effect).

Completely ignoring what happens at “insufficiently” large scales is the first step towards building a cosmological solution to General Relativity. Armed with these observational facts, we will impose homogeneity and isotropy on the metric and distribution of matter (i.e. the energy-impulsion tensor).

2Universe of maximum symmetry

Given the cosmological principle, we seek to determine the form that the metric of a Universe of maximum symmetry must take, i.e. a Universe whose properties are invariant by rotation and translation Weinberg, 1972[p. 379].

2.1Metric of an isotropic Universe

The space-time diagram of a homogeneous, isotropic Universe has a time axis orthogonal to the spatial base vectors, due to the isotropy of the Universe. Each hypersurface at time t is then a Universe of homogeneous density.

The space-time diagram of a homogeneous, isotropic Universe has a time axis orthogonal to the spatial base vectors, due to the isotropy of the Universe. Each hypersurface at time tt is then a Universe of homogeneous density.

First of all, if the Universe is isotropic, we can check that the crossed components gi0g_{i0} and gi0g_{i0} are zero. If this were not the case, we would have a privileged direction in the universe. We can convince ourselves of this by noticing that these components are non-zero if we perform a Lorentz transformation (6)-(7), precisely when we take a frame of reference in uniform translation with respect to another, thus moving in a chosen direction.

Another way of convincing ourselves is to consider a 2D spacetime. If the metric has the form :

g=(g00g01 g01g11)g=\begin{pmatrix} g_{00} & g_{01} \\\ g_{01} & g_{11} \end{pmatrix}

then the equation for light-like trajectories is :

ds2=0=g00c2dt2+2g01,cdxdt+g11dx2. \dd s^2=0=g_{00}c^2\dd t^2+2g_{01},c\,\dd x\,\dd t+g_{11}\dd x^2.

We can then check, by solving the second-degree equation in dt\dd t that if g010g_{01} \neq 0, two opposite dx\dd x give two different values of positive dt\dd t. In other words, an observer will receive at different times the light pulses emitted simultaneously by two sources located at the same distance in opposite ±dx\pm \dd x directions. This obviously breaks isotropy. The g0ig_{0i} and gi0g_{i0} terms of the metric are therefore zero. This means that the time vector e0\vec e_0 is orthogonal to the spatial basis vectors ei\vec e_i.

What’s more, if the Universe is homogeneous, then the g00(t,x)g_{00}(t,\vec x) component can only depend on time tt, so that clocks do not depend on position in space. So g00(t,x)=g00(t)g_{00}(t,\vec x) = g_{00}(t) Weinberg, 1972[p. 403]. If we call the parameter tt time, we see that we have a universal time at every point in space, called cosmological time. Since the Universe is homogeneous, this means that each date can be associated with a density of matter or energy that is identical for all observers, so that with a densimeter we can construct a clock common to all observers present in the Universe.

Combining the two previous results, the space-time interval can be written in the following form:

ds2=g00(t)c2dt2+dl2\dd s^2= g_{00}(t) c^2 \dd t^2 + \dd \vec l^2

where dl\dd \vec l is an elementary space vector. It is then possible to set g00g_{00} to -1, even if it means redefining the variable time[2]. The metric thus takes the form :

gμν=(10000γ11γ12γ130γ12γ22γ230γ13γ23γ33)g_{\mu\nu}=\begin{pmatrix} -1& 0 & 0 & 0\\ 0 & \gamma_{11} & \gamma_{12} & \gamma_{13} \\ 0&\gamma_{12} & \gamma_{22} & \gamma_{23} \\ 0&\gamma_{13} & \gamma_{23} & \gamma_{33} \end{pmatrix}

where γij\gamma_{ij} is the spatial metric, which can depend on time and position, and has 6 independent unknown components (a metric is a symmetrical tensor).

2.2Geometry of a maximally symmetrical Universe

Now let’s find an explicit form for dl2\dd \vec l^2. A maximally symmetrical Universe (homogeneous and isotropic) must have spatially constant curvature. This is fairly intuitive, but can also be demonstrated in General Relativity Weinberg, 1972 [p. 381]. Let aa be the associated radius of curvature, and let ξ=(ξ1,ξ2,ξ3)\vec \xi = (\xi^1, \xi^2, \xi^3) be a position vector in 3D space:

dl2=γijdξidξj,withi=1,2,3\dd \vec l^2 = \gamma_{ij} \dd \xi^i \dd \xi^j, \quad \text{with}\quad i=1,2,3

First of all, if this space has zero curvature, then the elementary distance dl\dd \vec l is simply written :

dl2=δijdξidξj,γij=δij\dd \vec l^2 = \delta_{ij} \dd \xi^i \dd \xi^j,\quad \gamma_{ij} = \delta_{ij}

Now let’s work on the case where the curvature is non-zero. To describe the curvature of a surface with the usual geometric notions, let’s study it in a space with an extra dimension. Let’s place this non-Euclidean (curved) 3D space in a 4D space of metric CABC_{AB} with Cartesian coordinates (x,y,z,w)(x, y, z, w). Let r2=x2+y2+z2r^2 = x^2 + y^2 + z^2 be the Euclidean distance in the 3D subspace. We then have two possible 3D hyper-surfaces, of constant Gaussian curvature 1/a21/a^2 Baumann, 2022 :

  • a 3-sphere of radius aa if immersed in a 4D Euclidean space:
CAB=diag(1,1,1,1),r2+w2=a2,dl2=dr2+dw2\quad C_{AB} = \mathrm{diag}(1,1,1,1), \quad r^2 + w^2= a^2,\quad \dd \vec l^2 = \dd r^2 + \dd w^2
  • a 3-hyperboloid of curvature aa if it is immersed in a 4D Lorentzian space:
CAB=diag(1,1,1,1),r2w2=a2,dl2=dr2dw2\quad C_{AB} = \mathrm{diag}(1,1,1,-1),\quad r^2 - w^2= -a^2,\quad \dd \vec l^2 = \dd r^2 - \dd w^2
2D spherical space of radius a represented in a 3D Euclidean space (x,y,w). At a coordinate \vec \xi = (\xi_1, \xi_2), we wish to express the length of an elementary vector tangent to the sphere \dd \vec l in both spaces.

2D spherical space of radius aa represented in a 3D Euclidean space (x,y,w)(x,y,w). At a coordinate ξ=(ξ1,ξ2)\vec \xi = (\xi_1, \xi_2), we wish to express the length of an elementary vector tangent to the sphere dl\dd \vec l in both spaces.

The last two cases of strictly positive or negative curvature are therefore defined by the constraint equation :

r2±w2=±a2(t)r^2 \pm w^2= \pm a^2(t)

where here we allow the radius a(t)a(t) to depend on time, since a priori γij\gamma_{ij} can depend on time.

The infinitesimal distance dl2\dd \vec l^2 between two points of the hypersurface defined in curved 3D space of metric γij\gamma_{ij} must be identical to that measured in 4D space, so :

dl2=γijdξidξj=dr2±dw2\dd \vec l^2= \gamma_{ij} \dd \xi^i \dd \xi^j = \dd r^2 \pm \dd w^2

where the ++ case corresponds to spherical geometry, the - case to hyperbolic geometry Weinberg, 1972 [p. 390-391].

However, the closure equation (9) links r,wr, w and aa, so we can replace dw\dd w by an expression that is a function of rr and aa (i.e. without the fourth dimension). Differentiating equation (9) gives the relation

(rdr)±wdw=0, (\vec r \cdot \dd \vec r) \pm w\dd w=0,

so, injecting equation (9) again, we obtain :

(rdr)2=(wdw)2(dw)2=(rdr)2w2=(rdr)2a2(t)r2(\vec r \cdot \dd \vec r)^2=(w\dd w)^2 \Rightarrow (\dd w)^2= \frac{(\vec r \cdot \dd \vec r)^2}{w^2} = \frac{(\vec r \cdot \dd \vec r)^2}{a^2(t) \mp r^2}

The infinitesimal distance between 2 points in non-Euclidean 3D space of constant non-zero curvature a2a^{-2} is then :

dl2=dr2±(rdr)2a2(t)r2\dd \vec l^2= \dd \vec r^2 \pm \frac{(\vec r \cdot \dd \vec r)^2}{a^2(t)\mp r^2}

At this stage, we can now combine the result obtained for the two non-zero curvatures with the Euclidean case by introducing the parameter of curvature kk :

k={+13-sphere}0flat space-13-hyperboloidk = \left\lbrace \begin{array}{rl} +1 & \text{3-sphere} \} 0 & \text{flat space} \text -1 & \text{3-hyperboloid} \\ \end{array}\right.

For the three possible geometries[^flat] of a maximally symmetrical Universe, we have :

dl2=dr2+k(rdr)2a2(t)kr2\dd \vec l^2= \dd \vec r^2 + k\frac{ (\vec r \cdot \dd \vec r)^2}{a^2(t) - k r^2}

where in the case of a flat space we’d have r=ξ\vec r = \vec \xi.

Finally, let’s introduce the rescaled variable σ=r/a(t)\vec\sigma=\vec r/a(t), and we get a new expression:

dl2=a2(t)(dσ2+k(σdσ)21kσ2)\dd \vec{l}^2= a^2(t) \left(\dd \vec{\sigma}^2 + k\frac{(\vec{\sigma} \cdot \dd \vec{\sigma})^2}{1 - k \sigma^2} \right)

The Friedmann-Lemaître-Robertson-Walker metric describing a homogeneous, isotropic Universe is finally written :

ds2=c2dt2+a2(t)(dσ2+k(σdσ)21kσ2)\dd s^2=-c^2\dd t^2 + a^2(t) \left(\dd \sigma^2 + k\frac{(\sigma \cdot \dd \sigma)^2}{1 - k \sigma^2} \right)

The Friedmann-Lemaître-Robertson-Walker (FLRW) metric is the basic framework of the Standard Cosmological Model. The assumptions of homogeneity and isotropy led directly to a metric describing a universe with only three possible geometries (flat, 3-sphere, 3-hyperboloid) and a scaling factor a(t)a(t) affecting distances. Note that, thanks to the imposition of homogeneity and isotropy symmetries, we have reduced the writing of the metric gμνg_{\mu\nu} (which is a symmetric tensor) consisting a priori of 10 unknown independent components to a tensor with a single unknown function a(t)a(t).

It’s important to understand the physical significance of the expansion factor a(t)a(t). First of all, according to equation (17), this factor relates the physical distance r\vec r and the coordinate distance σ by r=a(t)σ\vec r=a(t)\vec \sigma. A particle with fixed spatial coordinates σ\vec \sigma will see its physical distance from an observer in σ=0\vec \sigma=\vec 0 increase (or decrease) with time. This variation in distance is achieved at apparent speed:

drdt=da(t)σdt=a˙σ+aσ˙=a˙ar \frac{\dd \vec r}{\dd t} = \frac{\dd a(t)\vec \sigma}{\dd t} = \dot a \vec \sigma + a \dot{\vec \sigma} = \frac{\dot a}{a} \vec r

because σ˙=0\dot{\vec \sigma} = \vec 0 if the particle has no motion of its own, with the point   ˙  \;\dot{}\; expressing a derivative with respect to time tt. We thus obtain a direct relationship between distance from a central observer and apparent velocity: this is Hubble’s law. The distance rate is given by the Hubble parameter, which quantifies the rate of change of the scale factor:

H(t)=a˙(t)a(t)\boxed{\displaystyle H(t) = \frac{\dot a(t)}{a(t)}}

For a spherical universe, the scale factor a(t)a(t) also represents its radius of curvature. A dynamic spherical universe therefore corresponds to a universe with a time-varying radius of curvature. A flat space has no characteristic scale, so the value of a(t)a(t) is not a physical observable. The physically meaningful quantity for such a universe is the Hubble parameter H(t)H(t).

To simplify notation, the time dependence of the parameters is not always made explicit, so a(t)=aa(t)=a. Parameters evaluated at the present time t0t_0 are designated by the subscript or exponent 0, so that a(t0)=a0a(t_0)=a_0. In the following, we’ll work in the system where a0a_0 isn’t fixed at 1.

Table 1:Representation of the two-dimensional equivalents of the three solution spaces of cosmological principles: the 2-sphere, the plane, the 2-hyperboloid.

sphereplanhyperboloid

3Comobiles coordinates

It’s important to note that not all observers see the Universe as isotropic, but only so-called mobile observers, who are locally at rest with most of the matter in their vicinity. We, for example, are not mobile observers: when we observe the temperature of the CMB, the first feature we see is a large temperature dipole (warmer on one side, colder on the opposite side), which is the result of the particular motion of our solar system in the galaxy, and of our galaxy in the Universe (and of our group of galaxies). If we were to subtract this own motion from the CMB frame of reference, then we would be comoving observers. Thus, we can define a coordinate system associated with observers without proper motion, whose relative proper distances increase only with the scaling factor a(t)a(t) Weinberg, 1972[p. 409].

In the FLRW metric, where the expansion of the Universe is factorized by a scale factor a(t)a(t), the spatial coordinates σ\vec \sigma are called comoving coordinates. There is considerable freedom in the choice of comoving coordinates.

3.1Spherical coordinates

We often prefer spherical coordinates (ct,σ,θ,ϕ)(ct, \sigma, \theta, \phi) with the observer (ourselves) at the origin, such as :

σ1=σsinθcosϕσ2=σsinθsinϕσ3=σcosθ\begin{aligned} \sigma_1 &= \sigma \sin \theta \cos \phi \\ \sigma_2 &= \sigma \sin \theta \sin \phi \\ \sigma_3 &= \sigma \cos \theta \end{aligned}

After some simple but lengthy calculations (see notebook and here), in all three cases of curvature the FLRW metric is written :

ds2=c2dt2+a2(t)(11kσ2dσ2+σ2dθ2+σ2sin2θdϕ2)\dd s^2=-c^2\dd t^2 + a^2(t) \left( {1 \over 1-k\sigma^2}\dd \sigma^2 + \sigma^2 \dd \theta^2 + \sigma^2 \sin^2 \theta \dd \phi^2\right)

3.2Cartesian coordinates

The flat-universe case greatly simplifies the calculations that follow. Since the assumption of zero curvature is compatible with the ever-stricter constraints of cosmological observations, we will henceforth concentrate our analytical developments on the flat universe, mentioning results for the general case where necessary. In the case of zero curvature, it may be convenient to use Cartesian comoving coordinates (ct,x,y,z)(ct, x, y, z), such as :

σ1=x,σ2=y,σ3=z,σ2=x2+y2+z2\sigma_1 = x,\quad \sigma_2 = y,\quad \sigma_3 = z,\quad \sigma^2 = x^2 + y^2 + z^2

The FLRW metric is written in a flat universe:

gpmatrix=(10000a2(t)0000a2(t)0000a2(t))g_{pmatrix} = \begin{pmatrix} -1 & 0 & 0 & 0 \\ 0 & a^2(t) & 0 & 0\\ 0 & 0 & a^2(t) & 0 \\ 0 & 0 & 0 & a^2(t) \end{pmatrix}

4Geodesics in FLRW metric

What is the trajectory of a free-falling particle in an FLRW metric? From General Relativity, we know that such a particle moves along a geodesic xμ(s)x^\mu(s) whose equation is as follows:

d2xμds2+Γνκμdxνdsdxκds=0,{\dd^2 x^\mu \over \dd s^2} +\Gamma^{\mu}_{\,\,\nu\kappa} {\dd x^\nu \over \dd s}{\dd x^\kappa \over \dd s}=0,

where ss is any parameter describing the position along the geodesic (e.g. proper time). Another form of the geodesic equation will help us here, obtained from the definition of the covariant derivative (28):

d2xμds2Γμκνdxνdsdxκds=0.{\dd^2 x_\mu \over \dd s^2} -\Gamma^{\nu}_{\,\,\mu\kappa} {\dd x_\nu \over \dd s}{\dd x^\kappa \over \dd s}=0.

Let’s define the quadri-velocity UμU^\mu along a Universe line by Uμ=dxμ/dsU^\mu = \dd x^\mu / \dd s. Then :

dUμds=ΓμκνUνUκ=12gαβxμUαUβ.{\dd U_\mu \over \dd s} = \Gamma^{\nu}_{\,\,\mu\kappa} U_\nu U^\kappa= \frac{1}{2} \frac{\partial g_{\alpha\beta}}{\partial x^\mu} U^\alpha U^\beta.

From this form of the geodesic equation (see Hobson et al., 2006[p. 81] for a demonstration), let’s calculate the form that the contravariant vector UνU^\nu must take in a FLRW metric for a particle in free fall.

Let’s start with the case μ=3\mu=3 and use the spherical coordinates (σ,θ,ϕ)(\sigma,\theta,\phi). Since the FLRW metric does not depend on ϕ, then :

dU3ds=0 \frac{\dd U_3}{\dd s} = 0

so U3U_3 is a motion constant. Furthermore:

U3=g33U3=a2(t)(σsinθ)2U3 U_3 = g_{33} U^3 = a^2(t) (\sigma \sin \theta)^2 U^3

whose expression cancels at the origin at σ=0\sigma=0, where we observe. Since the U3U_3 component is constant, it is identically zero along the trajectory. We deduce:

U3=dϕds=0ϕ=constantU^3 = \frac{\dd \phi }{ \dd s} = 0 \Rightarrow \phi = \text{constant}

.

Let’s move on to the μ=2\mu=2 case. The only component of the metric depending on θ is g33g_{33} but U3U_3 is identically zero, so :

fracdU2ds=12gαβx2UαUβ=12g33x2U3U3=0.frac{\dd U_2}{\dd s} = \frac{1}{2} \frac{\partial g_{\alpha\beta}}{\partial x^2} U^\alpha U^\beta = \frac{1}{2} \frac{\partial g_{33}}{\partial x^2} U^3 U^3 = 0.

In the same way, we also have :

U2=g22U2=a2(t)σ2U2U_2 = g_{22} U^2 = a^2(t) \sigma^2 U^2

which cancels out at σ=0\sigma=0, so U2U^2 is zero all along the trajectory. From this we deduce:

U2=dθds=0θ=constantU^2 = \frac{\dd \theta }{ \dd s} = 0 \Rightarrow \theta = \text{constant}

5Spectral shift, or redshift

Notations for calculating redshift and cosmological distances in comoving coordinates.

Figure 6:Notations for calculating redshift and cosmological distances in comoving coordinates.

To measure the expansion history of the Universe, we need access to the scale parameter a(t)a(t). This is made possible by measuring the spectral shift of light coming from distant sources. In the FLRW metric, let’s place ourselves by convention at the center (σ=0\sigma=0), and consider an object located at comoving coordinates (σE,θE,ϕE)\left(\sigma_E,\theta_E,\phi_E\right), emitting a photon at time tEt_E (see Figure 6). For this photon, traveling at the speed of light, in the FLRW metric we have, at any instant:

ds2=0=c2dt2+a2(t)1kσ2dσ2.\dd s^2=0=-c^2 \dd t^2+\frac{a^2(t)}{1-k\sigma^2}\dd \sigma^2.

because along its geodesic θ and ϕ are constant (dθ=dϕ=0\dd \theta = \dd \phi=0). Let’s assume t0t_0 is the instant of reception of this wave at σ=0\sigma=0. Then, thanks to the previous equation, we have the relation :

tEt0cdta(t)=t0tEcdta(t)=0σEcdσsqrt1kσ2={σE si k=+1arcsinσE si k=0arcsh,σE si k=1.\int_{t_E}^{t_0} \frac{c\dd t}{a(t)} = -\int_{t_0}^{t_E} \frac{c\dd t}{a(t)} = \int_0^{\sigma_E}\frac{c\dd\sigma}{sqrt{1-k\sigma^2}} = \left\lbrace \begin{array}{cl} \sigma_E & \text{ si } k=+1 \arcsin \sigma_E & \text{ si } k=0 \\ \text{arcsh},\sigma_E & \text{ si } k=-1 \end{array} \right. .

with dσ<0\dd \sigma < 0 for dt>0\dd t > 0 considering a photon going from the source to the observer in 0.

For an electromagnetic wave with period TT, the expression (38) being valid at any instant, we can calculate the same integral for the wave emitted at instant tE+TEt_E+T_E and received at instant t0+T0t_0+T_0 (we assume that the period TT will vary with time):

tE+TEt0+T0cdta(t)=0σEdσ1kσ2.\int_{t_E+T_E}^{t_0+T_0} \frac{c \dd t}{a(t)}= \int_0^{\sigma_E}\frac{\dd \sigma}{1-k\sigma^2}.

Equating the expressions (39) and (40), since the period TT is small compared to the variations in the scale factor a(t)a(t) for usual electromagnetic waves, we obtain:

tE+TEt0+T0cdta(t)=tEt0cdta(t)tE+TEtEcdta(t)=t0+T0t0cdta(t)cT0a(t0)=cTEa(tE)λ0λE=a(t0)a(tE)\begin{aligned} \int_{t_E+T_E}^{t_0+T_0} \frac{c\dd t}{a(t)} & =\int_{t_E}^{t_0} \frac{c\dd t}{a(t)} \Leftrightarrow \int_{t_E+T_E}^{t_E} \frac{c\dd t}{a(t)} =\int_{t_0+T_0}^{t_0} \frac{c\dd t}{a(t)} \\ \Leftrightarrow \frac{cT_0}{a(t_0)} & = \frac{c T_E}{a(t_E)} \Leftrightarrow \frac{\lambda_0}{\lambda_E} = \frac{a(t_0)}{a(t_E)} \end{aligned}

Directly, if space is expanding, then a(tE)<a(t0)a(t_E) < a(t_0) and the received wavelength λ0\lambda_0 is therefore greater than the transmitted wavelength λE\lambda_E. We then define the spectral shift, commonly called redshift due to the fact that almost all observed galaxy spectra are redshifted, as :

z=λ0λEλE1+z=a0a(tE). \fbox{$ \displaystyle{z = \frac{\lambda_0-\lambda_E}{\lambda_E} \Leftrightarrow 1+z = \frac{a_0}{a(t_E)}} $}.

The spectral shift is both directly related to the scaling parameter a(t)a(t), and to an experimental quantity that can be directly measured on the emission spectrum of distant objects. Indeed, by looking at the position of the absorption and emission lines of distant objects, we can deduce their spectral shifts relative to the same chemical elements located on Earth, at rest. This experimental data is therefore often associated with the definition of distances in cosmology.

Solution to Exercise 1

For the first galaxy spectrum, the HβH\beta line is emitted at λE=4861\lambda_E = 4861\,Å and measured at λ05100\lambda_0\approx 5100\,Å, so the galaxy’s redshift is :

z=λ0λEλE=510048614861=0.049z = \frac{\lambda_0-\lambda_E}{\lambda_E} = \frac{5100-4861}{4861} = 0.049

For the second galaxy spectrum, the HβH\beta line is measured at λ05000\lambda_0\approx 5000\,Å so the redshift of the galaxy is :

z=λ0λEλE=500048614861=0.028z = \frac{\lambda_0-\lambda_E}{\lambda_E} = \frac{5000-4861}{4861} = 0.028

6Proper distance and comoving distance

Proper distance defines the physical distance between two objects at a time tt. Let a transmitting object be located at comoving coordinates (σE,θE,ϕE)(\sigma_E, \theta_E, \phi_E). By definition, the proper distance between this object and an observer located at the origin is along a curve at constant θ and ϕ and is equal to :

Dp(σE,t)=0σEgσσdσ=0σEa(t)dσ1kσ2=a(t)χ(σE)D_p(\sigma_E, t) = \int_0^{\sigma_E} \sqrt{g_{\sigma\sigma}}\dd\sigma' = \int_0^{\sigma_E}\frac{a(t)\dd\sigma'}{\sqrt{1-k\sigma'^2}} = a(t) \chi(\sigma_E)

where χ(σE)\chi(\sigma_E) is the distance comoving between this object and the observer:

χ(σE)=0σEdσ1kσ2={arcsinσE si k=+1σE si k=0arcsh,σE si k=1\chi(\sigma_E) = \int_0^{\sigma_E}\frac{\dd\sigma'}{\sqrt{1-k\sigma'^2}} = \left\lbrace\begin{array}{cl} \arcsin\sigma_E & \text{ si } k=+1 \\ \sigma_E & \text{ si } k=0 \\ \text{arcsh},\sigma_E & \text{ si } k=-1 \end{array} \right.

We can see that the proper distance DpD_p has the unit of a length, whereas the comoving distance is dimensionless. The latter represents distance in coordinate space and is independent of the expansion of the Universe. The proper distance, on the other hand, evolves over time with the scale factor.

Reciprocally, we define fk(χ)f_k(\chi) as follows:

σ=fk(χ)={sinχ si k=+1χ si k=0sinhχ si k=1\sigma = f_k(\chi) = \left\lbrace\begin{array}{cl} \sin\chi & \text{ si } k=+1 \\ \chi & \text{ si } k=0 \\ \sinh\chi & \text{ si } k=-1 \end{array} \right.

Now let’s imagine that this distance can be perceived through the journey of a photon. Light travels along a geodesic, so in the FLRW metric we have dθ=dϕ=0\dd \theta=\dd \phi=0 and :

ds2=0=c2dt2+a2(t)1kσ2dσ2.\dd s^2=0=-c^2 \dd t^2+\frac{a^2(t)}{1-k\sigma^2}\dd\sigma^2.

Therefore:

dσ1kσ2=cdta(t) \frac{\dd\sigma}{\sqrt{1-k\sigma^2}} = - \frac{c \dd t}{a(t)}

with the photon traveling along the direction dσ<0\dd \sigma<0 for dt>0\dd t > 0. The comoving distance is rewritten as

χ(σE)=0σEdσ1kσ2=t0tEcdta(t)=tEt0cdta(t)=χ(tE)\chi(\sigma_E) = \int_0^{\sigma_E}\frac{\dd\sigma}{\sqrt{1-k\sigma^2}} = \int_{t_0}^{t_E} -\frac{c\dd t'}{a(t')}= \int_{t_E}^{t_0} \frac{c\dd t'}{a(t')} = \chi(t_E)

Let’s express this distance in terms of the redshift zz, which, as we all know, is an experimental observable. Integrals can be transformed between the variables tt, aa by defining the expansion rate H(t)H(t), and between the variables aa and zz by defining the redshift:

H=1adadtdt=1aHdaH= \frac{1}{a}\frac{\dd a}{\dd t} \Rightarrow \dd t = \frac{1}{aH} \dd a
a=a01+zda=a0dz(1+z)2daa=dz1+za = \frac{a_0}{1+z} \Rightarrow \dd a = -a_0 \frac{\dd z}{(1+z)^2}\Rightarrow \frac{\dd a}{a} = -\frac{\dd z}{1+z}

Hence the comoving distance in terms of time tt, scale factor aa and redshift zz :

χ(σE)=χ(tE)=tEt0cdta(t)=aEa0cdaa2H(a)=z01+za0cdz(1+z)H(z)=1a00zcdzH(z)=χ(z)\begin{align} \chi(\sigma_E) & = \chi(t_E) = \int_{t_E}^{t_0} \frac{c\dd t'}{a(t')} = \int_{a_E}^{a_0} \frac{c\dd a}{a^2 H(a)} \\ & = - \int_z^0 \frac{1+z}{a_0}\frac{c\dd z}{(1+z)H(z)} = \frac{1}{a_0}\int_0^z\frac{c\dd z}{H(z)} = \chi(z) \end{align}

Table Table 3 summarizes how the parameters aa, tt and zz are converted at different times in the chronology of the Universe.

Table 3:Conversion of parameters aa, tt and zz. Note tUt_U, the age of the Universe today.

atz
reception (today)a0a_0tUt_U0
emission (past)aEa_EtEt_Ezz
beginning of the Universe00++\infty
end of time++\infty++\infty-1

The proper distance is the distance that could actually be measured between two objects at time tt. If we choose an object located at the comoving coordinate σE\sigma_E and a comoving observer in 0, then the proper distance today at t0t_0 is simply written for the three curvature cases:

Dp(z)=a0χ(z)=0zcdzH(z)D_p(z) = a_0 \chi(z) =\int_0^z\frac{c \dd z}{H(z)}

and is expressed in units of length. The notion of proper distance is illustrated Figure 7.

Proper distance between the Earth and a distant galaxy with no apparent proper velocity. (a) Today, the measured distance between the Earth and this galaxy is a_0 \sigma light-years in flat space. (b) At another date t, this distance evolves to a(t) \sigma. (c) Proper distance in spherical space.

Figure 7:Proper distance between the Earth and a distant galaxy with no apparent proper velocity. (a) Today, the measured distance between the Earth and this galaxy is a0σa_0 \sigma light-years in flat space. (b) At another date tt, this distance evolves to a(t)σa(t) \sigma. (c) Proper distance in spherical space.

Solution to Exercise 2 #
  1. The comoving coordinates χ is explicitly the polar angle between the observer and the galaxy. The proper distance is the arc length, which is Dp=a(t)χD_p = a(t) \chi as in usual geometry. The coordinate σ is related to the length bb of the chord at this angle χ by :
σ=sinχ=ba(t)\sigma = \sin \chi = \frac{b}{a(t)}

Note also that bb is the radius of the circle of colatitude χ: b=a(t)sinχb = a(t) \sin \chi.

  1. Using polar projection Figure 8, we see that θ is the angle delimiting an arc of size ll and radius bb, the radius of χ’s latitude circle. Thus :
θ=lb=la(t)sinχ=la(t)σ\theta = \frac{l}{b} = \frac{l}{a(t) \sin \chi} = \frac{l}{a(t) \sigma}

Thanks to the parameterization σ=fk(χ)\sigma=f_k(\chi), we can see that this equality works in both flat and curved space.

If we prefer to reason mathematically rather than graphically, using the FLRW metric, we can also demonstrate :

l=0ddl=0θgθθa(t)dθ=a(t)0θfk(χ)dθ=a(t)θsinχ=a(t)θσl = \int_0^d \dd l' = \int_0^\theta \sqrt{g_{\theta\theta}} a(t)\dd\theta' = a(t) \int_0^\theta f_k(\chi) \dd\theta' = a(t) \theta \sin\chi = a(t) \theta \sigma

So again we get θ=l/a(t)σ\theta = l / a(t) \sigma.

It’s easy to check that we have the same expressions in the flat case, and we’ll admit them in hyperbolic geometry. The σ coordinate is therefore useful for making calculations and drawings in the flat case and translating these results into the curved case (which isn’t so easy with χ).

Geometry in a spherical universe.

Figure 8:Geometry in a spherical universe.

7Cosmic time and conformal time

Time deserves a special mention. In our ideal Universe, with no over- or under-densities of matter, all clocks that follow expansion (i.e. with no motion of their own) beat the second at the same rate. With infinite time at our disposal to set all the clocks in the Universe, we can propose a universal convention for synchronizing our clocks. For example, when the temperature of the CMB reaches a given value, all the civilizations in the Universe can decide that this corresponds to a certain date. It is therefore possible to define a cosmic time, common to all comoving observers Weinberg, 1972 [p. 409].

Footnotes
  1. 1 parsec (pc) =3.262= 3.262 light-years =3.086×1016= 3.086\times 10^{16}\,m. 100Mpc3×108  100\,\text{Mpc}\approx 3\times 10^8\; light-years.

  2. We can introduce a new time variable tt' such that dt=g00dt\dd t' = \sqrt{\vert g_{00}\vert }\dd t.

References
  1. Scrimgeour, M. I., Davis, T., Blake, C., James, J. B., Poole, G. B., Staveley-Smith, L., Brough, S., Colless, M., Contreras, C., Couch, W., Croom, S., Croton, D., Drinkwater, M. J., Forster, K., Gilbank, D., Gladders, M., Glazebrook, K., Jelliffe, B., Jurek, R. J., … Yee, H. K. C. (2012). The WiggleZ Dark Energy Survey: the transition to large-scale cosmic homogeneity. Monthly Notices of the Royal Astronomical Society, 425(1), 116–134. 10.1111/j.1365-2966.2012.21402.x
  2. Mather, J. C., Fixsen, D. J., Shafer, R. A., Mosier, C., & Wilkinson, D. T. (1999). Calibrator Design for the COBE Far-Infrared Absolute Spectrophotometer (FIRAS). The Astrophysical Journal, 512(2), 511–520. 10.1086/306805
  3. The Planck Collaboration. (2013). Planck intermediate results. XIII. Constraints on peculiar velocities. Astronomy & Astrophysics, 561, A97. 10.1051/0004-6361/201321299
  4. Weinberg, S. (1972). Gravitation and cosmology: principles and applications of the general theory of relativity. http://www.lavoisier.fr/livre/notice.asp?ouvrage=1382255
  5. Baumann, D. (2022). Cosmology. Cambridge University Press. https://cmb.wintherscoming.no/pdfs/baumann.pdf