Cosmologie/Version imprimable

Ceci est la version imprimable de Cosmologie.
  • Si vous imprimez cette page, choisissez « Aperçu avant impression » dans votre navigateur, ou cliquez sur le lien Version imprimable dans la boîte à outils, vous verrez cette page sans ce message, ni éléments de navigation sur la gauche ou en haut.
  • Cliquez sur Rafraîchir cette page pour obtenir la dernière version du wikilivre.
  • Pour plus d'informations sur les version imprimables, y compris la manière d'obtenir une version PDF, vous pouvez lire l'article Versions imprimables.


Cosmologie

Une version à jour et éditable de ce livre est disponible sur Wikilivres,
une bibliothèque de livres pédagogiques, à l'URL :
https://fr.wikibooks.org/wiki/Cosmologie

Vous avez la permission de copier, distribuer et/ou modifier ce document selon les termes de la Licence de documentation libre GNU, version 1.2 ou plus récente publiée par la Free Software Foundation ; sans sections inaltérables, sans texte de première page de couverture et sans Texte de dernière page de couverture. Une copie de cette licence est incluse dans l'annexe nommée « Licence de documentation libre GNU ».

Le paradoxe d'Olbers

Nous allons commencer ce cours avec une question simple, qui est paradoxalement une excellente introduction à la cosmologie :

Pourquoi la nuit est-elle noire ?

Cette question peut sembler étrange, tant nous sommes habitués à voir les étoiles sur un fond noir. Bien que la lumière artificielle nous cache les étoiles et éclaire le ciel nocturne, nous savons que cette lumière cache le noir de la nuit, ses étoiles et la Lune. Mais pour un scientifique, même les choses évidentes sont dignes d'intérêt et demandent à être expliquées.

Les astronomes de l'antiquité n'étaient eux-mêmes pas surpris par l'obscurité de la nuit. Pour simplifier, ils pensaient les étoiles attachées à une sphère céleste, une sorte de sphère/coupole géante qui recouvre le ciel. S'ils supposaient que la Terre tourne sur elle-même, ils supposaient la sphère céleste fixe, expliquant le mouvement des étoiles au cours de la nuit et des saisons. Pour être plus précis, leur pensée utilisait plusieurs sphères : une pour le Soleil, une pour la Lune, une pour les étoiles proches, une pour les étoiles fixes des constellations (la plus lointaine) et bien d'autres. L'explication au noir de la nuit était que les étoiles étaient en nombre limitées et à des endroits bien précis, elles ne recouvraient pas entièrement le ciel. Beaucoup de savants avaient compris que ce modèle de sphères concentriques était un artifice de calcul utile, mais fondamentalement faux. Néanmoins, l'explication fondamentale au noir de la nuit restait valide : il n'y avait pas assez d'étoiles pour illuminer le ciel.

Terre dans la sphère céleste.
Sphères célestes emboîtées d'Anaximandre.

Avec les progrès de l'astronomie, et notamment l'arrivée des modèles géocentriques et héliocentriques, ce modèle de sphères emboîtées fût remis en cause. Pour les astronomes du moyen-âge et de la renaissance, les étoiles et astres n'étaient pas fixés sur des sphères proches, mais peuvent emplir tout l'espace. Leur vision de l'univers évolua vers un modèle d'univers infini, immobile, éternel, homogène. Par immobile, on veut dire que les effets de la gravité des astres s'annulent mutuellement sur de grandes distances. Chaque galaxie s'éloigne ou se rapproche de la nôtre, mais dans l'ensemble, ces différents mouvements se compensent et la moyenne des vitesses des galaxies est nulle. Par infini, on veut dire que l'univers a un volume infini : il n'a pas de début ni de fin, pas de frontières ou de bords, etc. Par éternel, on veut dire qu'il a un âge infini, qu'il n'a pas de début et de fin. Et par homogène, on veut dire que la répartition des étoiles est supposée globalement uniforme à grande échelle.

Mais les astronomes du moyen-âge et de la renaissance comprirent rapidement que le noir de la nuit était incompatible avec ce qu'ils savaient de l'univers.

Le paradoxe d'Olbers modifier

Le modèle d'un univers infini, immobile, éternel, homogène est assez logique, limite intuitive, mais elle entre en conflit avec le noir de la nuit. Plusieurs savant ont en effet montré qu'un tel univers devrait avoir une nuit extrêmement lumineuse, et certainement pas une nuit noire. Cela peut paraitre étonnant, mais plusieurs astronomes l'ont déduit mathématiquement ou sur la base d'arguments logiques. C'est un paradoxe, appelé paradoxe d'Olbers, nommé ainsi en l'honneur d'Olbers , un des astronomes à avoir évoqué ce paradoxe. D'autres astronomes avaient autrefois évoqué ce paradoxe, comme Thomas Digges en 15764, Johannes Kepler en 1605 ou Edmond Halley au XVIIIe siècle, mais pas aussi explicitement qu'Olbers.

Nombre d'étoiles visibles depuis la Terre en fonction de leur distance.

Logiquement, on s'attend à ce que la nuit soit noire car les étoiles lointaines nous envoient directement moins de lumière que les étoiles proches. Après tout, la lumière émise par un objet est répartie dans toutes les directions, ce qui fait que plus un objet est lointain, moins on reçoit sa lumière. La luminosité perçue depuis la Terre diminue avec la distance, et plus précisément avec le carré de la distance. Mais il se trouve que le diamètre apparent fait de même. Le résultat est que pour une portion du ciel d'une surface bien précise, on peut mettre soit peu d'étoiles proches, soit beaucoup d'étoiles lointaines. Il se trouve que sous certaines conditions, respectées dans l'univers infini, éternel et homogène, les deux phénomènes se compensent et que la luminosité d'une portion du ciel est globalement constante.

La lumière des étoiles se déplaçant dans le vide intersidéral en ligne droite, les étoiles sont visibles depuis la Terre tant qu'il n'y a pas d'obstacle entre elles et la Terre. Or, si l'univers est infini, éternel et homogène, alors on est certain qu'en tout point du ciel, on tombera sur au moins une étoile. Si l’univers est infini, il contient une infinité d'étoiles. Et s'il existe depuis un temps infini, la lumière de toutes ces étoiles a eu le temps d'atteindre la Terre. En conséquence, le ciel devrait être éclairé par une infinité d'étoiles. La répartition homogène des étoiles assure que ces étoiles soient réparties équitablement sur la sphère céleste, ce qui fait que peu importe où on regarde de la voûte céleste, on devrait y voir une étoile à cet endroit.

Avec ces conditions, la baisse de quantité de lumière reçue avec la distance est compensée par le plus grand nombre d'étoiles observables. Le schéma ci-contre illustre ce fait. Pour résumer, le ciel devrait être empli de lumière, cachant totalement le Soleil. C'était pour eux un paradoxe, appelé paradoxe d'Olbers, du nom de l'astronome qui le formalisa (d'autres astronomes, comme Kepler, avaient cependant mentionné ce paradoxe dans certains écrits, mais pas aussi explicitement qu'Olbers).

Pour mieux comprendre le problème, nous allons reprendre un développement mathématique assez connu qui formalise l'argument d'Olbers.

Modélisation de la répartition des étoiles modifier

En premier lieu, on découpe l'univers en sphères concentriques. Chaque couche, chaque sphère, a une épaisseur égale à la distance moyenne entre deux étoiles. De plus, chaque couche est à une distance de la Terre.

Illustration d'une couche sphérique dans la démonstration d'Olbers.

Calcul du nombre d'étoile dans chaque couche modifier

En second lieu, on calcule le nombre d'étoiles présentes dans chaque "sphère", dans chaque coque.

Pour cela, on utilise l'hypothèse suivante :

H1 : Les étoiles sont uniformément réparties.

Cette hypothèse dit que la densité d'étoile, à savoir le nombre d'étoiles par unité de volume, est constante. Dans ce qui suit, nous allons la noter .

Chaque "sphère", chaque coque, contient un nombre d'étoile égal à son volume multiplié par  :

Le volume d'une coque est approximativement égal à sa surface S multipliée par son épaisseur (la distance moyenne entre deux étoiles, notée D), ce qui donne :

En combinant les deux équations précédentes, on a :

Calcul du flux de lumière émis par les étoiles d'une couche modifier

En troisième lieu, on calcule le flux de lumière émis par une couche.

Pour cela, on part de l'hypothèse suivante :

H2 : Toutes les étoiles ont la même luminosité L.

On va utiliser cette hypothèse pour calculer le flux de lumière de chaque étoile. La physique du rayonnement nous dit que ce flux est égal à :

Le flux émis par toutes les étoiles d'une couche est la somme des flux de chaque étoile de la couche :

On utilise alors l'équation , démontrée plus haut :

On simplifie :

Cette équation nous dit que toutes les couches émettent la même quantité de lumière, ce qui n'est pas intuitif...

Calcul de la luminosité du ciel modifier

En quatrième lieu, on combine la luminosité de toutes les couches pour trouver la luminosité totale du ciel.

Pour cela, on utilise les hypothèses suivantes :

H3 : Toutes les étoiles sont visibles depuis la Terre, il n'y a pas d'obstacle entre une étoile et la Terre.

Cela signifie qu'il faut prendre en compte toutes les couches dans le calculs.

H4 : L'univers est infini.

Le fait que l'univers est infini nous dit qu'il existe une infinité de couches concentriques.

H4 : L'univers est éternel.

Le fait que l'univers est éternel nous dit que la lumière de toutes les étoiles a eu le temps d'atteindre la Terre.

Illustration du paradoxe d'Olbers, en ajoutant progressivement chaque couche d'étoiles.

Il faut donc additionner la luminosité de toutes les couches existantes pour trouver le flux de lumière visible dans le ciel, en faisant une intégrale. Nous n'allons pas faire le calcul, car il se trouve qu'on obtient un résultat infini. Intuitivement, la troisième étape suffit à comprendre pourquoi : chaque couche a une luminosité finie et identique à celle des autres couches, et il y a une infinité de couches. Le ciel devrait être infiniment lumineux !

Les "solutions" du paradoxe d'Olbers modifier

Divers savants se sont écharpés sur le paradoxe d'Olbers et de nombreuses réponses y ont été apportées. La démonstration d'Olbers n'ayant pas de problèmes mathématiques particuliers, il fallu se rendre à l'évidence : certaines hypothèses utilisées dans la démonstration sont fausses. Reste à trouver lesquelles.

On peut remettre en cause l'hypothèse de la luminosité fixe des étoiles, mais cela ne mène à rien. Même en supposant que certaines étoiles émettent peu ou pas de lumière, on se retrouve quand même avec une somme infinie.

Remettre en cause la répartition des étoiles a été envisagé, notamment dans des théories qui supposaient que la répartition des étoiles/galaxies était une distribution fractale. Ces théories sont cependant restées au stade de théories irréalistes, dont l'utilité est de montrer que telle piste de recherche est envisageable. La raison à cela est que les résultats empiriques ne suivent pas : dans ces théories, on devrait avoir des endroits du ciel qui seraient extrêmement lumineux, bien plus que le Soleil. Les observations astronomiques montrent de plus que si l'univers est assez hétérogène à petite distance, il est fortement homogène à grande distance.

Une première possibilité "crédible" était que les étoiles lointaines ne sont pas visibles depuis la Terre. Leur lumière est bien émise mais n'arrive pas à destination, reste à en trouver la raison. Une première explication fût que la lumière était absorbée par des nuages de gaz, de poussière, ou tout autre obstacle entre les étoiles et la Terre. Mais cette explication était incorrecte. L'obstacle en question absorbe le rayonnement de l'étoile, ce qui le chauffe. Or, tout corps chauffé émet un rayonnement dit de corps noir, dont l'intensité augmente avec la température. À force de chauffer, l'obstacle atteint une température d'équilibre, où le rayonnement absorbé par l'obstacle est intégralement réémis vers la Terre. L'obstacle est alors aussi lumineux que l'étoile qui le chauffe. Retour à la case départ.

La seule solution est que l'univers n'est pas infini et/ou qu'il a un âge fini. Les deux solutions font que l'on ne doit additionner que les couches les plus proches de la Terre, situées en-dessous d'un rayon maximal .

  • Si l'univers a un volume fini, il a de de facto un rayon maximal .
  • S'il a un âge fini , la lumière des étoiles très éloignées n'a pas eu le temps d'arriver sur Terre. Au-delà de la distance ( est la vitesse de la lumière), la lumière n'a pas encore pu arriver jusqu’à la Terre.



L'expansion de l'univers

Expansion de l'univers, image d'artiste.

Dès 1918, les astronomes ont entrepris de mesurer la vitesse de galaxies à partir de la lumière qu'elles émettent. Les observations actuelles utilisent soit des étoiles variables (des céphéides), soit des étoiles qui explosent : les supernovas. Plus précisément, les astronomes utilisent une classe bien précise de supernovas, qui ont pour particularité de générer systématiquement la même luminosité : les supernovas de type Ia. La luminosité perçue depuis la Terre de ces supernovas est proportionnelle à la distance. Et les observations ne collent pas du tout avec cette hypothèse ! Dans les grandes lignes, on observe que les galaxies s'éloignent plus qu'elles ne se rapprochent. La conclusion est claire : l'univers s'étend, gonfle.

La loi de Hubble modifier

Illustration de la relation entre la vitesse de fuite d'une galaxie et sa distance.

Hubble a été le premier astronome à mettre en équation ce comportement, dans son article daté de 1929. Il étudia un grand nombre d'observations provenant de ses collègues, ainsi que les observations qu'il avait effectuées lui-même. De ces observations, il induit une loi statistique, du nom de loi de Hubble. Cette loi dit que la vitesse d'éloignement d'une galaxie est proportionnelle à sa distance. Dit autrement, cette loi est identique à la formule qui suit, dans laquelle v est la vitesse de fuite, D la distance de la galaxie et H un facteur de proportionnalité nommée paramètre de Hubble.

Vous avez peut-être déjà entendu parler de la constante de Hubble pour désigner le paramètre de Hubble. C'était le terme utilisé au tout début de la cosmologie moderne, à une époque où on pensait qu'il s'agissait d'un paramètre constant. Mais il est aujourd'hui acquis que le paramètre de Hubble varie dans le temps. Ce n'est donc pas une constante, d'où le fait que nous parlerons de paramètre de Hubble dans de cours.

La valeur du paramètre de Hubble modifier

Mesures récentes du paramètre de Hubble.

Formellement, le paramètre de Hubble est égal à une vitesse divisée par une distance, ce qui fait qu'elle est l'inverse d'un temps. Le paramètre de Hubble H a donc pour unité , soit l'inverse d'une seconde. Mais dans les faits, les astronomes l'expriment en kilomètres par seconde par mégaparsecs (km/s/Mpc ou . Pour rappel, le mégaparsec est une unité astronomique qui vaut environ 3,26 millions d'années-lumière. Avec de telles unités, le paramètre de Hubble actuel, mesuré aujourd'hui, vaudrait approximativement :

Les mesures du paramètre de Hubble sont assez nombreuses, mais elles semblent converger vers la valeur de 70 mentionnée plus haut, avec cependant des incertitudes de mesure non-négligeables. Pour vous donner quelques exemples, les graphiques ci-dessous et ci-contre vous donnent les valeurs mesurées pour plusieurs compagnes d'observation assez anciennes. Vous voyez que la valeur n'est pas connue avec certitude.

Mesures du paramètre de Hubble par plusieurs campagnes scientifiques.

À l'heure actuelle, les cosmologistes ne savent pas expliquer pourquoi les mesures du paramètre de Hubble sont aussi différentes. Ce qui est sûr, c'est que le résultat dépend fortement de la méthode de mesure. Sans rentrer dans les détails techniques, il existe plusieurs méthodes indirectes pour estimer le paramètre de Hubble, qui estiment les distances des objets lointains en analysant la luminosité des supernovæ, des galaxies, des céphéides (des étoiles pulsatiles), et d'autres objets astronomiques. Et suivant la méthode utilisée ou l'objet observé, le paramètre de Hubble n'a pas la même valeur. Peut-être que les mesures sont entachées d'un biais systématique qui dépend de la mesure, peut-être que la valeur du paramètre de Hubble varie suivant la distance des objets considérés, peut-être que l'hypothèse d'un univers isotrope et homogène doit être abandonné, peut-être qu'une nouvelle physique se cache derrière ces résultats disparates, personne ne le sait.

Le temps de Hubble modifier

Comme dit plus haut, le paramètre de Hubble est l'inverse d'un temps. Ce temps en question est appelé le temps de Hubble et nous le noterons . Par définition, il vaut :

Ce qui se reformule en :

Il vaut environ 14,4 milliards d'années et cette valeur est très proche de celle actuellement admise par les scientifiques pour l'âge de l'univers (13,6 milliards d'années environ). Précisons que si les estimations actuelles nous disent que les deux sont assez proches, mais rien ne nous dit qu'ils sont exactement identiques. Les mesures et estimations sont en effet entachées d'une marge d'erreur assez large, ce qui fait que les deux mesures peuvent sembler se confondre alors qu'elles sont peut-être légèrement différentes.

Le temps de Hubble n'est pas l'âge de l'univers, attention à ne pas confondre ! Un bon moyen de s'en rendre compte est de regarder ce qui se passe dans un univers où serait constant, mais positif et non-nul. Le temps de Hubble serait alors constant, alors que l'univers serait en expansion.

La théorie nous dit que temps de Hubble et âge de l'univers n'ont pas de raison de coïncider. Dans la plupart des théories mathématiques de la cosmologie, et âge de l'univers ne coïncident qu'en un seul instant bien précis. À tous les autres instants, ces deux valeurs sont différentes. Dans de nombreuses théories cosmologiques, le temps de Hubble est plusieurs fois inférieur ou supérieur à l'âge de l'univers. La seule exception est un modèle théorique appelé le modèle (oui, le nom de cette théorie est bien une formule mathématique...), dans lequel le temps de Hubble et l'âge de l'univers se confondent à tout instant. Mais c'est l'exception qui confirme la règle. À l'heure actuelle, on ne sait pas si ce modèle décrit correctement l'univers actuel. Aussi, les scientifiques ne savent pas si la coïncidence actuelle entre temps de Hubble et âge de l'univers en est une ou est valide en permanence.

Le rayon de Hubble modifier

Les observations plus récentes ont toutes confirmé la loi de Hubble. Mieux, elles ont même montré que les galaxies vraiment lointaines semblent s'éloigner plus vite que la lumière. En effet, au-delà d'une certaine distance, le produit dépasse la vitesse de la lumière, si l'on en croit les observations. Cela peut sembler paradoxal, mais on verra dans ce chapitre qu'il n'en est rien. La distance à partir de laquelle les objets semblent aller plus vite que la lumière porte un nom : c'est le rayon de Hubble, noté . Par définition, on a :

On peut calculer le rayon de Hubble en utilisant la loi de Hubble . Pour cela, injectons celle-ci dans l'équation précédente :

Cette formule peut se réécrire comme ceci :

Ou encore :

Notons que cette dernière formule est une tautologie sans grande importance, sans sens physique. Elle vient de la manière dont on définit le temps de Hubble et le rayon de Hubble, mais ne dit rien sur l'interprétation physique du temps ou du rayon de Hubble. Dans la plupart des modèles cosmologiques, temps de Hubble et rayon associé n'ont pas de sens physique et ne sont que des conventions. Par exemple, le rayon de Hubble n'est en rien la limite au-delà de laquelle la lumière n'a pas eu le temps de nous parvenir (en termes scientifiques, ce n'est pas un horizon). On peut parfaitement voir des objets qui sont situés au-delà du rayon de Hubble, et les astronomes l'ont déjà fait. Bref, temps et rayon de Hubble sont des conventions. Un bon moyen de s'en rendre compte est, encore une fois, de regarder ce qui se passe avec un paramètre de Hubble constant, mais positif et non-nul. Le rayon de Hubble serait alors constant lui aussi, alors que l'univers serait en expansion.

L'expansion de l'univers et le facteur d'échelle modifier

Analogie de l'expansion de l'univers.

Aussi bizarre que cela puisse paraître, les scientifiques n'ont pas été étonnés de la découverte de la loi de Hubble. Il faut dire que la relativité générale était déjà bien avancée, et que les modèles d'univers en expansion ou en contraction étaient déjà étudiés à l'époque. De plus, cette équation avait été découverte quelques années auparavant par Georges Lemaître, un abbé féru de sciences, qui avait déduit cette relation des équations de la relativité générale. Les équations de la relativité expliquent la loi de Hubble avec le concept d'expansion de l'univers : les corps matériels de l'univers s'éloignent les uns des autres au fil du temps. Les interprétations de la relativité disent que l'expansion de l'univers ne provient pas d'un mouvement des objets dans l'espace, mais d'une modification de la manière de calculer les distances avec le temps. L'image qui est souvent donnée dans la vulgarisation scientifique compare l'univers avec un gâteau au raisin qui gonfle progressivement, les raisins étant les galaxies.

Expansion de l'univers.

Avec la loi de Hubble, il est évident que l'univers devait être plus "petit" par le passé (plus précisément, la portion de l'univers qui correspond aujourd'hui à l'univers observable). En renversant l'écoulement du temps, l'univers se contracte progressivement, et on peut facilement imaginer qu'après un certain temps, tout le contenu de l'univers soit rassemblé en un seul point : la singularité initiale. L'univers serait alors né d'une dilatation de cette singularité initiale, dilatation qui porte le nom de big-bang. Mais cette vue de l'esprit pose de nombreux problèmes mathématiques. En effet, cette singularité implique que de nombreux calculs dépendant des distances donnent des divisions par zéro. Par exemple, le calcul de la pression, de la température, de la densité, ou d'autres paramètres physiques ne sont pas calculables. Tout ce que peuvent faire les scientifiques, c'est étudier ce qu'il s'est passé quelques secondes ou minutes après le temps qui correspond à cette singularité hypothétique. Les calculs actuels ne donnent plus de résultats crédibles au-delà d'une certaine durée, la durée de Planck. Celle-ci vaut environ secondes.

Le facteur d'échelle modifier

Mettre en équation le phénomène d'expansion de l'univers est assez trivial. Du fait de l'expansion, toute distance entre deux points sera multipliée par un facteur multiplicatif après une durée t. Pour calculer ce facteur multiplicatif, les physiciens font intervenir ce qu'on appelle le facteur d'échelle, noté . Dans ce qui va suivre, nous allons supposer que nous prenons toutes les mesures dans un référentiel d'origine O, et que nous suivons la distance d'un objet matériel en fonction du temps. Nous allons comparer les distances entre un instant et un instant ultérieur. L'augmentation des distances à cause de l'expansion de l'univers se calcule comme suit :

Le rapport des facteurs d'échelle est un coefficient multiplicateur qui dit par combien les distances ont été multipliées entre l'époque actuelle et l'instant . Dit autrement, le facteur d'échelle est ce par quoi il faut diviser les distances actuelles pour obtenir les distances à l'instant .

Notons que le facteur d'échelle est sans dimensions (il n'a pas d'unité).

Pour simplifier les calculs, on considère souvent que le facteur d'échelle vaut 1 à un instant idéal pour simplifier les calculs. Sous cette hypothèse, l'équation précédente se simplifie :

L’interprétation de cette équation est assez simple : si le facteur d'échelle augmente de X %, les distances font de même.

Les distances propres et comobiles modifier

Différence entre distance comobile et distance propre. Sur ce schéma, les deux points sont immobiles, hors effet de l'expansion de l'univers. On voit que l'expansion éloigne les deux points en distance propre, si on mesure la distance avec une règle, sans tenir compte de la grille. Mais si on mesure la distance comobile, c'est à dire la distance mesurée avec les unités de la grille, la distance reste la même.

Reprenons la définition du facteur d'échelle vue précédemment :

Il est possible de réécrire la formule précédente comme suit :

Les distances et sont égale à une distance corrigée de l'influence des facteurs d'échelle. et . Elle peut s'interpréter comme la distance qu'auraient deux objets s'il n'y avait pas d'expansion. Elle est appelée la distance comobile. À contrario, la distance tient compte de l'expansion, qui augmente les distances entre deux objets. Les distances qui tiennent compte de l'expansion, opposées aux distances comobiles, sont appelées des distances propres. Par définition, les distances propres sont celles que l'on peut mesurer, qui ne sont pas corrigées de l'influence du facteur d'échelle.

Il existe une définition alternative de la distance comobile. C'est la distance mesurée à l'instant où .

Expansion de l'Univers illustré avec deux galaxies : on voit que la distance comobile montrée par la grille ne change pas, mais que la distance propre augmente.

L'augmentation des distances liée au facteur d'échelle n'est pas très intuitive. Dans le monde réel, une expansion a un centre, un point central d'où s'éloignent les autres. Mais avec l'expansion de l'univers, ce n'est pas du tout le cas. Peu importe où l'on soit dans l'univers, tous les autres objets semblent s'éloigner de nous. Un habitant de la Terre verra toutes les galaxies lointaines s'éloigner de la Terre, mais un habitant de la galaxie d'Andromède verra lui aussi l'ensemble des galaxies s'éloigner de lui et non de la Terre ! C'est cette particularité qui fait que l'on doit recourir à un facteur d'échelle pour expliquer l'expansion.

Schéma de l'expansion de l'univers. Les deux premiers schémas illustrent l'effet de l'expansion sur un ensemble de points, dont un point bleu et un point vert. Les deux schémas du bas montrent cette expansion du point de vue d'un observateur situé respectivement au point bleu, puis au point vert. On voit que chaque observateur voit les autres points s'éloigner de lui.

Le lien entre vitesse et facteur d'échelle modifier

La distinction entre distance comobile et propre peut aussi se faire pour les vitesses, volumes, surfaces et autres. Par exemple, il est possible de calculer une vitesse propre en dérivant la distance propre. Pour cela, on pourrait partir de la définition de la distance propre et en calculer la dérivée.

Mais les calculs seraient alors un petit peu longs, bien que pas difficiles. Pour éviter ce petit désagrément, nous allons ruser. À la place, nous allons calculer la dérivée de la distance comobile. La dérivée de la distance comobile est naturellement une vitesse, appelée la vitesse comobile. Elle correspond à la vitesse qu'aurait un objet s'il n'y avait pas d'expansion. Elle traduit le fait que les objets s'éloignent ou se rapprochent même sans expansion. Sachant que la distance comobile est définie par , on a :

On utilise la formule de la dérivée d'un quotient :

On simplifie par  :

Les deux termes : et sont techniquement des vitesses. Mais elles sont écrites en coordonnées comobiles, dans le référentiel où . Ce sont donc des vitesses que l'on ne peut pas mesurer. Pour passer dans le référentiel général, on doit multiplier par le facteur d'échelle :

Le terme n'est autre que la vitesse propre, la vitesse totale de l'objet qui incorpore les effets de l'expansion.

Cette équation peut se reformuler comme suit :

Le terme est ce qu'on appelle la vitesse locale. C'est la vitesse qu'a l'objet quand on retire l'effet de l'expansion, mais qu'on prend quand même en compte le facteur d'échelle. En effet, la vitesse comobile est mesurée dans un référentiel particulier, où . Et ce référentiel est situé arbitrairement dans le temps, pas forcément dans le temps présent. Pour obtenir la vitesse indépendante de l'expansion, on doit multiplier la vitesse comobile par , afin de tenir compte de la multiplication des distances au cours du temps.

On voit que la vitesse propre est la somme de deux vitesses. Une vitesse indépendante du facteur d'échelle et une autre qui dépend du facteur d'échelle. En clair, une vitesse indépendante de l'expansion et une qui y est proportionnelle. La première est une vitesse locale indépendante de l'expansion, alors que le second terme a pour origine l'expansion, d'où son nom de vitesse d'expansion.

Il faut noter que l'équation précédente nous explique pourquoi certaines galaxies très lointaines semblent s'éloigner de nous plus vite que la lumière. Si la vitesse locale ne peut dépasser la vitesse de la lumière, la vitesse de l'expansion n'est pas contrainte par . Ainsi, la vitesse supraluminique des galaxies lointaines provient de la vitesse de l'expansion et ne reflète pas une véritable vitesse supraluminique.

Le lien entre facteur d'échelle et paramètre de Hubble modifier

Pour les objets éloignés, la vitesse locale est négligeable par rapport à la vitesse d'expansion. Il est donc utile, de supposer la vitesse locale nulle. Sous cette hypothèse, l'équation précédente se simplifie comme suit :

Or, avec ces hypothèses, la vitesse se réduit à la vitesse de l'expansion, qui se calcule avec la loi de Hubble.

En combinant les deux équations, on trouve :

Comme on le voit, la démonstration précédente nous donne une nouvelle interprétation du facteur de Hubble : c'est la dérivée logarithmique du facteur d'échelle.

Pour rappel, la dérivée logarithmique d'une fonction se note et est définie par :

Elle porte ce nom car elle est égale à la dérivée du logarithme de la fonction initiale :

Cela nous permet de faire une remarque importante dans la suite du cours : l'intégrale d'une dérivée logarithmique est tout simplement le logarithme de la fonction. Nous ferons ce raccourci très souvent dans le cours.

Pour le dire plus clairement, le facteur de Hubble est le taux, du pourcentage auquel l'expansion a lieu. Pour information, la dérivée s’interprète comme la vitesse de l'expansion de l'univers, la vitesse à laquelle croît le facteur d'échelle. Plus la vitesse de l'expansion est grande, plus l'univers grandit vite et s'étend rapidement. De même, la dérivée seconde de est l'accélération de l'expansion de l'univers : plus elle est grande, plus l'expansion devient de plus en plus rapide avec le temps. Le facteur de Hubble est donc la vitesse de l’expansion divisée par le facteur d'échelle, soit le taux de variation du facteur d'échelle. Intuitivement, il indique approximativement si le facteur d'échelle augmente de 5 %, 10 % ou 20 % par unité de temps. Si H vaut 0,015, cela signifie que les distances augmentent de 1,5 % par seconde.

Faites attention à ne pas confondre la vitesse de l'expansion avec la vitesse d'expansion qui est la vitesse d'un objet acquiert à cause de l'expansion.

Maintenant, partons de l'équation précédente et intégrons-la :

On prend l'exponentielle et on réorganise l'équation :

Cette équation permet de calculer le facteur d'échelle quand on connait le facteur de Hubble. L'utilité de cette équation est qu'estimer le facteur de Hubble au cours du temps est possible, même si c'est par des moyens indirects. Les observations astronomiques permettent d'avoir une estimation précise de la valeur actuelle du facteur de Hubble, ainsi que de ses valeurs anciennes. Alors que l'évolution du facteur d'échelle ne l'est pas.

Le lien entre expansion et facteur de Hubble modifier

L'expansion, de part son action sur les distances, entraîne naturellement une variation des surfaces, volumes et densités. Prenons par exemple une sphère de rayon et de volume  : son rayon augmentant avec le facteur d'échelle, son volume fera de même. Quelques calculs triviaux nous disent que son volume évolue avec le facteur d'échelle selon la formule suivante, avec le volume de la sphère à l'instant .

Cette équation nous permet de déduire le rapport entre le facteur de Hubble et l'expansion des volumes. Pour cela, commençons par calculer la dérivée du volume :

Divisons ensuite par V, ce qui revient à diviser par :  :

On simplifie par et par  :

On utilise ensuite l'identité : pour simplifier le terme de droite, ce qui donne :

, qui peut aussi s'écrire de manière moins compacte comme ceci : .

Cette équation sera réutilisée plus tard dans le cours, quand nous démontrerons l'équation du fluide de Friedmann.

L'accélération de l'expansion de l'univers modifier

Il est intéressant de savoir si l'expansion se fait à vitesse constante, ou si l'expansion accélère/décélère. Pour simplifier les calculs, nous allons omettre la vitesse de la lumière et nous concentrer sur la vitesse de l'expansion de l'univers. Cela ne change rien aux résultats que nous allons obtenir vu que la vitesse de la lumière est constante : sa dérivée est donc nulle, ce qui la rend inutile dans les calculs de dérivée qui vont suivre.

L'accélération de l'expansion en fonction du paramètre de Hubble modifier

Une première étape pour savoir si l'expansion ralentit ou accélère, est de calculer l'accélération de l'expansion de l'univers. Cette accélération est simplement égale à la dérivée de la vitesse de l'expansion.

On peut alors appliquer la loi de Hubble pour déterminer la vitesse de l'expansion. En injectant dans l'équation précédente, on a :

Or, est simplement la vitesse d'expansion .

La vitesse de l'expansion peut se calculer avec la loi de Hubble, ce qui donne :

Le tout se simplifie en :

L'accélération de l'expansion en fonction du facteur d'échelle modifier

Une autre manière de quantifier l'accélération de l'expansion est de calculer la dérivée du facteur de Hubble. Le raisonnement derrière cette définition est assez simple. Le facteur de Hubble dit à quel taux les distances dans l'univers augmentent avec le temps. Si la dérivée est nulle, le facteur de Hubble est constant : l'expansion n’accélère pas plus qu'elle ne décélère. Inversement, une dérivée non-nulle nous dit comment le facteur de Hubble varie, et donc comment évolue l'expansion. Si la dérivée est positive, l’expansion accélère, et elle décélère pour le cas négatif.

Sachant que , la dérivée du facteur de Hubble est la suivante :

On utilise alors la formule

Ce qui se simplifie, en utilisant le facteur de Hubble :

On peut réorganiser les termes pour obtenir l'équation suivante :

On peut en profiter pour identifier cette équation avec l'équation , ce qui donne :

Le facteur de décélération modifier

Il est courant que les cosmologistes utilisent ce qu'on appelle le facteur de décélération, un nombre calculé à partir du facteur de Hubble. Celui-ci est positif si l'expansion est décélérée, négatif si elle accélère, et reste constant si la vitesse d'expansion reste constante. Par définition, ce facteur de décélération q vaut :

En utilisant la formule , la formule précédente se reformule comme suit :

On utilise alors la formule  :

Quelques manipulations algébriques donnent alors :

Les formules précédentes peuvent se réécrire sous la forme ci-dessous. Cette formule sera importante dans la suite du cours, car les deux termes et peuvent se calculer assez facilement. Dans les chapitres sur les équations de Friedmann, nous verrons que le premier terme se calcule à partir de la première équation de Friedmann et le second terme à partir de la seconde équation de Friedmann. Autant dire que l'étude des modèles cosmologiques est fortement facilitée quand on connaît cette relation.

Le lien entre temps de Hubble et facteur de décélération modifier

Une autre formule pour le facteur de décélération est la suivante :


Démonstration

La formule se démontre facilement si on calcule la dérivée du temps de Hubble :

En combinant avec les équations précédentes qui donnent le facteur de décélération en fonction de , on retrouve l'équation finale du facteur de décélération.

L'interprétation de cette équation est assez simple, si on regarde comment évolue la dérivée . Le cas où correspond au cas où le temps de Hubble augmente au même rythme que l'âge de l'univers, ce qui signifie que les deux sont égaux. Or, nous verrons dans quelques chapitres que ce n'est possible que si l'expansion de l'univers se fait à vitesse constante. Plus précisément, cela implique que le facteur d'échelle augmente linéairement avec le temps (), donc que sa dérivée première soit constante et sa dérivée seconde nulle. Les cas où et correspondent alors à une expansion supra- et infra-linéaire. Or, le facteur de décélération est définit de manière à valoir 0 pour un univers en expansion linéaire (à vitesse constante), positif pour une expansion supra-linéaire et négatif pour une expansion infra-linéaire. On voit que pour passer de la dérivée au facteur de décélération, il faut retrancher 1.

Une autre manière de réécrire cette formule, qui sera utile dans la suite du cours, est la suivante :

Le facteur de décélération moyen modifier

Rien n'implique que le facteur de décélération soit une constante. Il peut très bien varier dans le temps, si le facteur de Hubble varie lui aussi. Rien n’empêche d'avoir un univers dont l'expansion accélère, puis stoppe et décéléré, par exemple. Il est alors utile de calculer un facteur de décélération moyen, qui est définit par :

, avec l'âge de l'univers.

On utilise alors l'équation  :

L'intégrale d'une somme est la somme des intégrales :

On calcule la dernière intégrale :

On développe :

L'intégrale et la dérivée s'annulent si on néglige les constantes d'intégration, ce qui donne :

On voit que le facteur de décélération moyen dépend de l'âge de l'univers et du facteur de Hubble actuel, rien de plus. C'est un résultat très intéressant, qui permet soit de calculer le facteur de décélération moyen à partir de l'âge de l’univers, ou de faire l'inverse. L'âge de l'univers vaut donc :

Plus haut, nous avons dit que l'âge de l'univers qui fait actuellement consensus est proche du temps de Hubble. Si on en croit l'équation précédente le seul moyen d'avoir est que . En clair, l'univers a eu une expansion approximativement constante, en moyenne. Mais attention, cela ne signifie pas que l'expansion s'est faite de manière régulière tout le temps. Le consensus actuel est que l'univers a alterné entre décélération et d'expansion. Mais nous reparlerons de cela plus tard dans le cours, quand nous parlerons des modèles cosmologiques, de l'accélération de l'expansion de l'univers et de l'inflation.Toute la difficulté de la cosmologie est d'établir comment s'est déroulée l'expansion.


Le décalage vers le rouge (redshift)

Illustration du redshift. La lumière émise par l'objet est à gauche, alors que la lumière reçue par l'observateur est à droite. On voit que la lumière émise est perçue comme décalée, du point de vue du spectre, entre émetteur et observateur.

Si vous regardez une galaxie ou une étoile au loin, sa lumière semble être décalée vers les basses fréquences (le rouge), comparée à sa couleur d'émission. Le fait est que, comme on le verra juste après, il existe une relation entre la vitesse d'un objet et le décalage de la fréquence de sa lumière. La lumière est émise par l'objet à une certaine fréquence, mais la fréquence reçue par l'observateur n'est pas la même. C'est ce décalage vers le rouge qui était utilisé pour mesurer la vitesse des galaxies par Hubble et ses collègues. Ce phénomène s'appelle le décalage vers le rouge, noté .

Quantifier ce phénomène est assez facile. Pour cela, les physiciens utilisent le rapport entre le décalage des longueurs d'onde causé par l'expansion, et la longueur d'onde d'émission.

, avec la longueur d'onde actuelle, mesurée lors d'une observation, et la longueur d'onde lors de l'émission du rayonnement.

Cette quantité, notée , est appelée le décalage vers le rouge, ou encore le redshift. Dans ce qui va suivre, nous allons voir quelle est la relation entre la vitesse d'une galaxie et son décalage vers le rouge. Nous verrons ensuite que le redshift a un lien très fort avec l'expansion de l'univers.

Par souci de lisibilité, nous noterons parfois la dérivée première d'une variable , à savoir , comme ceci : . Le remplacement ne sera pas systématique, la notation étant plus courante et donc plus claire. La notation sera utilisée quand la notation est trop lourde, par exemple pour simplifier les formules de ce style : en .

La relation entre redshift et vitesse modifier

Les étudiants en physique apprennent que le décalage vers le rouge peut être causé par le mouvement d'un objet dans l'espace. Quand un objet s'éloigne de nous, à une certaine vitesse propre, la fréquence de la lumière qu'il émet baisse. Et réciproquement, un objet qui s'approche émettra une lumière plus bleuie, dont la fréquence a légèrement augmenté. C'est l'effet Doppler-Fizeau.

Illustration de l'effet Doppler-Fizeau.

En utilisant la physique de l'effet Doppler, on peut trouver une relation entre la vitesse de l'objet et le redshift. Suivant que l'on travaille en physique newtonienne, en relativité restreinte ou en relativité générale, les formules ne sont cependant pas les mêmes.

Le calcul du redshift en fonction de la vitesse modifier

La physique classique donne la formule suivante, entre la vitesse de l'objet en mouvement et le décalage vers le rouge qui en découle. Cette relation est cependant une approximation, dérivée des équations de Newton, valide uniquement pour des objets dont la vitesse d'éloignement est faible. En clair, cette formule ne vaut que pour des objets cosmologiques "proches". Les astres lointains vont avoir une vitesse élevée du fait de , ce qui fait que la formule n'est alors plus utilisable.

La relativité restreinte donne une formule encore plus précise, histoire d'obtenir des résultats corrects pour des objets lointains. Voici la formule en question :

, ou encore

Enfin la relativité générale nous donne une formule encore plus générale, qui est cependant rarement applicable telle quelle.

Le calcul de la vitesse en fonction du redshift modifier

On peut réécrire ces formules de manière à obtenir la vitesse de l'objet à partir de son redshift. La formule de la physique classique devient alors :

La relativité restreinte nous donne la formule suivante :

Enfin, la relativité générale nous donne la formule suivante :

La relation entre redshift et expansion de l'univers modifier

Relation théorique entre vitesse et redshift.

Appliquées de manière naïve à des galaxies lointaines, les formules précédentes donnent des vitesses supérieures à la vitesse de la lumière au-delà du rayon de Hubble. L’interprétation en termes d'effet Doppler impliquerait donc que les galaxies lointaines se déplacent plus vite que la lumière, en totale contradiction avec la relativité. La seule interprétation correcte de ce décalage vers le rouge cosmologique n'est donc pas une vitesse de déplacement dans l'espace, mais une modification du facteur d'échelle, valable quel que soit le référentiel. En clair, une partie du redshift est causée par l'effet Doppler proprement dit, mais une autre partie est liée à l'expansion de lumière. Cette dernière est un redshift un peu spécial dans le sens où il ne correspond pas à une vitesse propre, mais à une vitesse de fuite liée à l'expansion. Nous avons déjà vu dans les chapitres précédents que cette vitesse de fuite peut dépasser la vitesse de la lumière, vu que ce n'est pas une vraie vitesse, mais une modification de la géométrie de l'espace.

Le décalage vers le rouge cosmologique est la somme de deux effets : un effet Doppler et un redshift lié à l'expansion.

Une explication claire de ce processus demande d'utiliser le facteur d'échelle. Pour cela, on doit combiner étudier l'effet de l'expansion sur la longueur d'onde de la lumière.

L'effet de l'expansion sur la lumière modifier

Effect of the stretching of light on the light wavefront.

L'effet de l'expansion influe non seulement sur les distances entre corps matériels, mais aussi sur la lumière. La longueur d'onde de la lumière est une distance comme une autre, qui est modifiée par l'expansion de l'univers. Si une onde lumineuse est émise avec la longueur d'onde à un instant , sa longueur d'onde à un instant sera égale à :

Cela nous permet de calculer la fréquence d'une onde lumineuse en fonction du facteur d'échelle. En effet, il existe une relation entre la longueur d'onde et la fréquence pour la lumière (comme pour toute onde), les deux étant inversement proportionnelles. De cette relation, on peut déduire la relation suivante :

On peut utiliser cette relation entre longueur d'onde de la lumière et facteur d'échelle, pour calculer sa variation en fonction de l'expansion. Le fait est que le facteur d'échelle augmentant avec le temps, la longueur d'onde de la lumière augmente.

La relation entre redshift et facteur d'échelle modifier

Il est possible de démontrer une relation entre le facteur d'échelle et le redshift. Pour cela, partons de l'équation suivante :

On utilise alors l'équation , réécrite comme ceci :

En combinant les deux équations précédentes, on déduit la valeur du décalage vers le rouge en fonction du facteur d'échelle. Dans ce qui suit, on suppose que est l'instant d'émission de la lumière, alors que l'observation a lieu à l'instant .

Posons que le facteur d'échelle actuel vaut 1. L'équation obtenue est alors la suivante :

La relation entre redshift et paramètre de Hubble modifier

Quelques manipulations algébriques à partir des équations précédentes permettent d'exprimer le facteur de Hubble en fonction du redshift. Pour cela, partons de l'équation suivante :

Prenons la dérivée par rapport au temps :

Le calcul de la dérivée donne :

On applique la formule  :

On multiplie les deux côtés par

Maintenant, utilisons l'équation  :


L'univers observable

L'univers observable est la portion de l'univers que nous pouvons observer, compte tenu de la limite de la vitesse de la lumière. Des objets situés très loin ne peuvent pas être vus pour une raison très simple : la lumière qu'ils émettent n'a pas eu le temps de nous parvenir. La distance maximale à laquelle nous pouvons voir des objets (sans tenir compte d'éventuelles limitations techniques) dépend de l'âge de l'univers. S'il faut un temps supérieur à l'âge de l'univers pour nous parvenir, il nous est actuellement impossible de les voir, ce qui n'est pas le cas pour des objets situés plus près. Cette distance maximale est donc le rayon de l'univers observable. L'ensemble des points situés à la distance maximale des objets observables, à savoir la surface de l'univers observable, porte un nom : c'est l'horizon cosmologique.

Le rayon de Hubble modifier

On pourrait croire que le rayon de Hubble est une bonne estimation de la taille de l'univers observable. Pour rappel, ce rayon de Hubble est définit à partir du facteur de Hubble comme suit :

, avec c la vitesse de la lumière et le temps de Hubble définit par .

Il s'agit de la distance parcourue par la lumière pendant le temps de Hubble, lui-même définit comme l'inverse du facteur de Hubble. Mais le rayon de Hubble n'est malheureusement pas une bonne estimation du rayon de l'univers observable, pas plus que le temps de Hubble ne donne l'âge de l'univers. Pour calculer ces deux valeurs, il faut utiliser des développements mathématiques différents. La raison tient dans la définition du rayon de Hubble : c'est la distance à partir de la laquelle tous les objets s'éloignent de nous plus vite que la lumière. Dit autrement, c'est la distance au-delà de laquelle la vitesse propre tenant compte de l'expansion dépasse c, mais la vitesse comobile reste évidemment en-dessous. Le rayon de Hubble n'est en rien la limite de l'univers observable : on voit des objets sont situés au-delà et on peut calculer leur redshift ou tout autre mesure pertinente. Dans la plupart des modèles cosmologiques, temps et rayon de Hubble n'ont pas de sens physique.

L'évolution du rayon de Hubble modifier

On peut s'amuser à calculer dans quelles circonstances le rayon de Hubble est constant. Pour cela, prenons sa dérivée :

On utilise alors la formule qui relie la dérivée du temps de Hubble au facteur de décélération :

En combinant les deux équations précédentes, on trouve :

Cette dérivée s'annule quand . En analysant la dérivée, on remarque que le rayon de Hubble se contracte pour , augmente quand , et reste stationnaire pour .

Pour notre univers, tout semble indiquer qu'il augmente avec le temps. Et cela explique pourquoi on peut voir des objets qui se déplacent "plus vite que la lumière", dans le sens où l'expansion les éloigne de nous plus vite que la lumière. La raison est que, vu que le rayon de Hubble grandit, certains photons hors du rayon de Hubble finissent par y rentrer. Ils sont en quelque sorte rattrapés par le rayon de Hubble et nous finiront donc par les voir une fois arrivés sur Terre. Par contre, les objets qui ont émis cette lumière s'éloignent bien de nous plus vite que la lumière et restent au-delà du rayon de Hubble.

Prenons des objets situés au-delà du rayon de Hubble : la lumière émise par ces objets se déplace à la vitesse c, mais ces objets s'éloignent de nous plus vite que lumière du fait de l'expansion à une vitesse qui se calcule avec la loi de Hubble. La vitesse des photons venant vers nous, en coordonnées non-comobiles, est donc de . Pour les objets au-delà du rayon de Hubble, , ce qui fait que leur lumière ne devrait pas nous rejoindre si le rayon de Hubble restait statique, car elle est supérieure à c. Mais il faut maintenant prendre en compte que le rayon de Hubble grandit à la vitesse qu'on vient de calculer plus haut. Ce qui fait que ces photons dont se retrouvent quand même dans la sphère de Hubble, s'ils respectent la condition suivante :

Le rayon comobile de Hubble modifier

Notons que l'on peut calculer le rayon de Hubble comobile, c'est à dire en supprimant l'effet de l'expansion. Il est définit comme toute distance comobile, en divisant par le facteur d'échelle. Cela donne :

En utilisant la formule et en simplifiant, on trouve :

Ce rayon n'a pas plus de sens physique que le rayon de Hubble normal.

Le rayon de l'univers observable modifier

Le calcul du rayon cosmologique actuel est assez simple sur le principe, mais compliqué en pratique. Une méthode assez simple se base sur la vitesse d'éloignement de l'horizon, à savoir la vitesse à laquelle l'horizon cosmologique s'éloigne de nous. Une fois cette vitesse connue, il suffit de l'intégrer sur l'âge de l'univers.

La vitesse d'éloignement du rayon observable est définie par :

, avec le rayon de l'univers à l'instant t.

On réutilise l'équation du chapitre précédent, qui donne la vitesse en fonction de la distance: . Dans le cas de l'horizon cosmologique, la distance est égale au rayon de l'univers et la vitesse locale est égale à la vitesse de la lumière, ce qui donne :

On peut reformuler le tout en divisant par , ce qui donne :

On peut alors intégrer cette expression sur l'âge de l'univers pour obtenir le rayon de l'horizon cosmologique, du moins dans certaines conditions. Faire les calculs demande de connaître la loi d'évolution du facteur de Hubble . Mais dans les faits, elle nous est inconnue et on n'en a que quelques bribes. La théorie ne nous est pas d'un grand secours et le seul cas que l'on peut calculer simplement est celui où le facteur de Hubble est constant. En général, les modèles cosmologiques les plus simples supposent que le facteur d'échelle suit une loi de puissance, de type , ce qui fait que l'équation peut se résoudre avec quelques développements analytiques. Mais des modèles plus réalistes ne suivent pas vraiment une loi de puissance, ce qui complique les calculs. Pour nous éviter de longs calculs fastidieux, nous allons étudier le cas général, en utilisant quelques raisonnements astucieux.

Le rayon comobile de l'univers observable modifier

Une autre façon de faire les calculs est de passer par l'intermédiaire du rayon comobile. Pour rappel, ce rayon comobile est le rayon corrigé de l'influence du facteur d'échelle (et donc de l'expansion). Il vaut, par définition : .

La vitesse comobile est la dérivée de ce rayon comobile, qui est égale à :

On peut alors factoriser le rayon comobile  :

On simplifie :

En injectant l'équation dans la précédente, on a :

En développant, on trouve :

On voit que la dérivée est égale à ce qu'on appelle la vitesse comobile de la lumière. Par définition, la vitesse de la lumière est de , mais il s'agit d'une vitesse propre. On peut calculer sa vitesse comobile en divisant par le facteur d'échelle, ce qui n'est autre que le premier terme de l'équation précédente. . En intégrant l'équation précédente sur l'âge de l'univers, on a la distance comobile de l'horizon, celle à laquelle se situait l'horizon cosmologique quand la lumière de l'horizon a été émise. On se retrouve alors avec une équation très générale, qui marche même quand le facteur de Hubble est variable.

Seule la vitesse comobile de la lumière devant être prise en compte. On trouve alors que le rayon comobile se calcule avec la formule suivante, avec l'âge de l'univers :

Il est possible de simplifier fortement la formule précédente en faisant intervenir une variable particulière, appelée le temps conforme, définie par . On peut voir ce temps conforme comme l'équivalent de la distance comobile, mais pour les durées. Il n'a pas vraiment d'interprétation physique digne de ce nom et sert plus d'intermédiaire pour les calculs. Avec le temps conforme, la formule précédente devient alors :

On peut voir la formule précédente comme la généralisation de la formule , mais où le temps est remplacé par le temps conforme.

Le rayon actuel de l'univers observable modifier

On peut obtenir la distance propre par un calcul très simple, à partir de la distance comobile. Le passage de la distance comobile à la distance propre se fait simplement en multipliant par . On obtient alors la distance propre suivante :

Dans l'équation précédente, on peut factoriser la vitesse de la lumière, ce qui donne :

Comme pour le rayon comobile, on peut voir la formule précédente comme la généralisation de la formule , mais où la durée est remplacée par la valeur temporelle .

Les intégrales précédentes ne sont pas solubles telles quelles. Il faut préciser comment le facteur d'échelle évolue avec le temps, sans quoi les intégrales ne peuvent pas être calculées exactement. Manque de chance, la loi d'évolution du facteur d'échelle nous est inconnue. Nous sommes obligés de postuler des fonctions particulières, en espérant qu'elles collent au mieux aux observations. Nous ferons cela plus en détail dans le chapitre sur les modèles cosmologiques.


L'évolution de la matière

L'univers est peuplé de matière et de rayonnement. La matière est essentiellement composée de particules massives : baryons, quarks, électrons, etc. Pour plus de simplicité, on peut supposer que la matière de l'univers est un gaz. Cette hypothèse n'est pas si abusive vu l'état actuel de l'univers : 10% de la matière sert à fabriquer des étoiles, le reste étant localisé dans des nébuleuses et des nuages moléculaires dont la température ne dépasse pas la dizaine de degrés au-dessus du zéro absolu. L'univers est donc essentiellement composé de gaz. Évidemment, ce gaz de matière a une pression, un volume, une densité, une énergie, etc. Il reste de plus soumis aux lois de la thermodynamique. L'expansion va cependant faire varier continument sa densité (l'univers se dilue avec l'expansion), sa pression, sa température, etc. Dans ce chapitre, nous allons voir comment évolue ce gaz cosmologique en utilisant les relations de la thermodynamique, sous la contrainte de la loi de Hubble.

Rappels sur le gaz parfait modifier

Dans ce chapitre, et dans tous les chapitres suivants, nous allons supposer que la matière qui remplit l'univers est ce qu'on appelle un gaz parfait. Pour rappel, un gaz parfait est un gaz qui respecte la loi suivante :

, avec P la pression, V le volume, N le nombre de particules, T la température et la constante de Boltzmann.

La densité de particules modifier

On peut simplifier cette équation en utilisant la densité de particule, à savoir la quantité de particule présente par unité de volume, qui sera notée . Si on note V un volume et N le nombre de particules qu'il y a dedans, la densité de particules est définie par la formule suivante :

En combinant les deux équations précédentes, on a :

, avec la densité de particules.

L'expansion affecte surtout les termes de droite, à savoir la densité de particule et la température. Cette influence se répercute ensuite sur la pression. Notons que cette équation fonctionne pour un gaz d'atomes, de molécules, mais aussi pour décrire le rayonnement, la lumière, qui n'est autre qu'un gaz parfait de photons ! Mais avant d'étudier le cas du rayonnement, ce qui sera le fait du prochain chapitre, nous allons étudier le cas d'un gaz parfait composé de matière, de particules avec une masse, soumis à l'expansion.

L'énergie moyenne par particule d'un gaz parfait modifier

La température d'un gaz parfait est la somme des énergies cinétiques de chaque particule. La physique statistique nous donnent des formules pour relier la température et l'énergie cinétique moyenne d'une particule du gaz. Les formules en question dépendent du gaz, mais elles prennent toutes la forme d'une équation de la forme :

Le coefficient est appelé la capacité calorifique spécifique à volume constant. C'est l'énergie qu'il faut pour augmenter d'un degré la température du gaz, en considérant que son volume et sa masse sont constantes. La capacité calorifique est considérée comme constante avec un gaz parfait. Mais dans le monde réel, elle n'est pas totalement constante et peut varier suivant la température et la pression. Les variations les plus soudaines de la capacité calorifique impliquant souvent des transitions de phase elles-mêmes liées à des baisses/hausses de température.

A partir de l'équation précédente, on peut calculer la densité d'énergie et l'énergie moyenne par particule. L'énergie moyenne par particule se calcule en divisant l'équation précédente par N, ce qui donne :

En supposant que la capacité calorifique est constante, on peut reformuler cette équation comme suit :

On peut aussi calculer l'énergie moyenne par particule à partir de l'équation . Il suffit de la diviser par N, ce qui donne :

Le terme le plus à droite n'est autre que l'inverse de la densité de particule. En faisant le remplacement, on a :

La densité d'énergie d'un gaz parfait modifier

Partons maintenant de la formule suivante :

La formule précédente, une fois divisée par le volume, donne ceci :

En clair, la densité d'énergie, l'énergie moyenne par unité de volume, est proportionnelle à la pression. On peut alleer plus loin en reformulant l'équation précédente en fonction de la température. En combinant l'équation précédente avec l'équation , on a :

Le cas d'une expansion isotherme modifier

Dans ce qui suit, nous allons supposer que le nombre de particules est conservé. En clair, le nombre de particules matérielles dans l'univers est une constante : de nouvelles particules ne peuvent pas apparaître, pas plus que des particules existantes ne peuvent disparaître. Dans les faits, la constance du nombre de particules est quelque peu fausse : les lois de l’infiniment petit permettent l'apparition ou la disparition de particules, des conversions entre particules, et bien d'autres réactions. Et si ces réactions conservent la masse et l'énergie, la conservation du nombre de particules est loin d'être acquise. Par exemple, les réactions nucléaires peuvent faire disparaître des particules (tout en conservant la masse et l'énergie). Ou encore, les photons peuvent être absorbés ou émis entre interagissant avec la matière. Mais cela ne signifie pas que la constance des particules n'est pas une bonne approximation. Au global, on suppose que les réactions qui font apparaitre des particules s'équilibrent a peu-près avec celles qui en font disparaitre.

Vu que le volume de l'univers varie selon , la densité de particule varie donc comme :

La densité de particule est donc reliée au facteur d'échelle par l'équation

De même, on suppose que l'énergie est constante. L'énergie est conservée, ce qui parait raisonnable pour un univers remplit uniquement de matière. Cette hypothèse est clairement plus raisonnable que la constance du nombre de particule, la conservation de l'énergie étant vraiment importante en physique, là où faire apparaitre ou disparaitre des particules est fait régulièrement dans l'importe quel accélérateur de particule.

L'hypothèse de conservation de l'énergie, combinée à la conservation du nombre de particule, a des conséquences immédiates. En effet, cela impacte l'énergie par particule. Vu que l'énergie est conservée en même temps que le nombre de particule, pas de changement. L'énergie par particule n'est pas censée varier, vu que énergie et nombre de particules sont constants. Et vu que l’énergie par particule est proportionnelle à la température, alors la température est elle aussi constante. Dans un univers dont l'énergie est conservée, la température moyenne reste constante ! Vu que toutes les équations précédentes dépendent uniquement de la densité d'énergie et de la température, cela signifie qu'il nous reste à trouver le comportement de la densité de particules pour déterminer comment évolue le gaz parfait.

On peut maintenant déterminer la pression, maintenant que l'on connait le comportement de la température et de la densité de particule. Pour rappel, on a . Sachant que la température est constante, mais que la densité de particule diminue, alors la pression diminue avec l'expansion. On a alors :

Sachant que la densité d'énergie est proportionnelle à la pression, elle doit évoluer en elle aussi.

Pour résumer, on a :

  • une énergie par particule constante ;
  • une température constante ;
  • une densité d'énergie qui diminue comme  ;
  • une pression qui diminue comme .

Le cas d'une expansion adiabatique, où la chaleur est conservée modifier

Dans ce qui précédait, nous sommes parti du fait que l'énergie était conservée, et avons étudié l'évolution de la pression et de la température. Dans cette section, nous allons supposer qu'elle ne l'est pas, et nous allons établir des équations plus générales sur les gaz parfaits en expansion cosmologique. Nous allons étudier le cas d'une expansion adiabatique, à savoir dans laquelle il n'y a pas de conversion d'énergie en chaleur. La quantité d'énergie thermique de l'univers reste la même, mais pas forcément sa température. Il s'agit d'une modélisation simple, mais réaliste de l’expansion. En théorie, l'univers ne s'étend dans rien et l'expansion du gaz parfait qu'il contient ne pousse aucune paroi. Dans ces conditions, la pression du gaz n'exerce aucun travail susceptible de transformer sa chaleur en énergie cinétique. Son énergie thermique, sa chaleur, est donc globalement conservée si on omet les réactions chimiques.

Partons de l'équation précédente :

, avec un coefficient constant.

Si on suppose que N est constant, la dérivée de l'équation précédente donne :

Or, la thermodynamique nous dit qu'un gaz parfait en expansion adiabatique respecte la formule suivante :

En identifiant les deux équations précédentes, on trouve :

On peut reformuler la loi des gaz parfaits comme suit, en isolant la température :

En combinant les deux équations précédentes, on trouve :

Cette expression dit que la température et le volume évoluent globalement de la même manière.

L'évolution de la température lors de l'expansion modifier

On utilise alors la relation  :

On divise des deux côtés par k :

On prend alors l'intégrale, en utilisant la formule  :

On note et on réorganise les termes constants :

On applique la formule  :

On prend l'exponentielle des deux côtés :

L'exponentielle d'une somme est le produit des exponentielles () :

On pose :

Ou encore :

Pour un gaz monoatomique, on a , ce qui fait que la température de la matière varie donc comme l'inverse du carré du facteur d'échelle :

La relation entre pression et facteur d'échelle modifier

On peut combiner l'équation précédente avec l'équation vue plus haut qui donne la pression en fonction de la température et du facteur d'échelle :

On peut combiner l'équation précédente avec l'équation qui donne la pression d'un gaz parfait en fonction de sa température, vue précédemment :

Le tout se simplifie en :

L'évolution de la densité d'énergie lors d'une détente adiabatique modifier

Partons de l'équation de la conversation de l'énergie lors d'une détente adiabatique :

, avec E son énergie, P sa pression et V son volume.

On divise alors par dt :

Par définition, , avec la densité d'énergie. De plus, on a vu que . En faisant le remplacement, on a :

On utilise alors la formule du produit d'une dérivée sur le terme de gauche :

On réorganise les termes :

On simplifie :

En divisant par , on a :

On utilise alors la relation  :

Le facteur de Hubble étant égal à , l'équation précédente se réécrit comme suit :

Intégrons des deux côtés. Vu que  :

Ce qui signifie que, d'après les règles des puissances :

De cette équation, on peut déduire que :



L'évolution du rayonnement

La majorité du rayonnement dans l'univers forme un gaz de photons homogène et isotrope (nous détaillerons cette remarque quand nous aborderons le rayonnement de fond diffus cosmologique). Autrement dit, c'est un rayonnement de corps noir. Ce gaz de photons possède diverses propriétés, qui sont affectées par l'expansion : il possède une densité d'énergie, une température, et bien d'autres paramètres. Formellement, le rayonnement forme un gaz parfait, ce qui fait qu'on peut réutiliser les équations du chapitre précédent. Pour un gaz de photons, on peut prouver que . Les équations du chapitre précédent donnent alors :

(énergie par particule)
(densité d'énergie)

Pour les équations qui dépendent du facteur d'échelle, on a :

Dans la suite de ce chapitre, nous allons expliquer plus en détail d'où proviennent ces équations, ce qu'elles signifient physiquement et en donner des démonstrations alternatives.

L'évolution de la température du rayonnement avec l'expansion modifier

L'équation principale de ce chapitre dit que la température du rayonnement diminue au même rythme que l'augmentation du facteur d'échelle. Ainsi, le rayonnement né dans les premiers instants de l'univers, refroidit progressivement au fur et à mesure que l'univers s'étend.

Malheureusement, l'équation est rarement utilisable telle quelle. En effet, les facteurs d'échelle et ne sont pas connus et ne peuvent pas se mesurer. L'idéal est de remplacer les facteurs d'échelle par une grandeur qui se mesure. Le redshift fonctionne à merveille, d'autant qu'on sait qu'il est relié au facteur d'échelle par les équations vues il y a quelques chapitres. Dans ce qui suit, on suppose que est l'instant d'émission de la lumière, alors que l'observation a lieu à l'instant . Avec cette convention, on sait que . En faisant le remplacement dans l'équation précédente, on trouve :

Précisons que cette équation vient de la convention . Mais avec la convention inverse, à savoir , on trouverait l'équation inverse, à savoir :

L'évolution de la densité d'énergie d'un gaz de photons avec l'expansion modifier

Illustration de la loi de Planck.

Un gaz de photons est formellement décrit par la fameuse équation de Planck, que nous ne détaillerons pas ici. Celle-ci donne tout simplement la densité d'énergie des photons qui possèdent une fréquence f, dans un gaz de photons de température T. Et cette équation contient un terme, la fréquence, sensible au facteur d'échelle. Voici cette équation de Planck :

La distribution des photons suivant leur fréquence est illustrée par le schéma de droite. Celui-ci montre que le pic d'intensité, à savoir la fréquence où le nombre de photons est maximal. Or, ce pic est proportionnel à la température moyenne du gaz de photon, d'après la relation suivante : .

En intégrant l'équation précédente sur toutes les fréquences, on trouve la fameuse loi de Stephan, qui donne la densité d'énergie d'un gaz de photons en fonction de sa température. Voici cette loi, avec une constante, la constante de Stephan, et la température :

On peut la reformuler comme suit :

Sachant que l'on a , on a :

Dans ce qui suit, nous noterons la quantité comme suit : . L'équation précédente s'écrit alors, d'une manière plus formelle :

Au passage, on aurait pu avoir l'intuition de ce résultat dès le début du chapitre, quand j'ai dit que pour le rayonnement, on a : . Vu que la pression est proportionnelle à la densité d'énergie, on aurait du en déduire que celle-ci évoluait en .

L'évolution du nombre de particules d'un gaz de photons modifier

Des équations précédentes, il est possible de déterminer le nombre de photons par unité de volume, appelée la densité de photons par analogie avec la densité de matière, et la manière dont elle évolue avec le facteur d'échelle. Pour la calculer, il suffit de diviser la densité d'énergie par l'énergie moyenne d'un photon. On obtient alors :

En simplifiant, on a :

Ce qui se reformule comme suit :

Sachant que l'on a , on a :

L'interprétation physique de la variation de densité énergétique du rayonnement modifier

On a vu que la densité d'énergie du rayonnement varie selon la puissance quatrième du facteur d'échelle :

On peut donner un sens physique à cette équation. Premièrement, l'énergie du rayonnement est diluée dans un volume plus grand, égal à la puissance troisième du volume initial. La densité est donc divisée par la puissance troisième. À cela, il faut ajouter la diminution de la longueur d'onde causée par le facteur d'échelle. La somme de ces deux contributions donne la formule précédente. Pour nous en rendre compte, on peut partir de la définition de la densité d'énergie du rayonnement :

La variation de la densité d'énergie provient de deux sources : l'une est la variation de et l'autre est la diminution de . La première varie comme , car le nombre de photons reste globalement constant et que l'expansion fait augmenter le volume en . Cela sera justifié dans la section suivante, mais le résultat est le même que celui vu dans le chapitre précédent. Si l'énergie par particule demeurait constante lors de l'expansion, on aurait l'équation . On le voit, il manque un facteur pour obtenir l'équation finale. On peut facilement deviner son origine quand on sait que l'énergie moyenne d'un photon dans un gaz de photons est approximativement de :

On applique alors l'équation

En posant , on a :

En clair, la densité d'énergie diminue avec le facteur d'échelle pour deux raisons : l'expansion dilue le rayonnement dans un volume plus élevé, et l'expansion réduit la température du rayonnement.

Un autre argument qualitatif (et peu rigoureux) nous permet de justifier pourquoi l'expansion réduit l'énergie par particule. Rappelons qu'un photon de fréquence f a une énergie égale à , avec h la constante de Planck. Or, on a vu que la fréquence varie avec l'inverse du facteur d'échelle, l'énergie par photon doit faire de même, ce qui donne : math>{E \over N} \propto f \propto a^{-1}</math>. En clair, l'expansion étire la longueur d'onde des photons, ce qui leur fait perdre de l'énergie. Cependant, cette dérivation n'est pas parfaite, vu qu'on mélange la fréquence d'un photon unique avec la température d'un gaz de plusieurs photons. Ce qui nuit à la généralité de l'argument.

Cela a une conséquence assez importante : l'énergie de l'univers ne se conserve pas, mais diminue avec le temps ! Et ce n'est pas un problème qui serait réglé en relativité générale : il y a réellement une perte d'énergie quel que soit le modèle utilisé. À l'heure actuelle, on ne sait pas comment résoudre ce problème (si tant est que ce soit vraiment un problème).


Les processus de baryogenèse et nucléosynthèse

Au tout début de sa formation, l'univers était clairement chaud et dense. Les températures quelques microsecondes après le big-bang dépassaient le million voire le milliard de degrés. L'univers était en première approximation un gaz parfait de particules très différentes : neutrons, protons, neutrinos, électrons, photons, quarks, et autres. Les températures étaient tellement fortes que toutes les populations de particules réagissaient entre elles. N'importe quelle particule pouvait échanger de l'énergie avec n'importe quelle autre, homogénéisant les températures. Le mélange était tel que l'on pouvait définir une température moyenne valable pour tous les types de particules : les neutrons avaient une température moyenne similaire à celle des protons, elle-même identiques à celle des quarks, etc. On dit que l'équilibre thermique est respecté.

Puis, l'univers s'est refroidi en raison de l'expansion. En se diluant avec l'expansion, sa densité a diminué, les particules se sont espacées entre elles. Elles se sont éloignées au point que leurs interactions sont devenues plus rares. Leur proximité rendait leurs interactions faciles, chaque particule ayant rapidement accès à une voisine pour échanger de la quantité de mouvement. Mais avec la dilution, les particules se sont progressivement isolées les unes des autres, rendant les échanges de plus en plus difficiles. De nombreuses interactions ont disparues, perturbant les processus de mélange thermique. De plus, la température a aussi fortement baissé. Rappelez-vous le résultat du chapitre précédent : la température du rayonnement diminue comme l'inverse du facteur d'échelle (). Le gaz de particules élémentaires s'est alors progressivement condensé, donnant naissance à des particules composites. D'une soupe de particules se sont ainsi formés les nucléons, puis les noyaux et enfin les atomes. Ce chapitre vise à expliquer comment s'est produite cette condensation.

La baryogenèse modifier

Au tout début, on pouvait voir l'univers comme un mélange de plusieurs gaz composés de particules élémentaires. Du temps des fortes températures, quelques micro-secondes avant le big-bang, les particules composites ne pouvaient pas se former à partir de quarks. La température trop intense faisait que les particules composites étaient brisées par le chaos ambiant quelques microsecondes après leur formation. C'était essentiellement les photons et neutrinos qui réagissaient avec la matière et brisaient les structures ainsi formées. Il a fallu attendre que la température du rayonnement baisse pour que les quarks puissent s'assembler en protons et neutrons sans interagir avec un photon qui passe sur le chemin. Ce processus de formation des protons et neutrons s'appelle la baryogenèse, ce qui signifie formation des baryons (les protons et neutrons sont des exemples de baryons, d'où le nom).

Le rapport protons/neutrons modifier

La théorie du big-bang nous permet de déterminer comment s'est produit ce processus. Une réussite de la théorie tient dans le fait qu'elle prédit le rapport entre le nombre de protons et de neutrons dans l’univers. Celui-ci peut se calculer à partir du raisonnement suivant. Avant que les noyaux se forment, les protons et neutrons étaient libres et formaient un plasma de nucléons. La température de ce plasma a diminué progressivement avec l'expansion. Peu avant la formation des noyaux, la température était faible comparée à la masse des protons et neutrons (). Dans ces conditions, le gaz peut être décrit par ce qu'on appelle la distribution de Maxwell-Boltzmann. Celle-ci dit que la quantité de particules d'énergie par unité de volume est de :

Ainsi, on peut calculer le rapport entre protons et neutrons. Il suffit de faire le calcul de la densité de protons, et de la densité de neutrons séparément, et de diviser le premier par le second :

Les protons et neutrons forment un plasma tant que protons et neutrons peuvent interagir. Il arrive notamment que des protons se transforment en neutrons et réciproquement. Ces transformations, des réactions nucléaires, ont une probabilité d’occurrence qui dépend de la température. Quand le produit descend en-dessous de 0.8 Mev, ces réactions deviennent de plus en plus rares, au point que l'équilibre thermique du plasma est brisé. Les quantités de protons et de neutrons sont alors figées, de même que le rapport de leurs densités volumiques. Les calculs donnent 6 protons pour 1 neutron : . Dit autrement, ème de la matière baryonique est sous la forme de neutrons, alors que ème sont des protons. Cela correspond à environ 12% de neutrons pour 88% de protons.

Par la suite, ce rapport va cependant évoluer à cause de désintégrations de neutrons en protons (désintégration bêta). Précisons que ces désintégrations n'ont lieu qu'en-dehors des noyaux atomiques, et ne peut pas toucher l’abondance des éléments chimiques. Par contre, elle change l'abondance des protons et neutrons libres, isolés des atomes. Ces désintégrations suivent la fameuse loi de désintégration radioactive , avec égal à 880,3 s. On a alors :

La nucléosynthèse primordiale modifier

Une fois les protons et neutrons formés, l'univers était rempli de protons, de neutrons, d'électrons et de neutrinos, qui formaient un gaz à haute température. Durant un temps assez court, protons et neutrons ne pouvaient pas s'assembler pour former des noyaux, la température brisant les noyaux qui avaient l'occasion de se former. Mais, la température diminuant, cela ne dura guère. Après un certain temps, protons et neutrons ont pu s'assembler pour former des noyaux, quand la température a atteint un certain seuil.

La formation des premiers noyaux porte le nom de nucléosynthèse primordiale. Un nom barbare assez simple à comprendre : nucléosynthèse veut dire "synthèse de noyaux atomiques", et primordiale pour dire qu'elle a eu lieu peu après le big-bang. Ce terme sert à la distinguer de la nucléosynthèse qui a lieu actuellement au cours des étoiles, la nucléosynthèse stellaire. Les différences entre les deux sont assez nombreuses. Déjà, la nucléosynthèse primordiale s'est faite sur un temps très court, d'à peine quelques secondes grand maximum, alors que la nucléosynthèse des étoiles est un processus continu qui dure durant plusieurs milliards d'années. Ensuite, la nucléosynthèse primordiale a majoritairement créé des éléments chimiques légers, mais guère plus. Elle a permis de fabriquer de l'hydrogène, de l'hélium, du béryllium et du lithium, mais pas plus. Les autres éléments chimiques ont été fabriqués ultérieurement par la nucléosynthèse stellaire, qui a donné naissance à du carbone, de l'oxygène, de l'azote, et d'autres noyaux lourds.

Isotopes de l'hydrogène.

Suite à la nucléosynthèse primordiale, la quasi-totalité de la matière est composée d’hydrogène et d'hélium : environ 3/4 d'hydrogène et 1/4 d'hélium, le reste étant présent en quantités négligeables. C'est pour cela que la quasi-totalité de la matière des étoiles et planètes est sous la forme d'hélium et d'hydrogène, des particules formées par l'assemblage de neutrons et de protons. Pour rappel, Un noyau d'hydrogène est formé d'un simple proton, le nombre de neutrons variant de zéro à deux neutrons. La plupart de l'hydrogène ne contient pas de neutron, cette forme d'hydrogène étant appelé du protium. L'isotope avec un neutron est appelé le deutérium, alors que celui avec deux neutrons est appelé le tritium. Le protium est de loin la forme d’hydrogène dominante, les autres formes n'étant présentes que dans les étoiles, rarement dans le milieu interstellaire. Quant à l'hélium, il possède deux protons, avec un nombre variable de neutrons. Ses deux isotopes les plus fréquents possèdent deux neutrons pour l'hélium-4, un seul pour l'hélium-3.

Les réactions de la nucléosynthèse primordiale modifier

Résumé des réactions nucléaires principales de la nucléosynthèse primordiale, sous forme de liste.

La nucléosynthèse commence avec la fabrication de deutérium, un des isotopes de l'hydrogène. C'est à partir du deutérium que peuvent s'enclencher les réactions qui donnent naissance au tritium, à l'hélium, au lithium et au béryllium.

Dans ce qui suit, on notera D un noyau de deutérium, T un noyau de tritium, p un proton et n un neutron.
Illustration des réactions nucléaires donnant naissance au deutérium, au tritium et à l'hélium (3 et 4).

Le deutérium se forme en fusionnant un proton avec un neutron. La formation du deutérium ne s'est produite qu'une fois la température suffisamment basse. Au-dessus de cette température, les noyaux de deutérium ne survivent pas bien longtemps, à cause de la photodissociation. Les photons énergétiques brisent ces noyaux en quelques microsecondes, ne laissant que des protons et des neutrons. Mais une fois que la température descend sous la température critique du deutérium, les photons ne sont plus assez énergétiques pour briser les noyaux de deutérium, qui survivent.

Une fois le deutérium formé, des réactions donnent naissance soit à de l'hélium-3, soit à du tritium. L'hélium-3 peut se former de deux manières : soit par addition d'un proton, soit par fusion de deux noyaux de deutérium.

Le tritium peut lui aussi se former de deux manières différentes. Dans le premier cas, il est formé par la fusion de deux noyaux de deutérium. Dans le second cas, il est formé à partir d'un noyau d'hélium-3, dans lequel on remplace un proton par un neutron (par désintégration bêta, ou par capture/émission de nucléons).

Une fois le tritium ou l'hélium-3 formé, l'hélium-4 peut enfin apparaître. Il se forme soit à partir d'hélium-3, soit à partir de tritium. Et les méthodes pour ce faire assez diverses. Dans les deux cas, il se forme en ajoutant un noyau de deutérium, suivi par l'émission d'un nucléon. Le nucléon émit est un proton pour la fusion avec l'hélium-3, un neutron pour la fusion avec le tritium. Une autre manière consiste à ajouter un proton à du tritium, ou du neutron à de l'hélium-3.

Enfin, les autres éléments légers peuvent se former à partir de l'hélium-4. En fusionnant de l'hélium-4 avec soit du tritium, soit de l'hélium-3, on obtient respectivement du lithium et du béryllium. Le lithium peut aussi se former à partir du béryllium, par remplacement d'un neutron en proton. De plus, le lithium peut fusionner avec un proton pour donner deux noyaux d'hélium-4.

L'ensemble de ces réactions est résumé dans le schéma ci-dessous.

Réactions nucléaires principales de la nucléosynthèse primordiale, en schéma.

Le calcul de l'abondance de l'Hélium modifier

Des calculs théoriques poussés, basés sur la physique nucléaire, nous disent que la concentration en éléments chimiques a évolué rapidement au cours de la nucléosynthèse primordiale, avant de stabiliser. Le résultat est que les deux éléments majoritaires sont l'hydrogène (le protium) et l'hélium-4. Le deutérium et l'hélium-3 sont plus rares et ont une concentration assez similaire. Le tritium et les neutrons libres sont eux encore plus rares. Enfin, le lithium et le béryllium ferment la marche et sont les éléments les plus rares. Pour simplifier, on peut dire que l'univers est rempli presque exclusivement d'hydrogène et d'hélium-4. Les autres éléments sont tellement rares qu'ils sont presque négligeables.

Évolution de la concentration de chaque élément chimique dans l'univers, au cours de la nucléosynthèse primordiale.

Sans recourir à ces calculs compliqués, on peut calculer l'abondance des éléments principaux. L'idéal est de se concentrer sur l'hydrogène (le protium) et l'hélium-4 uniquement. Négliger les autres éléments n'est pas un problème tant ils sont rares. En faisant cela, on doit considérer que tous les neutrons ont été capturés dans les noyaux d'hélium-4, vu qu'il n'y en a pas dans les noyaux de protium.

Rappelons que le rapport protons/neutrons est de 1/7. Cela veut dire que sur 16 baryons, 2 sont des neutrons et 14 sont des protons (ce qui est équivalent à dire que sur 8 baryons, 1 est un neutron et 7 sont des protons). Avec ces 16 baryons, on peut créer un atome d'hélium-4 avec 2 neutrons et 2 protons, ce qui laisse 12 protons. On a donc 12 noyaux d’hydrogènes pour un noyau d'hélium-4. Vous avez peut-être vu d'autres chiffres dans la vulgarisation, notamment un rapport de 75% d'hydrogène contre 25% d'hélium. Mais ces pourcentages sont exprimés en termes de masse, non de nombre d'atomes. Pour retrouver ce rapport 75%/25%, il faut prendre en compte la masse relative des noyaux d'hélium-4 et d'hydrogène. Un noyau d'hélium-4 étant approximativement 4 fois plus massif qu'un noyau d'hydrogène, il faut diviser le nombre de noyaux d'hydrogène par 4 pour utiliser la même unité de masse. En faisant cela, on trouve alors que 75% de la masse de l'univers est composée d'hydrogène et les 25% restants d'hélium-4.

Abondance des éléments chimiques suite à la nucléosynthèse primordiale.

Il est possible de retrouver ces résultats par le calcul, ce que nous allons faire de ce pas. Dans les calculs qui vont suivre, nous noterons le nombre d'atomes d'hélium alors que le nombre d'atomes de protium sera noté . Les nombres de neutrons et de protons seront notés et . L'hypothèse comme quoi tous les neutrons sont capturés dans les noyaux d'hélium-4 signifie que :

, car il y a deux neutrons dans un noyau d'hélium-4.

La première étape est de calculer la quantité totale de baryons enfermés dans les noyaux d'hélium-4. Un atome d'hélium-4 contient 4 baryons, deux neutrons et deux protons. En multipliant par le nombre de noyaux, on obtient cette quantité totale, égale à :

Pour la seconde étape, on a besoin du nombre total de baryons dans l'univers. Par définition, il est égal à la somme (la somme du nombre de protons et de neutrons). En combinant les deux équations précédentes, on obtient le rapport entre le nombre de baryons dans les atomes d'hélium et le nombre total de noyaux, que nous noterons .

On utilise alors l'équation  :

On peut réécrire cette équation en utilisant uniquement le rapport protons/neutrons calculé dans la section précédente. Il suffit pour cela de diviser l'équation précédente par le nombre de protons. On a alors :

On a vu dans la section précédente que le le rapport protons/neutrons suite au big-bang est de . En utilisant cette valeur, on trouve que , ce qui veut dire que 75% de la masse de l'univers est composé d'hydrogène et 25% d'hélium-4. Cette valeur est très proche de la valeur observée. À l'heure actuelle, 74% de la masse de l'univers est composée d'hydrogène et 25% d'hélium-4, le reste se partageant le 1 % restant. Précisons qu'il s'agit d'un pourcentage en masse, non en nombre d'atomes.



Le découplage des photons et neutrinos

Dans les chapitres précédents, nous avons vu quelque chose de très important : la température du rayonnement diminue avec l'expansion. Un point important est que l'univers, à ses débuts, est composé d'un plasma, à savoir un gaz de particules chargées. L'univers ne comprend au départ par d'atomes, mais seulement des électrons, des protons et neutrons, etc. La matière était composée d'un plasma d'électrons libres, de baryons et de photons. A cette époque, le rayonnement interagissait fortement avec les électrons libres, les noyaux d'atomes, les électrons liés aux atomes, bref : avec la matière au sens général. Vu que le rayonnement et la matière interagissaient beaucoup, les deux étaient à la même température. Matière et rayonnement étaient en équilibre thermique, à la même température. Si la matière chauffait trop, elle redistribuait l'excès au rayonnement, et réciproquement. Ce plasma formait un fluide unique, avec une pression, une température, une densité d'énergie, etc.

Les photons interagissaient fortement avec les électrons, par diverses processus (diffusion Compton, et autres), ce qui redistribuait la température. Les photons chauffaient les électrons et réciproquement. Un équilibre thermique s'était ainsi installé entre photons et électrons libres, les deux ayant la même température/énergie cinétique moyenne. Cet équilibre incluait aussi les baryons, bien que les baryons n'interagissaient pas directement avec les photons. Les baryons interagissaient fortement avec les électrons, qui servaient d'intermédiaires avec les photons. Ce plasma avait naturellement des propriétés thermodynamiques simples : une pression, une température, un volume, etc. Sa pression était essentiellement causée par la pression de radiation des photons, avec une participation mineure de la pression des électrons libres et des baryons.

La matière neutre est transparente, ce qui induit un découplage des photons.

Mais la température baissant avec l'expansion, elle a finit par atteindre une température suffisamment basse pour que les atomes se forment. La matière est alors passée de l'état de plasma à un gaz d'atomes. Cela s'est produit 380 000 ans après le big bang, à une température d'environ 3000 degrés. Le résultat est que les interactions entre rayonnement et matière se sont faites plus rares, et surtout que la redistribution d'énergie entre matière et rayonnement s'est faite plus rare, plus lente.On nomme découplage de telles situations où deux populations de particules n'interagissent plus à la suite d'une baisse de température. Et ici, il s'agit du découplage des photons.

Après le découplage, le rayonnement et la matière font leur vie chacun de leur côté et cela s'est ressenti sur leur température et leur évolution. La première conséquence du découplage est que l'équilibre thermique est rompu : les deux populations de particules s'isolent thermiquement et leur températures deviennent indépendantes. La température de la matière est devenue différente de la température des photons. De plus, elles ne se sont pas refroidies de la même manière. La matière garda sa température constante lors de l'expansion, alors que celle du rayonnement chuta avec l'expansion. Le rayonnement s'est refroidi plus vite que la matière, à cause de la diminution de fréquence des photons, abordées il y a quelques chapitres.

Le fond diffus cosmologique modifier

Pour résumer, le plasma s'est divisé en deux gaz indépendants : un gaz de matière et un gaz de photon. Le gaz de matière s'est condensé pour donner des galaxies et autres structures, alors que le gaz de photons a subsisté jusqu’à aujourd'hui sous la forme d'un ensemble de photons de faible température, que l'on peut capter avec certains instruments. Ce gaz est appelé le fond diffus cosmologique, aussi appelé le rayonnement de fond diffus cosmologique, ou encore le CMB (Cosmic Microwave Background).

La même chose a eu lieu pour les neutrinos et anti-neutrinos qui se sont découplés de la matière et des photons un peu avant les photons. Ce fond diffus de neutrinos est malheureusement nettement moins étudié que le fond diffus cosmologique, car les neutrinos n'interagissent pas beaucoup avec la matière, et qu'ils sont donc difficiles à détecter. Nous n'en parlerons donc pas dans ce cours, par manque d'informations à leur sujet. Pour le moment, concentrons-nous sur le découplage des photons.
COBE monopole dipole and primordial perturbations

Le CMB a été théorisé avant d'être découvert. Dans un article de 1948, Alpher et ses collègues théorisèrent l'existence du CMB à partir d'un modèle de big-bang usuel. Mais il fallut attendre 1965 pour que ce signal soit observé pour la première fois, par Penzias et Wilson. Ceux-ci utilisaient une antenne de grandes dimensions, pour tester la fiabilité des communications entre satellites, et étudiaient des interférences radio qui apparaissaient à haute fréquence. Leurs investigations leur ont permis de capter un signal dans la bande de 4Ghz, qui avait des caractéristiques étranges : isotrope, non-polarisé et libre de toute variation saisonnière. L'origine de ce signal est restée inconnue durant quelques années, mais les scientifiques (dont Penzias et Wilson) avaient éliminé toute origine terrestre. Il fallu que Dicke et ses collaborateurs fassent le lien avec l'article d'Alpher. Par la suite, diverses campagnes d'observation ont permis d'obtenir une carte assez détaillée du fond diffus. De nombreux projets d'observations scientifiques ont ainsi observé le fond diffus cosmologique avec une précision de plus en plus grande : COBE, puis WMAP, et enfin la mission PLANCK.

Les observations de Penzias et Wilson montraient un CMB relativement uniforme. Par la suite, les observations du satellite COBE ont montré que le CMB a l'air d'avoir une structure en forme de dipôle, à savoir qu'il a un pôle chaud opposé à un pôle froid. Les observations plus récentes éliminent cet effet Doppler par divers traitements informatiques, et montrent un CMB sans dipôle, mais avec quelques inhomogénéités. On observe notamment une zone plus chaude au niveau de l'équateur, liée à la présence de la voie lactée (notre galaxie), qui réchauffe quelque peu le CMB de par son rayonnement.

La surface de dernière diffusion modifier

Surface de dernière diffusion.

Le fond diffus capté à l'heure actuelle correspond aux photons émis lors du découplage, environ 380 000 ans après le big-bang. Du point de vue de la Terre, les photons du CMB captés par les instruments de mesure sont répartis sur une sphère centrée sur la Terre, appelée la surface de dernière diffusion, dont le diamètre correspond à la distance parcourue par ces photons depuis le découplage.

Pour comprendre pourquoi on parle de surface de dernière diffusion, il faut faire un petit rappel sur les interactions entre photons et matière. Les photons ont en effet tendance se cogner sur les particules qui croisent leur chemin, et à rebondir dessus. Typiquement, ils rebondissent sur des électrons dans un processus appelé la diffusion Compton. Lors de ce choc entre électron et photon,une partie de l'énergie du photon est transférée à l'électron, ou inversement.

Diffusion Compton

Les photons du CMB sont soit des photons qui ont été émis au moment du découplage, soit des photons qui existaient avant, mais ont rebondit une ou plusieurs fois sur des particules de matière. Les diffusions Compton étaient très fréquentes avant le découplage, car les électrons libres étaient très nombreux. Les photons ne pouvaient pas faire un pas sans interagir avec un électron libre, leur temps de libre parcours moyen (temps entre deux collisions avec un électron libre) était très faible. Le mouvement d'un photon dans le plasma primordial était chaotique : les photons ne pouvaient pas aller en ligne droite, mais n'avaient de cesse de se cogner sur les atomes au point où leur trajectoire ressemblait à marche aléatoire, à un mouvement brownien. En conséquence, le plasma primordial avant le CMB était opaque, dans le sens où il ne laissait pas traverser la lumière facilement.

Mais après le découplage, les électrons se sont liés aux atomes, ne laissant presque aucun électron libre pour participer à des diffusions Compton. Les photons du CMB, autrefois diffusés fréquemment, n'ont ensuite plus été diffusé par diffusion Compton, ou presque. Leur libre parcours moyen, la distance parcoure entre deux collisions/diffusions est devenus beaucoup plus grande. Les photons ont alors pu se déplacer en ligne droite entre quelques rares collisions, et donc le plasma primordial est devenu transparent. Le découplage est le moment où ils ont diffusé pour la dernière fois, d'où le fait que l'on parle de surface de dernière diffusion.

Précisons cependant que le découplage ne s'est pas fait instantanément, mais est un processus progressif qui s'est étalé sur un temps non-négligeable. Et cela a un effet sur la surface de dernière diffusion. Il s'agirait bien d'une surface sphérique si le découplage était instantané, mais la réalité est que la surface de dernière diffusion a une certaine épaisseur. Son épaisseur est égale à la distance parcourue par un photon durant la durée du découplage. Cela a des effets mesurables, notamment dans le spectre de puissance du fond diffus, mais nous en reparlerons dans quelques chapitres.

Le CMB est un rayonnement de corps noir modifier

Les observations montrent que, peu importe à quel endroit on le regarde, le CMB est un rayonnement de corps noir quasiment parfait, ce qui est en accord avec la théorie. Le fait que le CMB soit un rayonnement de corps noir signifie que l'on peut lui attribuer une température, une température en chaque point du ciel pour être précis.

Comparaison du CMB avec un rayonnement de corps noir.

Le dipôle du CMB modifier

CMB vu par COBE, sans traitement (en haut), en retirant le dipôle (milieu) et en retirant la luminosité de la voie lactée (en bas).

Le fait que le CMB aie une forme approximativement dipolaire est en soi un gros problème, que l'on arrive pas à bien expliquer théoriquement. On pourrait croire que cela réfute l'idée d'un univers isotrope, mais ce n'est pas forcément le cas.

L'explication usuelle quant à l'existence de ce dipôle est qu'il serait d'origine cinématique, lié au mouvement de la Terre par rapport au CMB. L'explication détaillée part du principe qu'il y aurait un référentiel dans lequel le CMB serait isotrope et sans aucun dipôle. Cependant, la Terre se déplacerait à une certaine vitesse notée v par rapport à ce référentiel. Rappelons que le CMB est avant tout un rayonnement électromagnétique avec une certaine fréquence, fréquence qui dépend de sa température. Or, la fréquence du CMB est alors déformée par une sorte d'effet Doppler, ce qui fait que la température dépend de la direction dans laquelle on se déplace, par rapport au référentiel où le CMB est isotrope. Les zones du CMB qui s'éloignent de nous sont vues comme refroidies, alors que les zones qui s'approchent (opposées, donc) sont vues comme plus chaudes.

Pour rappel, à de faibles vitesse, le CMB subit un décalage vers le rouge de :

La relativité restreinte donne une formule encore plus précise :

Une explication beaucoup plus détaillée est cependant encore inconnue. On fait face à un cas assez mal compris théoriquement : établir la température d'un objet alors qu'on se déplace par rapport à lui à grande vitesse. Et outre la température, on peut aussi analyser la forme du rayonnement de corps noir observée. Température et corps noirs sont deux concepts de la thermodynamique, une branche de physique parfaitement maitrisée tant qu'on reste dans le cadre de la physique classique de Newton. Mais le cas qui nous intéresse, à savoir à quoi ressemble le fond diffus quand on se déplace à grande vitesse, demande qu'on utilise la relativité restreinte d'Einstein. Et c'est là que le bas blesse : on n'a pas de théorie validée pour la thermodynamique relativiste. La thermodynamique relativiste vise à reformuler les lois de thermodynamique dans le cadre de la relativité restreinte d'Einstein. A l'heure actuelle, les physiciens ont plusieurs théories pour cela, mais elles donnent des résultats complètement contradictoires !

Et pour comprendre pourquoi, faisons quelques rappels de relativité restreinte. La relativité nous dit que quand on se déplace par rapport à un référentiel immobile, les distances et les durées sont modifiées par un facteur appelé le facteur de Lorentz. Ce dernier vaut :

, avec est la vitesse de déplacement, et est la vitesse de la lumière dans le vide.

Concrètement, si je mesure une distance L dans un référentiel immobile, la même distance mesurée en allant à la vitesse v sera de . La distance mesurée est plus petite que dans un référentiel immobile et elle est d'autant plus petite que l'on va vite. C'est l'effet de contraction des longueurs. L'effet est inverse pour les durées, à savoir que la durée d'un évènement/phénoméne n'est pas la même selon la vitesse du référentiel de mesure. La durée mesurée dans un référentiel en mouvement est plus longue que celle mesurée dans un référentiel immobile. Cette fois-ci, une durée D devient . C'est l'effet de dilatation des durées.

Maintenant, prenons le cas où on cherche à calculer la température observée du CMB, ou même de tout autre rayonnement de corps noir observé depuis un objet en mouvement. Le CMB, dans son référentiel où il est isotrope (immobile dirait-on), a une température notée . L'observateur se déplace à une vitesse dans une direction bien précise. En théorie, plusieurs situations sont possibles.

  • La température observée est égale à .
  • La température observée est égale à .
  • La température observée est égale à .

A l'heure actuelle, nous ne savons pas quelle est la bonne solution. Les trois possibilités sont impossibles à distinguer expérimentalement, et les trois sont défendues par des arguments théoriques solides. Une revue de ces arguments est disponible dans le document nommé "Black Body Radiation in Moving Frames", par Kamran Derakhshani. Mais l'explication actuelle du dipole du CMB part du principe que l'on a une variante de la seconde solution. En 1968, Bracewell et Conklin ont sorti une étude où ils montraient que, si leur démonstration était correcte, alors un observateur en mouvement dans le CMB mesurait une température anisotrope, à savoir dépendante de la direction. La température mesurée en un point du CMB était la suivante :

, avec l'angle par rapport à la direction de déplacement.

L'équation nous dit que le CMB doit semble plus chaud à l'avant, et plus froid, par rapport à la direction de déplacement. On retrouve donc le dipole observé dans le CMB, à condition que la vitesse de déplacement ait la bonne valeur. Elle serait d'environ 300 kilomètres par secondes. C'est l'explication actuelle du dipole du CMB.

L'explication précédente du dipôle du CMB est actuellement contestée. A vrai dire, il s'agit d'une hypothèse parmi beaucoup d'autres et elle est juste un consensus assez mou, que les astrophysiciens cherchent à vérifier ou infirmer depuis longtemps. Diverses expériences ont tenté de vérifier l'origine cinétique du dipôle, de plusieurs façons. La première est de vérifier si ce dipôle apparait non seulement dans le CMB, mais aussi dans les observations des quasars et des sources radio lointaines. En théorie, si l'explication cinétique est bonne, alors on devrait retrouver un effet similaire en observant la positions des galaxies très lointaines. Diverses études ont tenté de vérifier cela, mais les résultats ne sont pas encore très fiables. De plus, l'explication cinétique implique qu'il y ait un référentiel dans lequel le CMB serait quasi-uniforme et sans dipôle. Cela pourrait en faire un référentiel privilégié et briserait le principe de relativité avec l'invariance galiléenne/lorentzienne, mais rien n'est moins sûr.

Outre l'existence de ce dipôle, d'autres observations semblent indiquer que le CMB aurait un axe privilégié, nommé l'axe du mal (axis of evil en anglais), qui impacterait le CMB, mais aussi l'expansion de l'univers. Mais le fait que l'axe soit aligné avec le plan de écliptique (le plan sur lequel la Terre tourne autour du Soleil) met le doute quant à son origine. Une origine cinématique locale est en tout cas plus probable qu'une origine cosmologique, sauf coïncidences. Mais là encore, rien n'est moins sûr et il se pourrait que les données soient biaisées, d'autres études n'ayant pas retrouvé cet axe du mal dans leurs données.

L'axe du mal du CMB.

Le découplage et le fond diffus modifier

Au moment du découplage, on sait que le gaz de photons devait avoir la même température que le plasma. Sans expansion, cette température serait égale à la température du plasma au moment du découplage, qui a été conservée par le gaz de photons. Mais l'expansion a décalé ce rayonnement de corps noir vers le rouge, diminuant sa température. La température du fond diffus a donc diminué en conséquence. De nos jours, les mesures donnent une température d'environ 2,735 Kelvin. Mais on est en droit de se demander quelle était sa température au moment de sa formation. On peut aussi se demander combien de temps de refroidissement du CMB a duré.

L'âge de la recombinaison modifier

La température du fond diffus au moment du découplage est estimée à 3000 degrés Kelvin, température de condensation d'un plasma en atomes. À partir de la température mesurée actuellement, et de la valeur théorique de formation d'un plasma, on peut déduire l'âge qui s'est écoulé depuis la formation du CMB. En théorie, on peut déduire la température du CMB en utilisant la formule suivante, établie dans le chapitre "L'évolution du rayonnement" :

Mais utiliser cette formule présuppose de connaître et . Cependant, on peut ruser en remplaçant le facteur d'échelle par le redshift. Pour cela, on utilise la formule vue dans le chapitre "L'évolution du rayonnement". Le décalage vers le rouge est mesuré entre l'époque actuelle et la recombinaison, ce qui fait qu'on le notera .

Rappelons que cette formule utilise la convention . Ici, le temps d'émission est l'instant où a eu lieu la recombinaison. En clair, la température est la température de découplage .

La température au moment du découplage était d'environ 3000 degrés Kelvin, alors que la température actuelle du fond diffus est d'environ 2,735 degrés Kelvin.

Le redshift calculé ainsi est de :

.

En utilisant un modèle cosmologique, on peut déduire une relation âge-redshift, et donc calculer combien d'années se sont écoulées entre le big-bang et la formation du CMB.

La température de recombinaison modifier

La recombinaison a eu lieu quand l'univers a atteint une certaine température, qu'il est important de connaître. En effet, grâce à elle, on peut calculer quand a eu lieu le découplage, et donc dater le CMB. Autant dire que calculer celle-ci est d'une importance primordiale. En théorie, la température du CMB est la température à laquelle un plasma se condense en atomes quand on le refroidit. Dit autrement, c'est la température d'une transition de phase. Vous avez peut-être déjà entendu que cette température est d'environ 3000 degrés Kelvin, ce qui est la température mesurée sur Terre. Reste qu'il vaut mieux la calculer et en rendre compte théoriquement. Les mesures réalisées sur Terre ne sont peut-être pas représentatives des conditions de l'univers primordial : la pression est plus élevée, la densité différente et j'en passe. Dans cette section, nous allons calculer la température théorique à laquelle le découplage a eu lieu.

L'approximation par l'énergie photonique moyenne modifier

Une première méthode est de comparer l'énergie d'ionisation de l'hydrogène avec l'énergie d'un photon. Dans un gaz de photons de température , l'énergie moyenne d'un photon est de . L'énergie d'ionisation d'un atome d'hydrogène (la plus faible de toutes) est de 13,6 électrons-Volts (l'EV est une unité d'énergie). Si l'énergie moyenne des photons est supérieure à l'énergie d'ionisation de l'hydrogène, alors la matière restera ionisée. Si ce n'est pas le cas, les photons ne sont pas assez énergétiques pour ioniser la matière, qui se condense. On peut alors calculer une approximation de la température de découplage avec le quotient suivant :

On voit que la température obtenue est diablement haute, comparée aux valeurs réelles : plus de 10 fois la valeur réelle. Cela vient d'un phénomène simple : l'énergie moyenne n'est qu'une moyenne, qui cache le fait que certains photons sont plus énergétiques que la moyenne. Même si l'énergie moyenne d'un photon est de 13,6 eV, de nombreux photons ont une énergie suffisante pour ioniser un atome dans le gaz de photon.

On peut obtenir une approximation plus crédible en étudiant plus en détail la dispersion des énergies des photons. Dans un gaz de photons, tous les photons n'ont pas la même énergie. Et l'énergie n'est pas répartie équitablement entre les photons, mais suit une loi assez compliquée. Mais pour simplifier, on peut estimer que le nombre de photons suit la loi de Boltzmann. Pour résumer, le nombre de photons qui a au moins une énergie est proportionnel à l'exponentielle de leur énergie. Mis sous forme de formule, cela donne :

.

Supposons maintenant qu'il faille un photon par atome à ioniser pour que l'ionisation se fasse. Sachant qu'il y a environ photons par baryon.atome, on peut trouver la température de découplage suivante :

L'approximation par l'équation de Saha modifier

Il est possible d'obtenir une approximation plus précise avec l'équation de Saha. Celle-ci permet de déduire le degré d'ionisation d'un gaz. Le gaz en question correspond à un gaz d'hydrogène, composant principal de l'univers, qui s'est justement formé lors du découplage. Avant le découplage, on peut considérer que l'univers était rempli d'un plasma formé par ionisation du gaz d'hydrogène, à savoir un gaz qui mélangeait protons et électrons. L'équation de Saha nous dit que, si on note :

  • , et la concentration en électrons, protons et hydrogène ;
  • la masse de l'électron ;
  • l'énergie d'ionisation d'un atome d'hydrogène.
Vous remarquerez que l'équation de Saha ressemble beaucoup à la distribution de Boltzmann.

Pour diverses raisons techniques, les physiciens décrivent souvent l'ionisation d'un gaz en utilisant la fraction d'électrons libres. Elle correspond au rapport en nombre d'électrons libres et nombre de protons (libres ou appartenant à un atome d'hydrogène). Elle vaut donc :

Le plasma primordial est neutre électriquement. Or, la neutralité électrique de la matière signifie que . Avec cette contrainte, l'équation de Saha se réécrit comme suit :

L'équation nous dit que le découplage n'est pas un évènement qui a eu lieu à une température bien précise, mais un processus continu dans lequel l'ionisation a baissé lentement. À 5000 Kelvins, la matière est ionisée à près de 99%. Le taux d'ionisation chute ensuite progressivement avec la température, pour atteindre 50% à 4000 Kelvins, puis 1% à 3000 Kelvins. Les scientifiques estiment, par pure convention, que le découplage a eu lieu quand le degré d'ionisation descend en-dessous de 1%, c'est à dire à une température d'environ 3700 Kelvin.

Les approximations plus précises modifier

Les observations ne sont pas complètement compatibles avec cette approximation, bien que le modèle de Saha colle à-peu-près. Les mesures semblent indiquer non seulement que la température calculée n'est pas tout à fait exacte, mais qu'en plus, le découplage a pris plus de temps, s'est déroulé plus lentement. Pour obtenir des résultats plus précis, divers modèles ont été inventés par les physiciens. Le plus simple de ces modèles est le modèle de Peebles, aussi connu sous le nom de modèle d'atome à trois étages. Il a été complété par de nombreux modèles, qui sont devenus de plus en plus complexes avec le temps, incluant de plus en plus d'acquis théoriques provenant de la physique atomique. Nous ne parlerons pas de ces modèles, qui sont assez compliqués pour ce cours et qui font notamment appel à quelques concepts de physique quantique.



Une introduction aux modèles cosmologiques

Les chapitres précédents nous ont appris beaucoup de choses. Ils ont fourni une description de l'expansion de l'univers basée sur le facteur d'échelle, de quoi calculer le rayon de l'univers observable, et ont décrit comment la matière et le rayonnement réagissent à l'expansion. Maintenant, nous allons voir une portion bien plus intéressante de la cosmologie. Nous allons étudier les modèles cosmologiques, des modèles théoriques qui décrivent comment l'univers lui-même a évolué au cours du temps. Les modèles cosmologiques que nous allons voir dans ce cours donnent au minimum les résultats suivants : l'évolution du facteur d'échelle au cours du temps (la fonction , le rayon observable de l'univers, son âge (s'il en a un), une formule pour le facteur de Hubble, et parfois d'autres données supplémentaires (T(t) : température à tout instant, ...).

Dans ce chapitre, nous allons voir quelques modèles cosmologiques ad hoc. Par ad hoc, on veut dire qu'on fait des hypothèses purement arbitraires pour démontrer le modèle. De tels modèles sont à contraster avec des modèles démontrés à partir de principes physiques fondamentaux. Les modèles démontrés à partir de lois physiques plus fondamentales sont en effet un peu plus compliqués et seront vus dans les chapitres qui suivent. Le prochain chapitre porte d'ailleurs sur les modèles de Friedmann, qui sont les modèles les plus utilisés en cosmologie physique. Dans ce chapitre, nous allons voir des modèles plus simples. La démarche de ce chapitre consistera simplement à postuler une loi d'évolution pour le facteur d'échelle, à postuler la fonction . Une fois cela fait, nous pourrons établir une formule pour le facteur de Hubble, le rayon cosmologique, etc. Une telle démarche demande de choisir une fonction en fonction de ses propriétés mathématiques, de sa facilité pour faire les calculs, etc. Dans les modèles cosmologiques physiques, cette loi est démontrée à partir de principes physiques plus fondamentaux.

Dans ce chapitre, nous allons voir plusieurs modèles. Le premier est celui pour lequel la loi est linéaire. Il s'agit d'un modèle très simple, mais qui n'est pas réaliste du tout. Cependant, ce modèle permet d'introduire deux grandeurs fondamentales : le temps et le rayon de Hubble. De plus, il sert de point de comparaison avec le second type de modèles : ceux pour lesquels la loi est une loi de puissance. Enfin, nous verrons ce qui se passe quand le facteur d'échelle croit de façon exponentielle. En soit, ces trois modèles font que l'univers est en expansion, mais de manière différente. Une croissance linéaire n'a rien à voir avec une croissance exponentielle, ni avec une croissance en loi de puissance. Nous verrons cependant que les modèles cosmologiques physiques, ceux de Friedmann, donnent une croissance soit exponentielle, soit en loi de puissance. Nous pourrons donc réutiliser les résultats de ce chapitre dans les chapitres ultérieurs sur les modèles physiques. Les modèles ad hoc que nous allons voir vont nous donner des résultats qui faciliteront l'étude de modèles plus physiques, plus compliqués.

Le modèle à croissance linéaire modifier

Pour commencer, nous allons étudier un premier modèle : le modèle à croissance linéaire. Il est aussi appelé modèle , en raison d'une des équations du modèle[1]. Les cosmologistes ont étudié ce modèle à partir de l'année 2012, car de nombreux problèmes des autres modèles disparaissent dans ce modèle, mais il semble qu'il ne colle pas trop aux observations actuelles[2]. Quoi qu’il en soit, ce modèle est intéressant à étudier et il sert de bonne introduction, avant de passer à des modèles plus compliqués.

Avec ce modèle, on suppose que l'univers gonfle de manière relativement constante avec le temps. Mathématiquement, cela se traduit par la loi d'évolution suivante pour le facteur d'échelle :

, avec le temps, une constante quelconque.

Une autre formulation, souvent utilisée, est la suivante :

, avec l'origine de temps choisie arbitrairement et le facteur d'échelle à ce même instant.

Notons que les deux sont équivalentes, sous réserve que . Dans ce qui va suivre, nous utiliserons la première formulation, sauf pour ce qui est l'étude du décalage vers le rouge. Avec ces formules, nous allons calculer quelle est la valeur du facteur de Hubble qui correspond, quel est l'âge de l'univers, quel est son rayon, la vitesse de l'expansion de l'univers, etc.

Avant de poursuivre, on peut calculer la dérivée première et seconde du facteur d'échelle, deux résultats qui seront utiles pour la suite. Le calcul de la dérivée première donne :

Celle-ci étant une constante, la dérivée seconde est donc nulle.

Rappelons que la dérivée seconde du facteur d'échelle nous dit si l'expansion de l'univers accélère, ralentit ou reste constante. Pour une valeur nulle, l'expansion se fait à rythme constant.

Le calcul du facteur de Hubble et de l'âge de l'univers modifier

Pour calculer le facteur de Hubble, nous allons partir de la dérivée du facteur d'échelle, qui est égale à :

En divisant par a(t), le terme de gauche devient le facteur de Hubble. On a alors :

On peut comparer l'équation précédente avec l'identité vue dans le second chapitre. On obtient alors :

L'âge de l'univers est égal au temps de Hubble dans ce modèle.

L'évolution du facteur de Hubble dans le temps modifier

De l'équation , on peut déduire comment le facteur de Hubble H évolue au cours du temps. Pour cela, il suffit d'étudier ce qui se passe quand l'âge de l'univers augmente. Au fur et à mesure que le temps passe, l'âge de l'univers augmente. Ce faisant, le dénominateur de croit et le facteur de Hubble diminue donc. Pour le dire en termes mathématiques, le facteur de Hubble est une fonction décroissante du temps. Un bon moyen de s'en rendre compte est de calculer la dérivée du facteur de Hubble :

La dérivée étant négative, on en déduit que le facteur de Hubble décroit au cours du temps. Et cette décroissance est d'autant plus rapide que l'univers est âgé.

Le calcul du paramètre de décélération modifier

À partir des résultats précédents, on peut calculer la valeur du paramètre de décélération. Pour rappel, celui-ci est un nombre qui dit à quelle vitesse l'expansion de l'univers accélère ou ralentit. L'expansion est décélérée s'il est positif, accélère s'il est négatif et reste constante sinon. Il se calcule en utilisant la formule suivante, vue dans le chapitre sur l'univers observable :

On injecte alors les équations et , ce qui donne :

En clair, l'expansion garde un rythme constant, elle n’accélère pas et ne ralentit pas, comme dit précédemment.

Le rayon de Hubble modifier

Par définition, le rayon de Hubble est égal à :

En combinant avec l'équation , on a :

Le modèle à expansion linéaire est aussi appelé modèle , en raison de ce résultat. D'ailleurs, c'est le seul modèle dans lequel cette égalité est tout le temps vraie. Si l'on a pas une expansion linéaire, alors le rayon de Hubble n'est égal à que pour une valeur de bien précise et fausse pour toutes les autres.

Le rayon comobile de Hubble est définit par la formule :

Par définition, on a , qui est une constante.

On remarque que toute la partie de droite ne contient que des termes constants (c et k sont tous deux constants). Ce qui signifie que le rayon comobile de l'univers est lui aussi constant. Certes, le rayon de Hubble grandit, mais il le fait au rythme que le facteur d'échelle. Les deux se compensent, donnant un rayon comobile constant.

Le calcul du rayon de l'univers observable modifier

Dans ce modèle, calculer le rayon de l'univers observable n'est pas une mince affaire. On a vu dans le chapitre sur l'univers observable qu'il faudrait idéalement utiliser cette formule :

Malheureusement, si on se met à faire les calculs à la main, on s’aperçoit rapidement que l'intégrale ne converge pas ! Le modèle nous dit donc que le rayon observable de l'univers est infini. Cela est un avantage, car cela résout un problème majeur, appelé problème de l'horizon, qui touche les autres modèles cosmologiques. Nous en reparlerons vers la fin du cours plus en détail.

Le calcul du décalage vers le rouge modifier

Dans ce qui suit, nous allons supposer qu'un objet à une distance émet de la lumière à un instant , lumière qui est reçue par l'observateur à un instant . Nous allons calculer le décalage vers le rouge de cette lumière, le redshift. Pour gérer ce cas, nous allons devoir utiliser la formule suivante, vue dans le chapitre sur le rayonnement, qui donne le redshift en fonction du facteur d'échelle :

On injecte alors l'équation  :

En simplifiant, on trouve :

Il est aussi possible d'écrire le facteur de Hubble en fonction du décalage vers le rouge et réciproquement. Par exemple, imaginons qu'un objet émette de la lumière à un instant , quand le facteur de Hubble est égal à et le facteur d'échelle à . Sa lumière est reçue par l'observateur à l'instant , quand le paramètre de Hubble et le facteur d'échelle valent et . On peut alors démontrer l'équation suivante :

Qui peut se reformuler comme suit :


Démonstration

Pour démontrer la formule précédente, partons de la relation :

Or, on sait que . On a donc et

En faisant le remplacement dans l'équation précédente, on trouve :

En réarrangeant les termes, on trouve la formule à démontrer, mais sous une forme légèrement différente :

Un modèle trop simplifié pour être réaliste modifier

Les scientifiques ont estimé le rayon de l'univers observable, et leurs résultats semblent relativement proches du rayon de Hubble. Là encore, comme pour la coïncidence entre âge de l'univers estimé et temps de Hubble, personne ne sait si c'est une coïncidence ou quelque chose de plus important. En théorie, l'égalité est respectée en permanence dans le modèle linéaire. Par contre, elle n'est respectée que pour un instant bien précis dans les autres et violée dans tous les autres cas. Le fait que l'égalité soit respectée dans notre monde est donc soit une grosse coïncidence, soit elle trahit quelque chose de fondamental qui est encore mal compris. Par contre, les résultats empiriques concernant les mesures de redshift, ou via d'autres mesures indirectes, ne collent pas du tout avec ce modèle. Autant l'égalité est respectée, autant les autres prédictions du modèle ne le sont pas.

Et les mesures ne sont pas seules à nous dire que ce modèle est trop simpliste pour décrire fidèlement l'univers. La théorie nous dit que ce modèle n'est valide que pour un univers vide, c'est à dire en absence de matière. Toute présence de matière biaise les résultats du modèle[3]. Et aux dernières nouvelles, la matière existe...

Le modèle à croissance en loi de puissance modifier

Maintenant, nous allons étudier le cas d'une expansion en loi de puissance. En clair, nous allons étudier le cas où le facteur d'échelle est égal à :

, avec le temps, une constante quelconque et une puissance quelconque.

Une autre formulation, souvent utilisée, est la suivante :

, avec l'origine de temps choisie arbitrairement et le facteur d'échelle à ce même instant.

Les deux formulations sont équivalentes si l'on postule . Dans ce qui suit, nous utiliserons la première formulation.

On peut d'or et déjà remarquer que l'expansion linéaire n'est qu'un cas particulier d'expansion en loi de puissance, où . Ne vous étonnez donc pas si les calculs se ressemblent quelque peu.

Le calcul du facteur de Hubble et des paramètres associés modifier

Pour calculer le facteur de Hubble, nous allons partir de la dérivée du facteur d'échelle, qui est égale à :

On divise alors par a(t) pour obtenir le facteur de Hubble :

Pour résumer :

L'évolution du facteur de Hubble modifier

Pour savoir comment évolue le facteur de Hubble dans le temps, le mieux est encore une fois de calculer sa dérivée. La dérivée du facteur de Hubble est égale à :

On pourrait croire qu'elle est négative, en raison du signe -, mais tout dépend de la valeur de n. C’est le cas si n est positif ou nul, mais il peut être négatif.

Le calcul du paramètre de décélération modifier

À partir des résultats précédents, on peut calculer la valeur du paramètre de décélération. Pour le calculer, on applique la formule suivante, vue dans le chapitre sur l'univers observable :

On injecte l'équation dans l'équation précédente, ce qui donne :

On injecte alors l'équation , ce qui donne :

Ce qui se simplifie en :

On peut réécrire le tout comme suit :

En clair, l'expansion dépend du coefficient .

  • Si , alors on retombe sur le cas de l'expansion linéaire. L'expansion a un rythme constant : elle n’accélère pas et ne ralentit pas.
  • Si , l'expansion accélère.
  • Enfin, si , l'expansion de l'univers ralentit.

Le calcul de l'âge de l'univers modifier

L'équation : permet de calculer l'âge de l'univers à partir du facteur de Hubble. En isolant dans l'équation précédente, on trouve en effet :

Reformulons en utilisant le temps de Hubble :

On voit que l'âge de l'univers est un multiple du temps de Hubble, le coefficient de proportionnalité n'étant autre que l'exposant de la loi de puissance.

Le calcul du rayon de Hubble modifier

Il est possible, à partir des résultats précédents, de calculer le rayon de Hubble assez facilement. Pour rappel, celui-ci est la distance à partir de laquelle les objets semblent aller plus vite que la lumière à cause de l'expansion. Pour rappel, il vaut :

On injecte alors l'équation :  :

Dans ce modèle, le rayon de Hubble augmente avec le temps.

Le calcul du rayon de l'univers observable modifier

À partir des résultats précédents, nous pouvons trouver le rayon de l'univers observable, qui n'est pas le même que le rayon de Hubble. Pour cela, il faut cependant se rappeler ce qu'on a vu dans le chapitre sur l'univers observable. Pour calculer le rayon de l'univers observable, il faut procéder en deux étapes. Premièrement, il faut calculer le rayon comobile, le rayon calculé en supprimant l'effet de l'expansion. Ensuite, on corrige ce rayon comobile pour ajouter l'expansion de l'univers. Si on note le rayon de l'univers à l'instant t, alors le rayon comobile est égal à . Le passage du rayon comobile au rayon observé se fait donc simplement en multipliant par le facteur d'échelle . La deuxième étape est donc triviale. Par contre, la première étape, le calcul du rayon comobile, fait intervenir une intégrale assez compliquée. Pour rappel, le rayon comobile se calcule avec la formule suivante :

Pour calculer le rayon comobile, nous allons introduire l'équation dans la précédente, ce qui donne :

La constante k peut être sortie de l'intégrale, ce qui donne :

On utilise alors la formule  :

On utilise alors la formule  :

On regroupe les termes comme suit :

Par définition, le terme est égal à , ce qui donne :

On peut donc en déduire le rayon propre de l'univers observable en multipliant des deux côtés par  :

Si on prend , on retombe sur les résultats du modèle vus précédemment. Et encore une fois, on voit que le rayon de l'univers diverge : l'application de la formule précédente donne une division par zéro.

La relation entre rayon de Hubble et rayon de l'univers observable modifier

Les équations précédentes montrent que le rayon de l'univers observable n'est pas égal au rayon de Hubble ! Il y a une différence entre et le rayon de l'univers observable . La relation entre les deux variables se calcule comme suit :

Multiplions au numérateur et au dénominateur par 1 - n :

En réorganisant les termes, on trouve :

Le terme le plus à droite n'est autre que le rayon de l'univers observable :

On peut aussi inverser la relation, ce qui donne :

Les deux rayons sont égaux seulement dans le cas où , ce qui donne .

On peut aussi remarquer que l'on a :

Le calcul du décalage vers le rouge modifier

Voyons maintenant le calcul du redshift. Dans ce qui suit, nous allons supposer qu'un objet à une distance émet de la lumière à un instant , lumière qui est reçue par l'observateur à un instant . Pour calculer le décalage vers le rouge de cette lumière, le redshift, on utilise la formule suivante, vue dans le chapitre sur le rayonnement :

, avec l'origine de temps choisie arbitrairement et le facteur d'échelle à ce même instant.

On injecte alors l'équation  :

On simplifie :

Sachant que , on a aussi :

En remarquant que est une constante, les équations précédentes se simplifient en :

Le modèle à croissance exponentielle modifier

Pour terminer, nous allons étudier le cas d'une expansion exponentielle, décrite par :

, avec le temps, et deux constantes quelconques.

Précisons que nous prenons , sans quoi la fonction serait décroissante, ainsi que , pour garantir que l'univers n'a pas de volume négatif. Voir le graphique ci-dessous pour comprendre pourquoi.

Illustration de la courbe de la fonction exponentielle.

Précisons que les propriétés mathématiques de l'exponentielle ont une conséquence assez intéressante. En effet, une exponentielle ne peut pas s'annuler, ce qui fait qu'il n'existe pas de temps t où le facteur d'échelle s'annule. En conséquence, on n'a pas d'instant où tout l'univers est rassemblé en un unique point, on n'a pas de singularité initiale. L'univers n'a donc pas d'âge, il existe depuis un temps infini. Ce modèle est donc un petit peu particulier, dans le sens où il n'y a pas de big-bang dedans. Cela peu sembler bizarre que les cosmologistes utilisent un tel modèle, mais nous verrons dans certains chapitres que ce modèle est plus réaliste qu'il n'y parait. Et à bien y penser, l'absence de singularité est plus un avantage qu'un inconvénient.

Le calcul du paramètre de Hubble modifier

Une démonstration simple part du principe que le facteur de Hubble est la dérivée logarithmique du facteur d'échelle, c'est à dire la dérivée du logarithme. Pour en faire le calcul, on commence par prendre le logarithme du facteur d'échelle, ce qui colle parfaitement avec la formule exponentielle.

Le logarithme d'un produit est égal à la somme des logarithmes des opérandes :

Le logarithme et l'exponentielle s'annulent, ce qui donne :

On prend ensuite la dérivée, ce qui donne :


Démonstration

Une démonstration alternative part de la dérivée du facteur d'échelle :

On divise par , ce qui donne :

Le paramètre de Hubble est donc la constante r de l'exponentielle. On peut donc reformuler la loi d'expansion comme suit :

Le calcul du temps et du rayon de Hubble modifier

À partir du facteur de Hubble calculé juste avant, on peut calculer le temps de Hubble :

Même chose pour le rayon de Hubble :

On voit que vu que le paramètre de Hubble est constant, le temps de Hubble et le rayon de Hubble le sont aussi. Cela peut paraître bizarre que le rayon de Hubble soit constant, mais rappelez-vous qu'il ne s'agit que du rayon au-delà duquel la vitesse de fuite d'un objet est égale à . Il n'a pas de signification physique.

Le calcul du rayon de l'univers observable modifier

Pour trouver le rayon de l'univers observable, on utilise encore une fois l'équation suivante :

On fait le remplacement avec  :

Vu que k est une constante, on la sort de l'intégrale :

On simplifie :

Le calcul de l'intégrale donne :

On reformule :

On simplifie :

Le terme n'est autre que le rayon de Hubble  :

On voit que le rayon de l'univers observable est plus grand que le rayon de Hubble.

Le calcul du paramètre de décélération modifier

À partir des résultats précédents, on peut calculer la valeur du paramètre de décélération. Pour cela, on utilise la formule vue dans le chapitre sur l'univers observable :

Le facteur de Hubble est constant, ce qui fait que sa dérivée est nulle. En injectant dans l'équation précédente, on trouve :

Ce résultat nous dit que l'expansion de l'univers accélère.

Le calcul du décalage vers le rouge modifier

Passons maintenant au calcul du redshift. Pour le calculer, on utilise la formule suivante, vue dans le chapitre sur le rayonnement :

, avec l'origine de temps choisie arbitrairement et le facteur d'échelle à ce même instant.

On injecte alors l'équation  :

On simplifie :

On voit que le redshift augmente exponentiellement.

References modifier

  1. The rh=ct universe
  2. Whe don't live in the Rh = ct universe
  3. Matter matters : unphysical properties of the Rh=ct model


Les équations de Friedmann newtoniennes

Armé de la loi de Hubble et du facteur d'échelle, nous pouvons aborder des équations très importantes de la cosmologie actuelle : les équations de Friedmann. Elles décrivent un univers de densité homogène, sans préférence au niveau de la direction (isotrope), et permettent de calculer le facteur de Hubble en fonction de la densité de l'univers et du facteur d'échelle. Elles sont au nombre de trois et portent les noms de "première équation de Friedmann", "seconde équation de Friedmann" et "équation du fluide cosmologique".

Démontrer ces équations demande d'utiliser la relativité générale, ce qui est affreusement compliqué. Heureusement, il existe une autre manière de déduire cette équation, nettement plus simple et intuitive, en utilisant la physique newtonienne. Les deux approches donnent des résultats similaires, avec cependant quelques différences que nous expliquerons plus tard. Dans ce chapitre, nous allons démontrer ces équations dans un cadre newtonien. Nous n'allons pas faire usage de la relativité et nous utiliserons seulement les outils de la physique classique. Ce qu'on verra dans le chapitre est donc assez simple et accessible.

Dans ce qui suit, les notations sont les suivantes : est la densité de matière, est la constante de gravitation et est un terme dont nous verrons la signification dans ce qui suit.

Les trois équations de Friedmann, en version non-relativiste, sont les suivantes :

Première équation de Friedmann
Seconde équation de Friedmann
Équation du fluide cosmologique de Friedmann

La première équation de Friedmann modifier

Pour démontrer la première équation de Friedmann, nous allons considérer que l'univers observable est une sphère de rayon R, de masse M, de volume V, de masse volumique homogène ρ, qui grossit à une vitesse v (la vitesse de l'expansion de l'univers). Le centre de cette sphère est situé là où se trouve l'observateur (la Terre, donc). Prenons maintenant une galaxie située sur le bord de l'univers observable, et posons que son énergie mécanique est notée . Ses énergies potentielle et cinétique sont donc respectivement égales à et , ce qui donne :

Modifions l'ordre des termes :

Nous pouvons alors utiliser la loi de Hubble , pour remplacer la vitesse de la galaxie.

Isolons maintenant en divisant des deux membres par .

Ensuite, rappelons-nous que par définition, . En posant qu'à l'instant t0, le facteur d'échelle est égal à 1, nous avons : . Nous n'allons cependant faire le remplacement que dans le second terme et pas dans le premier. Quelques simplifications ultérieures permettront de se débarrasser du rayon dans le premier terme.

Nous allons maintenant poser que est une constante, que nous allons noter . Prenez garde au signe de cette formule !

Exprimons maintenant la masse de l'univers comme étant égale à sa masse volumique multipliée par le volume : . On a alors :

Après quelques simplifications algébriques, nous obtenons la première équation de Friedmann.

La seconde équation de Friedmann modifier

Par souci de lisibilité, nous noterons la dérivée première d'une variable comme ceci : . Même chose pour la dérivée seconde, notée . De plus, toutes les dérivées sont par rapport au temps. En clair, et .

La seconde équation de Friedmann est la suivante.

Celle-ci nous dit comment varie la dérivée seconde du facteur d'échelle, qui elle-même dit si l'expansion de l'univers accélère, ralentit ou reste constante. Il s'agit donc d'une équation très importante pour comprendre la dynamique de l'expansion et comment celle-ci évolue avec le temps.

On peut la démontrer comme suit.


Démonstration

Supposons que nous soyons au centre de l'univers observable et prenons un point situé à sa surface, à la limite de l'univers observable. Ce point s'éloigne de nous à une vitesse égale à , avec le rayon de l'univers observable. Maintenant, essayons de calculer son accélération. Partons du principe que l'accélération de ce point est intégralement liée à la gravité, c'est à dire qu'il est attiré vers le centre de l'univers observable (on néglige les masses qui ne sont pas dans l'univers observable, bien que ce soit un peu de la triche). L'accélération de ce point se calcule avec la formule de Newton et vaut :

, avec M la masse totale de l'univers.

Divisons par R(t) :

La masse de l'univers observable est égale au produit de sa densité par son volume :

L’univers observable est une sphère, donc son volume est de . En injectant cette formule dans la précédente, on trouve :

En simplifiant, on trouve :

On applique ensuite la formule , valable pour toutes les distances :

En simplifiant, on a :

Une autre démonstration, équivalente, fait appel à l'équation de Poisson pour le champ gravitationnel.


Démonstration

Pour rappel, l'équation de Poisson relie l'accélération de la pesanteur à la densité du corps qui génère la gravité. L'équation en question utilise l'opérateur divergence, défini par :

Dans le cas qui nous intéresse, l'accélération n'est autre que l'accélération de l'expansion de l'univers. L'équation de Poisson est donc la suivante :

Nous avons vu il y a quelques chapitres l'équation suivante, qui donne l'accélération de l'expansion en fonction du facteur de Hubble et du rayon de l'univers r :

En combinant les deux équations, on trouve :

L'opérateur est une dérivée, ce qui fait qu'on peut utiliser la formule .

Maintenant, appliquons la définition de  :

Par définition, on a , ce qui donne :

On réorganise les termes :

L'équation du fluide modifier

Enfin, on peut aussi démontrer une troisième équation, appelée l'équation du fluide de Friedmann. Pour la démontrer, partons du principe que l'univers a une masse M égale au produit de son volume par sa densité :

Le volume de l'univers dépend du cube du facteur d'échelle, comme nous l'avons vu dans les chapitres précédents :

Divisons des deux côtés par  :

Prenons la dérivée :

Maintenant, supposons que la masse de l'univers est conservée. Vu que est une constante, l'est alors aussi, ce qui fait que sa dérivée est nulle. L'équation précédente devient donc :

Le calcul de la dérivée donne :

La dérivée d'un cube vaut  :

On divise par  :

Sachant que , on a :



Résoudre les équations de Friedmann newtoniennes

Nous avons établi les trois équations newtoniennes de Friedmann dans le chapitre précédent. Pour rappel, les voici :

Maintenant que nous avons ces trois équations, nous allons tenter de les résoudre. La première étape pour cela passe par l'équation du fluide de Friedmann. Avec cette équation, on peut trouver une équation qui relie la densité avec le facteur d'échelle. On peut ensuite injecter cette équation dans la première et la seconde équation de Friedmann, pour mieux les résoudre. Et cela est plus difficile à dire qu'à faire. Une solution générale existe, mais nous allons aussi étudier le cas simple où le terme de courbure est considéré comme nul.

La résolution de l'équation du fluide modifier

Commençons donc par la première étape, en partant de l'équation du fluide, écrite comme suit :

Divisons par  :

Il faut alors se rappeler que le facteur de Hubble dépend du facteur d'échelle, qui lui-même dépend du temps, d'après l'équation suivante, vue dans les chapitres précédents :

On utilise alors la formule  :

, avec et deux constantes d'intégration.

On réorganise les termes, ce qui donne :

, avec .

On prend l'exponentielle :

On utilise alors la formule  :

Pour conserver la cohérence des unités, correspond à une densité. Et plus précisément, l'équation n'a de sens que s'il s'agit de la densité quand . Notons-la . On a donc :

Cette équation peut se reformuler comme ceci :

Cette équation n'est autre que l'équation de conservation de la masse de l'univers.

De plus, cette équation est valable quelle que soit la courbure, ce qui est important à remarquer.

La reformulation de la première équation de Friedmann modifier

Maintenant que l'on a obtenu l'équation précédente, on peut l'injecter dans la première équation de Friedmann et dans la seconde. Commençons par la première. On a alors :

Ce qui s'écrit aussi comme suit :

Maintenant, remplaçons H par son expression dépendant du facteur d'échelle, vue plus haut :

Simplifions le terme de gauche :

On multiplie par  :

On prend la racine carrée :

L’univers dominé par la matière modifier

L'équation précédente est difficile à résoudre en raison du terme de courbure, qui pose quelques problèmes. Mais en postulant une courbure nulle, on peut aller plus loin. Le modèle cosmologique obtenu sous ces hypothèses est appelé l'univers dominé par la matière, ce qui trahit le fait que l'univers en question ne contient que de la matière, sans courbure.

La détermination de la loi d'évolution du facteur d'échelle modifier

Avec l'hypothèse d'une courbure nulle, l'équation précédente devient la suivante en réorganisant les termes :

On utilise alors la formule  :

Prendre la racine carrée est équivalent à élever à la puissance 1/2, ce qui fait qu'on peut faire le remplacement :

Le tout se reformule comme suit :

, avec .

On réorganise les termes :

On intègre :

On utilise la formule  :

La réorganisation des termes donne :

On peut simplifier le tout en :

L'équation précédente est une loi de puissance. Or, nous avons vu les modèles cosmologiques basés sur des lois de puissance dans le chapitre "Une introduction aux modèles cosmologiques", ce qui fait que l'on pourrait réutiliser les résultats de ce chapitre directement. On a en effet établi des formules pour le facteur de Hubble, le décalage vers le rouge, le rayon de l'univers observable et bien d'autres pour de tels modèles. En appliquant ces formules, on trouve les résultats du tableau suivant :

Paramètre Formule du cas général Formule dans le cas de l'univers dominé par la matière
Facteur de Hubble
Temps de Hubble
Rayon de Hubble
Rayon de l'univers observable
Age de l'univers
Paramètre de décélération

Ce qu'on peut tirer de ces formules, est surtout que l'expansion de l'univers ralentit avec le temps. En effet, le facteur de décélération est positif.

La relation entre densité et âge de l'univers modifier

On peut aussi déterminer une relation entre l'âge de l'univers et la densité dans ce modèle. Pour cela, partons de l'équation suivante :

En élevant l'équation au carré, on a :

On peut alors combiner cette équation avec l'équation suivante :

En égalisant les deux équations, on trouve :

En divisant par , on trouve :

Pour simplifier, on a :

On voit que la densité diminue avec le carré du temps et que la constante de proportionnalité est assez simple.

L'univers avec matière et courbure modifier

Après avoir vu le cas d'un univers de courbure nulle, il est temps de voir le cas général, celui où la courbure n'est pas nulle. Le cas général n'a pas de solution analytique (sous forme de formule simple), sauf dans les cas où la courbure est respectivement égale à -1, 0 et 1. Cependant, on peut comparer un univers à courbure non-nulle avec l'univers de courbure nulle et en tirer quelques conséquences. Pour faire cette comparaison, nous allons devoir parler du concept de densité critique, qui n'est autre que la densité compatible avec le facteur de Hubble pour un univers à courbure nulle.

Pour rappel, dans un univers à courbure nulle, la première équation de Friedmann est la suivante :

On peut alors isoler la densité, ce qui donne :

La densité d'un tel univers de courbure nulle s'appelle la densité critique. Vous remarquerez qu'il existe une valeur de densité différente pour chaque valeur de la constante de Hubble et qu'elle évolue donc avec le temps.

Les cosmologistes utilisent souvent le rapport entre la densité mesurée expérimentalement et la densité critique, ce rapport étant appelé le paramètre de densité. Celui-ci vaut, par définition :

À partir de cette équation, on peut montrer que est la déviation par rapport la densité critique, exprimée en proportion de la densité critique.

Rappelons que la densité varie avec le temps comme , et que cela vaut aussi pour la densité critique. En clair, le terme fait pareil. On a donc :

.

La reformulation de la première équation de Friedmann avec la densité critique modifier

Armé de ce concept de densité critique, on peut reformuler la première équation de Friedmann avec celle-ci, même dans le cas d'un univers avec courbure. Pour cela, on part de la première équation de Friedmann sous cette forme :

Par définition, on a : , ce qui donne :

Divisons par  :

On utilise alors l'équation  :

Le tout peut se reformuler comme suit :

Le destin de l'univers et le scénario du Big-Crunch modifier

Fort de l'équation précédente, trouvons les conditions pour lesquelles le facteur de Hubble s'annule à un moment bien précis. Pour cela, on pose simplement . On a alors :

On a alors :

On multiplie par  :

Maintenant, étudions pour quelles valeurs de cela peut se produire. On voit que cette équation n'a de sens que si , sans quoi le facteur d'échelle est soit nul, soit négatif.

  • Si on pose que , le facteur de Hubble s'annule pour un facteur d'échelle positif. Une analyse simple nous dit que le facteur de Hubble est positif mais décroit progressivement, avant de s'annuler quand , avant de décroitre jusqu’à ce que l'univers atteindre un volume nul.
  • Si on pose que , c'est à que l'univers a une courbure négative, les seules solutions sont un facteur d'échelle négatif, ce qui n'a aucun sens. Il n'est donc pas possible que l'expansion de l'univers s'arrête dans un univers à courbure négative.
  • Si on pose que , c'est à dire quand l'univers a une courbure nulle, l'équation précédente est une division par zéro. On peut interpréter ce résultat comme le fait que le facteur de Hubble tend vers zéro quand le facteur d'échelle est infini.
Les trois scénarios pour le destin de l'univers en fonction de la courbure k.

On voit donc que trois scénarios sont possibles, suivant la courbure de l'univers :

  • Dans le premier scénario, l'expansion de l'univers finit par cesser et s'inverse, l'univers se contracte et le volume de l'univers observable diminue. En clair, l'univers s'effondre sur lui-même dans un grand big-crunch.
  • Dans le second cas, l'expansion de l'univers ne s’arrête jamais et ralentit, mais sans pour autant que le facteur de Hubble s'annule, sauf après un temps infini. C'est le scénario du big-rip.
  • Et enfin, dans le dernier scénario, l'expansion de l'univers ne s’arrête jamais, mais celle-ci ralentit progressivement. L'univers commence par s'étendre, mais son rythme de croissance diminue peu à peu, jusqu’à s'annuler après un temps infini. Dans ce scénario, l'univers ne grossit pas indéfiniment et verra son volume tendre progressivement vers un volume maximum. Ce scénario est appelé le big chill.
Courbure Destin de l'univers Comportement de l'expansion
Positive Big-Crunch L'expansion ralentit, s'annule et s'inverse.
Nulle Big-Chill L'expansion ralentit et tend vers zéro avec le temps.
Négative Big-Rip L'expansion ralentit et mais tend vers une limite non-nulle avec le temps.
Univers à courbure positive Univers de courbure nulle Univers de courbure négative

À l'heure actuelle, il semblerait que la courbure soit nulle, ou tout du moins tellement faible qu'on peut la considérer comme nulle. Toutes les mesures, réalisées par les satellites WMAP et Planck donnent bien une valeur quasiment nulle, aux imprécisions expérimentales près. Les mesures les plus récentes, provenant du satellite Planck, nous disent qu'il y a 95% de chances pour que le paramètre de densité soit compris entre 1.0008 et −1.0029.



La masse et la relativité générale

Avant de passer au prochain chapitre, nous devons parler rapidement de la masse. Vous avez peut-être l'intuition de ce qu'est la masse d'un objet, en lien avec son poids. Vous savez aussi que la masse est ce que l'on trouve dans la formule de Newton, à savoir ce qui relie l'accélération d'un objet à la force auquel on le soumet. Mais en réalité, il faut être plus précis que cela et distinguer plusieurs sortes de masses. De plus, nous allons voir que masse et énergie ont des liens assez profonds.

Les différents types de masse modifier

Comme dit juste auparavant, il existe plusieurs sortes de masse : la masse inertielle, la masse propre et les deux masses graves (active et passive). La masse inertielle est celle de la loi de Newton, alors que les masses graves sont liées à la gravité. Ces distinctions, que nous allons voir en détail, auront une importance capitale dans de futurs chapitres.

La masse propre modifier

La masse propre est définie par la loi de Newton, à savoir par la relation :

, avec la force et l'accélération.

Dans le cadre de la relativité restreinte, la formule de Newton devient :

, avec .

L'équation précédente peut se reformuler de deux manières différentes. La première est la suivante, où l'on regroupe le terme et le terme . C'est de loin la version la plus correcte du point de vue physique, vu que les termes et le terme sont tous deux des termes cinématiques (liés au déplacement de l'objet) :

, avec .

La masse inertielle modifier

La masse inertielle est définie à partir de l'impulsion, la quantité de mouvement.

Dans la théorie de Newton, l'impulsion n'est autre que le produit suivant :

Dans la relativité restreinte, elle est égale à :

Ces deux définitions sont intéressantes, car la masse inertielle semble strictement identique à la masse propre. En effet, la force est par définition la dérivée de l'impulsion par rapport au temps. On a donc, par définition :

En supposant la masse inertielle et la masse propre constantes, les deux masses s'égalisent. Mais un point de la définition de la masse inertielle est qu'elle peut se définir à partir de la loi de la conservation de l'impulsion. Mais il existe quelques circonstances où les deux masses ne s'égalisent pas et où l'on a :

La masse grave modifier

La masse grave se décline en deux types : la masse grave passive et la masse grave active.

La masse grave passive est la masse sur laquelle la gravité agit. Rappelons que dans la théorie de Newton, la gravité agit sur les objets qui ont une masse et seulement ceux-ci. La gravité est considérée comme étant une force parmi tant d'autres. Si un objet est soumis à une force de gravité , il subira une accélération dans la direction de la force de gravité. L'accélération en question est appelée l'accélération de la pesanteur. La loi de Newton nous dit qu'il existe une masse m telle que :

La masse grave active est la source de la gravité. Un objet ayant une masse grave active attire les autres objets vers lui, sous l'effet de sa gravitation. Prenons la loi de gravitation de Newton, qui donne la force qu'un objet de masse M impose à un autre objet de masse m. La masse M est la masse grave active, alors que la masse grave passive est m. La formule est la suivante :

, avec G la constante de gravitation de Newton et r la distance entre les deux objets.

L'objet attiré subit une accélération égale à :

On voit que la masse grave active est la source de la force de gravité, et de l'accélération que les objets alentours subiront.

De l'équivalence des différentes formes de masse modifier

Dans la théorie de Newton, toutes les formes de masses sont équivalentes. Mais dans la relativité restreinte, c'est autre chose. La raison à cela est que Einstein a montré que la masse et l'énergie sont reliées par des équations assez simples.

Dans ce qui suit, la masse inertielle et la masse propre seront confondues, par simplicité. De plus, nous confondrons aussi les masses grave, passive et active. Masses inertielle et propre seront toutes deux notées , alors que la masse grave sera notée

La masse en relativité restreinte modifier

Relation masse-énergie en relativité restreinte.

Pour simplifier, quitte à dire des choses fausses, la masse et l'énergie sont deux choses presque équivalentes. L'équivalence masse-énergie d'Einstein dit que l'énergie d'un corps (ici l'univers) est reliée à sa masse. Si l'objet est en mouvement, on a :

, avec E l'énergie, M la masse et c la vitesse de la lumière et p son impulsion.

On voit que la masse n'est pas tout à fait confondue avec l'énergie. Cependant, quand la vitesse de l'objet est faible, son impulsion est négligeable. L'équation précédente devient alors :

Dans le cas d'un objet immobile, l'impulsion est nulle et l'approximation précédente devient une égalité exacte. On a alors :

On voit que l'énergie est proportionnelle à la masse, pour un objet immobile.

Maintenant, divisons l'équation précédente par le volume.

Par définition, le terme de gauche n'est autre que la densité d'énergie, alors que le premier terme de droite est la densité de masse.

On vient donc de trouver une relation entre densité de masse et densité de matière, dans le cas d'un système immobile. Ce système peut être un fluide, un solide ou tout autre corps. Mais précisons que cette relation n'est qu'approximative, dans le sens où nous avons négligé l'impulsion. Dans le cas d'un système immobile, ce n'est pas un problème, tant qu'on reste dans le cadre de la relativité restreinte. Mais Dans le cadre de la relativité générale, même un fluide immobile contient de l'impulsion, ce qui change les calculs, comme nous allons le voir dans ce qui suit.

La masse en relativité générale modifier

La définition de la masse grave est presque impossible en relativité générale. Pour comprendre pourquoi, partons de la formule newtonienne suivante, qui relie le potentiel gravitationnel à la densité de masse grave active :

Le terme de gauche regroupe tout ce qui a trait au potentiel gravitationnel. Le terme de droite regroupe tout ce qui a trait à la masse grave active : il est composé de la masse grave active, multipliée par un coefficient constant. Dans la relativité générale, la formule précédente n'est pas valable, mais il existe une formule équivalente. Elle relie un équivalent du potentiel gravitationnel, à un équivalent de la masse grave. Dans la théorie d'Einstein, il y a plusieurs potentiels gravitationnels, qui sont regroupés dans ce qu'on appelle le tenseur d'Einstein, noté . L'équivalent de la masse grave active est ce qu'on appelle le tenseur énergie-impulsion . C'est un objet mathématique qui regroupe divers termes liés à l'énergie et l'impulsion de objet gravitant. L'équation de champ d'Einstein s'écrit comme suit :

En soi, cette équation ne permet pas de définir une masse grave active, du fait de l'absence de force de gravité dans les équations d'Einstein. Mais on peut, avec quelques approximations, calculer une approximation semi-classique de la force de gravité. Ce faisant, les équations nous donnent une relation équivalente à la loi de Newton, mais tirée de la relativité générale. Le tenseur énergie-impulsion regroupe plusieurs termes, dont l'énergie, l'impulsion et la pression (pour un fluide). Chacun de ces termes joue un rôle dans la masse grave active.

En premier lieu, la masse grave active se confond avec l'énergie. L'énergie a un effet gravitationnel, peu importe que cette énergie corresponde à de la masse grave newtonienne, à de l'énergie cinétique, ou à de l'énergie potentielle. Pour simplifier, on a donc :

, avec la densité de masse grave active.

L'équation précédente n'est cependant qu'une approximation, qu'on retrouve dans beaucoup de contextes, mais qui échoue dans certains cas. En plus de l'énergie proprement dite, il faudrait aussi rajouter la pression comme sources de masse. On a donc :

, avec la pression.

En supposant que les relations sont linéaires, on devrait trouver quelque chose du genre :

, avec la densité de masse grave active.

En relativité générale, on n'est pas certain que masse inertielle, masse grave passive et masse grave active soient identiques. Elles sont censées être identiques localement, pour un objet, en vertu du principe d'équivalence. Mais les calculs ne sont pas vraiment d'accord. Dans certains cas particuliers, les trois masses sont complètement différentes. Par exemple, prenons l'exemple de la lumière, d'un photon pour être précis. La lumière n'a pas de masse inertielle, que ce soit dans la théorie de la relativité ou dans la physique de Newton. Par contre, elle dispose d'une énergie et a donc une masse grave active non-nulle (même si l'effet n'a jamais été mesuré). Elle a aussi une masse grave passive, vu que la lumière est déviée par la gravité ! On a donc un cas où masse grave et masse inertielle sont différentes. Un autre exemple est le cas des fluides, qui sont abondamment utilisés dans la cosmologie relativiste, qu'il nous faut aborder pour comprendre les chapitres qui vont suivre.

La densité de masse grave d'un fluide relativiste : démonstration intuitive modifier

La densité d'énergie d'un fluide n'est pas si difficile à comprendre. En théorie, l'énergie d'un fluide est la somme de plusieurs contributions : l'énergie de masse du fluide, son énergie cinétique et son énergie potentielle. L'énergie potentielle est souvent négligée, et ce sera le cas dans les développements qui vont suivre. Nous allons supposer que l'énergie du fluide est la somme de son énergie de masse et de son énergie cinétique. Cette simplification est peu rigoureuse, mais elle permet de comprendre ce qui va suivre sans recourir à des développements rigoureux. On a donc :

On divise alors par le volume :

Le terme de gauche est la densité d'énergie, par définition, alors que le premier terme à droite est la densité de matière.

En mécanique classique, on peut prouver que la densité d'énergie cinétique d'un fluide est proportionnelle à la pression (même si les deux ne sont pas exactement la même chose), le coefficient de proportionnalité dépendant du fluide. Par exemple, pour un gaz de photons, la pression est égale au tiers de la densité d'énergie. Pour un gaz parfait matériel, le coefficient de proportionnalité est de . Et pour de la matière solide, il est de zéro. En clair, on a :

Pour obtenir la densité de masse grave active, on peut, bien que ce ne soit pas rigoureux, appliquer bêtement la formule . Il suffit donc de diviser l'équation précédente par  :

La densité de masse grave d'un fluide relativiste : explications rigoureuses modifier

Le développement précédent n'est cependant pas rigoureux. La dérivation relativiste exacte de cette formule est assez compliquée, le concept même de masse étant ambigu dans la relativité générale. Cependant, dans certains cas, on peut définir une masse particulière, appelée la masse de Komar. Dans le cas d'un fluide idéal, la masse de Komar du fluide est égale à :

, avec E l'énergie du fluide, P sa pression et V son volume.

La densité de masse grave devrait être donc égale à :

Pour aller plus loin, les calculs relativistes donnent aussi des valeurs différentes pour la masse grave active et la masse grave passive. Les calculs en hydrodynamique relativistes donnent :

, où est la densité de matière inertielle.
, avec la densité de matière inertielle.

On voit que masse inertielle, masse grave active et masse grave passive sont différentes. On devra donc les utiliser à bon escient dans les chapitres suivants.

Quoi qu'il en soit, cela nous permet de comprendre pourquoi la démonstration de la section précédente n'était pas rigoureuse. Elle est imparfaite du fait de plusieurs erreurs majeures :

  • Elle intervertit la place de et dans les équations.
  • Elle confond et , alors que l'équation précédente nous dit que les deux valeurs ne sont pas égales.
  • Elle intègre le terme dans la densité d'énergie volumique, alors que la pression est avant tout un phénomène de surface.

La pression et la masse grave active modifier

Les équations précédentes nous montrent que la pression est une forme particulière de masse grave active. Néanmoins, la gravité liée à la pression nous pose un problème. Imaginons que l'on ait une sphère (une étoile, par exemple), remplie d'un gaz quelconque, autour de laquelle gravite un satellite. La masse grave active de cette sphère est égale à l'intégrale sur son volume de sa densité de masse grave active.

Cette dernière vaut, comme on l'a vu, . On injecte donc dans l'intégrale :

Dans le cas d'une matière sans pression, on trouve :

Mais dans le cas du rayonnement, la masse est nulle et la pression vaut le tiers de la densité d'énergie. En clair, . En calculant l'intégrale, on trouve :

On voit que la masse grave active n'est pas la même pour une densité d'énergie identique. La même quantité d'énergie n'a pas la même masse grave active, selon qu'elle est sous la forme de matière ou de rayonnement. Et cela pose un grave paradoxe. Imaginons qu'il se passe quelque chose dans la sphère, et que cette transformation transforme toute la matière en rayonnement. La masse grave active double du fait de cette transformation. Deux conséquences contradictoires sont possibles : soit la conservation de l'énergie est violée lors de la transformation, soit énergie et masse grave ne sont pas synonyme et la relativité générale est fausse. À l'heure actuelle, personne ne sait réellement comment résoudre correctement ce paradoxe, appelé le paradoxe de Tolman. Il y a bien quelques pistes, mais rien de bien folichon. Peut-être que le tenseur énergie-impulsion est incorrectement défini, peut-être que les équations mentionnées dans le chapitre sont fausses, peut-être que la théorie des fluides relativistes est à revoir, ou peut-être que tout simplement que les démonstrations actuelles oublient de prendre en compte des détails techniques subtils.


Les équations de Friedmann relativistes

Dans un chapitre précédent, nous avons vu les équations de Friedmann, mais avions fourni des versions newtoniennes de ces équations. Dans ce qui va suivre, nous allons voir les versions relativistes de ces équations. Nous allons voir comment passer des versions newtoniennes aux équations relativistes, en appliquant les concepts vus dans le chapitre précédent.

Dans ce qui suit, les notations sont les suivantes : est la densité d'énergie, est la densité de matière, est la pression et est un terme de courbure, qui décrit la géométrie de l'univers (nous reviendrons sur l'interprétation de ce terme dans ce qui suit).

Les trois équations de Friedmann relativistes sont les suivantes :

Première équation de Friedmann
Équation du fluide cosmologique de Friedmann
Seconde équation de Friedmann

Les équations de Friedmann relativistes ressemblent à leurs cousines newtoniennes, si ce n'est que la densité masse est remplacée par les équations du chapitre précédent. Mais il y a un problème : quelle est la masse utilisée dans les équations de Friedmann ? En toute généralité, les densités de masse inertielle, masse grave active et masse grave passive ne sont pas égales. Pire : leur valeur dépend de la situation considérée, ce qui fait qu'il faut faire le remplacement intelligemment. Il existe un cas où l'on sait calculer exactement ces trois densités de masse : le cas d'un fluide parfait. Aussi, on va supposer que l'univers observable est rempli d'un fluide qui regroupe la matière, le rayonnement, les galaxies, et de manière générale tout le contenant de l'univers. On postule que ce fluide est homogène, c'est à dire que sa température est la même partout, de même que la densité, la pression, et toutes les variables du même genre. On peut alors utiliser les équations vues au chapitre précédent pour faire le remplacement.

La première équation de Friedmann modifier

La première équation de Friedmann, en version non-relativiste, est la suivante :

La relation entre densité de masse et densité d'énergie se déduit de l'équivalence masse-énergie d'Einstein, même si cette démonstration est loin d'être rigoureuse. Pour convertir l'équation de Friedmann dans sa version relativiste, il suffit de remplacer la masse volumique par . Pour le terme , il faut faire l'inverse, à savoir la multiplier par , pour des raisons techniques assez compliquées.

L'interprétation de cette formule est assez simple à comprendre. Elle nous dit que l'expansion de l'univers dépend de deux forces opposées. Le premier terme correspond à l'influence gravitationnelle de la matière, alors que l'autre traduit l'effet de la courbure de l'univers.

Avant de passer à la suite, je tiens à faire remarquer une chose assez subtile, mais très importante. Dans le chapitre sur les équations de Friedmann newtoniennes, nous avions démontré la première équation de Friedmann, à partir d'un raisonnement newtonien simple. Dans cette démonstration, la masse utilisée pour calculer la densité de masse était la masse grave active, et non la masse inertielle. Le raisonnement tenu juste auparavant est donc faux. Il aurait fallu utiliser la relation . Mais la démonstration dans le cadre de la relativité générale donne bien l'équation obtenue plus haut, en remplaçant la densité de masse par la densité de masse inertielle. Cela est sans doute lié au fait que les trois types de masse sont égaux dans la théorie newtonienne.

Le paramètre de courbure modifier

La géométrie de l'univers peut prendre trois formes, selon que le paramètre de densité est nul, positif ou négatif.

Dans la démonstration newtonienne, nous avions posé pour simplifier les calculs. Mais ce paramètre a une interprétation totalement différente dans la relativité générale. Dans la relativité générale, ce terme s'appelle le paramètre de courbure. Il décrit la courbure, un paramètre qui décrit la géométrie de l'espace-temps :

  • un espace de courbure nulle a une géométrie euclidienne ;
  • un espace de courbure positive a une géométrie qui équivaut en trois dimensions à la surface d'une sphère en deux dimensions ;
  • un espace de courbure négative a une géométrie qui équivaut en trois dimensions à la surface d'une « selle de cheval » en deux dimensions (cette selle de cheval étant à une hyperbole ce que la sphère est au cercle).

Pour simplifier, la courbure décrit comment varie le facteur d'échelle dans l'espace et dans le temps. Avec un espace-temps courbe, le facteur d'échelle varie en fonction de la position et du temps. Pour simplifier les explications, imaginons que l'espace seul soit courbé. Un espace de courbure nulle a un facteur d'échelle identique partout et a donc une géométrie euclidienne. Avec un espace courbé, le facteur d'échelle varie en fonction de la position et les calculs de distance euclidiens ne marchent plus. La courbure décrit avec quelle intensité le facteur d'échelle varie en fonction de la direction, à quel point les calculs de distance varient des calculs euclidiens. Cependant, dans la relativité générale, c'est l'espace-temps qui est courbé, et pas seulement l'espace. Ce qui fait que le facteur d'échelle varie en fonction de la position, mais aussi du temps considéré. On verra dans quelques chapitres qu'une courbure positive fait décroître linéairement le facteur d'échelle avec le temps, alors qu'une courbure négative entraîne une croissance linéaire du facteur d'échelle. Mais passons.

Précisons que le paramètre de courbure décrit la courbure qu'aurait l'univers s'il était vide. Il s'agit d'un paramètre global, qui décrit l’univers tout entier. Mais il peut y avoir des variations locales, liées à la gravité. La gravité n'est autre qu'une variation de la courbure à proximité d'une masse. Toute masse grave active déforme l’espace-temps et le courbe, ce qui fait que le facteur d'échelle varie à proximité d'une masse. Le facteur d'échelle à proximité d'une masse décroît continuellement, ce qui fait que les distances se réduisent. Ainsi, un objet placé à proximité d'une masse se rapprochera progressivement de la masse considérée, non pas parce qu'elle bouge dans l'espace, mais parce que l'espace ne cesse de se contracter à proximité d'une masse grave active. On pourrait décrire la gravité comme l'expansion, avec une variation du facteur d'échelle avec le temps et la distance, ce que fait la relativité générale d’Einstein, mais cela demande une artillerie mathématique des plus élaborées. Pour résumer, on a :

Masse (grave active) Courbure locale Variation spatiale et temporelle du facteur d'échelle.

Le rayon de courbure modifier

En relativité générale, la formule pour le paramètre de courbure est différente. Déjà, les formules ne sont pas identiques et on a l'équation suivante :

, avec K la pseudo-courbure newtonienne et k la courbure telle que décrite dans la théorie de la relativité générale.

À comparer avec la formule newtonienne :

En combinant les deux équations précédentes, on trouve :

, avec le rayon de l'univers observable dans la démonstration newtonienne.

En utilisant la formule , on trouve alors :

Ce qui se simplifie en :

, avec R un rayon dont l'interprétation est donnée plus bas.

On voit que la courbure est proportionnelle à l'inverse d'un rayon. Dans la démonstration newtonienne, le rayon en question est celui de l'univers. Mais en relativité générale, ce rayon n'est pas le rayon de l'univers, mais un rayon dit "rayon de courbure".

Cercle osculateur.
Rayon de courbure.

Pour comprendre ce qu'est ce rayon, prenons une courbe en une dimension, dessinée sur un plan. En ce point, on peut décrire la courbure de cette courbe en faisant une analogie avec la courbure d'un cercle. Un cercle est en effet une ligne dont la courbure est constante. Qui plus est, plus un cercle est grand, plus sa courbure est faible (plus on courbe une ligne uniformément, plus elle forme un cercle petit). Pour la courbe prise en exemple, et sur le point choisi, on peut définir un cercle tangent en ce point, appelé le cercle osculateur. Le rayon de ce cercle. La courbure en ce point est alors définie comme la courbure du cercle tangent, qui est elle-même inversement proportionnelle au rayon du cercle osculateur. Cela permet de définir la courbure comme l'inverse du rayon du cercle osculateur.

Courbure locale de la courbe Courbure du cercle osculateur ( 1 / Rayon du cercle osculateur )

En deux, trois, quatre dimensions, on peut définir la même chose, ce qui permet de décrire un cercle osculateur et le lier avec la courbure locale/globale. Sauf que le coefficient de proportionnalité n'est pas le même, sans compter que les puissances changent. Pour résumer, la courbure de l'univers est proportionnelle à l'inverse de son rayon de courbure (plus précisément de son carré).

L'interprétation de la densité dans la version relativiste modifier

Formellement, la densité d'énergie prend en compte aussi bien la densité de matière que la densité d'énergie liée au rayonnement (lumière) , et potentiellement d'autres formes d'énergie. Pour simplifier, la densité d'énergie est la somme de la densité de matière, la densité de rayonnement, et quelques autres termes.

Il est possible aussi de rendre compte de la courbure via une densité d'énergie de courbure. En effet, nous avons vu que celle-ci est proportionnelle à une énergie dans la démonstration newtonienne. En appliquant l'équivalence masse-énergie, nous nous retrouvons donc avec une densité d'énergie de courbure dans l'équation de Friedmann. En injectant le tout dans l'équation de Friedmann, nous trouvons :

, avec

Si on pose , on peut reformuler cette équation comme ceci :

On verra dans les chapitres suivants que les densités de rayonnement et de matière dépendent du facteur d'échelle.

  • La densité de matière varie selon l'équation : .
  • Pour la densité d'énergie du rayonnement, l'équation a été vue dans le chapitre sur le rayonnement : .

La première équation de Friedmann se reformule ainsi, quand on injecte les relations précédentes dans l'équation :

Il est aussi possible de reformuler la première équation de Friedmann avec le paramètre de densité, à savoir le rapport entre la densité réelle et la densité critique obtenue quand on considére un univers de courbure nulle. Cependant, cela demande de fournir plusieurs paramètres de densité. En effet, il ne faut pas oublier l'influence différentielle du facteur d'échelle sur la matière, l'énergie de rayonnement et la courbure. Pour cela, il faut utiliser différents paramètres de densité : un pour la matière, un autre pour le rayonnement, et un autre pour la courbure. Le premier est égal à la densité de matière divisée par la densité critique et est noté . Le second est égal à la densité de rayonnement divisée par la densité critique et est noté . Même principe pour le rapport entre densité de courbure et densité critique .

L'équation du fluide de Friedmann modifier

L'équation du fluide de Friedmann est, dans sa version relativiste, reformulée de la même manière que la première équation. On remplace la densité de masse par sa valeur compatible avec la relativité, obtenue dans le chapitre précédent.

L'équation du fluide de Friedmann newtonienne est la suivante :

Dans la démonstration de l'équation de Friedmann d'il y a deux chapitres, on peut remarquer que la masse utilisée est indéterminée. Mais on retrouve l'équation relativiste si on suppose qu'il s'agit de la masse grave passive. Or, on a vu dans le chapitre précédent que celle-ci est égale, pour un fluide relativiste, à :

, où est la densité de matière inertielle.

En faisant le remplacement, on a :


Démonstration

Dans cette section, nous allons donner une dérivation de la version relativiste, qui n'est cependant pas très rigoureuse. En faire une véritable dérivation demanderait de passer par les mathématiques de la relativité générale, ce que nous ne ferons pas dans ce cours.

Nous allons considérer que l'univers est un système isolé, à savoir qu'il n'y a pas d'échange de matière ou d'énergie avec l'extérieur (il est en effet difficile de donner un sens à l'« extérieur de l'univers »). Dans ces conditions, la première loi de la thermodynamique s'applique. Celle-ci dit que toute variation de l'énergie interne de l'univers provient des variations de chaleur et de volume/pression. De plus, on va supposer que l'expansion est adiabatique : il n'y a pas d'échange de chaleur entre l'univers observable et un éventuel extérieur. Dans ces conditions, toute variation de l'énergie interne de l'univers se calcule avec la formule suivante, avec P la pression et V le volume de l'univers.

Utilisons ensuite la relation  :

Le facteur étant constant, nous pouvons le sortir de la dérivée.

Divisons par  :

Reformulons la masse de l'univers en fonction du volume et de la densité de matière

Appliquons la formule sur la dérivée d'un produit.

Factorisons maintenant le terme .

Divisons maintenant par V.

Puis divisons par dt.

Appliquons maintenant l'égalité suivante, vue dans le premier chapitre : .

On peut la récrire de la manière suivante, sous certaines hypothèses :

, où est la densité d'énergie.


Démonstration

Démonstration de la deuxième version. Partons de la première équation de Friedmann :

Et la deuxième :

Prenant la dérivée temporelle de la première des deux équations :

Que l'on peut mettre sous la forme :

Et donc :

En identifiant avec la seconde équation de Friedmann :

On bouge un peu les termes et on identifie le dernier terme avec la première équation de Friedmann :

Il reste à simplifier comme :

On trouve donc enfin :

La seconde équation de Friedmann modifier

La seconde équation de Friedmann est la suivante, dans sa forme newtonienne.

Dans la démonstration newtonienne, la masse utilisée était la masse grave active, ce qui fait qu'on doit utiliser l'équation suivante :

, avec la densité de masse inertielle et P la pression.

En combinant les deux équations précédentes, on retrouve la seconde équation de Friedmann :

Cette seconde équation de Friedmann nous dit quelque chose d'assez contre-intuitif concernant la pression : celle-ci lutte contre l'expansion. L'intuition nous dirait pourtant le contraire, mais celle-ci est trompeuse. La pression de la matière ou du rayonnement ne pousse pas sur les bords de l'univers, et ne peut donc guider l'expansion. Rappelons que l'horizon cosmologique n'est pas une barrière matérielle sur laquelle la pression pourrait pousser, mais une simple limite liée à la vitesse de la lumière. L'effet de la pression sur l'expansion ne peut donc provenir de ce mécanisme intuitif.

Pour comprendre l'effet contre-expansionniste de la pression, il faut se rappeler que la pression est proportionnelle à la densité d'énergie. Par exemple, on peut facilement démontrer que la pression d'un gaz parfait est égale aux deux tiers de sa densité d'énergie interne. Même chose pour le rayonnement, quoique le coefficient de proportionnalité soit différent. Or, rappelez-vous que l'énergie et la masse sont reliées dans la théorie de la relativité, l'équation en étant la plus simple expression. Pour simplifier, l'énergie a un poids, un effet gravitationnel qui va lutter contre l'expansion. Plus la pression et donc la densité d'énergie sont grandes, plus cet effet gravitationnel sera fort, plus l'expansion sera contrecarrée.

On peut démontrer la seconde équation de Friedmann (relativiste) à partir de la première équation et de l'équation du fluide.


Démonstration

Partons de la première équation de Friedmann, écrite comme ceci :

Multiplions des deux côtés par  :

Dérivons maintenant par rapport au temps :

Divisons par 2 les deux membres de l'équation :

Calculons la dérivée dans le second terme :

Divisons par  :

On utilise alors la relation pour simplifier le second terme :

Simplifions par  :

Divisons par  :

Or, l'équation du fluide nous dit que . Nous pouvons donc faire la substitution dans l'équation précédente, ce qui donne :

En simplifiant, nous obtenons :



L'énergie noire et la constante cosmologique

La seconde équation de Friedmann, vue au chapitre précédent, nous donne l'accélération de l'expansion de l'univers. Elle nous dit si l'expansion accélère ou ralentit, en fonction de la densité de la pression de l'univers.

Quant à la première équation de Friedmann, elle donne la vitesse de l'expansion.

Mais un terme manque dans ces équations. Et pour comprendre pourquoi, il nous faire un peu d'histoire.

L'univers statique d'Einstein modifier

Les deux équations nous disent une chose assez importante : l'univers est en expansion. Le seul moyen pour avoir un univers sans expansion avec ces équations est que la pression et la densité soient nulles, que la matière n'existe pas. Rien de surprenant en soi, mais à l'époque où ces équations ont été déduites à partir de la relativité générale, c'était une surprise. L'expansion de l'univers n'était pas encore connue et la loi de Hubble n'avait pas encore été découverte. Le consensus commun de l'époque était que l'univers était statique et un univers en expansion paraissait contre-intuitif. Il fallait réconcilier un univers statique avec ces équations, ce qu'a fait Einstein.

L'univers statique sans constante cosmologique modifier

Pour avoir un univers statique, il faut que le facteur d'échelle soit constant, ce qui implique que , et à tout instant t. En injectant ces conditions dans la première équation de Friedmann, on trouve :

Ce qui se simplifie en :

Cette condition dit qu'un univers statique est possible, à condition que la courbure compense l'effet de la matière et du rayonnement. Ce qui est assez intuitif : la présence de matière/rayonnement fait s'effondrer l'univers sur lui-même sous l'effet de la gravité. Pour compenser cela, on doit avoir une courbure qui crée une expansion exactement opposée.

En mettant les conditions , et dans la seconde équation de Friedmann, on trouve :

Ce qui est équivalent à :

L'équation nous dit que si matière il y a, sa densité doit être compensée par sa pression. La matière doit donc avoir une pression négative, ce qui n'a pas de sens physique.

On peut raisonnablement penser que la densité doit aussi rester constante, la matière n'étant pas diluée par l'expansion. On a alors , ce qui simplifie l'équation du fluide en :

.

Cette équation est incompatible avec la précédente, sauf dans le cas où et P sont toutes deux nulles.

L'impossibilité d'un univers statique avec ces équations n'est pas surprenante. En effet, les équations de Friedmann nous disent que l'expansion dépend de la courbure et de la gravité. La gravité est censée lutter contre l'expansion, en attirant toute la matière de l'univers vers un même point central, dans un cadre newtonien. La première équation de Friedmann nous dit qu'un univers statique n'est possible que si la courbure compense exactement la gravité. Mais les deux autres équations nous disent que c'est impossible, sauf à placer des contraintes non-physiques sur la matière.

L’invention de la constante cosmologique modifier

Quand Einstein vit que les équations de Friedmann impliquaient un univers en expansion, il leur ajouta un terme censé corriger ce "problème", rendre statique l'univers des équations de Friedmann. Pour obtenir un univers stationnaire (stable, sans expansion), il pris les équations de la relativité générale et ajouta un terme au bon endroit (dans le tenseur énergie-impulsion, pour les connaisseurs). Les équations de Friedmann étaient alors modifiées, et sont devenues ceci :

Le terme , ajouté aux équations de Friedmann, a reçu le nom de constante cosmologique.

L'ajout de la constante cosmologique modifie la dynamique des équations. Dans la première équation de Friedmann, la constante sert à compenser l'effet de la gravité. L'idée est que la somme courbure + constante cosmologique compense exactement la densité, ce qui annule l'expansion. En soi, on pouvait faire la même chose avec la courbure, mais avec un défaut : la courbure n'est pas présente dans la seconde équation de Friedmann, ce qui ne réglait qu'une partie du problème. Mais la constante cosmologique n'a pas ce problème. Elle est bien présente dans la seconde équation de Friedmann et peut compenser l'effet de la densité avec la valeur appropriée. On a alors :

Le calcul de la valeur de la constante cosmologique pour un univers statique modifier

La seconde équation se reformule alors comme suit :

On multiplie par  :

Le calcul de la courbure pour un univers statique modifier

Partons de la seconde équation de Friedmann avec constante cosmologique :

En injectant dans la première équation de Friedmann, on trouve :

On développe :

On regroupe les termes et on simplifie :

On factorise et on réorganise les termes :

Ce qui donne :

Le lien entre constante cosmologique et courbure dans le modèle d'Einstein modifier

Les deux sections précédentes nous ont démontré les deux équations suivantes :

Si on suppose que la pression est négligeable, on peut combiner les deux équations précédentes pour trouver :

En clair, la constante cosmologique implique que l'univers soit courbe. L'obligation d'avoir un univers courbe est un des défauts des univers statiques avec constante cosmologique, mais il est loin d'être le seul. Mais dans les grandes lignes, tous ces défauts pointent vers un même problème : pour obtenir un univers statique, la valeur de la constante cosmologique et de la courbure doivent être minutieusement calibrées pour obtenir un univers statique. L'ensemble de ces défauts font que la constante cosmologique est une réponse imparfaite pour ceux qui veulent un univers statique avec les équations de Friedmann. Aussi, Einstein a abandonné son idée de constante cosmologique pour obtenir un univers statique.

L'accélération de l'expansion modifier

Si la constante cosmologique a été mise de côté suite à la découverte de l'expansion de l'univers, les physiciens la réutilisèrent pour expliquer un autre phénomène : l'accélération de l'expansion. Il a longtemps été cru que l'expansion de l'univers devait décélérer avec le temps. Pour comprendre pourquoi, rappelons que l'accélération de l'expansion est donnée par la dérivée seconde du facteur d'échelle. Or, la seconde équation de Friedmann sans constante cosmologique nous donne justement cette dérivée seconde, qui vaut :

La densité et la pression sont forcément positives pour la matière et le rayonnement, le terme entre parenthèses est donc positif, ce qui fait que la dérivée seconde est négative en raison du signe moins. Et un négatif signifie que l'expansion décélère.

Mais quelques observations ont remis cette affirmation en cause. Deux expériences ont fait voler en éclat ce consensus, en faisant état d'une accélération de l'expansion : le Supernova Cosmology Project et le High-Z supernovae search team. Ces observations ont été confirmées par des observations faites avec le télescope Hubble, à l'heure où j'écris ces lignes (fin 2016). Autant dire que la découverte de l'accélération de l'univers remettait en question les équations de Friedmann telles qu'elles sont formulées dans les chapitres précédents.

Illustration schématique de l'accélération de l'expansion de l'univers.

Pour résoudre le problème, les physiciens ont envisagé des formes de matière de densité ou de pression négatives. Si l'hypothèse d'une densité/masse négative n'a pas été retenue, ce n'est pas le cas de l'autre possibilité. Puis, les physiciens se sont rendus compte que l'introduction d'une constante cosmologique permet de rendre compte de l'accélération de l'expansion, sous certaines conditions. En effet, dans la seconde équation de Friedmann, la constante cosmologique apparaît comme un terme positif, qui permet de faire passer la dérivée seconde au-dessus de zéro.

La dérivée seconde est positive sous la condition suivante :

L'énergie noire modifier

Contenu de l'univers.

La constante cosmologique est encore un mystère. On ne sait pas à quel phénomène physique elle correspond, quelle est son origine. Pour l'expliquer, les cosmologistes supposent qu'il existe, dans l'univers, quelque chose qui a une énergie positive, mais dont la pression est négative. Ce quelque chose, ils l'ont appelé l'énergie noire. Cette énergie noire assez mal connue serait quelque chose de physique, qui contient une certaine densité d'énergie, sous la forme de masse et de pression. Mais les caractéristiques de cette énergie noire font qu'elle ne se comporte pas comme de la matière ordinaire ou du rayonnement. Les observations semblent indiquer que 30% de l'univers est composé de matière, tandis que le reste de la densité d'énergie est intégralement composé d'énergie noire.

Les propriétés de l'énergie noire modifier

L'énergie noire a une masse positive, mais a aussi une pression négative. Or, la relativité générale nous dit que la pression est une forme d'énergie, qui a donc un effet gravitationnel. Dans cette théorie, toute énergie a un poids, une sorte de masse fictive, et la pression ne fait pas exception. Mais pour la matière ordinaire, la pression est positive, c'est à dire qu'elle a un effet gravitationnel normal, similaire à une masse positive. Elle a donc tendance à ralentir l'expansion, dans une certaine mesure. Mais l'énergie noire a une pression négative, ce qui est particulièrement original (même si on observe des phénomènes similaires dans certains cristaux ou lors de changements de phases). Du point de vue de la relativité générale, une pression négative implique un effet anti-gravitationnel, une attraction gravitaire dirigée vers l'extérieur. En clair, elle agit comme une force répulsive qui contrecarre la gravité et augmente l'expansion au lieu de la ralentir.

L'énergie noire se caractériserait par une densité constante au cours du temps, même en tenant compte de l'expansion de l'univers. La constance de la densité d'énergie noire a une conséquence assez intéressante : la densité étant constante, le facteur de Hubble est constant lui aussi. On peut remarquer que cela implique que l'énergie de l'univers ne se conserve pas, vu que cette densité constante est couplée à une augmentation du volume de l'univers par l'expansion.

Les scientifiques ont tenté d'expliquer cette énergie noire de plusieurs manières. La tentative la plus connue se base sur une prédiction assez singulière de la théorie quantique des champs., à savoir que le vide est censé être rempli d'une énergie appelé énergie du vide. C'était le candidat idéal pour expliquer l'énergie noire, jusqu'à ce que les calculs montrent que l'expansion mesurée implique une énergie du vide près de fois plus grande qu'observée. Ce résultat est aujourd'hui considéré comme la "pire prédiction théorique de l'histoire de la physique".

Énergie noire et équations de Friedmann modifier

La constante de la densité d'énergie noire permet de simplifier la première équation de Friedmann, en utilisant la densité d'énergie noire .

Rappelons que la densité de matière, de rayonnement et de courbure suivent les relations suivantes : , et . En substituant dans l'équation de Friedmann, on trouve :

Par définition, . En faisant le remplacement dans l'équation précédente et en faisant quelques manipulations algébriques, on trouve :

On multiplie par  :

À l'heure actuelle, la constante cosmologique n'est qu'un paramètre ad-hoc qui permet de faire coller prédictions théoriques et observations. On ne connaît pas encore son origine, même si les physiciens ont quelques pistes. L'interprétation de ce terme varie suivant les physiciens.

Il est possible de reformuler la première équation de Friedmann avec le paramètre de densité, à savoir le rapport entre la densité réelle et la densité critique obtenue quand on considère un univers de courbure nulle.



Les équations de Friedmann pour un univers vide

Grâce aux chapitres précédents, nous sommes enfin armés des équations de Friedmann, qui décrivent l'expansion de l'univers en fonction de sa masse, de sa courbure et de la constante cosmologique. Pour rappel, ces équations sont les suivantes :

Il est possible de résoudre ces équations dans des cas particuliers, le cas général étant trop compliqué. Autant on peut faire des simulations informatiques avec les équations générales, autant dériver des formules exactes n'est pas à notre portée. Pour se simplifier la vie, il est possible de négliger certains termes dans les équations. On peut par exemple négliger la constante cosmologique, ou négliger la courbure, voire négliger la pression. C'est ce que nous allons faire dans ce chapitre, ainsi que dans le suivant.

Dans ce chapitre, nous allons voir ce qu'il en est pour les univers vides, c'est à dire sans matière dedans, mais avec une constante cosmologique et/ou une courbure. Dans les équations, il suffit de poser que la densité et la pression sont nulles. Sous ces conditions, l'équation du fluide disparaît et les deux autres équations de Friedmann deviennent alors :

La disparition de l'équation du fluide est assez simple à comprendre. Cette équation décrit comment la matière et le rayonnement sont dilués par l'expansion. Mais sans matière et rayonnement, il n'y a rien à diluer avec l'expansion, et donc cette équation devient inutile. La dilution de la matière et du rayonnement avec l'expansion n'existant pas, on devine que cela a des conséquences sur la dynamique de l'expansion.

Résoudre ces équations dans le cas général, avec une courbure non-nulle et une constante cosmologique, est assez difficile et donne des équations relativement complexes. Mais il est beaucoup plus facile d'obtenir des résultats simples dans le cas où la courbure vaut soit 1, soit 0, soit -1. On obtient alors les résultats suivants :

Dans ce qui va suivre, nous allons étudier des cas particuliers, pour lesquels obtenir des résultats est assez facile. Dans ce qui va suivre, nous allons étudier des cas particuliers qui sont faciles à traiter, à savoir :

  • un univers sans courbure et sans constante cosmologique ;
  • un univers avec courbure, mais sans constante cosmologique ;
  • et un univers sans courbure, mais avec une constante cosmologique.

On peut se demander pourquoi étudier de tels univers hypothétiques, alors que l'univers n'est évidemment pas vide. Mais leur étude montre que les univers vides ont des propriétés intéressantes, qui ne sont pas respectées par les univers remplis de matière et/ou de rayonnement. Par exemple, les univers avec une courbure ont une croissance linéaire et sont les seuls dans ce cas ! Bref, commençons.

Le modèle cosmologique dominé par la courbure modifier

Le premier modèle que nous allons voir est celui d'un univers vide, sans constante cosmologique. Les équations de Friedmann se simplifient alors fortement. L'équation du fluide et la deuxième équation disparaissent complètement, et il ne reste que la première. Voici le résultat :

On peut alors utiliser la formule  :

On multiplie par  :

On peut alors résoudre cette équation pour trouver une formule du type . Cependant, la résolution est légèrement différente selon que l'on ait une courbure positive ou une courbure négative. Il faut dire que résoudre l'équation précédente demande de prendre une racine carrée, ce qui se marie assez mal avec les nombres négatifs. Dans ce qui va suivre, nous allons d'abord voir le cas à courbure positive et nulle, avant de voir le cas à courbure négative.

L'univers dominé par une courbure nulle modifier

Si on suppose une courbure nulle, on a :

En intégrant, on trouve que :

, avec la constante .

Traduit en langage commun, cela signifie que l'univers est statique, sans expansion. Mais ce cas est celui d'un univers complètement vide, sans courbure, ni constante cosmologique, ni matière, ni rayonnement, ni quoique ce soit. La présence de matière ou de tout autre composant rend la solution instable et brise la stabilité de l'univers, sauf coïncidence extraordinaire et/ou choix bien précis de paramètres. Un univers statique n'existe donc que si l'effet de la matière est compensé par quelque chose ayant un effet inverse de même ampleur. Raison pour laquelle Einstein avait inventé la constante cosmologique, mais les calculs ne collaient pas vraiment...

L'univers dominé par une courbure positive modifier

Pour étudier le cas d'une courbure positive mais non-nulle, repartons de l'équation précédente :

Prenons la racine carrée de l'équation précédente :

Pour simplifier les calculs, on utilise la constante K telle que  :

En intégrant, on trouve que :

L'univers est ici en décroissance linéaire.

L'univers dominé par une courbure négative modifier

Pour étudier un univers à courbure négative, nous allons partir de l'équation suivante :

En supposant que la courbure est négative et est égale à , on a alors :

Ce qui se simplifie en :

En prenant la racine carrée, on a :

Là encore, on pose  :

En intégrant, on trouve :

En négligeant la constante d'intégration , on a alors :

Traduit en langage commun, cela veut dire que l'univers est en expansion linéaire, à rythme constant. Le modèle obtenu n'est autre que le modèle à expansion linéaire que nous avions étudié dans le chapitre "Introduction aux modèles cosmologiques".

Le modèle cosmologique dominé par l'énergie noire modifier

Le cas suivant est celui d'un univers dominé par l'énergie noire, c'est à dire un univers avec de l'énergie noire, mais sans matière, ni rayonnement, ni courbure, ni quoique ce soit d'autre. Un tel univers hypothétique est appelé un univers de De Sitter. Sous ces hypothèses, les équations de Friedmann deviennent celles-ci :

Tous les paramètres de cette équation, , c et les autres, sont tous constants, ce qui fait que le facteur de Hubble l'est aussi. Le fait que le facteur de Hubble est constant est à la base des développements mathématiques qui vont suivre.

Le calcul du facteur d'échelle modifier

La constance du facteur de Hubble a des conséquences assez intéressantes, qui peuvent se déduire de la définition du facteur de Hubble. Pour rappel, la définition du facteur de Hubble en fonction du facteur d'échelle donne :

Multiplions par des deux côtés :

Intégrons sur  :

Dans le premier terme, H est une constante, ce qui donne :

Dans le second terme, l'intégrale n'est autre que l'âge de l'univers. Dans ce qui suit, nous noterons l'âge de l'univers . Le calcul de l'intégrale de droite donne :

Prenons l'exponentielle des deux côtés :

En injectant l'équation , on a alors :

Remarques sur ce modèle modifier

Dans ce modèle, le facteur d'échelle suit une relation exponentielle avec le temps. Il s'agit donc d'un cas particulier de modèle à croissance exponentielle, ces derniers ayant déjà été abordés dans le chapitre "Une introduction aux modèles cosmologiques". Tout ce qui a été dit sur ces modèles peut donc s'adapter à ce modèle d'univers dominé par l'énergie noire. Par exemple, leur facteur de décélération est toujours positif, égal à -1. Un tel univers va donc une croissance exponentiellement croissante, ce qui est beaucoup plus rapide que dans les deux prochains modèles, où la croissance est plus lente. Ce plus, on sait aussi qu'un tel l'univers a un âge infini : vu qu'une exponentielle ne peut être nulle, le facteur d'échelle et le volume de l'univers ne peuvent pas être nuls. Dit autrement : la singularité initiale est totalement évitée !



Les équations de Friedmann pour un univers plat et sans constante cosmologique

Par souci de lisibilité, nous noterons la dérivée première d'une variable comme ceci : . Même chose pour la dérivée seconde, notée . De plus, toutes les dérivées sont par rapport au temps. En clair, et .

Après avoir vu le cas d'un univers vide, sans matière, nous allons maintenant étudier le cas d'un univers remplit de matière et/ou de rayonnement. Cependant, nous allons supposer que l'univers n'a pas de courbure, afin de simplifier les calculs. De plus, nous allons supposer que la constante cosmologique est nulle. En clair, nous allons étudier le cas d'un univers plat (c.a.d sans courbure et sans constante cosmologique. Avec un tel univers, la première équation de Friedmann se simplifie et devient :

Les deux autres équations de Friedmann ne sont pas touchées, mais il reste à les résoudre. Une solution pour simplifier les calculs consiste à réduire le nombre d'inconnues à un seule. Au lieu de travailler avec la pression et la densité, il est possible de ne travailler qu'avec la densité. Cela demande de postuler une relation entre la densité d'énergie et la pression, à savoir une équation d'état. Généralement, on postule que la pression est proportionnelle à la densité d'énergie. La relation entre pression et densité est donc de la forme :

La valeur de dépend selon que l'on considère la matière, le rayonnement ou l'énergie noire.

  • Pour la matière, on part du principe que celle-ci est un gaz parfait qui emplit l'espace. La pression d'un tel gaz est alors définie par la formule , avec la densité du gaz et v la vitesse moyenne de ses particules. Cependant, par souci de simplification, il est supposé que la vitesse des particules du gaz est très faible, au point qu'on peut la supposer nulle (ce qui marche bien pour de la matière qui va à faible vitesse). Dans ces conditions, la pression s'annule, peu importe la densité. On a alors : .
  • Pour le rayonnement, il est établi par la physique du rayonnement que .
  • Pour l'énergie noire, on sait que la pression est négative pour une densité positive. On sait donc que , mais guère plus. Pour simplifier les calculs, on suppose que . Empiriquement, les mesures réalisées par le satellite Planck semblent compatibles avec la valeur w = −1.028 ± 0.032.

Le cas général, avec w indéterminé modifier

Pour commencer, nous allons étudier le cas général, dans lequel la valeur de n'est pas précisée. Nous allons étudier comment les équations de Friedmann se reformulent quand on injecte l'équation à l'intérieur. Les trois équations obtenues sont donc :

Remarquons que la première équation de Friedmann ne change pas, vu qu'elle n'a pas de terme de pression.

Pour obtenir cette version de la seconde équation de Friedmann, partons de la seconde équation de Friedmann originelle :

Injectons la formule  :

On simplifie par et on factorise  :

La détermination de l'équation de la densité en fonction du facteur d'échelle modifier

Pour la troisième équation de Friedmann, on part de l'équation relativiste :

On fait le remplacement :

On simplifie par , et on factorise  :

On divise par et on réorganise les termes :

On utilise alors la relation , qui relie facteur de Hubble et facteur d'échelle :

On utilise alors la formule  :

, avec et deux constantes d'intégration.

On réorganise les termes, ce qui donne :

, avec .

On prend l'exponentielle :

On utilise alors la formule  :

Pour conserver la cohérence des unités, correspond à une densité. Et plus précisément, l'équation n'a de sens que s'il s'agit de la densité quand . Il s'agit d'une densité de base que l'on notera . On a donc :

De plus, cette équation est valable quelle que soit la courbure, ce qui est important à remarquer.

La détermination de l'équation du facteur d'échelle modifier

Maintenant, essayons de trouver une équation pour le facteur d'échelle en fonction du temps. Pour cela, partons de l'équation précédente, reformulée comme suit :

Maintenant, prenons la première équation de Friedmann écrite comme suit :

En combinant les deux, on a :

On prend alors la racine carrée, ce qui est équivalent à élever à la puissance 1/2. Cela donne l'équation suivante :

On utilise alors la formule  :

On multiplie par a(t), ce qui revient à ajouter 1 à l'exposant pour le terme de droite :

On simplifie l'exposant :

Résoudre cette équation différentielle dans le cas général est quelque peu complexe. Toujours est-il que le résultat est le suivant :

Les équations générales du modèle modifier

L'équation précédente nous dit que le facteur d'échelle suit une loi de puissance, ce qui fait qu'on peut utiliser les formules du chapitre "Une introduction aux modèles cosmologiques". On trouve donc les résultats suivants :

La relation entre densité et âge de l'univers modifier

On peut aussi déterminer une relation entre l'âge de l'univers et la densité. Pour cela, partons de l'équation suivante :

En élevant l'équation au carré, on a :

On peut alors combiner cette équation avec l'équation suivante :

En égalisant les deux équations, on trouve :

En divisant par , on trouve :

Pour simplifier, on a :

On voit que la densité diminue avec le carré du temps, mais que la constante de proportionnalité dépend du paramètre de l'équation d'état.

Le facteur de décélération et les équations de Friedmann modifier

Maintenant, étudions le comportement de l'accélération de l'univers, à partir du facteur de décélération. Nous avons vu plus haut quelle était la valeur du facteur de décélération dans le cas général, mais sachez qu'il est possible de démontrer la valeur du facteur de décélération autrement.


Démonstration

Pour rappel, la définition du facteur de décélération est la suivante :

Le terme de gauche se calcule à partir de la seconde équation de Friedmann, alors que se calcule à partir de la première. Dans le cas qui nous intéresse, on a alors :

, vu que la courbure est considérée comme nulle.
, du fait de l'hypothèse .

On a alors :

On simplifie par  :

D'où on dérive :

Pour rappel, ce dernier nous dit sous quelles condition l'expansion ralentit, accélère ou reste stable. Pour cela, calculons la valeur de pour laquelle . On a donc :

On trouve alors que :

En étudiant un petit peu la relation , on en déduit le tableau suivant.

L'expansion accélère L'expansion est constante (modèle à croissance linéaire) L'expansion ralentit

L'accélération de l'univers ne peut donc s'expliquer que si l'énergie noire a une équation d'état de la forme :

Faisons maintenant le rapprochement avec les valeurs de obtenues pour le rayonnement, la matière et l'énergie noire. Pour un univers composé uniquement d'énergie noire, on déduit que l'expansion accélère sans cesse, vu que l'on a . Mais pour la matière et le rayonnement, on a , ce qui fait qu'un univers composé intégralement de matière et de rayonnement doit voir son expansion décélérer. Pour un univers qui contient un mélange d'énergie noire et de matière/rayonnement, tout dépend des parts respectives de chaque composant.

Le modèle cosmologique dominé par la matière modifier

Le cas de l'univers qui ne contient que de la matière, sans rayonnement, ni constante cosmologique est le second cas que nous allons aborder. Dans ce modèle, la matière est un gaz parfait dont les particules sont des galaxies ou des amas de galaxies. Cette hypothèse est crédible dans le sens où les amas de galaxies sont relativement éloignés et interagissent peu. Comme autre simplification, nous allons prendre le cas d'une matière froide, au zéro absolu. Cette autre hypothèse n'est pas si abusive vu l'état actuel de l'univers : seul 10% de la matière sert à fabriquer des étoiles, le reste étant localisé dans des nébuleuses et des nuages moléculaires dont la température ne dépasse pas la dizaine de degrés au-dessus du zéro absolu. L'univers est donc vraiment très froid ! En appliquant la loi des gaz parfaits avec une température au zéro absolu, on trouve que la pression est nulle quelle que soit la densité. Dit autrement, .

Avec cette hypothèse, les équation de Friedmann deviennent les suivantes :

On remarque qu'on retrouve les équations de Friedmann newtoniennes dans le cas où la courbure est nulle. Et nous avions vu que ce modèle permettait de démontrer les formules suivantes :

Paramètre Formule dans le cas de l'univers dominé par la matière
Densité en fonction du facteur d'échelle
Facteur d'échelle en fonction du temps
Facteur de Hubble
Temps de Hubble
Rayon de Hubble
Rayon de l'univers observable
Age de l'univers
Paramètre de décélération
Densité

Le modèle cosmologique dominé par le rayonnement modifier

Dans le cas où on considère un univers entièrement rempli de rayonnement, on, postule que le rayonnement est formé d'un gaz parfait de photons. Dans ces conditions, le comportement des photons fait que :

Avec cette valeur, les équations du cas général se simplifient comme suit :

Paramètre Formule du cas général Formule dans le cas de l'univers dominé par le rayonnement
Densité en fonction du facteur d'échelle
Facteur d'échelle en fonction du temps
Facteur de Hubble
Temps de Hubble
Rayon de Hubble
Rayon de l'univers observable
Age de l'univers
Paramètre de décélération
Densité

Les équations précédentes disent que l'univers a une expansion plus rapide qu'avec l'univers dominé par la matière. La raison est que le rayonnement se dilue plus vite que la matière, comme le dit l'équation , en raison de la dilution par l'expansion de l'univers, mais aussi de l'expansion des longueurs d'onde. La dilution étant plus grande, la gravité perd plus rapidement en importance et l'expansion est donc plus rapide.

On peut noter que le rayon de Hubble et le rayon de l'univers observable sont égaux dans l'univers dominé par le rayonnement.

Le paramètre de décélération étant positif, on en déduit que l'expansion de l'univers ralentit progressivement avec le temps.



Le destin de l'univers

Dans les chapitres précédents, nous avons établi les deux équations de Friedmann. Pour rappel, les voici :

On peut déduire bien des choses à partir de ces équations, notamment comment le facteur de Hubble évolue au cours du temps. Cela nous permet de déduire ce qu'il adviendra de l'univers. Continuera-t-il à s'étendre indéfiniment ? Ou au contraire, l'expansion cessera-t-elle au bout d'un certain temps ? L'univers finira-t-il par s'effondrer sur lui-même ? Il n'y a pas 36 possibilités et seuls trois scénarios sont possibles :

  • Dans le premier scénario, l'expansion de l'univers s'inverse après un certain temps, l'univers se contracte et s'effondre sur lui-même : c'est le big-crunch.
  • Dans le second cas, l'expansion de l'univers ne cesse jamais : c'est le scénario du big-rip.
  • Et enfin, dans le dernier scénario, l'expansion de l'univers ne s’arrête jamais, mais celle-ci ralentit progressivement jusqu’à s'annuler après un temps infini : c'est le scénario du big chill.
Trois possibilités pour l'évolution de l'univers.

Dans le second chapitre sur les équations de Friedmann newtonienne, nous avions vu que le destin de l'univers dépendait uniquement de sa courbure, mais pas de son contenu en matière. Et autant ce résultat marche avec seulement de la matière, autant cela change quand on ajoute la constante cosmologique. Nous allons refaire les calculs et l'analyse dans le cas général, avec les équations de Friedmann relativistes.

Facteur de Hubble, facteur d'échelle et destin de l'univers modifier

Le scénario qui se matérialise dépend du facteur de Hubble. Rappelons en effet que le facteur de Hubble est le taux d'expansion de l'univers, le taux auquel l'univers augmente de volume. Un H négatif signifie que l'univers se contracte, un H positif signifie que l'univers est en expansion, et un H nul signifie que l'univers est stationnaire. Pour l'étude du destin de l'univers, on se préoccupe du facteur de Hubble obtenu après un temps assez long, pour un âge de l'univers très important. Idéalement, on doit étudier la limite de H quand le temps tend vers l'infini : . On distingue les trois scénarios précédents selon que la limite de H est positive, négative ou nulle.

  • Positive : l'expansion de l'univers ne s’arrête jamais et c'est le big rip qui se matérialise.
  • Nulle : le big chill se matérialise si la limite tend vers 0.
  • Négative : l'expansion s'inverse si la limite devient négative et l'univers finit en big-crunch.

Il est aussi possible de voir les choses à partir du facteur d'échelle, en utilisant la dérivée du facteur d'échelle. Quand , le facteur d'échelle augmente au cours du temps, ce qui implique un univers en expansion. Par contre, implique un facteur d'échelle qui se réduit au cours du temps et donc un univers qui se contracte. Enfin, implique un univers stable, qui n'est ni en expansion ni en contraction.

En théorie, on devrait étudier le destin de l'univers, ce qui veut dire son état quand l'âge de l'univers est très avancé, voire à la limite . Mais les équations précédentes sont écrites avec le facteur d'échelle, pas avec le temps. Cependant, le facteur d'échelle dépend du temps et on peut l'utiliser en remplacement, sous certaines hypothèses.

  • Dans le cas du big crunch, l'univers gonfle avant de se rétracter. Au niveau du facteur d'échelle, cela signifie que le facteur d'échelle est initialement croissant, cesse de croître à un temps bien précis, et décroît au-delà. Au niveau du facteur de Hubble, cela signifie qu'il est d'abord positif (univers en expansion), puis s'annule, et enfin devient négatif. On est donc dans le cas où la dérivée du facteur d'échelle s'annule pour un temps fini.
  • Dans le cas du big chill, la dérivée s'annule, mais seulement pour la limite . L'univers ne cesse de croître et ce n'est qu'après un temps infini qu'il se stabilise.
  • Dans le cas du big rip, la dérivée reste positive quand on passe à la limite . L'univers ne cesse de croître et ne se stabilise pas, même après un temps infini.
Big Crunch Négative Nulle
Big Chill Nulle Infinie
Big Rip Positive

On peut calculer la dérivée du facteur d'échelle en utilisant la formule . Et donc, on peut la calculer indirectement à partir de la première équation de Friedmann. Celle-ci s'écrit en effet, comme on l'a vu dans le chapitre sur l'énergie noire :

Pour rappel, la densité d'énergie est égale à la densité de matière et la densité de rayonnement :

Pour rappel, la densité de matière et de rayonnement évoluent comme suit : et . En injectant dans l'équation précédente, on a :

On simplifie :

On prend alors la limite quand t tend vers l'infini :

En sortant les termes indépendants de t de la limite, on a :

Le cas du big crunch modifier

Dans le cas du big crunch, l'univers décroît au-delà d'un certain temps et se ratatine sur lui-même au point d'atteindre un volume nul. Traduit mathématiquement, cela se traduit par les équations suivantes :

On part de l'équation des sections précédentes :

On combine les deux équations précédentes :

On réorganise :

Or, on a vu que le big crunch impliquait que : . En injectant dans l'équation précédente, on trouve :

Cette inégalité nous dit que le big crunch n'est possible que si l'univers a une courbure positive et que celle-ci a une intensité suffisante. Il faut alors que l'effet de la constante cosmologique et de la gravité soient compensés quand l'âge de l'univers a atteint une valeur suffisante.

Le cas du big rip modifier

Pour le big-rip, on a, par définition du big rip :

.
.

On part de l'équation des sections précédentes :

On combine les deux équations précédentes :

On réorganise :

On utilise alors la condition . On en déduit que les facteurs et deviennent très petits, au point d'être négligeables, et vont s'annuler avec l'augmentation progressive de . On a alors :

L'équation nous dit que la présence d'une constante cosmologique non-nulle suffit à elle seule pour que le big-rip survienne. Mais dans le cas où elle est nulle, le big rip dépend de la valeur de la courbure. Supposons que la constante cosmologique soit nulle, ce qui donne :

Ce qui se simplifie en :

On voit que le big rip est possible si , mais à condition que la courbure soit négative.

Pour résumer, le big rip a lieu soit si la constante cosmologique est positive, soit si la courbure est négative (et même suffisamment négative pour compenser l'effet de la gravité/densité). Dit autrement, le destin de l'univers dépend uniquement de la valeur de la constante cosmologique et de la courbure.

Le cas du big chill modifier

Dans le cas du big chill, le facteur d'échelle tend vers l'infini et sa dérivée tend vers zéro :

.
.

On part de l'équation des sections précédentes :

On combine les deux équations précédentes :

On réorganise :

On ajoute alors la contrainte . Dans ce cas, on a et . On a alors :

Ce qui n'est possible que si la constante cosmologique est nulle, idem pour la courbure. Un univers en big chill ne peut atteindre un volume infini que si sa courbure et sa constante cosmologique sont toutes les deux nulles, en présence de matière et/ou de rayonnement.

Résumé modifier

Dans les paragraphes précédents, on sait quel est le destin de l'univers en fonction de son contenu. Tout dépend de la valeur de la constante cosmologique et de sa courbure, la présence de matière et de rayonnement important peu. On peut donc résumer le tout dans un tableau à deux entrées, qui donne le destin de l'univers en fonction du signe de la courbure et de la constante cosmologique.

Nul Positif
Courbure négative Big Rip Big Rip
Courbure nulle Big chill
Courbure positive Big Crunch Dépend de qui l'emporte



Le fluide cosmologique

On peut modéliser l'ensemble matière + rayonnement comme un fluide non relativiste (en clair, un fluide dont la vitesse des particules est petite devant la vitesse de la lumière). Un tel fluide est gouverné par les équations de la mécanique des fluides, et on pourrait penser qu'il suffit de réutiliser les équations de la mécanique des fluides telles quelles. Cependant, les équations de la mécanique des fluides usuelle ne s'appliquent pas sans modifications, du fait de l'expansion de l'univers.

Pour comprendre pourquoi, il faut savoir que les équations de la mécanique des fluides décrivent l'évolution d'un petit volume infinitésimal de fluide, appelé particule de fluide. Ce volume est pris très très petit, de manière à donner l'illusion que le fluide est un milieu continu. Dans le cas de l'univers, le très très petit signifie une taille petite comparée à l'horizon cosmologique, ce qui fait que les volumes des particules fluides peuvent être extrêmement grands comparés à la normale. Dans ces conditions, tenir compte du volume des particules fluides est absolument primordial. Or, les équations de la mécanique des fluides "usuelle" sont démontrées en partant du principe que ce volume infinitésimal de fluide reste constant. Ce serait le cas si l'expansion de l'univers n'existait pas, mais l'expansion fait que ce volume augmente avec l'expansion. Et il en est de même avec la surface des particules fluides, ainsi que des distances parcourues par le fluide... Autant dire qu'il faudra donc redémontrer les équations de la mécanique des fluides pour tenir compte de ce phénomène.

Rappels sur les opérateurs vectoriels
Le gradient est aux vecteurs ce que la dérivée est aux scalaires/fonctions. Dans le cas des nombres, les deux se confondent d'ailleurs.
La divergence est la somme des termes du gradient :

On peut remarquer que sa notation explicite bien la nature de la divergence. L'opérateur n'est autre que l'opérateur gradient, alors que l'opérateur signifie que l'on fait la somme de tous les termes d'un vecteur. Attention : l'opérateur se place après le dans l'écriture, mais il est en réalité appliqué après que le gradient ait été calculé.

La divergence est un opérateur dit linéaire, ce qui veut dire que :

  • La divergence d'une somme est la somme des divergences.
  • La divergence d'un multiple est le multiple de la divergence.

Si f est un scalaire constant et uniforme, la formule précédente se simplifie et devient alors :

Le laplacien est tout simplement la divergence du gradient. Ce qui revient à appliquer deux fois l'opérateur gradient, avant d'additionner les termes du vecteur obtenu.

Les équations de conservation "usuelles" modifier

Les trois équations de la mécanique des fluides sont des équations de conservation. Elles portent d'ailleurs les noms d'équation de conservation de la masse, d'équation de conservation de l'énergie et d'équation de conservation de la quantité de mouvement. Elles disent que la masse, l'énergie et la quantité de mouvement se conservent lors de l'évolution d'un fluide. Démontrer ces équations de conservation demande de simplement faire le bilan des variations, entrées et sorties de la particule fluide.

Variation interne + Entrées/sorties = ...

Le terme de variation interne est naturellement un terme de la forme :

Le flux entrant/sortant d'une quantité se calcule en multipliant la vitesse par un opérateur appelé divergence, qui donne l'intégrale du flux entrant à travers la surface de la particule fluide :

Pour résumer, toute équation de la mécanique des fluides prend la forme suivante. Elle décrit une quantité A, entraînée à la vitesse v, et dont une quantité S est produite à chaque instant  :

La conservation de la masse modifier

Par exemple, l'équation de conservation de la masse s'établit en faisant la somme des entrées/sorties de masse du volume (flux entrant/sortant) et de la variation de la masse de la particule fluide (variation interne). Vu que la masse se conserve, toute variation interne de masse de la particule doit provenir d'un flux entrant ou sortant : la masse de la particule diminue avec un flux sortant et augmente avec un flux entrant. La somme des deux termes est nulle, ce qui traduit le fait que toute variation de ces grandeurs provient d'un flux égal à travers la surface du volume.

Vu qu'il n'y a pas de source de masse, on a l'équation suivante :

Variation interne = Flux entrant/sortant, ou encore : Variation interne + Entrées/sorties = 0.

On a donc :

En divisant le tout par le volume infinitésimal, on a :

On peut sortir la densité des parenthèses :

La conservation de la quantité de mouvement modifier

Pour la quantité de mouvement, les choses sont plus compliquées. Le fait est que du point de vue de la particule de fluide, les forces sont des sources de quantité de mouvement, qui donnent l'illusion que de la quantité de mouvement est créée. La somme entre variation interne et flux ne peut pas être nulle en présence de forces, celles-ci s'ajoutant à la quantité conservée de quantité de mouvement. Dans ce cas, la variation de quantité de mouvement et son flux sont causées par des forces, ce qui donne :

Variation interne + Entrées/sorties = Forces.

On suppose l'existence de deux forces : la force causée par le gradient de pression et la gravité. L'accélération liée à la gravité est elle-même liée à un gradient : le gradient du potentiel gravitationnel . On a donc :

La masse étant constante, on peut simplifier par celle-ci :

Les équations au complet modifier

Il y a une équation de conservation pour la masse et la quantité de mouvement, à laquelle on pourrait ajouter une équation de conservation de l'énergie. À ces équations, il faut ajouter une équation qui décrit le champ de gravité.

En supposant que le volume de la particule fluide est constant, les équations obtenues sont alors les suivantes :

L'équation de conservation de la masse du fluide :

L'équation de bilan de la quantité de mouvement :

Les équations du fluide avec l'expansion modifier

On peut retrouver les équations de Friedmann à partir des équations précédentes : il suffit d'utiliser l'équation , et de faire quelques manipulations algébriques. Dit autrement, on ne prend en compte que l'effet de l'expansion et on néglige l'effet des vitesses locales. Les équations se simplifient alors.

Dérivation de l'équation du fluide de Friedmann modifier

Par exemple, prenons l'équation de conservation de la masse :

Injectons la loi de Hubble :

Appliquons maintenant la définition de la divergence : .

On a alors l'équation :

Cette équation ressemble à l'équation du fluide de Friedmann, bien qu'il manque la pression.

Dérivation de la seconde équation de Friedmann modifier

La seconde équation de Friedmann se dérive en prenant la divergence de l'équation d'Euler et en utilisant l'équation dite de Poisson pour calculer le terme gravitationnel.

Mais avant de faire cela, nous devons simplifier l'équation d'Euler. Il nous faut notamment déterminer la valeur du terme . Il se trouve que la première équation de conservation nous a permis de déduire que :

En injectant dans l'équation d'Euler, on trouve :

Il ne reste plus qu'à prendre la divergence. De laborieux calculs permettent alors de retrouver la seconde équation de Friedmann.

Les équations de conservation en coordonnées comobiles modifier

Pour décrire l'évolution du fluide dans un univers en expansion, il faut prendre en compte le fait que le volume de la particule fluide augmente avec le facteur d'échelle. Ce que les équations précédentes ne faisaient pas, supposant un volume constant. Les équations précédentes ne fonctionnent pas telles quelles. Par exemple, prenons l'équation de conservation de la masse et supposons qu'il n'y ait pas de flux entrant ou sortant de matière. On devrait avoir . Cependant, l'expansion va faire gonfler la particule fluide, ce qui diluera son contenu : sa masse restant constante, sa densité va diminuer avec le temps, ce qui contredit l'équation précédente ! Les équations précédentes ne donnent pas une version correcte de la conservation de la masse, de la quantité de mouvement et de l'énergie.

L'utilisation des coordonnées comobiles modifier

On peut cependant remarquer que l'on peut utiliser des distances, vitesses et volumes comobiles pour résoudre ce problème. Le volume comobile de la particule fluide ne change pas avec l'expansion, pas plus que sa surface comobile ou sa vitesse comobile. Intuitivement, la densité comobile (masse / volume comobile) se comporte comme la densité dans les équations précédentes. Même chose pour les autres grandeurs comobiles. Cependant, les opérateurs divergence et gradient doivent tenir compte du fait qu'on travaille en coordonnées comobiles, ce qui demande de les reformuler. Nous allons donc devoir reformuler les équations de manière à utiliser des coordonnées comobiles . La reformulation du gradient est assez simple :

La divergence n'est que la somme des termes du gradient, ce qui donne :

Enfin, la vitesse à utiliser doit être la vitesse locale, à savoir :

Sa dérivée est la suivante :

La reformulation des équations précédentes modifier

Maintenant, reformulons les équations précédentes avec les opérateurs définis dans la section précédente.

Prenons la première équation, celle de la conservation de la masse du fluide.

La seule modification à faire est de remplacer la divergence usuelle par celle en coordonnées comobiles, ce qui demande juste de diviser par le facteur d'échelle a(t). On a alors :

Modifions ensuite la seconde équation, celle de la conservation de l'impulsion, que voici :

Premièrement, on doit remplacer la dérivée, en utilisant l'équation . On a alors :

Ensuite, on divise par le facteur d'échelle a(t), ce qui permet de remplacer la vitesse usuelle par la vitesse comobile, sans compter que cela fait de même pour la divergence. On a alors :

Pour résumer, on obtient les équations suivantes :


Les perturbations cosmologiques

L'origine des galaxies et autres super-structures de plusieurs milliards d'étoiles est encore assez mal connue. Cependant, on sait que divers processus sont à leur origine. La gravité a naturellement joué un rôle prédominant, en forçant les étoiles à se rapprocher en amas de plusieurs millions/milliards d'étoiles. Mais la nature exacte et le déroulement de ce processus sont encore assez mal connus : les étoiles se sont-elles formées en premier avant de se regrouper en galaxies ou bien les galaxies se sont-elles formées avant de se fragmenter en amas puis en étoiles ? Personne ne le sait à l'heure actuelle. La seule chose qui est certaine est que ces structures se sont formées à partir de zones de surdensités, qui ont grossit de plus en plus sous l'effet de la gravité au point de donner des structures comme des galaxies. Il existe une théorie qui explique les grandes lignes de cette évolution des surdensités : c'est la théorie des perturbations cosmologiques. Nous allons l'aborder dans ce chapitre.

Les précurseurs des galaxies, amas de galaxies et autres "structures cosmologiques" sont des zones où la densité est plus grande qu'aux alentours. La densité n'était en effet pas uniforme et on observait des zones où la densité était plus faible ou plus grande que la moyenne. Ces petites variations de densité ont reçu le nom d'inhomogénéités. La densité dans une surdensité est modélisée comme une petite perturbation par rapport à la densité moyenne , ce qui donne :

Gaussian 2d

Les inhomogénéités sont supposées provenir d'agitations aléatoires de la matière et du rayonnement dans l'univers primordial. Les scientifiques ont de bonnes raisons de penser que ces fluctuations aléatoires étaient des fluctuations dites gaussiennes : la probabilité d'observer la perturbation en un endroit est définie par une fonction gaussienne, la fameuse courbe en cloche. Mais peu importe la manière dont elles sont générées, cela ne change rien à leur évolution future. Ces fluctuations de densité vont ensuite évoluer : certaines vont grossir, d'autres vont diminuer, etc. L'évolution des perturbations est gouvernée par divers processus que nous aborderons dans les chapitres suivants. Quoiqu'il en soit, la distribution initiale sera modifiée au cours du temps. On peut résumer cela en disant qu'on peut obtenir la distribution à un instant à partir de la distribution à un instant t en multipliant cette dernière par une fonction de transfert. Celle-ci va atténuer ou accentuer certaines perturbations initiales, donnant une nouvelle distribution statistique des perturbations.

Les équations du "fluide primordial" avec les inhomogénéités modifier

Pour rendre compte de l'évolution des inhomogénéités, on réutilise les équations en coordonnées comobiles du chapitre précédent. Celles-ci sont, pour rappel :

On doit injecter dans ces équations.

L'équation de conservation de la masse perturbée modifier

Prenons l'équation :

Injectons l’équation  :

Le terme étant constant, on peut le sortir de la dérivée du terme de gauche.

Divisons par  :

Simplifions la dérivée :

Les équations linéarisées modifier

Les équations obtenues sont les suivantes :

Ces équations ne peuvent cependant pas être résolues à l'heure actuelle, car ce ne sont pas des équations linéaires. On peut cependant supposer que la perturbation est petite : . Si l'on prend cette approximation, on a immédiatement : . Les équations deviennent alors :

De plus, cette approximation fait que les termes non-linéaires peuvent être négligés et supprimés des équations. Dit autrement, on ne conserve que les termes linéaires, en supprimant tout carré ou terme de puissance > 1. La conséquence est que le terme dans l'équation d'Euler disparaît. Les équations après linéarisation sont les suivantes :

La divergence de l'équation d'Euler modifier

Dans la suite des calculs, nous aurons à faire de nombreuses linéarisations, ce qui fait que les calculs se simplifieront progressivement. La première simplification, consiste à prendre la divergence de l'équation d'Euler. L'utilité de cette reformulation deviendra évidente sous peu. Prenons donc l'équation d'Euler et prenons sa divergence des deux côtés de l'équation :

La divergence d'une somme est la somme des divergences, ce qui simplifie le terme de gauche :

On utilise ensuite la formule , ce qui simplifie le terme de droite, ainsi que le terme : .

La divergence d'un gradient est identique au laplacien, ce qui donne :

On applique alors l'identité suivante :

Les relations entre densité et autres paramètres de l'équation d'Euler modifier

Les inhomogénéités vont naturellement influencer la pression et le potentiel gravitationnel, ainsi que la vitesse du fluide. Les zones plus denses que leur voisinage auront naturellement une gravité plus grande que leur entourage. De même, leur pression sera supérieure vu que la gravité va compresser la matière dans la surdensité, augmentant donc sa pression. Enfin, il en est de même pour la vitesse du fluide. Sous l'effet de la pression, la matière va tendre à fuir la surdensité où elle est comprimée : la matière va donc avoir une vitesse sortante supérieure à l'environnement. Mais cette vitesse est contrariée par la gravité, qui tend à faire rentrer la matière et donc à lui imposer une vitesse entrante non-nulle. Mais les équations précédentes ne permettent pas de rendre compte de ce phénomène : les termes de pression et de potentiel gravitationnel ne sont pas exprimés en fonction de la densité. On doit donc trouver des relations entre densité, pression et potentiel gravitationnel. Avec ces relations, on pourra reformuler les équations du fluide avec seulement la densité.

La relation entre densité et vitesse locale modifier

Le premier paramètre de l'équation d'Euler est la divergence de la vitesse locale . Celle-ci est en effet le seul constituant de son premier terme, avec le facteur de Hubble et le facteur d'échelle. On peut relier cette divergence à la densité en utilisant l'équation de conservation de la masse. Celle-ci donne en effet une relation entre la divergence de la vitesse et la densité. Ce n'est pas une relation directe avec la vitesse locale, mais c'est déjà un bon début.

On réorganise :

On multiplie par  :

On injecte alors cette équation dans les équations linéarisées :

Le calcul des dérivées donne :

On divise alors par pour simplifier :

On simplifie :

L'équation de Poisson modifier

Maintenant, nous allons traiter le terme lié à la gravité, à savoir le gradient du potentiel gravitationnel. Il est relié à la densité par la fameuse équation de Poisson, une forme alternative de la loi de la gravité de Newton que vous avez peut-être déjà étudié dans vos cours de physique. La voici :

En injectant la relation , on a:

Si on linéarise l'équation, seul le terme proportionnel à est conservé, ce qui donne :

On peut alors injecter cette équation dans l'équation finale de la section précédente, ce qui donne :

La relation entre pression et densité modifier

Il ne nous reste plus qu'à calculer le terme . Il se trouve que, par définition du laplacien , celle-ci est simplement la dérivée seconde de la pression. Avant de calculer cette dérivée seconde, nous allons calculer la dérivée première, avant d'en déduire la dérivée seconde. De manière générale, la pression dépend de la densité et uniquement de celle-ci dans un gaz parfait et on a donc :

Fait intéressant, le terme est, par définition, proportionnel au carré de la vitesse du son dans le fluide considéré. En effet, les formules nous disent que pour un gaz parfait, on a la relation suivante :

, avec le coefficient d'expansion adiabatique.

Dans ce qui suit, on négligera le coefficient d'expansion adiabatique , en supposant qu'il vaut 1. En notant le terme et la vitesse du son, on a :

En dérivant une seconde fois, on trouve la dérivée seconde, qui n'est autre que le terme  :

En injectant l'équation précédente dans l'équation d'Euler, on a :

L'équation finale modifier

En simplifiant l'équation précédente, on obtient l'équation finale :

On peut voir que cette équation fait intervenir un terme qui décrit à quelle vitesse l'inhomogénéité grandit. Les autres termes décrivent la manière dont la croissance de la perturbation va ensuite retentir sur sa pression et son champ de gravité.

  • Le terme est un terme de pression hydrostatique. On voit que plus la perturbation est importante, plus sa pression tendra à repousser les flux entrants de masse, et dont à contrecarrer la croissance de la perturbation.
  • Le terme : est un terme gravitationnel, qui dit que plus la perturbation est grande, plus elle sera massive et attirera de nouvelle matière. La perturbation va donc, sous l'influence de ce facteur seul, tendre à grossir de plus en plus, à une "accélération" constante.
  • Enfin, le terme : décrit l'effet de dilution de la perturbation suite à l'expansion. L'expansion tend à diluer la perturbation dans un espace plus grand et donc à réduire celle-ci de plus en plus vite. Ce phénomène est appelé l'entraînement de Hubble. Il tend à lutter contre la croissance de la perturbation, du fait de son signe négatif.

Si les termes de pression et l’entraînement de Hubble étaient nuls, toute perturbation ne ferait que grossir de plus en plus vite. La moindre surdensité aurait une gravité supérieure à son entourage, ce qui attire de la matière et fait grossir la perturbation progressivement. Ce faisant, elle aurait une masse encore plus grande, ce qui augmenterait sa gravité, et ainsi de suite. La croissance d'une telle inhomogénéité serait exponentielle. Mais la pression et l’entraînement de Hubble limitent ce phénomène, dans une certaine mesure. Pour décrire plus en détail l'évolution des inhomogénéités, il faut étudier les solutions de l'équation précédente, qui disent ce qu'il advient des inhomogénéités. On pourra ainsi savoir si celle-ci vont se contracter ( va augmenter), ou osciller. Tout le challenge tient dans le fait que nous avons étudié des perturbations ponctuelles, là où les grumeaux de matière cosmologiques ont une certaine taille, une extension spatiale. Il nous faut donc étudier ce qui se passe quand les grumeaux de matière/rayonnement ont une taille non-nulle.



L'évolution adiabatique et isotherme des perturbations

Dans les chapitres précédents, nous avons utilisé des équations où la densité était un terme unique. Celle-ci n'est autre que la densité d'énergie, qui est elle-même la somme de la densité de matière, du rayonnement et de la courbure. La courbure moyenne est considérée comme homogène, identique en tout point de l'univers. Les modifications locales de courbures sont causées par l'accumulation de matière et/ou de rayonnement, pas par une fluctuation locale de densité de courbure. Toute surdensité peut ainsi être une surdensité de matière et/ou de rayonnement.

Reste à expliquer comment les surdensités de matière sont liées avec les surdensités de rayonnement. Il existe plusieurs possibilités : soit les surdensités de matières sont liées à des surdensités de rayonnement, soit les surdensités de matières sont compensées par des sous-densités de rayonnement, soit les surdensités de matière et de rayonnement sont indépendantes. Ces trois possibilités sont étudiées dans des modèles théoriques différents, notamment dans l'étude de l'inflation ou de l'univers primordial. Dans ce qui va suivre, nous allons nous intéresser à deux types de perturbations, résumées dans le tableau ci-dessous.

Les types de perturbations
Nom Description Équations
Perturbations adiabatiques Surdensité à la fois de matière et de rayonnement.

Le rapport entre densité de rayonnement et de matière reste le même dans la surdensité et en-dehors. De même, l'entropie reste la même dans et en-dehors de la perturbation, ce qui fait qu'on donne parfois le nom d'isentropique donné à ce type de perturbation.

Perturbations isothermes (d'isocourbure) Les surdensités de matière sont compensées par une sous-densité de lumière.

Du fait de cette compensation, la densité totale reste constante partout. En conséquence, la gravité/courbure reste la même dans la perturbation et dans son environnement, d'où le nom d'isocourbure donné à ce type de perturbation.

Dans ce chapitre, nous allons voir que les perturbations adiabatiques rendent bien compte de l'état de la matière avant le découplage des photons, alors que le découplage a rendu les perturbations de matière et de rayonnement indépendantes.

Le rapport entre densités de matière et de rayonnement modifier

Dans la suite du cours, nous allons avoir besoin du rapport entre densité de matière et de rayonnement avant et après le découplage pour rendre compte de ces perturbations. Dans les quatre premiers chapitres du cours, nous avons établit que :

et

Les équations précédentes peuvent s'écrire comme suit :

et

Il est intéressant d'utiliser le rapport suivant, pour simplifier certaines équations qui vont suivre :

La vitesse du son dans l'univers primordial modifier

Pour calculer la vitesse du son, partons de sa définition :

Sachant que l'univers primordial contient du rayonnement et de la matière (l'énergie noire n'a pas encore un effet notable), on a :

La matière n'ayant pas de pression, la formule se simplifie en :

On peut ensuite utiliser la formule pour le rayonnement, ce qui donne :

On sait que , . En injectant dans l'équation précédente, on a :

Divisons par en haut et en bas de la fraction :

On utilise alors le rapport R :

L'horizon sonore modifier

Chaque zone de surdensité, supposée sphérique, a créé des ondes sonores dont le front d'onde a une forme sphérique. La taille maximale que peut prendre cette sphère dépend du temps écoulé depuis le big-bang et la recombinaison. La sphère de taille maximale correspond à une onde sonore créée lors du big-bang et figée lors de la recombinaison : c'est pour cela qu'on l'appelle l'horizon sonore. Le rayon de cette sphère se calcule en faisant l'intégrale de la vitesse du son dans le plasma primordial sur l'âge de l'univers au moment de la recombinaison.

La vitesse du son dans un gaz de photons est égale à . Elle augmente quelque peu du fait de la présence des baryons dans la réalité.

Il se trouve que l'horizon sonore correspond à la longueur de l'onde sonore de plus grande amplitude possible : la surdensité a eu le temps maximal pour se compresser.

Prenons le cas où deux régions de l'espace sont à une distance inférieure à l'horizon sonore. Dans ce cas, une onde sonore a pût se transmettre de la première région à la seconde, avant la recombinaison. Dans le cas inverse, où les deux régions sont plus éloignées que l'horizon sonore, c'est l'inverse. On est certain que les deux régions sont trop éloignées pour qu'une onde sonore ait pu passer d'une région à l'autre.

Avant le découplage : les perturbations adiabatiques modifier

Pour décrire le comportement de la matière et du rayonnement avant le découplage, il faut étudier les interactions entre matière ionisée et rayonnement. Avant le découplage des photons, les interactions entre photons et électrons libres étaient fréquentes, suffisamment pour redistribuer l'énergie des photons aux électrons (et réciproquement). Voyons dans le détail ces interactions, avant de voir leurs conséquences sur la densité de l'univers.

Le processus de diffusion Compton sans expansion modifier

Illustration du processus de diffusion Compton.

Avant le découplage, le processus d'interaction principal était la diffusion Compton, ainsi que le processus inverse, la diffusion Compton inverse. Celle-ci est une collision entre un photon et un électron libre, où les deux particules échangent de la quantité de mouvement et de l'énergie cinétique.

Les physiciens rendent compte de la fréquence de ces collisions par deux paramètres : le temps moyen entre deux collisions et la distance moyenne entre deux collisions. Le temps moyen entre deux collisions est appelé le temps de libre parcours moyen, tandis que la distance moyenne entre deux collisions est appelée le libre parcours moyen. Évidemment, les deux sont liés entre eux par la vitesse des particules du gaz : la distance de libre parcours moyen est ainsi égale au produit de la vitesse moyenne par le temps de libre parcours moyen. Ainsi, plus un gaz est froid, plus le temps et la distance entre deux collisions sont faibles. Il en est de même pour les gaz denses : ceux-ci ont un libre parcours moyen assez faible, vu que les particules sont plus proches.

Il est possible de calculer le libre parcours moyen en utilisant une relation très simple entre celui-ci et deux autres variables : le nombre de particules par unités de volume et la probabilité moyenne d'une collision. Le nombre de particules par unité de volume remplace en quelque sorte la densité : plus il y a de particules par unité de volume, plus celles-ci ont de chances d'entrer en collision. Quant à la probabilité d'une collision, elle se calcule à partir des équations de la diffusion Compton. Nous ne ferons pas ces calculs ici, d'autres articles sur la diffusion Compton le feront certainement mieux que moi. Elle se calcule avec la formule suivante, avec la constante de structure fine et la masse d'un électron.

Le libre parcours moyen est égal, par définition, à l'équation suivante :

, avec le libre parcours moyen, le nombre d'électrons par unité de volume et la probabilité d'interaction Compton entre deux particules.

Le temps de libre parcours moyen se calcule en divisant cette distance par la vitesse d'un photon, ce qui donne :

Avec expansion, les équations de la diffusion Compton change un petit peu. Avec l'expansion, le libre parcours moyen va naturellement augmenter, à cause de l'évolution du facteur d'échelle. On peut remarquer que l'expansion n'influence pas la section efficace d'interaction Compton : son influence se fait sur la densité des électrons, qui diminue d'un facteur avec l'expansion. Suite à l'expansion, il arrive un moment où le libre parcours moyen devient plus grand que l'horizon cosmologique, que l'univers observable. Lorsque cela arrive, la diffusion Compton entre un photon du CMB et un électron libre devient "impossible", ou tout du moins fortement improbable. Les photons cessent d'interagir complètement avec les électrons libres, qui sont alors libres de former des atomes avec les noyaux. C'est ainsi qu'à lieu le découplage. Et cela a des conséquences sur l'évolution des perturbations.

Avant le découplage, les perturbations étaient adiabatiques, isentropiques. La diffusion Compton redistribuait la quantité de mouvement entre matière et rayonnement, et le couplage fort liait les densités de photons et de baryons entre elles. Du fait des interactions entre rayonnement et matière, le rayonnement restait prisonnier des surdensités de matière. Là où il y avait de la matière, le rayonnement devait suivre. Ne pouvant s'échapper sous forme de rayonnement, la chaleur restait coincée dans les surdensités. La chaleur ne pouvant s'échapper, l'entropie reste constante dans ces perturbations adiabatiques. Ce n'est qu'après le découplage des photons que la matière et le rayonnement ont commencé à évoluer indépendamment l'une de l'autre. Le rayonnement a pu alors quitter les zones de surdensité, vu qu'il n'interagissaient plus avec elles, et échanger de la chaleur avec l'environnement. Le rayonnement s'est progressivement dilué dans l'espace, alors que les surdensités de matière sont restées, donnant des perturbations de matière pure. Les perturbations ne sont alors ni isothermes, ni adiabatiques.

L'approximation du couplage fort modifier

La diffusion Compton ne faisait pas que redistribuer la quantité de mouvement : elle couplait aussi fortement le nombre de photons et de baryons. Si on prenait un volume d'espace fini, le nombre de photons et de baryons ne différait pas sensiblement l'un de l'autre. Il y avait bien quelques variations ci et là, mais on peut parfaitement supposer qu'en moyenne, le nombre de particules par unité de volume était le même pour les photons et les baryons. Le résultat est ce qu'on appelle l'approximation du couplage fort : le nombre de baryons et de photons dans un volume d'espace fini est approximativement le même. La densité de baryon et de rayonnement étaient alors fortement liés par une relation que nous allons démontrer.

L'approximation du couplage fort nous dit que, avec et le nombre de photons et de baryons par unité de volume, on a :

On sait de plus que et que , ce qui donne :

Ces informations nous sont utiles pour calculer la vitesse du son. Pour rappel, la vitesse du son se calcule avec la formule suivante :

En combinant les deux équations précédentes, on trouve :

Après le découplage : les perturbations isothermes et leur évolution modifier

Après le découplage, la matière et le rayonnement ne sont plus couplées. Les densités de matière et de rayonnement deviennent indépendantes l'une de l'autre et évoluent chacun dans leur coin. Elles suivent alors les relations vues dans les premiers chapitres :

Une vitesse du son en chute libre modifier

Rappelons que la vitesse du son est la suivante :

, avec .

Calculons le rapport R :

En combinant les deux équations précédentes, on trouve :

L'équation précédente dit que la vitesse du son diminue avec l'expansion, du moins après le découplage.

On peut s'en rendre compte avec l'approximation , qui transforme l'équation précédente en :

Si on attend suffisamment longtemps après le découplage, le facteur d'échelle sera devenu suffisamment grand pour que l'on ait : . On obtient alors :

La réduction de la vitesse du son n'est pas sans conséquences sur la formation des perturbations et leur évolution. En effet, les ondes sonores sont en effet des ondes de pression (compression/décompression). Elles naissent quand la pression varie en un point d'un objet et transmettent cette variation aux alentours. La vitesse du son est aussi la vitesse de transmission des forces de pression. Or, la pression est ce qui permet de lutter contre la gravité dans les surdensités. Si la vitesse de transmission de la pression est trop faible, les forces de gravité peuvent l'emporter sur la pression dans certaines circonstances. Une surdensité peut ainsi s'effondrer si les ondes de pression sont trop lentes et ne permettent pas de lutter assez rapidement contre leur contraction gravitaire.



L'évolution des perturbations avant le découplage

Sans forcément utiliser le spectre de puissance du chapitre précédent, on peut utiliser directement la décomposition en séries de Fourier du champ de densité. Celle-ci décrit le champ de densité comme la somme/intégrale d'une infinité d'ondes périodiques simples. Si on note chaque onde élémentaire , le théorème de Fourier nous donne l'équation suivante :

, avec le vecteur d'onde, un vecteur de norme .

On peut alors combiner cette équation avec l'équation suivante, dérivée il y a quelques chapitres :

On peut combiner les deux équations. Le résultat, très difficile à obtenir, est pourtant surprenamment simple :

L'équation différentielle précédente ressemble formellement à une équation d'oscillateur harmonique amorti. On rappelle qu'un oscillateur harmonique est un système soumis à une force de rappel proportionnelle à la distance : . La force de rappel est ici la somme de la pression et de la gravité, qui poussent le système à osciller de manière périodique. Pour obtenir un système amorti, il faut ajouter un terme de friction, qui tend à réduire l'amplitude des oscillations au fil du temps. Ici, la force de friction est le poids des baryons, qui force le système à se compresser. La forme générale d'une équation d’oscillateur harmonique est la suivante :

On voit que les deux équations précédentes sont identiques si on pose : , et .

Les solutions de cette équation sont de la forme :

, avec et les solutions de l'équation du second degré .

Le comportement de la solution dépend selon la valeur des racines. On distingue trois cas, qui sont décrits dans le tableau ci-dessous.

Discriminant Interprétation dans l'équation du fluide cosmologique
Cas périodique amorti Le discriminant de l'équation du second degré est positif. La pression et la gravité se compensent l'une l'autre.
Cas apériodique Le discriminant de l'équation du second degré est négatif. La pression est plus forte que la gravité.
Cas critique Le discriminant de l'équation du second degré est nul. La gravité est plus forte que la pression.

Nous allons voir chaque cas l'un après l'autre.

Le cas stable : la distance de Jeans modifier

Lorsque la pression est égale à la gravité, on a :

Cette équation dit simplement que la perturbation est stable : elle s'accentue au même rythme que l'expansion la dilue (terme d’entraînement de Hubble). En clair : la perturbation est stable et ne change pas. Mais rappelons-le, cela n'arrive que pour une longueur d'onde précise, appelée rayon de Jeans. Cette longueur d'onde donne la taille de la surdensité pour laquelle la perturbation se stabilise. Cela fonctionne pour les perturbations sphériques, même si elles ne sont pas périodiques (la magie des fonctions de Fourier...).

Le rayon de Jeans modifier

Le rayon de Jeans se calcule comme suit :

On utilise alors l'équation  :

On développe le carré :

Simplifions par  :

Prenons l'inverse de l'équation précédente :

Prenons la racine carrée :

Simplifions :


Démonstration

On peut donner une autre démonstration du rayon de Jeans. Pour cela, imaginons une surdensité sphérique et homogène de rayon et de masse . Et prenons le cas où il n'y a aucune force de pression. Toute surdensité va alors s'effondrer sur elle-même en mettant un certain temps, que l'on va appeler temps d'effondrement. Les équations de la physique classique disent que ce temps est proportionnel à :

Maintenant, ajoutons les forces de pression, dont on suppose qu'elles contrecarrent totalement la gravité. La transmission des forces de pression se fait à la même vitesse que le son. Une surdensité ne peut pas s'effondrer si le temps d'effondrement est égal au temps que met la pression pour parcourir la surdensité. Dit autrement, une surdensité s'effondre sur elle-même si :

Pour une surdensité de rayon et une vitesse pour la pression, le temps de parcours vaut :

En combinant les deux équations précédentes, on a:

On peut alors calculer le rayon R, appelé rayon de Jeans. Celui-ci est simplement le rayon au-delà duquel la surdensité ne peut pas s'effondrer sur elle-même.

L'évolution du rayon de Jeans avant et après le découplage modifier

Les calculs précédents nous disent que le rayon de Jeans dépend de la vitesse du son. Or, rappelez-vous que la vitesse du son dépend du rapport . La relation exacte est la suivante :

Avant le découplage, le rapport R vaut 1, et ce durant toute la période avant le découplage. La vitesse du son est constante et le rayon de Jeans fait de même.

Mais après le découplage, le rapport R diminue progressivement et la vitesse du son évolue comme ceci :

Cela impacte le rayon de Jeans, qui évolue avec le facteur d'échelle et donc avec l'expansion. Celui-ci se calcule alors comme suit :

Ou encore :

Ce qui se simplifie en :

On voit que le rayon de Jeans diminue avec l'expansion.

Le cas où G > P : la croissance/décroissance infinie modifier

Prenons maintenant une perturbation dont la taille est supérieure à la longueur de Jeans, où la gravité l'emporte sur la pression. Il n'y a qu'une seule solution à cette équation différentielle, qui est de la forme :

On voit que cette solution est la somme de deux solutions particulières :

  • la première est une fonction , strictement croissante ;
  • la seconde est une fonction , strictement décroissante.

Dit autrement, toute structure d'une taille supérieure au rayon de Jeans va soit s'effondrer sur elle-même sous l'effet de sa propre gravité, soit gonfler à cause de l'expansion.

Le cas où P < G : les oscillations acoustiques de baryons modifier

Une perturbation dont la taille est inférieure au rayon de Jeans correspond au cas amorti. Dans celui-ci, la solution de l'équation est une fonction d'onde de la forme :

Illustration de l’oscillation de densité - équation d'un oscillateur harmonique amorti.

L'équation n'a que des solutions ondulatoires dont l'amplitude diminue au cours du temps. Dit autrement, la perturbation pulse, elle grossit avant de dégonfler et ainsi de suite, indéfiniment. On peut interpréter cette solution physiquement comme suit. Sous l'effet de la gravité, de petites zones de surdensités vont se former et vont se contracter sous l'effet de la gravité. Cette contraction va réchauffer la matière dans la zone de surdensité, ce qui en augmentera la pression de radiation (on néglige totalement la pression de la matière baryonique). Cette pression luttera contre la force de gravité, forçant la surdensité à gonfler et à se dilater. Cette dilatation fait chuter sa température, et sa pression fait alors de même. La gravité reprend alors ses droits, en devenant plus forte que la pression, et un nouveau cycle commence. Ce gonflement/dégonflement périodique de la surdensité va retentir sur son environnement : elle va former des ondes de compression/décompression qui se propagent dans le plasma primordial. Ces ondes ne sont autres que des ondes sonores ! En clair : l'évolution des inhomogénéités cause des ondes sonores primordiales, aussi appelées oscillations acoustiques de baryons.

Lors du découplage, ces ondes sonores se sont figées, les photons n'agissant plus sur la matière. La pression de radiation reste la même, mais celle-ci influence relativement peu la matière, laissant la gravité faire son travail. Ne pouvant plus contrecarrer la gravité, la pression de radiation cesse de donner naissance aux oscillations acoustiques de baryons. Les zones de surdensité cessent d’osciller et s'effondrent sur elles-mêmes, donnant naissance à diverses structures : de futures galaxies pour les surdensités de grande échelle, notamment. Cependant, on trouve encore une trace des oscillations acoustiques de baryons dans le fond diffus cosmologique. On verra que ces oscillations acoustiques de baryons ont laissé de petites perturbations de température dans le fond diffus. Nous en reparlerons dans le prochain chapitre.



L'évolution des perturbations après le découplage

Dans ce chapitre, nous allons étudier ce qui se passe après le découplage. Pour cela, on repart de l'équation habituelle :

On a vu il y a deux chapitres que, après le découplage, le terme dépendant de la vitesse du son disparaît. L'équation précédente devient alors :

La raison à cela est assez simple à comprendre. Vu que la matière n'a par définition pas de pression, seul le rayonnement est à l'origine du terme de pression dans l'équation originale. Après le découplage, le rayonnement se dilue avec l'expansion et la pression diminue donc. Au bout d'un moment, la pression disparaît. Par contre, même si la matière n'a pas de pression, elle a une masse qui influence l'évolution des perturbations, de par sa gravité. Le terme de gravité reste donc, alors que le terme de pression s'annule.

Dans ce chapitre, nous allons voir quelle est la dynamique de cette équation, ce qui nous enseignera comment les perturbations ont été affectées par l’expansion.

Une croissance/décroissance linéaire modifier

L'équation différentielle précédente n'a aucun coefficient dépendant des coordonnées spatiales, sans compter qu'il n'y a pas de dérivées par rapport à ces coordonnées. Seul le terme dépend des coordonnées spatiales. Dans ce cas, on est assuré par les lois mathématiques que la solution est égale au produit de delta par un terme indépendant des coordonnées spatiales. Par contre, ce terme dépend du temps, vu que les dérivées et les coefficients dépendent du temps. On a donc une solution de la forme :

L'équation a précisément deux solutions qui ont une croissance ou décroissance linéaire. Mais la solution décroissante peut être oubliée et l'on peut se concentrer sur la solution croissante. Celle-ci nous dit qu'après le découplage, les inhomogénéités gonflent au même rythme que l'expansion. En coordonnées comobile, leur forme est restée figée après le découplage, sans aucune modification autres que celles liées à l'expansion.

En injectant cette solution croissante, l'équation originale devient :

Le terme étant indépendant du temps, il est une constante du point de vue des dérivées. On peut donc le sortir des dérivées, ce qui donne :

On divise alors par  :

Ce qui peut d’écrire d'une manière plus succincte comme suit :

Les cas particuliers : univers dominé par la matière, le rayonnement ou l'énergie noire modifier

Formellement, la fonction correspond à la vitesse de croissance des structures. On peut démontrer qu'elle est proportionnelle au facteur d'échelle :

On peut alors comparer la vitesse de croissance entre un univers dominé par la matière, un autre dominé par le rayonnement et un dernier dominé par l'énergie noire. Le premier cas correspond à un univers après le découplage, alors que le second correspond à ce qu'on observe avant le découplage.

L'univers dominé par la matière modifier

Pour l'univers dominé par la matière, les équations vues dans le chapitre sur les modèles cosmologiques de Friedmann nous disent que :

Rappelons l'équation originale :

Injectons l'équation  :

Simplifions par  :

Puis, injectons l'équation  :

Simplifions :

En supposant que , on a :

En calculant les dérivées, on trouve :

On divise par  :

Les deux solutions possibles sont respectivement et . La première correspond à des fluctuations qui décroissent avec le temps, la seconde est plus intéressante. En injectant dans , on trouve :

L'univers dominé par le rayonnement modifier

Pour l'univers dominé par le rayonnement, les équations vues dans le chapitre sur les modèles cosmologiques de Friedmann nous disent que :

Rappelons l'équation originale :

Injectons l'équation et simplifions :

Négligeons le dernier terme :

Supposons que la solution soit une loi de puissance et injectons-la dans l'équation précédente :

On peut sortir la constante k des dérivées :

Simplifions par k :

Le calcul des dérivées donne :

Si on pose , on peut remarquer que c'est une solution, et même une double solution (elle annule les deux termes). La résolution de cette équation nous dit donc que la solution générale est de la forme :

Les perturbations croissent donc de manière logarithmique dans un univers dominé par le rayonnement. Au passage, cette situation représente bien ce qui s'est passé avant le découplage des photons.

L'effet Mezaros modifier

On vient de voir que le facteur vaut :

  • pour un univers dominé par la matière ;
  • pour un univers dominé par le rayonnement ;
  • et pour un univers dominé par l'énergie noire.

On voit que dans des univers dominés par la matière, les perturbations ne grandissent pas très vite mais le font en suivant une loi de puissance. Cette croissance est inférieure à celle observée dans un univers statique, pour lequel on aurait :  ! C'est avant le découplage, dans un univers dominé par le rayonnement, que la croissance des structures est la plus lente. Leur croissance a été plus rapide après le découplage, leur vitesse de croissance ayant quelque peu augmenté lors du passage à un univers dominé par la matière. Ce résultat porte le nom d'effet Meszaros.



Le spectre de puissance des perturbations

Dans le chapitres précédent, nous avions éliminé la pression des équations, afin de les rendre plus simples. Il faut dire que résoudre ces équations avec la pression est compliqué, mais c'est faisable en utilisant un formalisme mathématique précis. Celui-ci est le formalisme de la transformée de Fourier. Et quitte à introduire ce formalisme, autant parler plus en détail de la manière de modéliser le champ de densité, qui a un lien très fort avec le formalisme de Fourier. Faisons donc une petite parenthèse sur le sujet.

Nous sommes parti du principe que la distribution de la densité nous était inconnue. Nous supposions simplement qu'elle existe, mais sans en savoir plus sur ses propriétés. Et il est vrai que dans le cas général, on ne peut rien dire sur la distribution des perturbations. Par contre, on peut en donner quelques propriétés statistiques plus ou moins pertinentes. Ce champ de densité, peu importe sa forme exacte, a une densité moyenne, une certaine dispersion autour de cette moyenne, et ainsi de suite. Dans ce chapitre, nous allons voir diverses mesures statistiques du champ de densité et voir comment elles se marient avec les équations du chapitre précédent.

La fonction de corrélation modifier

Pour commencer, prenons un exemple assez simple, mais en apparence détaché des perturbations cosmologiques. Si je regarde un point de l'espace et un volume autour de ce point (supposé assez gros à l'échelle cosmologique). Il y a une probabilité que dans le volume , je trouve une galaxie. Si on suppose qu'il y a en moyenne galaxies par unité de volume, alors cette probabilité est de :

De même, en un point , cette probabilité est de :

Maintenant, on souhaite savoir quelle est la probabilité que j'observe une galaxie à la fois au point et au point . Intuitivement, on penserait avoir :

Mais cela ne vaut que si les positions des galaxies sont totalement indépendantes, ce qui n'est pas garanti. Dans les faits, il est possible qu'il y ait une relation dans la distribution des galaxies qui fait que si on observe une galaxie en x, alors sa présence en y est plus probable, en fonction des positions x et y. On doit donc tenir compte de telles corrélations. La formule exacte, qui en tient compte, est la suivante :

, où est la fonction de corrélation.

Par analogie, on peut utiliser le même formalisme, mais pour prédire non pas la présence d’une galaxie, mais l'intensité d'une fluctuation de densité. On a alors :

, où est la densité moyenne du champ de densité.
Il est possible de définir des fonctions de corrélation pour trois points, quatre points, voire beaucoup plus.

L'hypothèse d’homogénéité statistique modifier

Si on considère que l'univers est statistiquement homogène, alors la fonction de corrélation ne dépend que de la distance entre x et y, mais pas de la position exacte de x et y.

Ce qui peut s'écrire comme suit :

, avec r la distance entre les deux points x et y.

Sous cette hypothèse, le calcul de la densité moyenne est assez simple. Il suffit de prendre la moyenne spatiale de la densité. La corrélation moyenne entre deux points peut se calculer en prenant un grand nombre de points x et y et en calculant la corrélation pour chaque paire de points. Il suffit de faire la moyenne des corrélations obtenues, pour obtenir la corrélation moyenne. les mesures semblent montrer que la fonction de corrélation suit une loi de puissance de la forme suivante :

, avec .

La transformée de Fourier du champ de densité modifier

Outre la fonction de corrélation, on peut aussi utiliser le formalisme dit des séries de Fourier. Pour rappel, le terme est une fonction qui associe une perturbation de densité à tout endroit de l'espace et à chaque instant. On dit aussi que cette fonction décrit un champ de densité. Or, il existe un théorème qui nous dit que tout champ peut être décomposée en champs périodiques semblables à des cosinus ou sinus. Formellement, ces champs périodiques sont des ondes de forme cosinusoïdales ou sinusoïdales. Dans notre cas, la forme de ces ondes est l'équivalent en trois dimension d'un sinus/cosinus. En additionnant un certain nombre (potentiellement infini) de ces ondes de base, pondérées par un coefficient, on peut obtenir n'importe quel champ résultant. C'est ce qu'on appelle la transformée de Fourier des fonctions continues. Le champ de densité ne fait pas exception et on peut utiliser ce théorème pour décomposer le champ de densité en une somme d'ondes.

Illustration des séries de Fourier.

Le spectre de puissance des perturbations modifier

Décrire le champ de densité avec des séries de Fourier permet de définir le spectre de puissance du champ de densité. Il donne l'ensemble des amplitudes de chaque onde sinusoïdale et s'obtient avec la transformée de Fourier. Dit autrement, il donne l'amplitude de chaque perturbation en fonction de la taille (de sa longueur d'onde). Dans le cas qui nous intéresse, la longueur d'onde correspond à la taille d'une perturbation périodique. Le spectre de puissance donne donc l'intensité de la surdensité en fonction de sa taille et est donc une fonction du type :

, avec le nombre d'ondes qui est défini par .

Le spectre de puissance et la fonction de corrélation sont reliés l'un à l'autre, par la relation suivante :

.

Qui peut aussi s'écrire comme suit :

Les deux équations précédentes nous disent que le spectre de puissance et la fonction de corrélation (en fait son intégrale) sont proportionnels.

Les spectres en loi de puissance et le spectre de Harrison-Zeldovitch modifier

Dans le cas général, connaître le spectre de puissance n'est pas suffisant pour décrire complètement le champ de densité, du moins d'un point de vue statistique. Il en est de même avec la fonction de corrélation qui est elle aussi un résumé imparfait de la distribution. Cependant, il existe des distributions statistiques pour lesquelles la connaissance du spectre de puissance et/ou de la fonction de corrélation suffit à décrire totalement les propriétés statistiques du champ décrit. Ce sont les champs aléatoires gaussiens, pour lesquels la densité suit une distribution gaussienne (la fameuse courbe en cloche). Et ce sont ces gaussiennes qui sont utilisées pour modéliser le champ de densité cosmologique, faute de mieux. Le spectre de puissance de tels champs aléatoires gaussiens suit une loi de puissance de la forme :

Quand l'exposant vaut 1, l'amplitude des fluctuations ne dépend pas de l'échelle. Le spectre de puissance avec est appelé le spectre de Harrison-Zeldovitch, quand il est utilisé en cosmologie.

Précisons que si le spectre de puissance est une loi de puissance, alors la fonction de corrélation est aussi une loi de puissance de la forme :

, avec .

Un avantage de ce spectre de puissance est qu'il est invariant si les perturbations sont stables en coordonnées comobiles. Par exemple, c'est le cas après le découplage, où toutes les perturbations évoluent linéairement avec l'expansion, c'est à dire qu'elles grossissent au même rythme que l'expansion. Dit autrement, elles ont toujours la même taille en coordonnées comobiles et sont figées. Dans ce cas, on peut décrire de telles perturbations comme suit :

La fonction de corrélation associée devient alors :

Et le spectre de puissance est de :

En clair, avec l'expansion, le spectre de puissance évolue, mais cette évolution est prévisible. Il est simplement multiplié par une valeur qui dépend du carré du temps écoulé. Cela veut dire deux choses. Premièrement, le spectre de puissance est décalé avec le temps. Deuxièmement, les amplitudes sont réduites par un coefficient multiplicatif dépendant du temps. Du moins, c'est le cas pour les structures qui suivent cette loi d'évolution.

Les défauts des spectres en loi de puissance modifier

Les spectres en loi de puissance ont des propriétés intéressantes. Cependant, des arguments techniques nous font dire que le spectre des perturbations ne peut pas être en loi de puissance sur l'ensemble des longueurs d'onde possibles. Ce n'est pas compatible avec l’homogénéité du CMB, qui impose que :

si

La relation entre spectre de puissance et fonction de corrélation impose alors que :

En clair, si il y a des corrélations en loi de puissance pour un intervalle de k précis, alors il existe des anti-corrélations en-dehors de cette intervalle. Les corrélations et anti-corrélations se compensent et annulent l'intégrale précédente.



Les anisotropies du fond diffus

Carte du fond diffus cosmologique, qui illustre les anisotropies.

Le fond diffus capté à l'heure actuelle correspond aux photons émis 380 000 ans après le big-bang environ, sur une sphère centrée sur la Terre : cette sphère est appelée la surface de dernière diffusion. Si on regarde le fond diffus, on peut remarquer que celui-ci n'est pas totalement homogène, avec quelques variations de température assez aléatoires et dispersées appelées anisotropies. Les scientifiques définissent les anisotropies comme toute déviation de la température CMB par rapport à sa température moyenne. Dans une anisotropie, il existe une petite différence de température en pourcentage par rapport à la moyenne du fond diffus, qui est égale à :

L'origine des anisotropies du CMB modifier

Les anisotropies ont diverses origines, les scientifiques font ainsi la différence entre anisotropies primaires et secondaires. Les anisotropies primaires se sont formées lors du découplage, de petites variations de densité ayant laissé leur trace dans le fond diffus. Les anisotropies secondaires après le découplage, quand la matière des étoiles a interagi avec ce fond diffus (rappelez-vous que la matière est à l'état de plasma dans une étoile, plasma qui interagit fortement avec les photons). Les anisotropies secondaires ne sont pas importantes pour le cosmologiste qui étudie la théorie du big-bang. De plus, ces anisotropies secondaires sont de petite taille, peu visibles sur les résultats actuels (à l'exception des anisotropies formées par effet Sache-Wolfe). On peut donc les omettre sans perdre en généralité.

Les empreintes des perturbations de densité dans le CMB modifier

La formation des anisotropies primaires est liée aux perturbations cosmologiques, à la présence de régions de surdensité ou de sous-densité. Les perturbations de densité se traduisent par des perturbations de température, vu que la densité et la température sont reliées par l'équation suivante :

On verra à quoi correspond la variable x plus tard, quand nous aborderons l'effet Sache-Wolfe. Toujours est-il qu'on peut la négliger pour le moment, ce qui donne :


Démonstration

On a vu dans le chapitre sur l'évolution du rayonnement que pour le rayonnement :

Quelques manipulations algébriques donnent le résultat voulu.

L'effet Sache-Wolfe modifier

Les surdensités à l'origine des anisotropies sont des zones où la température est plus élevée qu'alentour, cette différence de température entraînant un surcroît de luminosité. Néanmoins, ce sont aussi des zones où la gravité est plus importante. Or, la relativité générale nous dit qu'un photon qui part d'une zone de forte gravité est décalé vers le rouge. Ce décalage gravitationnel est d'origine gravitationnelle, purement lié à la relativité générale. Cela vient du fait que la gravité, dans la relativité générale, est une déformation de la géométrie de l'espace. On peut simplifier en disant que les distances ne se calculent pas de la même manière dans un champ de gravité qu'en dehors. On peut rendre compte de ce phénomène d'une manière simplifiée (et même fausse) en disant que la gravité entraîne une modification stable du facteur d'échelle, par rapport aux environs (c'est une modification "locale"). En sortant d'un champ gravitationnel, un photon garde la même longueur d'onde comobile, ce qui se traduit par une modification de la longueur d'onde perçue (qui tient compte de la modification locale du facteur d'échelle). La modification fait que la longueur d'onde augmente en "sortant" du champ gravitationnel.

Dans le détail, la différence locale de facteur d'échelle est égale au potentiel gravitationnel .

Du fait de ce décalage vers le rouge, la surdensité semble avoir une température inférieure par rapport à celle qu'elle aurait sans ce décalage vers le rouge. C'est ce qu'on appelle l'effet Sache-Wolfe. Pour résumer, on peut décomposer la température perçue par un observateur lointain comme étant la somme du décalage vers le rouge avec la température locale (sans décalage vers le rouge). En notant et les températures perçues et locales, on a :

Dans leur article original, Sache et Wolfe ont postulé que la valeur de était proportionnelle à . Rien d'étonnant à cela : la température dépend de la densité, qui dépend elle-même du potentiel gravitationnel. Leurs calculs et approximations leur ont donné : . On peut en faire une démonstration relativement simple, basée sur quelques arguments physiques.

Par définition, on sait que :

On a vu dans le chapitre sur les modèles cosmologiques que dans un univers où se mélangent rayonnement et matière, on a :

On combine les équations précédentes, on a :

On peut alors se demander à quoi correspond le terme . Il faut savoir que dans un champ de gravitation, le temps s'écoule plus lentement : c'est le phénomène de dilatation des durées. La relativité générale nous dit que ce ralentissement du temps est approximativement proportionnel au potentiel gravitationnel. Autrement dit : . En combinant les deux équations précédentes, on a :

En supposant que la matière domine au fond diffus, on trouve :

En injectant dans l'équation précédente, on trouve :

Ce qui se simplifie en :

Le spectre de puissance du fond diffus modifier

La taille des anisotropies a beaucoup de choses à nous dire concernant la courbure de l'univers, le rapport entre masse visible et énergie noire, et ainsi de suite. Pour comprendre pourquoi, il faut faire appel à ce que l'on appelle le spectre de puissance du fond diffus. Pour rappel, la surface de dernière diffusion est une sphère. La localisation d'un point sur cette sphère demande d'utiliser un système de cordonnées avec : un méridien et un équateur : on peut alors déterminer une latitude et une longitude pour chaque point, celles-ci permettant de localiser le point sur la surface de la sphère. Tout point du fond diffus est donc identifié par une longitude et une latitude . En chaque point, l'intensité de la lumière qui compose le fond diffus varie, à cause de différences de température liées aux anisotropies. Les scientifiques cartographient le fond diffus en notant pour chaque point, la différence de température en pourcentage par rapport à la moyenne du fond diffus :

Illustration de quelques harmoniques sphériques.

Les différences de température entre chaque point sont causées par des différences de densité et de pression, elles-mêmes décrites par un champ de densité. Il y a quelques chapitres, nous avions vu que le champ de perturbation peut se décrire comme la somme de plusieurs champs de densité d'amplitudes différentes. En effet, un théorème très utilisé en physique, le théorème de Fourier, nous dit que toute fonction peut être décomposée en fonctions plus simples, de formes sinusoïdales. En additionnant un certain nombre (potentiellement infini) de ces ondes de base, pondérées par un coefficient, on peut obtenir n'importe quelle fonction résultante. Dans le cas du CMB, qui est de forme sphérique, ces ondes simples sont appelées des harmoniques sphériques.

Il est possible d'analyser l'amplitude de chaque harmonique sphérique du CMB en fonction de leur fréquence. Il suffit de reporter ces amplitudes sur un graphique, dont l'abscisse donne la fréquence de chaque harmonique, et l'ordonnée son amplitude. Dit autrement, ce graphique donne l'amplitude de chaque anisotropie en fonction de sa taille : telle anisotropie d'une taille de quelques kilomètres aura une amplitude de tant, telle anisotropie plus grosse aura une amplitude de température plus élevée, etc. Le tout est ce qu'on appelle le spectre de puissance du CMB. Le résultat donne la figure suivante, bien connue des cosmologistes :

Spectre de puissance du CMB.

Le spectre de puissance du CMB peut se déduire des développements des chapitres précédents. Pour résumer rapidement les chapitres précédents, les perturbations ont évolué de deux manières différentes en fonction de leur taille. Les grosses perturbations sont restées figées si on les regarde en coordonnées comobiles, ce qui signifie qu'elles ont grossi en même temps que l'expansion. Dans le chapitre où nous avions introduit le concept de spectre de puissance, nous avions vu que pour de telles perturbations, le spectre de puissance était décalé par l'expansion et que les amplitudes étaient multipliées par un coefficient dépendant du temps, mais le spectre de puissance n'était pas déformé. À l'inverse, pour les petites longueur d'onde, les oscillations acoustique modifient le spectre de puissance et forment des pics assez visibles.

Le plateau de Sachs-wolfe modifier

Comme dit plus haut, les grosses perturbations n'évoluent plus et ne sont affectées que par l'expansion, elles restent figées en coordonnées comobiles. Leur spectre de puissance est assez stable et ne montre pas de pics particuliers, vu qu'il n'y a pas de taille prépondérante pour de telles perturbations. Il n'y a pas de longueur caractéristique qui se démarque. Le spectre de puissance, dans cette portion, suit une loi de puissance.

La taille minimale de ces perturbations se dérive à partir du calcul de l'horizon sonore, la distance à partir de laquelle les oscillations acoustiques n'ont pas d'effet. Pour une distance angulaire inférieure à 100 (l < 100), on est certain que les deux points étudiés sont séparés par une distance supérieure à l’horizon sonore. Précisons que cette distance angulaire tient compte de l'effet Sachs-Wolfe. On est donc certain que les oscillations acoustiques de baryons n'ont pas pu agir dans cette région du spectre, pas plus que les autres phénomènes. Seules les fluctuations primordiales ont laissé leur trace dans cette région du spectre. Le spectre initial de ces fluctuations a été perturbé par l'effet Sachs-Wolfe, ainsi que par des parasites liés aux anisotropies secondaires. Sur le spectre, on voit que cette région est relativement plate et se trouve juste avant le premier pic. Ce qui lui a valu le nom de plateau de Sachs-Wolfe.

Les pics du CMB modifier

Lorsqu'on observe le spectre de puissance du CMB, on observe la présence de plusieurs pics à certains endroits. Les mathématiques nous disent que ces pics sont localisés à des fréquences bien précises, multiples d'une fréquence dite fondamentale. Les pics sont situés respectivement à la fréquence fondamentale pour le premier pic et aux harmoniques pour les suivants. On observe plusieurs de ces pics sur le CMB, dont l'amplitude nous renseigne sur certains paramètres, comme la courbure ou le rapport entre masse visible et matière noire. Les deux premiers sont déterminés par la taille des oscillations acoustiques de baryons, alors que ce n'est pas le cas des suivants.

Le premier pic donne des indications sur la fréquence fondamentale des ondes sonores, qui est liée à la courbure de l'univers K. Précisément, la fréquence fondamentale correspond à une surdensité dont l'extension était maximale au moment du big-bang, et qui s'est compressée avant le découplage. Au moment du découplage, celle-ci a atteint sa compression maximale et s’apprêtait à se détendre. La fréquence fondamentale ayant l'amplitude maximale, sa longueur d'onde est donc égale à l'horizon sonore. Depuis, cet horizon sonore a gonflé du fait de l’expansion de l'univers, ce qui se traduit par une diminution de sa taille. La fréquence fondamentale actuelle est tout simplement le rapport entre distance de la surface de dernière diffusion et rayon de l'horizon sonore. L'analyse du premier pic donne une courbure nulle.

Le second pic correspond au double de la fréquence fondamentale, qui correspond à une surdensité qui a eu le temps de faire une période complète entre le big-bang et le découplage. Il donne des indications sur la quantité de matière baryonique, à savoir formée de protons et de neutrons (en réalité, de baryons, des particules composites formées de quarks et gluons). Cela vient du fait que la valeur moyenne des ondes sonores vaut, comme on l'a vu plus haut : . Le rapport entre densité de matière baryonique et de rayonnement modifie la valeur de R, et donc la valeur moyenne. Pour comprendre le mécanisme exact, rappelons que la fréquence des oscillations dépend de la différence entre pression et gravité dans les surdensités. La pression dépend essentiellement de la pression de radiation, mais pas du tout de la pression de la matière baryonique. La gravité par contre, dépend autant de la matière baryonique que du poids du rayonnement. Sans matière baryonique, les amplitudes des compressions et décompressions devraient être égales, symétriques par rapport à la valeur moyenne. Si on ajoute de la matière baryonique, l'amplitude maximale (des décompressions) reste la même, vu qu'elle dépend surtout de la pression de radiation. Mais l'amplitude maximale des compressions est amplifiée par le poids de la matière baryonique. Les amplitudes sont alors inégales, plus amples pour les compressions. Ainsi, le rapport entre matière baryonique et rayonnement influence l'amplitude des ondes sonores, ce qui se remarque dans le spectre de puissance, sur le second pic.

Les pics après le troisième sont quelque peu atténués. Pour comprendre pourquoi, il faut se rappeler, les ondes sonores se propagent par des collisions entre particules. Avant et au moment de la recombinaison, les photons et autres particules avaient un temps de libre parcours moyen très élevé. Dans les unités de distance actuelle, le libre parcours moyen serait de 10 000 années lumières. Il va de soi que sur le CMB actuel, les ondes sonores dont la longueur d'onde est plus petite sont atténuées. Vu que seule une portion des particules a pu entrer en collision, les particules n'ont pas eu le temps de transmettre totalement l'onde sonore par collisions. Cela en réduit l'amplitude. La région du spectre dans laquelle se manifeste cet effet, localisée pour des longueurs l tels que l > 1000, est appelée la queue amortie.

En savoir plus modifier



Preuves de la théorie du big-bang

Dans les chapitres précédents, nous avons abordé en détail les bases théoriques de la cosmologie. En soit, la théorie cosmologique actuelle est un petit peu plus compliquée. Dans le détail, les observations astronomiques sont utilisées pour déterminer la valeur de nombreux paramètres, laissés libres dans les chapitres précédents. Le modèle actuel, le modèle LCDM, précise la valeur du facteur de Hubble, l'âge de l'univers, la courbure spatiale de l'univers, la densité de matière ordinaire, la densité de matière noire, la densité d'énergie noire et quelques autres paramètres dérivés. Ce modèle permet d'expliquer de nombreuses observations, mais il faut face à quelques problèmes. Des théories concurrentes existent, comme le modèle Rh=ct, et bien d'autres. Dans ce qui va suivre, nous allons voir quelles sont les preuves en faveur de la théorie LCDM et de la théorie de l'expansion en général. Puis, nous allons voir quels sont les problèmes rencontrés par le modèle LCDM.

Les observations en faveur de la théorie du big-bang modifier

On peut se demander quelles sont les preuves de l’expansion de l'univers (et donc, indirectement, du big-bang). On sait que la loi de Hubble et le phénomène de décalage vers le rouge de la lumière sont des prédictions de la théorie de l'expansion qui sont confirmées expérimentalement depuis belle lurette. Mais d'autres témoins sont demandés à la barre, pour que le jury puisse se faire son opinion !

La nucléosynthèse primordiale modifier

La première preuve tient dans les conséquences du découplage thermique sur la formation des noyaux et des baryons. Les équations permettent de calculer exactement le rapport entre le nombre de protons et de neutrons né de ce découplage, ainsi que le rapport du nombre d'atomes d'hydrogènes, d'hélium et de lithium. La théorie prédit que 92% des atomes de l'univers est composé d'hydrogène, 7% des atomes d’hélium, et le reste est composé d'autres éléments chimiques. Malheureusement, les étoiles ont produit de l’hélium et d'autres éléments chimiques, modifiant les rapports de proportion dans les galaxies récentes : les supernovas ont enrichi le milieu interstellaire en Carbone, Oxygène, Hélium, etc. C'est pour cela qu'on observe des rapports de l'ordre de 75% d'hydrogène, et 24% d'hélium. Pour éliminer le problème, il suffit d'observer la composition chimique des toutes premières galaxies formées, que l'on peut observer à la lisière de l'univers observable (les fameux quasars). Les observations montrent que dans ces galaxies, on observe des rapports exactement identiques à ceux prédits par la théorie du big-bang.

Les mêmes équations permettent de calculer le rapport nombre de protons sur nombre de neutrons. Là encore, celui-ci a été modifié par certaines désintégrations radioactives, mais on peut quantifier le phénomène et prédire les rapports actuels en fonction des rapports prédits lors du découplage. Et là encore, les conclusions sont exactement celles prédites par la théorie de l’expansion de l'univers.

Le fond diffus cosmologique modifier

La théorie du big-bang est la seule à prédire l'existence du fond diffus cosmologique de photons, abordé il y a quelques chapitres. Le fond diffus capté à l'heure actuelle correspond aux photons émis sur une sphère centrée sur la Terre, 380 000 ans après le big-bang environ : cette sphère est appelée la surface de dernière diffusion. Si on regarde le fond diffus, on peut remarquer que celui-ci n'est pas totalement homogène, avec quelques variations de température assez aléatoires et dispersées appelées anisotropies. On a vu il y a quelques chapitres que ces anisotropies sont les traces laissées par des ondes sonores qui se propageaient dans le plasma primordial. Diverses campagnes d'observation ont déterminé la taille des surdensités, des grumeaux de matière, et ont réussi ce qui a permis d'induire les valeurs des paramètres utilisés pour calculer la taille des surdensités.

On a vu il y a quelques chapitres que le spectre du fond diffus possède quelques pics, chacun ayant une signification relativement précise. En postulant un univers sans courbure, et avec quelques hypothèses supplémentaires, on arrive à avoir un accord parfait avec les données observationnelles, comme le montre la courbe ci-dessous. La théorie arrive ainsi à déterminer la taille des zones de surdensité en fonction de certains paramètres comme l'âge de l'univers, le facteur de Hubble, etc. C'est la première réussite de la théorie du big-bang : réussir à prédire l'état du fond diffus cosmologique, seulement à partir des équations vues plus haut. Les campagnes d'observation du fond diffus se succèdent les unes aux autres, pour tenter d'invalider ou de valider divers modèles d’univers, ou théories spéculatives sur le big-bang et les mécanismes thermodynamiques en place avec la recombinaison.

Les problèmes de la théorie cosmologique LCDM modifier

La théorie moderne de la cosmologie n'est pas sans défauts. En utilisant les paramètres mesurés par les astronomes, le modèle cosmologique marche bien pour prédire la plupart des variables cosmologiques pertinentes. Mais aucune théorie n'est parfaite et les scientifiques doivent chercher leurs défauts, en vue de les réfuter ou de les améliorer. À l'heure actuelle, l'existence de l'énergie noire, ainsi que de la matière noire, est un challenge que les théories ne peuvent pas encore expliquer correctement. Aucun mécanisme physique, aucune particule, aucun champ physique ne permet d'expliquer énergie et matière noire. Mais outre ces problèmes, le modèle standard de la cosmologie fait face à deux autres problèmes principaux : le problème de la platitude et le problème de l'horizon.

Le problème de la platitude modifier

Le premier problème de la théorie du big-bang est que l'univers est plat. On pourrait croire que ce n'est qu'une coïncidence, mais les physiciens ne sont pas vraiment satisfaits par cette explication. Beaucoup pensent que le fait que la courbure soit nulle n'est pas naturel et que c'est quelque chose qui doit être expliqué. C'est ce qu'on appelle le problème de la platitude : pourquoi la courbure semble-t-elle nulle ?

Le problème de l'horizon modifier

Problème de l'horizon.

Les observations montrent que le CMB a une température relativement uniforme, proche d'un équilibre thermique. Même des points opposés sur la surface de dernière diffusion ont des températures similaires et sont relativement homogènes entre eux. En soi, cette observation n'est pas un problème dans tous les modèles cosmologiques. Par exemple, un univers en expansion linéaire est parfaitement compatible avec un CMB en équilibre thermique total, de même qu'un univers en expansion exponentielle. Mais dans le modèle cosmologique LCDM, qui fait consensus parmi la communauté des cosmologistes, un tel équilibre thermique parfait ne peut pas exister. Des points éloignés n'ont pas pu interagir entre eux pour harmoniser leurs températures, même indirectement, sauf à supposer que l'on puisse dépasser la vitesse de la lumière. On voit ainsi apparaître le problème de l'horizon : pourquoi le fond diffus est-il si homogène ? C'est là un des points que la théorie LCDM ne peut pas expliquer.

Pour comprendre pourquoi, il faut remarquer que même en allant à la vitesse de la lumière, les photons n'ont pût parcourir qu'une certaine distance entre le big-bang et la recombinaison. On peut calculer cette distance assez simplement. Pour cela, on peut soustraire le rayon de l'univers observable actuel du rayon de l'univers observable lors de la recombinaison. Dit autrement, on souhaite calculer la différence suivante :

, avec l'âge actuel de l'univers et l'âge de l'univers lors de la recombinaison.

Utilisons alors la formule qui donne le rayon de l'univers observable :

En combinant ces deux équations, on trouve :

On factorise et on simplifie les intégrales, ce qui donne :

Après le découplage, l'univers était dominé par la matière, ce qui fait que l'on a :

En combinant les deux équations précédentes, on trouve que le tout vaut :

Le calcul donne une distance approximativement égale à :

Cette distance donne, quand elle est projetée sur la surface de dernière diffusion, un angle maximal bien précis estimé à 2°. Pour le dire plus simplement, sur ce qu'on voit du CMB depuis la Terre, deux points séparés de plus de 2° ne peuvent pas avoir communiqué entre eux et ne peuvent pas être en équilibre thermique. Pourtant, tout le CMB est à l'équilibre thermique et a une température relativement uniforme ! Même des points aux antipodes du CMB ont une température assez similaire et suivent la même loi du rayonnement de corps noir.

Rappelons que ce problème n’apparaît que pour un univers dont l'expansion ralentit dans le temps, c'est à dire un univers dominé par la matière, le rayonnement ou les deux. Ce qui colle avec le modèle LCDM actuel, qui a successivement été dominé par le rayonnement puis par la matière, avant que l'énergie noire commence à faire son effet sur une période récente. Les calculs disent que dans un univers en expansion constante (modèle linéaire) ou croissante, voire exponentielle (dominé par la constante cosmologique), le problème de l'horizon n’apparaît pas. Et cela a une importance capitale pour le prochain chapitre.



L'hypothèse de l'inflation

Inflation comparée à un modèle sans inflation.

Les problèmes du modèle LCDM peuvent être résolus par une hypothèse appelée l'inflation cosmologique. Celle-ci postule que l'univers a subi une phase d’expansion avant la recombinaison, avec une vitesse d'expansion supérieure à celle de la lumière. Avec l'inflation, la surface de dernière diffusion correspond à une portion d'espace suffisamment petite pour harmoniser les températures. Cette portion d'espace homogène a subi l'inflation cosmologique, jusqu'à devenir la surface de dernière diffusion. Le tout est illustré dans le schéma ci-contre, qui montre bien que la surface de dernière diffusion correspond à un espace assez grand sans inflation cosmologique, et à un espace bien plus petit avec elle.

Dans cette théorie, l'évolution de l'univers s'est déroulée en plusieurs étapes :

  • une étape antérieure à l'inflation, qu'on ne connaît pas très bien ;
  • une étape d'inflation exponentielle ;
  • une étape où l'univers est dominé par le rayonnement ;
  • une étape où l'univers est dominé par la matière.
Inflation

À l'heure actuelle (début 2016), il n'existe pas de preuve certaine de l'existence d'une inflation. À vrai dire, les scientifiques ne sont même pas certains du déroulement des évènements. On ne sait pas si cette inflation a eu lieu, ni quand. Par quand, on veut dire qu'elle peut avoir eu lieu n'importe quand avant la recombinaison. Dans le cas le plus simple, l'inflation aurait eu lieu avant la phase dominée par le rayonnement, mais il est aussi possible qu'elle ait eu lieu en plein milieu. Dans les faits, on ne peut pas encore discerner ces scénarios avec les observations dont on dispose et la théorie ne nous éclaire pas plus que cela... Au mieux, il s'agit d'une théorie qui fait des prédictions sur certains paramètres du fond diffus, mais ces derniers n'ont pas encore été mesurés. En 2014, la coopération internationale BICEP2 a effectuée des mesures du CMB, espérant trouver des traces de l'inflation dans le CMB, prévues par la théorie. Les résultats étaient initialement favorables et il semblait que l'inflation venait d'être validée scientifiquement, mais des analyses ultérieures ont montré que ce n'était pas le cas. Les analyses statistiques avaient omis divers effets parasites, comme l'influence de la poussière interstellaire/intergalactique et divers autres biais, réduisant les résultats à néant. Reste à voir ce que nous réserve l'avenir !

Le problème de la platitude et l'inflation modifier

Cette hypothèse résout le problème de l'horizon, mais aussi le problème de la platitude : la zone a tellement gonflé que les effets de la géométrie de l'univers seraient presque invisibles. Pour faire une analogie, imaginez que la zone gonflée par l'expansion est la surface de la Terre : vous êtes si petits par rapport à la Terre que vous ne vous rendez pas compte de sa rotondité et voyez un sol plat. Pour une explication plus mathématique, il faut recourir à la première équation de Friedmann, que voici :

Vous voyez que dans l'équation précédente, l'effet de la courbure sur l'expansion correspond au terme , où l'intensité de la courbure est modulée par le terme . Plus le facteur d'échelle augmente, plus le terme diminue et moins le terme de courbure a d'effet sur l'expansion. Sans inflation, on s'attend à ce que le facteur d'échelle suive une loi de puissance dont les exposants sont connus. D'abord on a une phase d'expansion dominée par le rayonnement avec un exposant proche de 1/2, qui dure environ 300 000 ans, puis une phase dominée par la matière de près de 13 milliards d'années où l'exposant est de 2/3. Les calculs nous disent que la courbure K devait être particulièrement basse pour coller avec les données observationnelles actuelles connues. Mais si l'inflation a bien eu lieu, a(t) aurait connu une hausse très importante pendant ou avant la phase dominée par le rayonnement, réduisant le terme à presque rien, en accord avec les données observationnelles.

Le mécanisme de l'inflation modifier

Le mécanisme supposé pour l'inflation est assez simple à comprendre à partir des équations de Friedmann. À vrai dire, une expansion accélérée s'explique déjà dans ce modèle par la constante cosmologique. Mais en supposant que la constante cosmologique permet d'expliquer l'inflation, on est face à un problème : l'inflation ne s’arrête jamais ! En effet, une constante cosmologique entraîne une expansion exponentielle de l'univers ad vitam æternam, alors que l'inflation a été un phénomène temporaire. La constante cosmologique ne peut donc pas expliquer l'inflation, et il faut trouver autre chose qui se comporte comme une constante cosmologique à un détail près : elle doit varier dans le temps, histoire d'être assez grand lors de la phase d'inflation, et disparaître ensuite.

Au niveau des équations, on peut rendre compte de cela de la même manière que pour la constante cosmologique : en rajoutant une densité d'énergie dans la première équation de Friedmann. Celle-ci s'écrit alors comme suit :

Il suffit juste que augmente pendant la phase d'inflation avant de se réduire à peau de chagrin ensuite et le tour est joué. Reste à inventer le mécanisme d'inflation qui se cache derrière . Et là, les théories sont assez nombreuses ! Il faut dire que n'importe quelle fonction mathématique ad hoc peut fonctionner tant qu'elle a quelques propriétés mathématiques peu contraignantes. La majorité des tentatives se basent sur la présence de ce qu'on appelle un champ scalaire. Un champ scalaire est quelque chose qui attribue un nombre à tout point dans l'espace à chaque instant, appelé le potentiel, noté . De plus, ce champ a une densité d'énergie en chaque point de l'espace, qui dépend du potentiel : celle-ci sera notée . Si on suppose que le champ a la même valeur en tout point, les équations nous donnent la densité d'énergie et la pression associée à ce champ.

La densité d'énergie causée par le champ dépend du potentiel et de sa variation dans le temps (qui est un "équivalent" de l'énergie cinétique pour un champ scalaire). La formule exacte est la suivante :

, avec .

La pression est de :

On peut alors en déduire , et , qui seront utilisées dans ce qui suit :

Injectons ces équations dans les équations de Friedmann, ce qui donne :

Maintenant, négligeons le terme , supposé bien plus petit que le potentiel et faisons de même pour la dérivée seconde. En clair, on suppose que , hypothèse souvent complémentée avec . Les équations précédentes deviennent ceci :

Les équations précédentes nous disent que l'univers est en expansion pour la première équation, que cette expansion est accélérée pour la seconde, et que le potentiel et la densité d'énergie restent constants malgré l'expansion pour l'équation du fluide. Les trois propriétés sont ce qu'on attend d'un univers dominé par une constante cosmologique, ou par de l'énergie noire de pression négative. Pour résumer, on obtient bien une expansion accélérée dans le cas où , hypothèse souvent complémentée avec .

Précisons cependant que si le potentiel a une valeur nulle, il ne cause aucune expansion : on peut alors modéliser l'inflation en supposant que diminue avec le temps, l'inflation s’arrêtant quand le champ est nul. Tout champ qui respecte ces conditions peut déclencher une inflation, ce qui laisse beaucoup de possibilités. Le champ lié au boson de Higgs est une solution possible, par exemple.

GFDL GFDL Vous avez la permission de copier, distribuer et/ou modifier ce document selon les termes de la licence de documentation libre GNU, version 1.2 ou plus récente publiée par la Free Software Foundation ; sans sections inaltérables, sans texte de première page de couverture et sans texte de dernière page de couverture.