Skip to article frontmatterSkip to article content

1Principe cosmologique

Figure 1:Distribution des galaxies compilée par l’étude eBOSS. Chaque point de ce diagramme « circulaire » est une galaxie, codée par type : vert pour les galaxies proches, magenta et rouge pour les vieilles galaxies rouges, bleu pour les jeunes galaxies, jaune et blanc pour les quasars. Crédit : A. Raichoor (EPFL) / A. Ross (Ohio State Univ.) / Collaboration SDSS

Carte de température du fond diffus cosmologique micro-onde (CMB) observée par le satellite Planck. L’écart relatif observé entre la  température des points chauds (rouges) ou froids (bleus) par rapport à la moyenne est de l’ordre de \delta \theta / \theta \approx 10^{-5}.

Figure 2:Carte de température du fond diffus cosmologique micro-onde (CMB) observée par le satellite Planck. L’écart relatif observé entre la température des points chauds (rouges) ou froids (bleus) par rapport à la moyenne est de l’ordre de δθ/θ105\delta \theta / \theta \approx 10^{-5}.

Pour être en mesure de construire un modèle de l’Univers, c’est-à-dire une construction théorique capable de décrire le contenu de l’Univers et son évolution, il faut parvenir à résoudre l’équation d’Einstein de la Relativité Générale. Cependant, jusqu’à quel niveau de détail est-il nécessaire d’aller pour le décrire suffisamment bien aux grandes échelles ? On se doute que paramétrer l’équation d’Einstein jusqu’à inclure l’échelle du système solaire est à la fois illusoire et non nécessaire. Quelle est la structure de l’Univers aux plus grandes échelles ? Ici la nature nous fait un beau cadeau, qui va simplifier considérablement l’écriture d’un modèle cosmologique à partir des équations de la Relativité Générale.

  1. l’Univers est homogène : la métrique ne dépend donc pas de la position d’un observateur dans l’espace, donc aucune position n’est particulière dans l’Univers. Cette affirmation, issue du principe Copernicien, n’est que statistiquement vraie car localement on observe bien que la matière a formé des grumeaux (planètes, étoiles, galaxies,...) au milieu de larges vides. Cependant l’observation de l’Univers à grande échelle montre que l’Univers est bien globalement homogène à des échelles plus grandes que 100100\,Mpc[1] (voir Figure 1 et par exemple Scrimgeour et al. (2012) pour une mesure de l’homogénéité de l’Univers par comptage de galaxies).

  2. l’Univers est isotrope: aucune direction n’est privilégiée. Ainsi, des observations effectuées dans deux directions différentes du ciel sont équivalentes. Ceci est bien vérifié par l’observation du fond diffus cosmologique micro-onde (CMB) dont la température est mesurée identique à 2.725±0.002K2.725\pm0.002\,\kelvin dans toutes les directions de l’espace Mather et al., 1999. Seules des fluctuations de température de l’ordre de 105K10^{-5}\,\kelvin sont détectées sur cette image de l’Univers jeune (voir Figure 2 et par exemple The Planck Collaboration (2013) pour une vérification du principe d’isotropie utilisant l’effet Sunyaev-Zeldovich).

Ignorer complètement ce qui se passe à des échelles “insuffisamment” grandes est la première étape pour construire une solution cosmologique à la Relativité Générale. Muni de ces faits observationnels, nous imposerons l’homogénéité et l’isotropie à la métrique et à la distribution de la matière (c’est-à-dire au tenseur énergie-impulsion).

2Univers de symétrie maximale

Étant donné le principe cosmologique, on cherche à déterminer la forme que doit prendre la métrique d’un Univers de symétrie maximale, c’est-à-dire un Univers dont les propriétés sont invariantes par rotation et translation Weinberg, 1972[p. 379].

Métrique d’un Univers isotrope

Le diagramme d’espace-temps d’un Univers homogène et isotrope possède un axe temporel orthogonal aux vecteurs de base spatiaux à cause de l’isotropie de l’Univers. Chaque hypersurface à un instant t est alors un Univers de densité homogène.

Le diagramme d’espace-temps d’un Univers homogène et isotrope possède un axe temporel orthogonal aux vecteurs de base spatiaux à cause de l’isotropie de l’Univers. Chaque hypersurface à un instant tt est alors un Univers de densité homogène.

Tout d’abord, si l’univers est isotrope, on peut vérifier que les composantes croisées gi0g_{i0} et g0ig_{0i} sont nulles. Si tel n’était pas le cas, on aurait une direction privilégiée dans l’univers. On peut s’en convaincre en remarquant que ces composantes sont non nulles si on réalise une transformation de Lorentz (6)-(7), justement lorsqu’on prend un référentiel en translation uniforme par rapport à un autre, donc se déplaçant dans une direction choisie.

Autre façon de se convaincre, prenons un espace-temps 2D. Si la métrique a la forme :

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

alors l’équation des trajectoires de type lumière est :

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

On peut alors vérifier, en résolvant l’équation du second degré en dt\dd t que si g010g_{01} \neq 0, deux dx\dd x opposés donnent deux valeurs différentes de dt\dd t positif. C’est-à-dire qu’un observateur recevra à des moments différents les impulsions lumineuses émises simultanément par deux sources situées à la même distance dans des directions opposées ±dx\pm \dd x. Cela rompt évidemment l’isotropie. Les termes g0ig_{0i} et gi0g_{i0} de la métrique sont donc nuls. Cela signifie que le vecteur temporel e0\vec e_0 est orthogonal aux vecteurs de base spatiaux ei\vec e_i.

De plus, si l’Univers est homogène alors la composante g00(t,x)g_{00}(t,\vec x) ne peut dépendre que du temps tt, de sorte que le battement des horloges ne dépende pas de la position dans l’espace. Donc g00(t,x)=g00(t)g_{00}(t,\vec x) = g_{00}(t) Weinberg, 1972[p. 403]. Si nous appelons le paramètre tt le temps, nous voyons que nous avons un temps universel en tout point de l’espace, appelé temps cosmologique. Comme l’univers est homogène, alors cela signifie qu’à chaque date on peut associer une densité de matière ou d’énergie identique pour tous les observateurs, donc qu’avec un densimètre on peut construire une horloge commune à tous les observateurs présents dans l’Univers.

En combinant les deux résultats précédents, l’intervalle espace-temps peut être écrit sous la forme suivante :

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

dl\dd \vec l est un vecteur élémentaire spatial. Il est alors possible de fixer g00g_{00} à -1 quitte à redéfinir la variable temps[2]. La métrique prend donc la forme :

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}

γij\gamma_{ij} est la métrique spatiale, qui peut dépendre du temps et de la position, et comportant 6 composantes indépendantes inconnues (une métrique est un tenseur symétrique).

Géométrie d’un Univers maximallement symétrique

Trouvons maintenant une forme explicite pour dl2\dd \vec l^2. Un Univers de symétrie maximale (homogène et isotrope) doit posséder une courbure constante spatialement. Cela se comprend assez intuitivement mais aussi se démontre en Relativité Générale Weinberg, 1972 [p. 381]. Notons aa le rayon de courbure associé, et soit ξ=(ξ1,ξ2,ξ3)\vec \xi = (\xi^1, \xi^2, \xi^3) un vecteur position dans l’espace 3D :

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

Tout d’abord, si cet espace possède une courbure nulle, alors la distance élémentaire dl\dd \vec l s’écrit simplement :

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

Travaillons maintenant sur le cas où la courbure est non nulle. Pour décrire la courbure d’une surface avec des notions de géométries habituelles, étudions-la dans un espace avec une dimension supplémentaire. Plaçons cet espace 3D non euclidien (courbé) dans un espace 4D de métrique CABC_{AB} avec des coordonnées cartésiennes (x,y,z,w)(x, y, z, w). Posons r2=x2+y2+z2r^2 = x^2 + y^2 + z^2 la distance euclidienne dans le sous-espace 3D. On a alors deux hyper-surfaces 3D possibles, de courbure de Gauss constante 1/a21/a^2 Baumann, 2022 :

  • une 3-sphère de rayon aa si elle est plongée dans un espace Euclidien 4D :
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
  • une 3-hyperboloïde de courbure aa si elle est plongée dans un espace Lorentzien 4D :
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
Espace sphérique 2D de rayon a représenté dans un espace euclidien 3D (x,y,w). À une coordonnée \vec \xi = (\xi_1, \xi_2), on souhaite exprimé la longueur d’un vecteur élémentaire tangent à la sphère \dd \vec l dans les deux espaces.

Espace sphérique 2D de rayon aa représenté dans un espace euclidien 3D (x,y,w)(x,y,w). À une coordonnée ξ=(ξ1,ξ2)\vec \xi = (\xi_1, \xi_2), on souhaite exprimé la longueur d’un vecteur élémentaire tangent à la sphère dl\dd \vec l dans les deux espaces.

Les deux derniers cas de courbures strictement positive ou négative sont donc définis par l’équation de contrainte :

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

où on autorise ici le rayon a(t)a(t) à dépendre du temps, car a priori γij\gamma_{ij} peut dépendre du temps.

La distance infinitésimale dl2\dd \vec l^2 entre deux points de l’hypersurface définie dans l’espace 3D courbe de métrique γij\gamma_{ij} doit être identique à celle mesurée dans l’espace 4D, donc :

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

où le cas ++ correspond à une géométrie sphérique, le cas - à une géométrie hyperbolique Weinberg, 1972 [p. 390-391].

Or, l’équation de contrainte (9) lie r,wr, w et aa donc on peut remplacer dw\dd w par une expression fonction de rr et aa (donc sans la quatrième dimension). La différentiation de l’équation (9) donne la relation

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

donc, en injectant de nouveau l’équation (9), on obtient :

(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}

La distance infinitésimale entre 2 points de l’espace 3D non euclidien de courbure non nulle constante a2a^{-2} est alors :

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}

A cette étape, nous pouvons maintenant combiner le résultat obtenu pour les deux courbures non nulles avec le cas euclidien en introduisant le paramètre de courbure kk :

k={+13-spheˋre0espace plat13-hyperboloı¨dek = \left\lbrace \begin{array}{rl} +1 & \text{3-sphère} \\ 0 & \text{espace\ plat} \\ -1 & \text{3-hyperboloïde} \\ \end{array}\right.

On a ainsi pour les trois géométries possibles d’un Univers maximallement symétrique :

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}

où dans le cas d’un espace plat on aurait r=ξ\vec r = \vec \xi.

Enfin, introduisons la variable rééchelonnée σ=r/a(t)\vec\sigma=\vec r/a(t), et nous obtenons une nouvelle 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)

La métrique de Friedmann-Lemaître-Robertson-Walker décrivant un Univers homogène et isotrope s’écrit finalement :

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

La métrique de Friedmann-Lemaître-Robertson-Walker (FLRW) constitue le cadre de base du modèle cosmologique standard. Les hypothèses d’homogénéité et d’isotropie ont directement conduit à une métrique décrivant un univers avec seulement trois géométries possibles (plat, 3-sphère, 3-hyperboloid) et un facteur d’échelle a(t)a(t) affectant les distances. Notez que, grâce à l’imposition des symétries d’homogénéité et d’isotropie, nous avons réduit l’écriture de la métrique gμνg_{\mu\nu} (qui est un tenseur symétrique) constituée a priori de 10 composantes indépendantes inconnues à un tenseur possédant une seule fonction inconnue a(t)a(t).

Il est important de comprendre la signification physique du facteur d’expansion a(t)a(t). Tout d’abord, d’après l’équation (17), ce facteur relie la distance physique r\vec r et la distance de coordonnées σ par r=a(t)σ\vec r=a(t)\vec \sigma. Une particule dont les coordonnées spatiales σ\vec \sigma sont fixes verra sa distance physique avec un observateur en σ=0\vec \sigma=\vec 0 augmenter (ou diminuer) avec le temps. Cette variation de la distance se réalise à la vitesse apparente :

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

car σ˙=0\dot{\vec \sigma} = \vec 0 si la particule n’a pas de mouvement propre, avec le point   ˙  \;\dot{}\; exprimant une dérivée par rapport au temps tt. On obtient donc une relation directe entre la distance à un observateur central et la vitesse apparente : c’est la loi de Hubble. Le taux d’éloignement est donné par le paramètre de Hubble, qui quantifie la vitesse de variation du facteur d’échelle :

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

Pour un univers sphérique, le facteur d’échelle a(t)a(t) représente également son rayon de courbure. Un univers sphérique dynamique correspond donc à un univers possédant un rayon de courbure variable dans le temps. Un espace plat ne possède pas d’échelle caractéristique, la valeur de a(t)a(t) n’est donc pas une observable physique. La quantité ayant un sens physique pour un tel univers est le paramètre de Hubble H(t)H(t).

Dans le but d’alléger les notations, la dépendance temporelle des paramètres ne sera pas toujours explicitée, de sorte que a(t)=aa(t)=a. On désignera les paramètres évalués au temps présent t0t_0 par l’indice ou l’exposant 0 si bien que a(t0)=a0a(t_0)=a_0. Dans la suite, nous travaillerons dans le système où a0a_0 n’est pas fixé à 1.

Table 1:Représentation des équivalents à deux dimensions des trois espaces solution des principes cosmologiques: la 2-sphère, le plan, le 2-hyperboloïde.

sphereplanhyperboloid

3Coordonnées comobiles

Il est important de noter que tous les observateurs ne voient pas l’Univers comme isotrope, mais seulement les observateurs dits comobiles, qui sont localement au repos avec l’essentiel de la matière dans leur voisinage. Nous, par exemple, ne sommes pas des observateurs mobiles : lorsque nous observons la température du CMB, la première caractéristique que nous voyons est un grand dipôle de température (plus chaud d’un côté, plus froid du côté opposé), qui est le résultat du mouvement particulier de notre système solaire dans la galaxie, et de notre galaxie dans l’Univers (et de notre groupe de galaxies). Si on retranche ce mouvement propre par rapport au référentiel du CMB, alors cela ferait de nous des observateurs comobiles. Ainsi, on peut définir un système de coordonnées associés à des observateurs sans mouvements propres, dont les distances propres relatives n’augmentent qu’avec le facteur d’échelle a(t)a(t) Weinberg, 1972[p. 409].

Dans la métrique FLRW où l’expansion de l’Univers est factorisée par un facteur d’échelle a(t)a(t), les coordonnées spatiales σ\vec \sigma sont appelées coordonnées comobiles. Il existe une grande liberté dans le choix des coordonnées comobiles.

Coordonnées sphériques

On privilégie souvent les coordonnées sphériques (ct,σ,θ,ϕ)(ct, \sigma, \theta, \phi) avec l’observateur (nous-mêmes) à l’origine, telles que :

σ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}

Après des calculs simples mais longs (voir notebook et ici), dans les trois cas de courbures la métrique FLRW s’écrit :

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)

Coordonnées cartésiennes

Le cas de l’univers plat simplifie beaucoup les calculs qui suivront. Étant donné que l’hypothèse d’une courbure nulle est compatible avec les contraintes de plus en plus fortes des observations cosmologiques, nous concentrerons désormais nos développements analytiques sur l’univers plat, en mentionnant des résultats pour le cas général lorsque cela est nécessaire. Dans le cas de courbure nulle, il peut être pratique d’utiliser les coordonnées comobiles cartésiennes (ct,x,y,z)(ct, x, y, z), telles que :

σ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

La métrique FLRW s’écrit dans un univers plat :

gμν=(10000a2(t)0000a2(t)0000a2(t))g_{\mu\nu} = \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}

4Géodésiques dans la métrique FLRW

Quelle est la trajectoire d’une particule en chute libre dans une métrique FLRW ? D’après la Relativité Générale, nous savons qu’une telle particule se déplace le long d’une géodésique xμ(s)x^\mu(s) dont l’équation est la suivante :

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,

ss est un paramètre quelconque décrivant la position le long de la géodésique (le temps propre par exemple). Une autre forme de l’équation des géodésiques va ici nous aider, obtenue à partir de la définition de la dérivée covariante (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.

Définissons la quadri-vitesse UμU^\mu le long d’une ligne d’Univers par Uμ=dxμ/dsU^\mu = \dd x^\mu / \dd s. Alors :

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.

A partir de cette forme de l’équation des géodésiques (voir Hobson et al., 2006[p. 81] pour une démonstration), calculons la forme que doit prendre le vecteur contravariant UνU^\nu dans une métrique FLRW pour une particule en chute libre.

Commençons par le cas μ=3\mu=3 et utilisons les coordonnées sphériques (σ,θ,ϕ)(\sigma,\theta,\phi). Comme la métrique FLRW ne dépend pas de ϕ, alors :

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

donc U3U_3 est une constante du mouvement. Or, on a par ailleurs :

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

dont l’expression s’annule à l’origine en σ=0\sigma=0, là nous observons. Comme la composante U3U_3 est constante alors elle est identiquement nulle le long de la trajectoire. On en déduit :

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

Passons au cas μ=2\mu=2. La seule composante de la métrique dépendant de θ est g33g_{33} mais U3U_3 est identiquement nul donc :

dU2ds=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.

Or de même on a par ailleurs :

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

qui s’annule en σ=0\sigma=0 donc U2U^2 est nul tout le long de la trajectoire. On en déduit :

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

5Le décalage spectral, ou redshift

Notations pour le calcul du redshift et des distances cosmologiques en coordonnées comobiles.

Figure 6:Notations pour le calcul du redshift et des distances cosmologiques en coordonnées comobiles.

Pour mesurer l’histoire de l’expansion de l’Univers, il faut avoir accès au paramètre d’échelle a(t)a(t). Ceci est possible par la mesure du décalage spectral de la lumière venant de sources distantes. Dans la métrique FLRW, plaçons-nous par convention au centre (σ=0\sigma=0), et considérons un objet situé aux coordonnées comobiles (σE,θE,ϕE)\left(\sigma_E,\theta_E,\phi_E\right), émettant un photon à l’instant tEt_E (voir Figure 6). Pour ce photon, voyageant à la vitesse de la lumière, dans la métrique FLRW on a, à tout 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.

car le long de sa géodésique θ et ϕ sont constants (dθ=dϕ=0\dd \theta = \dd \phi=0). Posons t0t_0 l’instant de la réception de cette onde en σ=0\sigma=0. Alors grâce à l’équation précédente on a la relation :

tEt0cdta(t)=t0tEcdta(t)=0σEdσ1kσ2={arcsinσE si k=+1σ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{\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. .

avec dσ<0\dd \sigma < 0 pour dt>0\dd t > 0 en considérant un photon allant de la source vers l’observateur en 0.

Pour une onde électromagnétique de période TT, l’expression (38) étant valable à tout instant, on peut calculer la même intégrale pour l’onde émise à l’instant tE+TEt_E+T_E et reçue à l’instant t0+T0t_0+T_0 (on suppose donc que la période TT va varier au cours du temps):

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}.

Par égalité des expressions (39) et (40), comme la période TT est petite devant les variations du facteur d’échelle a(t)a(t) pour les ondes électromagnétiques usuelles, on obtient:

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}

Directement, si l’espace est en expansion alors a(tE)<a(t0)a(t_E) < a(t_0) et la longueur d’onde reçue λ0\lambda_0 est donc supérieure à la longueur d’onde émise λE\lambda_E. On définit alors le décalage spectral, communément appelé redshift en raison du fait que la quasi-totalité des spectres des galaxies observées sont décalées vers le rouge, par :

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

Le décalage spectral est à la fois directement lié au paramètre d’échelle a(t)a(t), mais aussi à une grandeur expérimentale directement mesurable sur le spectre d’émission des objets distants. En effet, en regardant la position des raies d’absorption et d’émission des objets lointains, on peut en déduire leurs décalages spectraux par rapport aux mêmes éléments chimiques situés sur Terre, au repos. Cette donnée expérimentale est donc souvent associée à la définition des distances en cosmologie.

Solution to Exercise 1

Pour le premier spectre de galaxie, la raie HβH\beta est émise à λE=4861\lambda_E = 4861\,Å et mesurée à λ05100\lambda_0\approx 5100\,Å, le décalage vers le rouge de la galaxie est donc :

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

Pour le second spectre de galaxie, la raie HβH\beta est mesurée à λ05000\lambda_0\approx 5000\,Å donc le décalage vers le rouge de la galaxie est :

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

6Distance propre et distance comobile

La distance propre définit la distance physique entre deux objets à un instant tt. Soit un objet émetteur situé aux coordonnées comobiles (σE,θE,ϕE)(\sigma_E, \theta_E, \phi_E). Par définition, la distance propre entre cet objet et un observateur situé à l’origine est le long d’une courbe à θ et ϕ constants et vaut :

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)

où on fait apparaître χ(σE)\chi(\sigma_E) la distance comobile entre cet objet et l’observateur :

χ(σ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.

On voit que la distance propre DpD_p possède bien l’unité d’une longueur alors que la distance comobile est sans dimension. Cette dernière représente la distance dans l’espace des coordonnées et est indépendante de l’expansion de l’Univers. En revanche la distance propre évolue dans le temps avec le facteur d’échelle.

Réciproquement, on définit fk(χ)f_k(\chi) ainsi :

σ=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.

Imaginons maintenant que cette distance puisse être perçue par le voyage d’un photon. La lumière voyage en suivant une géodésique, donc dans la métrique FLRW on a bien dθ=dϕ=0\dd \theta=\dd \phi=0 et :

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.

Donc :

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

avec le photon voyageant le long de la direction dσ<0\dd \sigma<0 pour dt>0\dd t > 0. La distance comobile se réécrit :

χ(σ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)

Exprimons cette distance en fonction du redshift zz, qui, on le rappelle, est une observable expérimentale. Les intégrales peuvent se transformer entre les variables tt, aa grâce à la définition du taux d’expansion H(t)H(t), et entre les variables aa et zz grâce à la définition du 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}

D’où la distance comobile en terme de temps tt, facteur d’échelle aa et 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}

On résume dans le tableau Table 3 comment se convertissent les paramètres aa, tt et zz à différents moments dans la chronologie de l’Univers.

Table 3:Conversion des paramètres aa, tt et zz. On note tUt_U l’âge de l’Univers aujourd’hui.

atz
réception (aujourd’hui)a0a_0tUt_U0
émission (passé)aEa_EtEt_Ezz
début de l’Univers00++\infty
fin des temps++\infty++\infty-1

La distance propre est la distance que l’on pourrait mesurer effectivement à un instant tt entre deux objets. Si on choisit un objet situé à la coordonné comobile σE\sigma_E et un observateur comobile en 0, la distance propre aujourd’hui à t0t_0 s’écrit alors simplement pour les trois cas de courbure :

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

et s’exprime bien en unités de longueur. La notion de distance propre est illustrée Figure 7.

Distance propre entre la Terre et une galaxie lointaine sans vitesse propre apparente. (a) Aujourd’hui, la distance mesurée entre la Terre et cette galaxie est de a_0 \sigma années-lumière dans un espace plat. (b) A une autre date t, cette distance évolue et vaut a(t) \sigma. (c) Distance propre dans un espace sphérique.

Figure 7:Distance propre entre la Terre et une galaxie lointaine sans vitesse propre apparente. (a) Aujourd’hui, la distance mesurée entre la Terre et cette galaxie est de a0σa_0 \sigma années-lumière dans un espace plat. (b) A une autre date tt, cette distance évolue et vaut a(t)σa(t) \sigma. (c) Distance propre dans un espace sphérique.

Solution to Exercise 2 #
  1. La coordonnée comobile χ est explicitement l’angle polaire entre l’observateur et la galaxie. La distance propre est la longueur de l’arc qui est Dp=a(t)χD_p = a(t) \chi comme dans la géométrie habituelle. La coordonnée σ est liée à la longueur bb de la corde à cet angle χ par :
σ=sinχ=ba(t)\sigma = \sin \chi = \frac{b}{a(t)}

Notons également que bb est le rayon du cercle de colatitude χ : b=a(t)sinχb = a(t) \sin \chi.

  1. Avec la projection polaire Figure 8, on voit que θ est l’angle délimitant un arc de taille ll et de rayon bb, le rayon du cercle de latitude de χ. Ainsi :
θ=lb=la(t)sinχ=la(t)σ\theta = \frac{l}{b} = \frac{l}{a(t) \sin \chi} = \frac{l}{a(t) \sigma}

Grâce au paramétrage σ=fk(χ)\sigma=f_k(\chi), on voit que cette égalité fonctionne en espace plat et en espace courbe.

Si on préfère un raisonnement mathématique plutôt que graphique, en utilisant la métrique FLRW, nous pouvons démontrer également :

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

Donc on obtient de nouveau θ=l/a(t)σ\theta = l / a(t) \sigma.

Il est facile de vérifier que nous avons les mêmes expressions dans le cas plat, et nous les admettrons dans la géométrie hyperbolique. La coordonnée σ est donc utile pour faire des calculs et des dessins dans le cas plat et traduire ces résultats dans les cas courbes (ce qui n’est pas si facile avec χ).

Géométrie dans un univers sphérique.

Figure 8:Géométrie dans un univers sphérique.

7Temps cosmique et temps conforme

Le temps mérite une mention spéciale. Dans notre Univers idéal, sans surdensités ou sous-densités de matière, toutes les horloges qui suivent l’expansion (c’est-à-dire sans mouvement propre) battent la seconde à la même cadence. Avec un temps infini à notre disposition pour mettre à l’heure toutes les horloges de l’Univers, nous pouvons proposer une convention universelle pour synchroniser nos horloges. Par exemple, lorsque la température du CMB atteint une valeur donnée, toutes les civilisations de l’Univers peuvent décider que cela correspond à une certaine date. Il est donc possible de définir un temps cosmique, commun à tous les observateurs comobiles Weinberg, 1972 [p. 409].

Footnotes
  1. 1 parsec (pc) =3.262= 3.262 années-lumière =3.086×1016= 3.086\times 10^{16}\,m. 100Mpc3×108  100\,\text{Mpc}\approx 3\times 10^8\; années-lumière.

  2. On peut introduire une nouvelle variable temps tt' telle que 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
  6. 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