Physique

Les effets et les tests de la relativité générale

Une théorie vaut ce que valent ses prévisions. Celle d'Einstein annonce que la lumière se décale et se courbe, que les horloges dépendent de l'altitude, qu'une orbite tourne et qu'un signal radar prend du retard. Voici comment on l'a vérifié, chiffres à l'appui, jusque dans le GPS.

Tous les effets de cette page viennent de la même métrique, celle de Schwarzschild autour d'une masse. Ils tiennent dans un seul dessin. Les effets y sont très exagérés.

Autour du Soleil : l'orbite de Mercure tourne lentement sur elle-même, un rayon d'étoile est dévié, un signal radar ralentit en passant près du Soleil, et une horloge proche bat moins vite qu'une horloge lointaine.

Planche 1

Le décalage des fréquences et le ralentissement des horloges

Ce sont deux lectures d'un même fait. Dans un champ statique, tous les signaux mettent le même temps à monter : ceux qu'une horloge envoie en bas arrivent donc en haut au même rythme, lu avec le temps t de la métrique. Mais une horloge ne mesure pas t. Elle mesure dτ = √(1 − rs/r) dt. En bas, elle compte moins de battements. Le signal reçu en haut paraît donc plus lent : c'est le décalage vers le rouge gravitationnel.

Une onde monte, sa longueur d'onde s'allonge (très exagéré). À droite, les battements émis en bas et ceux reçus en haut, écart lui aussi exagéré.

Décalage Δν/ν, métrique exacte
Premier ordre gh/c²
Vitesse de source qui compenserait (Doppler)
Avance de l'horloge du haut, par an
Pound et Rebka (1960), sur 22,5 m, un sens(2,56 ± 0,25)·10⁻¹⁵
Pound et Rebka comparent, à Harvard, des rayons gamma qui montent et qui descendent une tour, grâce à l'effet Mössbauer : ils compensent le décalage en déplaçant la source à moins d'un micromètre par seconde. Gravity Probe A, en 1976, emporte une horloge à 10 000 km. En 2010, des horloges optiques au NIST voient l'effet pour une différence de hauteur de 33 cm.

Plus r se rapproche de rs = 2GM/c², plus l'horloge ralentit. Le rythme réel figure sur les deux cadrans : pour la Terre et le Soleil, il est indiscernable ; on ne voit l'écart qu'à l'approche d'un astre très compact.

À gauche, l'horloge lointaine (verte) et l'horloge locale (rouge), sans exagération. À droite, l'écart de rythme 1 − dτ/dt selon r/rs, en échelle logarithmique.

Rythme dτ/dt = √(1 − rs/r)
Écart 1 − dτ/dt
Premier ordre rs/(2r)
Retard accumulé en un an sur l'horloge lointaine
L'horloge locale ne ressent rien : pour elle, une seconde dure une seconde. Le ralentissement n'apparaît que dans la comparaison entre deux lieux. La naine blanche Sirius B fut l'un des premiers tests, et l'étoile S2 a montré en 2018 ce décalage près du trou noir de notre galaxie.

En octobre 1971, Hafele et Keating font le tour du monde en avion de ligne avec quatre horloges au césium. Deux effets s'y mêlent. L'altitude fait gagner du temps. La vitesse, vue d'un repère non tournant, change avec le sens de vol : vers l'est, l'avion s'ajoute à la rotation de la Terre et l'horloge perd ; vers l'ouest, elle gagne.

La Terre vue du pôle nord, dans un repère qui ne tourne pas. L'avion suit son parallèle. À droite, le modèle de cette page, la prévision publiée et la mesure, en nanosecondes.

Modèle : gain d'altitude gh/c² × T
Modèle : effet de vitesse
Modèle : total
Prévision publiée (gravité + vitesse)
Mesure publiée
Le modèle est volontairement simple : vol à vitesse, altitude et latitude constantes (30°). Le vrai calcul suivait les plans de vol. Il donne malgré tout le bon signe et le bon ordre de grandeur, et montre que les deux termes ont le même poids à ces vitesses.
Planche 2

La déviation de la lumière

En suivant la métrique de Schwarzschild, un rayon de lumière de paramètre d'impact b vérifie u″ + u = (3/2) rs u², avec u = 1/r. L'angle total est, pour b grand devant rs, α = 4GM/(c²b). Un corpuscule de lumière newtonien ne donnerait que la moitié : ce facteur 2, dû à la courbure de l'espace, est le vrai test.

À gauche, de vrais rayons calculés près d'un trou noir, sans exagération. Le rouge est celui que vous choisissez. À droite, l'angle selon b : calcul complet (rouge), formule (vert) et Newton (orange).

Déviation α de ce rayon (RK4)
Premier ordre 2rs/b
Newton, tan(α/2) = rs/(2b)
b critique de capture (dichotomie)
Formule (3√3/2) rs
Rayon frôlant le Soleil (RK4)
Sous le b critique, le rayon ne ressort plus : il est capturé. C'est ce bord que dessine l'ombre du trou noir imagée par l'Event Horizon Telescope. Près du Soleil, l'effet est mille fois plus petit qu'une pleine Lune, et c'est pourquoi il a fallu attendre une éclipse.

Pendant une éclipse, on photographie les étoiles proches du Soleil, puis on compare leur position à celle prise de nuit. L'écart suit Δ = L · R☉/r. Tout est dans l'ajustement de L, en présence d'erreurs de mesure. Reproduisez-le : vous choisissez la « vraie » théorie, le nombre d'étoiles et le bruit de mesure.

À gauche, le ciel autour du Soleil : déplacements des étoiles, très agrandis. À droite, les écarts mesurés selon R☉/r, avec la droite ajustée, la droite d'Einstein (vert) et celle de Newton (orange).

L ajusté au bord du Soleil (moindres carrés)
Incertitude de cet ajustement
Écart à Einstein, en σ
Écart à Newton, en σ
1919 : Sobral et Principe (écarts-types)1,98 ± 0,18″ ; 1,61 ± 0,45″
Un troisième instrument, à Sobral, donnait 0,93″ ; il fut écarté pour de grosses erreurs systématiques, ce qui a nourri une longue polémique. Une nouvelle analyse (Gilmore et Tausch-Pebody, 2020) conclut que ce rejet était justifié et ne change presque rien au résultat. Aujourd'hui, l'interférométrie radio donne γ = 0,9998 ± 0,0004, sans aucune éclipse.
Planche 3

Les lentilles gravitationnelles

Si une masse se trouve entre nous et une source lointaine, les rayons lui passent des deux côtés. On voit alors plusieurs images. Pour une masse ponctuelle, l'équation de la lentille est β = θ − θE²/θ, où β est la vraie position de la source, θ celle de l'image et θE le rayon d'Einstein.

À gauche, la source telle qu'elle est. À droite, ce qu'on voit : chaque pixel est calculé en remontant le rayon. Le cercle pointillé est le rayon d'Einstein. La croix bleue est la lentille.

Images comptées dans le calcul
θ₊ et θ₋ par dichotomie (γ = 0)
Formule (β ± √(β² + 4))/2
Aire des images / aire de la source (pixels, source étendue)
Source ponctuelle : (u² + 2)/(u√(u² + 4))
Essayez β petit : un anneau d'Einstein apparaît. Ajoutez une marée externe et réduisez β : quatre images se forment, comme dans la croix d'Einstein. La luminosité de surface se conserve : une lentille agrandit la source et ne chauffe rien. C'est pour cela que l'aire des images donne l'amplification.

Quand la lentille est une simple étoile, les images sont trop proches pour être séparées, et on ne voit qu'une augmentation de lumière. La source semble s'allumer puis s'éteindre, selon une courbe symétrique et indépendante de la couleur.

À gauche, le passage de la source derrière la lentille et ses deux images. À droite, l'amplification au cours du temps, en unités du temps de traversée.

Amplification actuelle : μ₊ + μ₋, dérivées numériques
Formule (u² + 2)/(u√(u² + 4))
Amplification maximale
Paczyński propose en 1986 de chercher ces événements parmi des millions d'étoiles. Les premiers sont annoncés en 1993. Une petite bosse sur la courbe trahit parfois une planète autour de la lentille : le microlentillage est aussi une méthode de détection d'exoplanètes.

Le rayon d'Einstein vaut θE = √(4GM DLS/(c² DL DS)). Il dépend de la masse totale de la lentille, visible ou non. Voilà pourquoi les lentilles pèsent la matière noire : les arcs ne se trompent pas sur la masse. En 1936, Einstein jugeait l'effet d'une étoile hors de portée. Un an plus tard, Zwicky proposait que des galaxies y parviendraient.

Source, lentille et observateur. Les distances ne sont pas à l'échelle et l'angle est très exagéré.

Rayon d'Einstein θE
Taille physique dans le plan de la lentille
Foyer du Soleil, R☉²c²/(4GM)
Refsdal a vu dès 1964 que le retard entre deux images donne la constante de Hubble. Les images d'un même quasar arrivent en effet à des dates différentes : une partie de ce délai est un effet Shapiro (planche 4). Le premier quasar double, Q0957+561, date de 1979, et le premier arc géant, découvert à Toulouse par Soucail et ses collègues, de 1987.
Planche 4

Le périhélie de Mercure et le retard de Shapiro

Le Système solaire offre un laboratoire d'une grande précision. Le premier test de la relativité générale y était connu cinquante ans avant elle. Le dernier grand test classique n'y a été imaginé qu'en 1964.

L'ellipse de Mercure tourne de 574″ par siècle. Les autres planètes en expliquent 531″. Le Verrier relève l'écart dès 1859 et imagine une planète invisible, Vulcain, qu'on ne trouvera jamais. En novembre 1915, Einstein obtient 43″ : aucun paramètre n'est ajusté.

À gauche, l'orbite calculée, l'effet amplifié de A. À droite, le bilan en secondes d'arc par siècle.

Avance par tour, orbite amplifiée (RK4)
A × 6πGM/(c²a(1 − e²))
Sans amplification (A = 1, RK4), par siècle
Observé, résidu (Lo et coll., 2013)42,98 ± 0,04″
Temps pour un tour complet à cette vitesse
L'effet est très petit : 43″ en cent ans, soit environ 0,1″ par tour. Il ne se voit que parce qu'il s'additionne sur des siècles d'observations. Ces 43″ tiennent en une formule, qui ne contient que G, M, c et l'orbite. Aux fortes amplifications, l'avance dépasse la formule, qui n'est qu'un premier ordre.

La même formule vaut ailleurs. Plus l'orbite est serrée et allongée, plus l'effet est grand. Il existe même des astres où il saute aux yeux : un pulsar double et une étoile qui frôle un trou noir.

CorpsaePar tourCumulRéférence

Avance par tour, calculée par la formule, en échelle logarithmique. Chaque graduation multiplie par 100.

Les valeurs « calculées » viennent de la période et de la masse, la distance a s'en déduisant par la loi de Kepler. Pour Mercure, 42,98″. Pour S2, 12,2′ par tour ici, 12,1′ selon GRAVITY, qui mesure 1,10 ± 0,19 fois cette valeur. Le pulsar PSR B1913+16 avance de plus de 4° par an, et sa perte d'énergie donne en plus la première preuve des ondes gravitationnelles.

Irwin Shapiro remarque en 1964 qu'un signal radar passant près du Soleil doit revenir plus tard. Avec la métrique, la vitesse dr/dt de la lumière, mesurée de loin, vaut environ c(1 − (1 + γ)GM/(rc²)). Pour un observateur local, elle reste c. L'aller-retour prend alors, en plus :

Δt = 2(1 + γ)(GM/c³) · ln[(r1 + r2 + R)/(r1 + r2 − R)]

À gauche, la Terre, le Soleil et la cible : un signal fait l'aller-retour, ralenti près du Soleil. À droite, le retard selon la position de la cible sur son orbite, avec le pic de la conjonction.

Retard aller-retour : intégrale le long du rayon (Simpson)
Formule ci-dessus
Plus petite distance du rayon au Soleil
Pic au ras du Soleil pour cette cible
Cassini (2003) : γ − 1(2,1 ± 2,3)·10⁻⁵
Les premières mesures sur Vénus et Mercure en 1966-1967 atteignaient environ 20 %. Les atterrisseurs Viking, sur Mars, descendent à 0,1 % en 1979. Cassini, en 2002, fait mieux en mesurant la variation de fréquence du signal lors de la conjonction : c'est le test le plus précis à ce jour de la courbure de l'espace.
Planche 5 · signature

Les corrections du GPS

Un récepteur GPS mesure l'heure d'arrivée de signaux envoyés par des satellites, et en tire sa position. Une erreur d'une nanoseconde fait 30 cm. Les horloges embarquées subissent justement les effets de la planche 1 : le GPS est un test de la relativité générale qui fonctionne en permanence.

Le satellite est plus haut : son horloge gagne. Il va vite : elle perd. Le bilan fait 38,6 µs par jour. On règle donc l'oscillateur en usine sur une fréquence très légèrement plus basse. Mais une orbite n'est jamais parfaitement circulaire : l'altitude et la vitesse varient, ce qui ajoute un terme périodique, que le récepteur corrige lui-même.

À gauche, les deux effets et le bilan. À droite, l'écart périodique de l'horloge sur un tour, dû à l'excentricité e : calcul pas à pas (rouge) et formule de Ashby (vert).

Altitude : gain par jour
Vitesse : perte par jour
Bilan par jour
Fréquence d'usine, avec Φ₀/c² du géoïde
Valeur publiée (Ashby)10,22999999543 MHz
Excentricité : amplitude (Simpson)
Formule 4,4428·10⁻¹⁰ e √a
Le tout premier test date de 1977 : le satellite NTS-2 mesure une avance de 442,5·10⁻¹², pour une prévision de 446,5·10⁻¹². Pour vérifier, son horloge a d'abord fonctionné sans la correction, puis on l'a activée. L'écart est inférieur à 1 %.

On lit partout que, sans relativité, le GPS se tromperait de 11 km par jour. C'est la dérive des horloges exprimée en distance, mais pas directement l'erreur de position. Le récepteur résout quatre inconnues, x, y, z et l'erreur de sa propre horloge. Une dérive commune à tous les satellites est donc absorbée par cette dernière. Ce qui nuit, c'est l'écart entre satellites. Résolvons le vrai système, avec six satellites au-dessus de Paris.

À gauche, le récepteur, les signaux et le décalage de chaque horloge. À droite, l'erreur de position trouvée selon le temps écoulé depuis la dernière remise à l'heure. Trait orange : la dérive de 38,6 µs par jour, en distance.

Erreur de position (moindres carrés)
Erreur d'horloge du récepteur, en distance
Dérive relativiste en distance, c × δ × t
Résidu de l'ajustement (cohérence du calcul)
Ce calcul simplifie : satellites fixes, pas d'erreurs atmosphériques, remises à l'heure imaginées et régulièrement espacées. Il montre que l'erreur suit l'écart entre satellites, et non le temps écoulé : environ 3 km par heure d'écart ici. La conclusion tient : sans le réglage relativiste, le temps GPS dériverait de 38 µs par jour par rapport au temps terrestre, et il suffirait que les satellites ne soient pas recalés ensemble pour que la position se dégrade. Il faudrait aussi transmettre et corriger beaucoup plus de dérive.

Toutes les corrections n'ont pas la même importance. Voici leur effet en distance, calculé pour un satellite GPS.

Effet en mètres, échelle logarithmique. Chaque graduation multiplie par 10. Les trois premières lignes s'entendent par jour.

Sagnac, maximum à l'équateur : Ω R a sinθ/c²
Valeur usuelle133 ns
Shapiro, satellite à 30° d'élévation, aller
Erreur d'une nanoseconde
Ce qu'il faut retenir. Horloges en altitude, vitesse orbitale, excentricité, rotation de la Terre : chaque terme est calculé et appliqué. Le système Galileo a de plus servi de test, grâce à deux satellites placés par erreur sur une orbite elliptique en 2014, qui ont vérifié le décalage gravitationnel mieux que Gravity Probe A.

Le bilan des tests

EffetExpériencePrévisionMesure
Décalage gravitationnelPound et Rebka, 19602,46·10⁻¹⁵ sur 22,5 m(2,56 ± 0,25)·10⁻¹⁵
Horloges et altitudeHafele et Keating, 1971, vol ouest+275 ± 21 ns+273 ± 7 ns
Déviation de la lumièreDyson, Eddington, Davidson, 19191,75″1,98 ± 0,18″ ; 1,61 ± 0,45″
Déviation (radio)VLBI, 2004γ = 10,9998 ± 0,0004
Périhélie de MercurePlanétologie, 201342,98″ par siècle42,98 ± 0,04″
Retard de ShapiroCassini, 2003γ − 1 = 0(2,1 ± 2,3)·10⁻⁵
Horloge de satelliteNTS-2, 1977446,5·10⁻¹²442,5·10⁻¹²

Le temps

Décalage et ralentissement sont un seul effet : √(1 − rs/r). On le mesure de la tour au trou noir.

L'espace

Déviation, lentilles et retard Shapiro font intervenir la courbure de l'espace. Elle est testée à γ = 1 à 2·10⁻⁵ près.

Les limites

Toutes ces mesures sont en champ faible ou modéré. Les tests en champ fort, ondes gravitationnelles et images de trous noirs, sont le terrain actuel.
Pour aller plus loin

Le paramètre γ. On évalue une théorie avec le formalisme PPN : γ dit combien d'espace courbe produit une unité de masse. La déviation et le retard sont proportionnels à (1 + γ)/2 ; le périhélie à (2 + 2γ − β)/3.

Que teste vraiment le décalage ? Il découle du principe d'équivalence : il teste moins les équations d'Einstein que l'idée de l'espace-temps courbe. La déviation complète, le périhélie et le retard testent davantage les équations.

Une prudence. Le test de Mercure dépend aussi de l'aplatissement du Soleil, qui ajoute environ 0,03″ : un faible aplatissement mesuré a clos le débat des années 1960.

Dans la même série. equations-einstein-v1.html pour d'où vient la métrique ; espace-temps-courbe-v1.html pour la géométrie ; principe-equivalence-v1.html pour les marées et le décalage en fusée.

Sources principales. Will, Living Reviews, 2014 ; Ashby, Living Reviews, 2003 ; Hafele et Keating, Science, 1972 ; Gilmore et Tausch-Pebody, 2020 ; GRAVITY Collaboration, Astronomy & Astrophysics, 2020 ; Everitt et coll. 2011.