# Équations de Friedmann

Les équations de Friedmann-Lemaître sont les équations de la relativité générale (appelées équations d'Einstein) écrites dans le contexte d'un modèle cosmologique homogène et isotrope, ce dernier étant représenté par une métrique de Robertson-Walker. Elles régissent donc l'évolution du taux d'expansion de l'Univers et par suite de la distance entre deux astres lointains (le facteur d'échelle) et en fonction du temps appelé dans ce contexte temps cosmique. L'évolution de ces quantités est déterminée par les propriétés du contenu matériel de l'Univers (rayonnement, atomes, matière noire, constante cosmologique, etc.), et éventuellement par la théorie de la gravitation considérée ; il est en effet possible de remplacer la relativité générale par une autre théorie relativiste de la gravitation. Il peut par exemple s'agir d'une théorie tenseur-scalaire. Une autre façon de changer le modèle est de considérer la relativité générale standard mais dans un univers possédant une ou plusieurs dimensions supplémentaires. C'est le cas des modèles de cosmologie branaire.
Les équations de Friedmann-Lemaître sont le fondement de la quasi-totalité des modèles cosmologiques, dont le Big Bang.

## Histoire

Les éponymes des équations de Friedmann-Lemaître sont le physicien russe Alexandre Friedmann (1888-1925) et le physicien belge Georges Lemaître (1884-1966). Friedmann est le premier des deux à les avoir proposées : d'abord dans un article$^,$ paru en 1922 et relatif au cas des espaces à courbure spatiale positive ; puis dans un article$^,$ paru en 1924 et plus général, incluant le cas d'une courbure spatiale négative. Peu de temps après Friedmann et indépendamment de celui-ci, Lemaître retrouve ces équations qu'il publie$^,$ en 1927$^,$ et interprète sans ambiguïté le phénomène d'expansion qu'elles impliquent et prédisant ou expliquant la loi de Hubble avant que celle-ci ne soit découverte en 1929. Auparavant, en 1917, Albert Einstein les a écrites dans le cas particulier d'un univers statique, ainsi que Willem de Sitter dans le cas d'un univers vide de matière, mais avec une constante cosmologique. La métrique de l'espace-temps et sa forme générale sont trouvées et présentées sous une forme unifiée par Howard P. Robertson en 1929, puis indépendamment par Arthur G. Walker en 1936. Pour toutes ces raisons, le type de modèle cosmologique décrit par ces équations est appelé univers de Friedmann-Lemaître-Robertson-Walker (abrégé en FLRW).
Bien que Friedmann et Lemaître soient considérés comme les codécouvreurs des équations, le nom de Lemaître est souvent omis. En particulier, d'une part, l'expression « équation de Friedmann » (au singulier) est employée pour désigner une des deux équations, à savoir : $H^{2}=\left({\dot {a}}/a\right)^{2}$ ; et, d'autre part, l'expression « équations de Friedmann » (au pluriel) est employée pour désigner les deux équations dans le cas $\Lambda =0$. L'ordre des noms est rarement inversé et celui de Friedmann est plus rarement omis que celui de Lemaître.

## Les deux équations de Friedmann-Lemaître

Les équations des Friedmann-Lemaître sont un ensemble de deux$^,$$^,$ équations exprimant l'évolution du facteur d'échelle cosmologique $a(t)$ en fonction du temps cosmologique $t$, et décrivant l'expansion de l'Univers. La première relie le taux d'expansion ou paramètre de Hubble H = ȧ/a, la courbure spatiale K = k/a et le facteur d'échelle $a$ à la densité volumique de masse ρ; la seconde relie la pression P à la dérivée temporelle du taux d'expansion. La variable de temps utilisée est le temps cosmique, qui correspond essentiellement au temps mesuré sur Terre. Dans le cadre de la relativité générale, ces deux équations s'écrivent :
    $3\left({\frac {H^{2}}{c^{2}}}+{\frac {K}{a^{2}}}\right)={\frac {8\pi G}{c^{2}}}\rho$,
    $-2{\frac {\dot {H}}{c^{2}}}-3{\frac {H^{2}}{c^{2}}}-{\frac {K}{a^{2}}}={\frac {8\pi G}{c^{4}}}P$,
où G est la constante de Newton et c la vitesse de la lumière (la quantité 8πG/c$^4$ est parfois appelée constante d'Einstein). Ces deux équations ne sont pas indépendantes : la seconde s'obtient en prenant la dérivée temporelle de la première, et en utilisant le principe des équations de conservation appliqué à la conservation de l'énergie, reliant la dérivée de la densité d'énergie à la pression. Quand on considère une extension de la relativité générale ou une autre théorie, ces équations sont modifiées. Pour cette raison, l'expression « équation de Friedmann » est parfois employée au singulier, auquel cas la seule équation considérée est la première.
Il peut arriver (par exemple dans les modèles à dimensions supplémentaires) qu'elles ne soient plus conséquences l'une de l'autre. Dans les autres cas, il est possible de combiner les deux équations de façon à modifier la seconde, soit pour ne faire apparaître qu'une certaine combinaison de la pression et de la densité d'énergie dans le membre de droite, soit pour ne faire apparaître que la dérivée seconde du facteur d'échelle dans le membre de gauche (voir paragraphe ci-dessous). Il est également possible d'effectuer un changement de variable pour utiliser le temps conforme plutôt que le temps cosmique. La résolution de ces équations s'effectue une fois la dépendance de la densité d'énergie et de la pression par rapport au temps ou au facteur d'échelle connus. Certaines solutions exactes sont connues.
Dans le cadre du modèle standard de la cosmologie (modèle Big Bang), l'allure générale du facteur d'échelle a(t) en fonction de la courbure spatiale K et de la constante cosmologique Λ est connue.

## Quelques solutions particulières

Quand le membre de droite ne comporte qu'une seule espèce, et en l'absence de courbure spatiale, les équations de Friedmann peuvent être résolues sans difficulté. En présence de plusieurs types de matière et/ou de courbure, des solutions analytiques exactes peuvent parfois être trouvées. Dans les autres cas, une résolution numérique des équations se fait sans difficulté.

### Univers de poussière

Quand l'Univers est empli de matière non relativiste (c'est-à-dire dont la pression est négligeable, d'où le terme de « poussière »), la densité d'énergie décroît uniquement du fait de la dilution due à l'expansion. On parle aussi d'ère de la matière. Pour un univers de densité critique, appelé univers d'Einstein-de Sitter, le facteur d'échelle évolue alors selon la loi
    $a(t)=a_{0}\left({\frac {t}{t_{0}}}\right)^{\frac {2}{3}}$,
les quantités a$_0$ et t$_0$ correspondant au facteur d'échelle et au temps évalués à une époque de référence donnée. De plus, l'âge de l'Univers à cette époque, t$_0$ se déduit du taux d'expansion à la même époque par $t_{0}={\frac {2}{3}}{\frac {1}{H_{0}}}$.
De fait de la nature de la matière, la densité d'énergie ρ$_m$ évolue selon la loi
    $\rho _{\rm {m}}=\rho _{\rm {m}}^{0}\left({\frac {a_{0}}{a}}\right)^{3}$,
l'indice 0 se référant aux quantités correspondantes évaluée à une époque de référence donnée. On a ainsi, en l'absence de courbure
    ${\frac {3}{c^{2}}}\left({\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}t}}\right)^{2}={\frac {8\pi G}{c^{4}}}\rho _{\rm {m}}^{0}\left({\frac {a_{0}}{a}}\right)^{3}$.
On peut évaluer cette équation à une époque donnée, ce qui donne
    ${\frac {3}{c^{2}}}H_{0}^{2}={\frac {8\pi G}{c^{4}}}\rho _{\rm {m}}^{0}$.
Si l'on note par x le « facteur d'échelle réduit », normalisé à 1 à l'époque de référence, la combinaison de ces deux équations peut s'écrire
    $\left({\frac {1}{x}}{\frac {{\rm {d}}x}{{\rm {d}}(H_{0}t)}}\right)^{2}={\frac {1}{x^{3}}}$.
La résolution de cette équation donne immédiatement
    $x(t)=\left({\frac {3}{2}}H_{0}t\right)^{\frac {2}{3}}$.
Une telle solution est appelée, pour des raisons historiques univers d'Einstein-de Sitter, bien que ces derniers ne soient pas les premiers à avoir exhibé cette solution.
Dans un tel modèle, l'univers est donc plus jeune que le temps de Hubble. Les premières déterminations de Hubble (quelque 520 kilomètres par seconde et par mégaparsec) étaient erronées par suite de mauvais étalonnage sur les Céphéides (il en existe de deux types !) ce qui conduisait à un âge de l'univers autour de 2 milliards d'années (!) bien plus jeunes que les espèces déjà vivant sur Terre. Cet âge étant significativement inférieur à l'âge estimé de nombreux objets astrophysiques (certaines étoiles, naines blanches, et la Voie lactée dans son ensemble), il est considéré qu'un tel univers de poussière ne peut correspondre à l'univers observable, ce qui est une des indications parmi d'autres de l'existence de l'énergie noire (voir ci-dessous).
On peut consulter utilement les récents articles qui résument l'évolution de la constante de Hubble, tenant compte de plusieurs générations d'expériences. Actuellement, d'après les résultats de WMAP (2003), la valeur mesurée de la constante de Hubble est donnée 70,1 ± 1,3 kilomètres par seconde et par mégaparsec, correspondant à un âge de 13,7 ± 1% milliards d'années. Il est aussi à noter que la dépendance du facteur d'échelle par rapport au temps est telle que sa dérivée seconde est négative, ce qui correspond à une phase d'expansion décélérée, compatible avec l'intuition que la nature attractive de la gravité a tendance à ralentir le mouvement d'expansion. Les observations, en particulier celle des supernovae de type Ia, suggèrent que la phase actuelle d'expansion est accélérée, ce qui est une autre indication de la nécessité de la présence d'énergie noire. Il est par contre certain que l'univers a connu par le passé une phase où son expansion était dominée par de la matière non relativiste, car c'est la seule époque où le mécanisme d'instabilité gravitationnelle peut se produire.

### Univers de radiation

Quand l'Univers est empli de matière relativiste ou de radiation (ère radiative), la densité d'énergie ρ$_r$ décroît plus vite que précédemment car en plus de la dilution due à l'expansion, l'énergie des particules individuelles décroît avec l'expansion (ce n'est rien d'autre que l'effet de décalage vers le rouge). La densité décroissant plus vite, le taux d'expansion décroît également plus vite, et l'expansion décélère plus rapidement que dans le cas d'un univers de poussière. Le facteur d'échelle évolue selon la loi
    $a=a_{0}\left({\frac {t}{t_{0}}}\right)^{\frac {1}{2}}$.
L'âge de l'univers dans un tel modèle s'écrit alors
    $t_{0}={\frac {1}{2}}{\frac {1}{H_{0}}}$.
De fait de la nature de la matière, la densité d'énergie ρ$_r$ évolue désormais selon la loi
    $\rho _{\rm {r}}=\rho _{\rm {r}}^{0}\left({\frac {a_{0}}{a}}\right)^{4}$,
l'indice 0 se référant à nouveau aux quantités correspondantes évaluée à une époque de référence donnée. On a ainsi, en l'absence de courbure
    ${\frac {3}{c^{2}}}\left({\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}t}}\right)^{2}={\frac {8\pi G}{c^{4}}}\rho _{\rm {r}}^{0}\left({\frac {a_{0}}{a}}\right)^{4}$.
On peut évaluer cette équation à une époque donnée, ce qui donne
    ${\frac {3}{c^{2}}}H_{0}^{2}={\frac {8\pi G}{c^{4}}}\rho _{\rm {r}}^{0}$.
Si l'on note par x le « facteur d'échelle réduit », normalisé à 1 à l'époque de référence, la combinaison de ces deux équations peut s'écrire
    $\left({\frac {1}{x}}{\frac {{\rm {d}}x}{{\rm {d}}(H_{0}t)}}\right)^{2}={\frac {1}{x^{4}}}$.
La résolution de cette équation donne immédiatement
    $x(t)=\left(2H_{0}t\right)^{\frac {1}{2}}$.
Comme dans le cas de l'Univers de poussière, l'expansion est décélérée et l'âge donné par ce modèle ne peut correspondre aux observations. Ce modèle correspond par contre bien aux époques reculées de l'histoire de l'Univers, quand la quasi-totalité de la matière était sous forme relativiste. En particulier l'époque de la nucléosynthèse s'est produite quand la densité d'énergie était majoritairement due à de la matière relativiste, hypothèse que l'on teste directement par la mesure de l'abondance des éléments légers. D'autres indications du fait que l'Univers a connu une période dominée par la radiation est l'existence du fond diffus cosmologique.

### Constante cosmologique

Une constante cosmologique peut s'interpréter comme une forme de matière de densité d'énergie constante et de pression exactement opposée à celle-ci. Quand la densité d'énergie est positive, les équations de Friedmann admettent, en l'absence de courbure spatiale, la solution exponentielle
    $a\propto e^{H_{0}t}$,
le paramètre de Hubble étant constant au cours du temps. La domination de la constante cosmologique est l'ère de l'énergie noire.
La première équation de Friedmann donne comme d'habitude
    $3{\frac {H^{2}}{c^{2}}}={\frac {8\pi G}{c^{4}}}\rho _{\Lambda }$,
la densité d'énergie de la constante cosmologique étant ici constante. Le paramètre de Hubble est donc constant, et l'on a donc
    ${\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}t}}=H_{0}={\rm {Constante}}$,
dont la solution est de façon évidente
    $a=a_{0}e^{H_{0}(t-t_{0})}$
Un tel modèle représente un univers éternel, sans commencement ni fin. C'est l'univers de de Sitter, obéissant au principe cosmologique parfait. Il existe une autre configuration donnant lieu à la même dynamique de l'expansion : les deux formes de matière sont alors de la matière ordinaire non relativiste, et une forme de matière extrêmement exotique, appelée champ C, responsable d'une création continue de matière compensant exactement la dilution due à l'expansion. C'est la base de la théorie de l'état stationnaire, aujourd'hui abandonnée en raison de l'absence totale de motivation théorique ou observationnelle du champ C. Enfin, en présence d'un champ scalaire, celui-ci peut, temporairement se comporter de façon extrêmement semblable à une constante cosmologique, comportement qu'il peut par la suite abandonner. Ainsi, il est possible que sous l'effet d'un tel champ scalaire, l'expansion de l'Univers connaisse une phase exponentielle. C'est la base des modèles d'inflation cosmique qui prédisent l'existence d'une telle phase à une époque très reculée de l'histoire de l'Univers.

### Équation d'état constante

Tant que la pression est proportionnelle à la densité, avec une constante de proportionnalité constante, on peut résoudre les équations de Friedmann. On note w le rapport de la pression à la densité,
    $w={\frac {P}{\rho }}$.
Pour toute forme de matière raisonnable, le facteur w est compris entre -1 et 1. Le cas où w est inférieur à -1 correspond à ce qui est appelé énergie fantôme, une forme d'énergie extrêmement spéculative donnant lieu à un scénario cosmologique étrange, le Big Rip. Les cas où w vaut 0, 1/3, -1 correspondent respectivement aux cas de la poussière, la radiation et de la constante cosmologique. Quand w est différent de -1, on obtient
    $a=a_{0}\left({\frac {3(1+w)}{2}}H_{0}|t|\right)^{\frac {2}{3(1+w)}}$.
Par les équations de conservation, l'évolution temporelle de la densité d'énergie est déterminée par
    ${\frac {{\rm {d}}\rho }{{\rm {d}}t}}=-3H(P+\rho )$.
En utilisant la relation entre densité d'énergie et pression et la définition de la constante de Hubble, cette équation se réécrit
    ${\frac {{\rm {d}}\rho }{{\rm {d}}t}}=-3{\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}t}}(1+w)\rho$,
soit
    $\rho \propto a^{-3(1+w)}$.
Les équations de Friedmann se réécrivent alors
    ${\frac {3}{c^{2}}}\left({\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}t}}\right)^{2}={\frac {8\pi G}{c^{4}}}\rho _{0}\left({\frac {a_{0}}{a}}\right)^{3(1+w)}$.
On peut évaluer cette équation à une époque donnée, ce qui donne
    ${\frac {3}{c^{2}}}H_{0}^{2}={\frac {8\pi G}{c^{4}}}\rho _{0}$.
Si l'on note par x le « facteur d'échelle réduit », normalisé à 1 à l'époque de référence, la combinaison de ces deux équations peut s'écrire
    $\left({\frac {1}{x}}{\frac {{\rm {d}}x}{{\rm {d}}(H_{0}t)}}\right)^{2}={\frac {1}{x^{3(1+w)}}}$.
Quand w = -1, on est dans le cas d'une constante cosmologique vu précédemment. Dans les autres cas, la résolution de cette équation donne immédiatement
    $x(t)=\left({\frac {3(1+w)}{2}}H_{0}|t|\right)^{\frac {2}{3(1+w)}}$.
Quand w est supérieur à -1, le facteur d'échelle tend vers 0 quand t tend vers 0. Il y a donc un Big Bang, et l'expansion se poursuit indéfiniment. À l'inverse, quand w est inférieur à -1, alors cette fois, le facteur d'échelle tend vers 0 quand t est négatif et tend vers l'infini, et il tend vers l'infini quand t tend vers 0. Dans ce cas, on a un univers existant de toute éternité, mais ayant une durée de vie finie dans le futur (voir ci-dessous).
Quand w est supérieur à -1, l'âge de l'univers s'écrit
    $t_{0}={\frac {2}{3(1+w)}}{\frac {1}{H_{0}}}$.
Les valeurs de w égales à 0 et 1/3 redonnent bien sûr les résultats précédents. Quand w est égal à -1, on a un univers sans commencement ni fin, et ne présentant pas d'évolution. Quand w est inférieur à -1, on a un univers sans commencement, mais atteignant une singularité gravitationnelle dans un temps fini, bien que restant toujours en expansion. En effet, quand le facteur d'échelle tend vers l'infini, alors la densité d'énergie diverge, et ce en un temps fini : c'est le Big Rip (litt. « grande déchirure »). Le temps qu'il reste à l'univers avant cette singularité est donné par
    $t_{f}={\frac {2}{3|1+w|}}{\frac {1}{H_{0}}}$.

## Solutions avec plusieurs espèces

Dans le cas où plusieurs espèces ou formes de matière emplissent l'Univers et participent notablement aux équations de Friedmann, il faut tenir compte de leurs équations d'état respective et de leur abondance relative. Dans le cas de deux espèces, par exemple, la première équation de Friedmann se réécrit
    $3\left({\frac {H^{2}}{c^{2}}}+{\frac {K}{a^{2}}}\right)={\frac {8\pi G}{c^{2}}}\left(\rho _{1}+\rho _{2}\right)$.
Si l'on appelle w$_1$ et w$_2$ les rapports de la pression à la densité d'énergie de chacune des deux espèces et qu'on les suppose constants au cours du temps, on peut réécrire cette équation à l'aide de diverses variables adimensionnées, sous la forme
    $\left({\frac {{\rm {d}}x}{xH_{0}{\rm {d}}t}}\right)^{2}=\left[\Omega _{1}^{0}x^{-3(1+w_{1})}+\Omega _{2}^{0}x^{-3(1+w_{2})}\right]$,
où l'on a défini, à partir d'un époque de référence notée avec l'indice ou l'exposant 0 la densité critique $\rho _{\rm {c}}^{0}$ par
    $\rho _{\rm {c}}^{0}={\frac {3c^{2}H_{0}^{2}}{8\pi G}}$,
ainsi que les quantités sans dimension appelées paramètres de densité $\Omega _{i}^{0}$ par
    $\Omega _{i}^{0}={\frac {\rho _{i}^{0}}{\rho _{\rm {c}}^{0}}}$,
et le facteur d'échelle normalisé, x, par $x=a/a_{0}$.
En prenant une époque de référence, que l'on note avec l'indice ou l'exposant 0, on a
    $3\left({\frac {H^{2}}{c^{2}}}+{\frac {K}{a^{2}}}\right)={\frac {8\pi G}{c^{4}}}\left[\rho _{1}^{0}\left({\frac {a_{0}}{a}}\right)^{3(1+w_{1})}+\rho _{2}^{0}\left({\frac {a_{0}}{a}}\right)^{3(1+w_{2})}\right]$.
On introduit la densité critique $\rho _{\rm {c}}^{0}$ par
    $\rho _{\rm {c}}^{0}={\frac {3c^{2}H_{0}^{2}}{8\pi G}}$,
ainsi que les quantités sans dimension appelées paramètres de densité $\Omega _{i}^{0}$ par
    $\Omega _{i}^{0}={\frac {\rho _{i}^{0}}{\rho _{\rm {c}}^{0}}}$.
Ainsi, on obtient, en introduisant le facteur d'échelle normalisé $x=a/a_{0}$,
    ${\frac {H^{2}}{H_{0}^{2}}}+{\frac {Kc^{2}}{H_{0}^{2}a_{0}^{2}}}{\frac {1}{x^{2}}}=\left[\Omega _{1}^{0}x^{-3(1+w_{1})}+\Omega _{2}^{0}x^{-3(1+w_{2})}\right]$.
En évaluant cette quantité à l'époque de référence, il vient alors
    $1+{\frac {Kc^{2}}{H_{0}^{2}a_{0}^{2}}}=\Omega _{1}^{0}+\Omega _{2}^{0}$,
d'où finalement
    $\left({\frac {{\rm {d}}x}{xH_{0}{\rm {d}}t}}\right)^{2}+\left(\Omega _{1}^{0}+\Omega _{2}^{0}-1\right){\frac {1}{x^{2}}}=\left[\Omega _{1}^{0}x^{-3(1+w_{1})}+\Omega _{2}^{0}x^{-3(1+w_{2})}\right]$.
Cette expression se simplifie quand la densité totale égale la densité critique, ou, de façon équivalente quand la somme des paramètres de densité vaut 1. On a alors
    $\left({\frac {{\rm {d}}x}{xH_{0}{\rm {d}}t}}\right)^{2}=\left[\Omega _{1}^{0}x^{-3(1+w_{1})}+\Omega _{2}^{0}x^{-3(1+w_{2})}\right]$.
Chacune de ces deux expressions se généralise à un nombre arbitraire de composantes. Des solutions exactes existent pour certaines valeurs des paramètres w$_1$ et w$_2$.
En particulier, dans le cadre du modèle standard de la cosmologie, l'Univers peut être décrit comme étant rempli de trois types d'espèce : de la matière relativiste (neutrinos et rayonnement), de la matière non relativiste (matière baryonique et matière noire), et de l'énergie noire, que l'on va ici approximer par une constante cosmologique. Les abondances relatives de ces espèces, qui varient avec le temps font qu'il n'arrive jamais qu'elles coexistent avec chacune des densités d'énergie significatives. À l'époque actuelle, on trouve essentiellement de l'énergie noire et de la matière non relativiste, alors qu'aux époques reculées, on trouvait de la matière relativiste et non relativiste, la constante cosmologique étant négligeable. Ceci vient du fait que les abondances de ces différents types de matière varient au cours du temps selon $1/x^{3(1+w)}$ : la densité d'énergie associée à la constante cosmologique reste constante, alors que celle de la matière relativiste ou non relativiste croît à mesure que l'on remonte vers le passé.

### Poussière et radiation

Quand on remonte vers le passé, les densités de matière non relativiste (« poussière ») ou relativiste (« radiation ») croissent comme x$^{-3}$ et x$^{-4}$ respectivement. On a donc
    $\left({\frac {{\rm {d}}x}{xH_{0}{\rm {d}}t}}\right)^{2}=\left[{\frac {\Omega _{\rm {m}}^{0}}{x^{3}}}+{\frac {\Omega _{\rm {r}}^{0}}{x^{4}}}\right]$,
avec
    $\Omega _{\rm {r}}^{0}=1-\Omega _{\rm {m}}^{0}$.
L'intégration donne
    ${\frac {4\left(1-\Omega _{\rm {m}}^{0}\right)^{\frac {3}{2}}}{3\left(\Omega _{\rm {m}}^{0}\right)^{2}}}\left[1-{\sqrt {1+{\frac {\Omega _{\rm {m}}^{0}x}{1-\Omega _{\rm {m}}^{0}}}}}\left(1-{\frac {1}{2}}{\frac {\Omega _{\rm {m}}^{0}x}{1-\Omega _{\rm {m}}^{0}}}\right)\right]=H_{0}t$.
On simplifie l'équation de départ en
    $\left({\frac {x{\rm {d}}x}{H_{0}{\rm {d}}t}}\right)^{2}=\left[\Omega _{\rm {m}}^{0}x+1-\Omega _{\rm {m}}^{0}\right]$,
soit
    ${\frac {x{\rm {d}}x}{\sqrt {\Omega _{\rm {m}}^{0}x+1-\Omega _{\rm {m}}^{0}}}}=H_{0}{\rm {d}}t$.
On effectue le changement de variable
    $u={\sqrt {\Omega _{\rm {m}}^{0}x+1-\Omega _{\rm {m}}^{0}}}$,
dont on tire
    $x={\frac {u^{2}-\left(1-\Omega _{\rm {m}}^{0}\right)}{\Omega _{\rm {m}}^{0}}}$,
    ${\rm {d}}x={\frac {2u{\rm {d}}u}{\Omega _{\rm {m}}^{0}}}$.
L'équation à résoudre devient
    ${\frac {2\left[u^{2}-\left(1-\Omega _{\rm {m}}^{0}\right)\right]{\rm {d}}u}{\left(\Omega _{\rm {m}}^{0}\right)^{2}}}=H_{0}{\rm {d}}t$,
qui s'intègre immédiatement en
    ${\frac {2}{\left(\Omega _{\rm {m}}^{0}\right)^{2}}}\left[{\frac {u^{3}}{3}}-\left(1-\Omega _{\rm {m}}^{0}\right)u\right]_{u_{i}}^{u}=H_{0}t$,
la quantité $u_{i}$ correspondant à la valeur $x=0$, soit ${\sqrt {1-\Omega _{\rm {m}}^{0}}}$. On trouve alors
    ${\frac {2}{\left(\Omega _{\rm {m}}^{0}\right)^{2}}}\left[{\sqrt {\Omega _{\rm {m}}^{0}x+1-\Omega _{\rm {m}}^{0}}}\left({\frac {\Omega _{\rm {m}}^{0}x+1-\Omega _{\rm {m}}^{0}}{3}}-\left(1-\Omega _{\rm {m}}^{0}\right)\right)+{\frac {2}{3}}\left(1-\Omega _{\rm {m}}^{0}\right)^{\frac {3}{2}}\right]=H_{0}t$,
soit
    ${\frac {4\left(1-\Omega _{\rm {m}}^{0}\right)^{\frac {3}{2}}}{3\left(\Omega _{\rm {m}}^{0}\right)^{2}}}\left[1-{\sqrt {1+{\frac {\Omega _{\rm {m}}^{0}x}{1-\Omega _{\rm {m}}^{0}}}}}\left(1-{\frac {1}{2}}{\frac {\Omega _{\rm {m}}^{0}x}{1-\Omega _{\rm {m}}^{0}}}\right)\right]=H_{0}t$.
Selon la valeur de x par rapport à ${\frac {1-\Omega _{\rm {m}}^{0}}{\Omega _{\rm {m}}^{0}}}$, qui correspond à l'époque de transition où densité de matière relativiste et non relativiste sont égales, on retrouve deux régimes précédemment étudiés. Pour les valeurs de x petites, on obtient
    ${\frac {1}{2}}x^{2}={\sqrt {\Omega _{\rm {r}}^{0}}}H_{0}t$,
ce qui correspond au résultat pour l'univers de radiation trouvé ci-dessus, la valeur de H$_0$ étant remplacée par ${\sqrt {\Omega _{\rm {r}}^{0}}}H_{0}$, correspond à la valeur qu'aurait la constante de Hubble en l'absence de matière non relativiste.
Dans l'autre régime, on trouve
    ${\frac {2}{3}}x^{\frac {3}{2}}={\sqrt {\Omega _{\rm {m}}^{0}}}H_{0}t$,
qui là encore redonne le résultat de l'univers de poussière à la correction sur la valeur de la constante de Hubble près.

### Poussière et constante cosmologique

Les seconds cas particuliers importants sont celui d'un univers dominé par de la matière non relativiste (w$_m$ = 0) et dont la constante cosmologique (w$_Λ$ = -1) devient dominante. C'est selon toute vraisemblance dans une situation similaire que se trouve notre Univers, l'accélération de l'expansion de l'univers attestant de l'existence d'une forme de matière se comportant de façon assez voisine à une constante cosmologique. Dans ce cas, l'intégration des équations de Friedmann donne, avec $\Omega _{\rm {m}}^{0}=1-\Omega _{\Lambda }^{0}$,
    ${\frac {a}{a_{0}}}=\left[{\sqrt {\frac {1-\Omega _{\Lambda }^{0}}{\Omega _{\Lambda }^{0}}}}\sinh \left({\frac {3{\sqrt {\Omega _{\Lambda }^{0}}}}{2}}H_{0}t\right)\right]^{\frac {2}{3}}$.
On reprend à nouveau la même équation avec les variables réduites. Dans un univers spatialement plat, on a
    $\left({\frac {1}{x}}{\frac {{\rm {d}}x}{{\rm {d}}(H_{0}t)}}\right)^{2}={\frac {1-\Omega _{\Lambda }^{0}}{x^{3}}}+\Omega _{\Lambda }^{0}$,
que l'on peut réécrire en
    $\left({\frac {{\rm {d}}x}{{\rm {d}}(H_{0}t)}}\right)^{2}={\frac {1-\Omega _{\Lambda }^{0}}{x}}+\Omega _{\Lambda }^{0}x^{2}$,
soit
    ${\frac {{\rm {d}}x}{\sqrt {{\frac {1-\Omega _{\Lambda }^{0}}{x}}+\Omega _{\Lambda }^{0}x^{2}}}}={\rm {d}}(H_{0}t)$,
ou bien
    ${\frac {x^{\frac {1}{2}}{\rm {d}}x}{\sqrt {(1-\Omega _{\Lambda }^{0})+\Omega _{\Lambda }^{0}x^{3}}}}={\rm {d}}(H_{0}t)$.
On pose
    $y=H_{0}t$,
    $u=x^{\frac {3}{2}}$,
ce qui permet de réécrire l'équation précédente en
    ${\frac {2}{3{\sqrt {1-\Omega _{\Lambda }^{0}}}}}{\frac {{\rm {d}}u}{\sqrt {1+{\frac {\Omega _{\Lambda }^{0}}{1-\Omega _{\Lambda }^{0}}}u^{2}}}}={\rm {d}}y$.
En posant
    $v={\sqrt {\frac {\Omega _{\Lambda }^{0}}{1-\Omega _{\Lambda }^{0}}}}u$,
d'où
    ${\frac {2}{3{\sqrt {\Omega _{\Lambda }^{0}}}}}{\frac {{\rm {d}}v}{\sqrt {1+v^{2}}}}={\rm {d}}y$.
Cette équation peut s'intégrer en
    ${\frac {2}{3{\sqrt {\Omega _{\Lambda }^{0}}}}}\operatorname {arsinh} (v)=y$,
que l'on inverse en
    $\sinh \left({\frac {3{\sqrt {\Omega _{\Lambda }^{0}}}}{2}}y\right)=v$,
soit
    $\sinh \left({\frac {3{\sqrt {\Omega _{\Lambda }^{0}}}}{2}}y\right)={\sqrt {\frac {\Omega _{\Lambda }^{0}}{1-\Omega _{\Lambda }^{0}}}}u$,
et finalement
    ${\frac {a}{a_{0}}}=\left[{\sqrt {\frac {1-\Omega _{\Lambda }^{0}}{\Omega _{\Lambda }^{0}}}}\sinh \left({\frac {3{\sqrt {\Omega _{\Lambda }^{0}}}}{2}}H_{0}t\right)\right]^{\frac {2}{3}}$.
Pour les temps petits, on peut effectuer un développement limité qui redonne le comportement en $t^{\frac {2}{3}}$ :
    ${\frac {a}{a_{0}}}\sim \left[{\frac {3}{2}}{\sqrt {\Omega _{\rm {m}}^{0}}}H_{0}t\right]^{\frac {2}{3}}$.
On retrouve exactement la formule déjà obtenue pour un univers de poussière, à ceci près que c'est ${\sqrt {\Omega _{\rm {m}}^{0}}}H_{0}$ qui apparaît au lieu de $H_{0}$. Ceci provient du fait qu'en l'absence de constante cosmologique, la constante de Hubble serait réduite d'un facteur ${\sqrt {\Omega _{\rm {m}}^{0}}}$. Le fait que l'on retrouve la formule de l'univers de poussière provient du fait qu'à ces époques, la constante cosmologique est négligeable par rapport à la densité de matière.
À l'inverse, pour les temps grands, on retrouve le comportement exponentiel du facteur d'échelle par rapport au temps :
    ${\frac {a}{a_{0}}}\sim \left({\frac {1-\Omega _{\Lambda }^{0}}{4\Omega _{\Lambda }^{0}}}\right)^{\frac {1}{3}}\exp \left({\sqrt {\Omega _{\Lambda }^{0}}}H_{0}t\right)$.
Dans ce cas, la constante de Hubble tend effectivement vers la valeur asymptotique ${\sqrt {\Omega _{\Lambda }^{0}}}H_{0}$ qui entre dans l'exposant. La transition entre ces deux régimes s'effectue au moment où l'argument du sinus hyperbolique est de l'ordre de 1. Quand cela se produit, la quantité $a/a_{0}$ est de l'ordre de $\left[\left(1-\Omega _{\Lambda }^{0}\right)/\Omega _{\Lambda }^{0}\right]^{\frac {1}{3}}$. Cela correspond au moment où la densité d'énergie associée à la constante cosmologique devient de l'ordre de celle de la matière. Avant cette époque c'est celle de la matière qui domine, et l'on retrouve les résultats de l'univers de poussière, après cette époque, c'est la constante cosmologique qui domine, et l'on retrouve la dynamique de l'univers de de Sitter.
Un point intéressant est de calculer la relation entre âge de l'univers et temps de Hubble $1/H_{0}$. On trouve ainsi
    $t_{0}={\frac {1}{H_{0}}}{\frac {2}{3}}{\frac {1}{\sqrt {\Omega _{\Lambda }^{0}}}}\operatorname {arsinh} {\sqrt {\frac {\Omega _{\Lambda }^{0}}{1-\Omega _{\Lambda }^{0}}}}$.
Pour des valeurs petites de $\Omega _{\Lambda }^{0}$, on retrouve la valeur habituelle de $t_{0}=2/3H_{0}$, alors qu'à mesure que $\Omega _{\Lambda }^{0}$ croît (en restant bien sûr inférieur à 1), l'âge de l'univers devient plus grand que le temps de Hubble.

## Autres écritures

Dans certains cas, on peut préférer exprimer les équations en fonction du facteur d'échelle et non du paramètre de Hubble. Il vient alors
    $3\left({\frac {1}{c^{2}}}{\frac {{\dot {a}}^{2}}{a^{2}}}+{\frac {K}{a^{2}}}\right)={\frac {8\pi G}{c^{2}}}\rho$,
    $-2{\frac {1}{c^{2}}}{\frac {\ddot {a}}{a}}-{\frac {1}{c^{2}}}{\frac {{\dot {a}}^{2}}{a^{2}}}-{\frac {K}{a^{2}}}={\frac {8\pi G}{c^{4}}}P$.
Le cas échéant, la seconde équation peut se réécrire en faisant disparaître le terme de courbure en utilisant la première, ce qui donne
    ${\frac {\ddot {a}}{a}}=-{\frac {4\pi G}{3}}\left(\rho +{\frac {3P}{c^{2}}}\right)$.
Cette dernière forme, qui donne l'accélération relative de deux objets distants du fait de l'expansion de l'univers est un cas particulier de l'équation de Raychaudhuri et pour cette raison parfois appelée ainsi. Plus généralement, toute écriture de la seconde équation de Friedmann faisant apparaître la dérivée seconde du facteur d'échelle (ou la dérivée première du taux d'expansion) peut être appelée ainsi.
Enfin, l'on peut également préférer extraire la constante cosmologique du contenu matériel de l'univers, lui donnant ainsi un rôle purement géométrique. Ce choix est aujourd'hui considéré comme peu opportun en cosmologie car la nature exacte de l'énergie noire est inconnue, mais correspond historiquement à celui d'Einstein et de Lemaître. On obtient alors, en considérant que la constante cosmologique est homogène à l'inverse du carré d'une longueur,
    $3\left({\frac {1}{c^{2}}}{\frac {{\dot {a}}^{2}}{a^{2}}}+{\frac {K}{a^{2}}}\right)-\Lambda ={\frac {8\pi G}{c^{2}}}\rho$,
    $-2{\frac {1}{c^{2}}}{\frac {\ddot {a}}{a}}-{\frac {1}{c^{2}}}{\frac {{\dot {a}}^{2}}{a^{2}}}-{\frac {K}{a^{2}}}+\Lambda ={\frac {8\pi G}{c^{4}}}P$.
Les écritures ci-dessus utilisent le temps cosmique. Il est possible et parfois utile de lui substituer le temps conforme, η, défini par la formule
    $a(\eta ){\rm {d}}\eta =c{\rm {d}}t\,$
ou encore
    $\eta (t_{0})=\int ^{t_{0}}c{\frac {{\rm {d}}t}{a(t)}}$
de sorte que l'élément de longueur devienne proportionnel à une métrique de Minkowski (on dit conformément minkowskien) :
    ${\rm {d}}s^{2}=a^{2}(\eta )\left({\rm {d}}\eta ^{2}-\gamma _{ij}{\rm {d}}x^{i}{\rm {d}}x^{j}\right)$.
Dans ce cas, on peut définir un « paramètre de Hubble conforme » ${\mathcal {H}}$ par
    ${\mathcal {H}}={\frac {1}{a}}{\frac {{\rm {d}}a}{{\rm {d}}\eta }}$.
Une telle quantité n'a pas d'interprétation physique immédiate, mais permet la réécriture des équations de Friedmann en termes du temps conforme :
    $3\left({\mathcal {H}}^{2}+K\right)={\frac {8\pi G}{c^{2}}}a^{2}\rho$,
    $-2\partial _{\eta }{\mathcal {H}}-{\mathcal {H}}^{2}-K={\frac {8\pi G}{c^{4}}}a^{2}P$.
La résolution de ces équations suit sensiblement les mêmes étapes que dans les cas précédents. En particulier, on trouve les dépendances suivantes du facteur d'échelle par rapport au temps conforme :
* $a\propto \eta$ dans un univers de radiation
* $a\propto \eta ^{2}$ dans un univers de poussière
* $a\propto -\eta ^{-1}$ avec une constante cosmologique
* $a\propto |\eta |^{\frac {2}{1+3w}}$ avec une forme de matière dont l'équation d'état est $P=w\rho$.
La relation
    ${\mathcal {H}}^{2}\propto a^{2}\rho$
a pour conséquence immédiate que
    ${\frac {{\rm {d}}a}{{\rm {d}}\eta }}\propto a^{2}{\sqrt {\rho }}$.
En utilisant de plus le fait que pour une équation d'état du type $P=w\rho$, la dépendance de ρ reste
    $\rho \propto a^{-3(1+w)}$,
il vient
    ${\frac {{\rm {d}}a}{{\rm {d}}\eta }}\propto a^{{\frac {1}{2}}-{\frac {3w}{2}}}$,
qui se résout en
    $a\propto |\eta |^{\frac {2}{1+3w}}$.
À noter que η varie de 0 à $+\infty$ selon une loi de puissance avec exposant positif tant que w est supérieur à -1/3, varie exponentiellement quand w vaut -1/3 et varie de $-\infty$ à 0 selon une loi de puissance avec exposant négatif quand w est inférieur à -1/3.
L'intérêt de résoudre les équations de Friedmann en termes du temps conforme vient du fait que le concept d'horizon des particules et d'horizon des événements est très étroitement relié à la relation $a(\eta )$, et notamment à son comportement pour les plus petites et plus grandes valeurs de η. De plus, les concepts de distance angulaire et de distance de luminosité dépendent eux aussi directement de cette même relation.

## Interprétation

Les équations de Friedmann en présence de matière non relativiste peuvent se retrouver (de façon heuristique) par un raisonnement purement newtonien. On peut en effet considérer l'évolution d'une sphère de matière, dont on suppose la densité constante à tout instant. Cette hypothèse est inexacte en général, mais la faire permet de se placer dans une situation assez semblable à celle d'un univers homogène et isotrope. auquel cas, le taux d'expansion de la sphère est relié à sa densité par la même formule que la première équation de Friedmann.
On considère donc une sphère de poussière (c'est-à-dire que l'on néglige les forces de pression) de masse M dont on étudie l'évolution du rayon R(t) en fonction du temps, sous l'effet de la gravité. Un point situé à la surface de la sphère est uniquement soumis aux forces de gravité, et obéit donc à l'équation
    ${\ddot {R}}=-{\frac {GM}{R^{2}}}$.
En multipliant par la dérivée de R, il vient
    ${\ddot {R}}{\dot {R}}=-GM{\frac {\dot {R}}{R^{2}}}$.
Les deux membres de l'équation correspondent à des dérivées par rapport au temps, que l'on peut intégrer en
    ${\frac {1}{2}}{\dot {R}}^{2}={\frac {GM}{R}}+E$,
la quantité E étant alors une constante d'intégration déterminée par les conditions initiales. Cette constante d'intégration n'est rien d'autre que l'énergie totale par unité de masse, égale à la somme de l'énergie cinétique et de l'énergie potentielle :
    $E={\frac {1}{2}}{\dot {R}}^{2}-{\frac {GM}{R}}$.
Si l'on remplace la masse de la sphère par le produit de son volume et de sa densité μ , on obtient
    ${\frac {1}{2}}({\dot {R}})^{2}={\frac {4\pi }{3}}G\mu R^{2}+E$,
que l'on peut réécrire en
    $3\left({\frac {{\dot {R}}^{2}}{R^{2}}}-{\frac {2E}{R^{2}}}\right)=8\pi G\mu$.
Enfin, en remplaçant la densité de masse μ par la densité d'énergie de masse ρ = μ c$^2$, on trouve
    $3\left({\frac {{\dot {R}}^{2}}{c^{2}R^{2}}}-{\frac {2E}{c^{2}R^{2}}}\right)={\frac {8\pi G}{c^{4}}}\rho$.
La différence cruciale entre ces deux approches vient du rôle joué par la courbure K (dans le modèle relativiste) qui s'apparente dans le modèle newtonien à une constante d'intégration sans signification géométrique. En relativité générale elle détermine les propriétés géométriques de l'espace.

## Dérivation des équations de Friedmann

Les équations de Friedmann ne sont rien d'autre que l'écriture des équations d'Einstein décrivant un univers homogène et isotrope. Leur dérivation ne pose pas de difficulté particulière, et elles représentent même une des solutions analytiques exactes les plus simples parmi celles connues à ces équations.
L'hypothèse de l'homogénéité et de l'isotropie des sections spatiales de l'univers permet d'écrire l'élément de longueur sous la forme
    ${\rm {d}}s^{2}=c^{2}{\rm {d}}t^{2}-a^{2}(t)\gamma _{ij}{\rm {d}}x^{i}{\rm {d}}x^{j}$,
où γ$_{ij}$ représente la métrique des sections spatiales, décrites par les coordonnées x$^i$. Un tel système de coordonnées est appelé coordonnées comobiles. Un observateur décrivant la trajectoire $x^{i}={\rm {Constante}}$ suit une géodésique. Il est appelé observateur fondamental. Son temps propre correspond ici exactement à la coordonnée t, appelée temps cosmique (le temps propre des observateurs fondamentaux). Une telle trajectoire correspond en première approximation aux trajectoires des galaxies, si l'on fait abstraction de leur mouvement propre. La distance entre deux observateurs fondamentaux augmente proportionnellement à la fonction a(t), qui est le facteur d'échelle. Une fois cette forme de la métrique supposée, les équations d'Einstein déterminent l'évolution temporelle du facteur d'échelle.
Les coefficients de la métrique s'écrivent
    $g_{00}=c^{2}$,
    $g_{0i}=0$,
    $g_{ij}=-a^{2}(t)\gamma _{ij}$.
La forme exacte de $\gamma _{ij}$ importe peu et dépend du type de coordonnées choisi (cartésiennes ou sphériques, par exemple, si on suppose les sections spatiales euclidiennes). Le seul résultat important est de savoir que les sections spatiales étant homogènes et isotropes forment un espace à symétrie maximale, et qu'un tel espace dont la métrique serait $\gamma _{ij}$ a un tenseur de Ricci donné par
    $^{3}R_{ij}=2K\gamma _{ij}$,
la quantité K étant la courbure spatiale associée.
L'inverse de la métrique a pour coordonnées :
    $g^{00}={\frac {1}{c^{2}}}$,
    $g^{0i}=0$,
    $g^{ij}=-{\frac {1}{a^{2}(t)}}\gamma ^{ij}$,
où $\gamma ^{ij}$ représente la métrique inverse de $\gamma _{ij}$.
Les symboles de Christoffel s'écrivent :
    $\Gamma _{00}^{0}=0$,
    $\Gamma _{0i}^{0}=0$,
    $\Gamma _{00}^{i}=0$,
    $\Gamma _{ij}^{0}=a^{2}{\frac {H}{c^{2}}}\gamma _{ij}$,
    $\Gamma _{0j}^{i}=H\delta _{j}^{i}$,
    $\Gamma _{jk}^{i}=^{(3)}\Gamma _{jk}^{i}$,
où H est le paramètre de Hubble, donne par $H={\rm {d}}a/a{\rm {d}}t$, et $^{(3)}\Gamma _{jk}^{i}$ correspond aux symboles de Christoffel associés à la métrique $\gamma _{ij}$.
Les coefficients du tenseur de Ricci s'écrivent alors
    $R_{00}=-3{\dot {H}}-3H^{2}$,
    $R_{0i}=0$,
    $R_{ij}=3a^{2}{\frac {H^{2}}{c^{2}}}\gamma _{ij}+a^{2}{\frac {\dot {H}}{c^{2}}}\gamma _{ij}+^{(3)}R_{ij}$,
La courbure scalaire s'écrit
    $R=-12{\frac {H^{2}}{c^{2}}}-6{\frac {\dot {H}}{c^{2}}}-6{\frac {K}{a^{2}}}$.
Le tenseur d'Einstein s'écrit enfin
    $G_{00}=3H^{2}+3{\frac {K}{a^{2}}}{c^{2}}$,
    $G_{0i}=0$,
    $G_{ij}=-\left(3{\frac {H^{2}}{c^{2}}}+2{\frac {\dot {H}}{c^{2}}}+{\frac {K}{a^{2}}}\right)a^{2}\gamma _{ij}$
Il faut maintenant déterminer le tenseur énergie-impulsion de la matière. Le seul tenseur énergie-impulsion compatible avec les hypothèses d'homogénéité et d'isotropie utilisées est celui d'un fluide parfait. Il est alors uniquement décrit par la densité d'énergie, la pression et la quadrivitesse du fluide selon la formule
    $T_{\mu \nu }=(P+\rho )u_{\mu }u_{\nu }-Pg_{\mu \nu }$.
Cette quadrivitesse correspond à celle des observateurs fondamentaux. La formule donnant la quadrivitesse étant
    $u^{\mu }={\frac {{\rm {d}}x^{\mu }}{{\rm {d}}\tau }}$,
et le temps propre étant donné par le temps coordonné, on a immédiatement
    $u^{0}=1$,
    $u^{i}=0$,
et les composantes contravariantes s'écrivent
    $u_{0}=c^{2}$,
    $u_{i}=0$.
Les composantes du tenseur énergie-impulsion sont donc
    $T_{00}=\rho c^{2}$,
    $T_{0i}=0$,
    $T_{ij}=Pa^{2}\gamma _{ij}$.
En utilisant la formule des équations d'Einstein, la composante 00 donne donc
    $3\left({\frac {H^{2}}{c^{2}}}+{\frac {K}{a^{2}}}\right)={\frac {8\pi G}{c^{2}}}\rho$,
et la composante ij donne
    $-3{\frac {H^{2}}{c^{2}}}-2{\frac {\dot {H}}{c^{2}}}-{\frac {K}{a^{2}}}={\frac {8\pi G}{c^{4}}}P$.