Physique

Les ondes gravitationnelles

Deux trous noirs fusionnent à plus d'un milliard d'années-lumière. Sur Terre, des bras de 4 km changent de longueur d'environ un deux-centième du rayon d'un proton. Comment une telle mesure est-elle possible, et qu'y apprend-on ?

Deux astres se rapprochent en spirale, de plus en plus vite. Leur mouvement émet une onde qui s'éloigne à la vitesse de la lumière. Elle n'est pas un mouvement dans l'espace : c'est l'espace-temps lui-même qui se déforme. Un anneau de grains, plus loin, en subit l'effet quand l'onde lui parvient.

Motif de l'onde calculé avec le retard de propagation : chaque point montre ce que la source faisait un peu plus tôt. Fréquence et amplitude montent jusqu'à la fusion, puis tout s'apaise. Amplitudes très exagérées.

Planche 1

Les déformations propagatives de l'espace-temps

Loin des sources, on écrit la métrique g = η + h, avec h très petit. Les équations d'Einstein deviennent alors une équation d'ondes : □h = 0. L'onde va à la vitesse c, et elle est transverse. Elle a deux états de polarisation, notés + et ×.

L'onde arrive perpendiculairement à l'écran. Quand un axe s'étire, l'autre se comprime d'autant. Le mouvement revient identique après une demi-rotation de 180° : c'est une onde de spin 2, à la différence de la lumière, de spin 1.

À gauche, l'anneau et ses grains libres (pointillés : position de repos). À droite, l'allongement relatif selon x au cours du temps, pour un grain. Le point rouge suit l'instant présent.

Plus grand allongement du grain, marée intégrée (RK4)
Formule ½(h₊x₀ + h×y₀)/x₀
Aire de l'anneau / aire de repos (polygone)
Formule 1 − (h₊² + hײ)/4
L'aire ne change pas au premier ordre : l'onde déforme sans comprimer. Le grain ne reçoit aucune force : il reste au repos dans ses coordonnées, et c'est la distance qui varie. Cette variation se mesure comme une marée, ξ̈ = ½ ḧ ξ. C'est la même courbure que dans la page sur le principe d'équivalence, mais qui voyage.

Une file d'anneaux, posés dans la direction de propagation. Chacun se déforme avec le retard z/c. Aucun ne se déplace en z : l'onde n'agit qu'en travers.

Anneaux vus en perspective. À droite, une impulsion d'onde qui obéit à ∂²h/∂t² = c²∂²h/∂z² : elle est calculée sur une grille.

Vitesse mesurée du sommet de l'impulsion (grille)
Vitesse attendue1 c
Déplacement selon z d'un grain0
Einstein prédit ces ondes en 1916, puis doute en 1936-1937 qu'elles existent, avant de se corriger. La question reste débattue jusqu'à la conférence de Chapel Hill (1957), où l'argument de la « tige à perles » de Feynman montre qu'elles transportent bien de l'énergie : si les perles glissent, elles chauffent la tige.

Pourquoi cela a-t-il pris un siècle ? Parce que l'espace-temps est très raide. Le signal typique arrivant sur Terre a h ≈ 10⁻²¹. Pourtant, il transporte beaucoup d'énergie, car le flux est F = (c³/16πG) ⟨ḣ²⟩, avec un énorme préfacteur.

Échelle des longueurs, en mètres. Le trait rouge est la variation ΔL = hL.

Variation ΔL = hL
Comparée au rayon du proton (≈ 8,4·10⁻¹⁶ m)
Flux d'énergie de cette onde
c³/(16πG)
À 100 Hz et h = 10⁻²¹, le flux dépasse 1 mW/m² : c'est du même ordre que la lumière de la pleine Lune. L'énergie est bien là, mais la matière n'y réagit presque pas. Elle traverse donc tout, sans être absorbée ni diffusée : c'est un messager très fiable.
Planche 2

Les sources binaires

Une onde gravitationnelle exige un quadrupôle de masse qui varie. Une boule qui gonfle ou une sphère qui tourne n'émet rien. Mais deux astres en orbite conviennent parfaitement. Einstein établit la formule en 1918 : la puissance vaut P = (G/5c⁵) ⟨Q⃛ijQ⃛ij⟩, où Q est le moment quadrupolaire.

Pour une orbite circulaire, Q revient identique après une demi-orbite. L'onde a donc la fréquence f = 2 forb. Et Peters obtient en 1963 P = (32/5) G⁴ μ² M³ / (c⁵ a⁵), avec M = m₁ + m₂ et μ = m₁m₂/M. Le 1/a⁵ dit tout : il faut des astres serrés.

À gauche, le couple (les tailles ne sont pas à l'échelle). À droite, Qxx au cours du temps : deux périodes de Q pour une orbite.

Puissance par la dérivée troisième de Q (numérique)
Formule de Peters
Fréquence de l'onde 2forb
Durée avant la fusion, 5c⁵a⁴/(256G³m₁m₂M)
La Terre et le Soleil rayonnent environ 200 W, soit quelques ampoules. Il leur faudrait 10²³ ans pour fusionner. Deux étoiles à neutrons à 1 000 km émettent, elles, 10¹⁴ fois plus que le Soleil en lumière, et fusionnent en une heure environ.

En perdant de l'énergie, les astres se rapprochent, tournent plus vite, et rayonnent encore plus. La fréquence monte comme un chirp, un gazouillis d'oiseau. Une seule combinaison des masses règle cette montée : la masse de chirp Mc = (m₁m₂)^{3/5}/M^{1/5}.

En haut, la déformation h(t) sur la dernière seconde ou moins. En bas, la fréquence selon le temps. Modèle à l'ordre le plus bas, valable jusqu'à la dernière orbite stable.

Masse de chirp, retrouvée depuis f et ḟ numériques
Formule (m₁m₂)^{3/5}/M^{1/5}
Durée depuis 20 Hz : RK4 sur a(t)
Formule (5/256)(GMc/c³)^{−5/3}(πf)^{−8/3}
Cycles depuis 20 Hz (somme) / formule
Fréquence à la dernière orbite stable, c³/(6^{3/2}πGM)
Le son que vous entendez a la vraie fréquence de l'onde, et dure le temps réel quand c'est possible : ce n'est pas un son réel dans l'espace, mais un signal de même forme. Les masses légères chantent longtemps et finissent haut. Les plus lourdes font un bref coup de sifflet grave.

La première preuve vient d'un pulsar découvert en 1974 par Hulse et Taylor. Son orbite excentrique rétrécit. Son passage au périastre arrive donc de plus en plus tôt, et l'écart cumulé grandit comme le carré du temps. Taylor et ses collègues le mesurent dès 1979 : prix Nobel en 1993.

Avance cumulée du passage au périastre depuis 1975, calculée avec la formule de Peters. Ce n'est pas le tracé des données publiées, mais la courbe qu'elles suivent.

Ṗ de la relativité générale (Peters, e = 0,617)
Ṗ mesuré, corrigé de la Galaxie (Weisberg et Huang)−(2,398 ± 0,004)·10⁻¹²
Facteur d'excentricité f(e)
Avance cumulée cette année
Durée avant la fusion (RK4 sur a, e)
Le rapport mesure/prévision vaut 0,9983 ± 0,0016. Ces ondes ont donc été détectées indirectement quarante ans avant LIGO. Les deux étoiles à neutrons finiront par fusionner, dans environ 300 millions d'années : il existe donc bien des binaires qui mènent à ce qu'observent les détecteurs.
Planche 3

La fusion des trous noirs et des étoiles à neutrons

Le signal a trois actes : la spirale (l'approche), la fusion (qu'on ne décrit que par simulation numérique, depuis 2005), et les oscillations de l'objet final. Les trous noirs ne laissent que de l'onde. Les étoiles à neutrons laissent aussi de la matière et de la lumière.

En haut, schéma du signal (spirale, fusion, oscillation amortie). En bas, ce qui se passe. Le raccord entre les phases est simplifié, et les valeurs de l'objet final viennent d'ajustements sur des simulations, pour des trous noirs sans rotation initiale : pour GW231123, dont les spins sont élevés, elles ne sont qu'indicatives.

Énergie rayonnée (ajustement)
Masse et spin de l'objet final
Fréquence de l'oscillation finale (formule de Berti)
Fréquence lue sur le signal (passages à zéro)
Durée d'amortissement
Aire finale / somme des aires initiales
Unité naturelle de puissance, c⁵/G
GW150914 a rayonné environ trois masses solaires en quelques dixièmes de seconde, avec un pic supérieur à la puissance lumineuse de tout l'Univers observable. L'aire finale dépasse la somme des aires, comme l'exige Hawking : c'est le test de GW250114. Selon la formule de Kerr, la fréquence finale dépend seulement de M et du spin : aucun « cheveu ».

Plus légères, les étoiles à neutrons restent bien plus longtemps dans la bande du détecteur. Elles se déforment par marée juste avant de se toucher, ce qui révèle leur structure. En fusionnant, elles éjectent de la matière très riche en neutrons, qui forme des éléments lourds et brille : une kilonova. C'est ce qui s'est passé le 17 août 2017.

Faites défiler le temps depuis la fusion, en échelle logarithmique. À gauche, un schéma ; à droite, l'arrivée successive des messagers de GW170817.

Temps dans la bande depuis 30 Hz : étoiles à neutrons
Le même pour GW150914
Fréquence quand les étoiles se touchent, √(GM/(2R)³)/π
Écart de vitesse |cg − c|/c, si le retard est tout dû à la vitesse
Limites publiées−3·10⁻¹⁵ à +7·10⁻¹⁶
Les rayons gamma sont arrivés 1,7 s après la fin du signal, après 130 millions d'années de voyage. La gravité voyage donc à la vitesse de la lumière à 10⁻¹⁵ près, ce qui a éliminé d'un coup de nombreuses théories alternatives. Les éléments comme l'or et le platine sont, au moins en partie, formés dans ces événements.
Planche 4

Le principe des détecteurs interférométriques

Un interféromètre compare la longueur de deux bras perpendiculaires. Quand l'onde passe, l'un s'allonge, l'autre raccourcit. Le laser se sépare en deux, va-et-vient, puis se recombine. On règle tout pour que les deux faisceaux s'annulent en sortie : la moindre variation laisse échapper de la lumière.

Pour h = 10⁻²¹ et L = 4 km, ΔL ≈ 4·10⁻¹⁸ m. Deux astuces gagnent énormément : les cavités de Fabry-Perot font faire à la lumière des centaines d'allers-retours, et le recyclage de puissance accumule la lumière dans les bras, jusqu'à 750 kW en régime nominal.

Vue de dessus, déformation très exagérée. Les miroirs bougent, le faisceau de sortie (photodiode, à droite) s'allume et s'éteint.

ΔL différentiel, hL
Longueur effective avec les allers-retours, 2FL/π
Déphasage Δφ = 4π h Leff/λ
Lumière en sortie, sin²(Δφ/2)
Approximation (Δφ/2)²
Bras optimal pour 100 Hz, c/(4f)
Limite quantique sur 10 ms, avec 750 kW (ordre de grandeur)
Il y a bien un paradoxe : si l'onde allonge aussi la lumière, comment la mesurer ? L'onde lumineuse est bien étirée, mais celle qui arrive ensuite vient de plus loin. Ce qui compte, c'est le temps d'aller-retour. Il change, et c'est ce qu'on lit à la sortie.

Un interféromètre n'a pas de lentille : il entend presque tout le ciel, mais avec une sensibilité qui dépend de la direction. Il faut plusieurs détecteurs et leurs temps d'arrivée pour localiser la source : Hanford et Livingston sont à 10 ms-lumière l'un de l'autre, ce qui dessine une couronne sur le ciel.

Sensibilité d'un détecteur en L selon la direction de la source (toutes directions, carte en longitude et latitude). Plus clair : plus sensible. Le point rouge marque la direction choisie.

Réponse √(F₊² + Fײ) pour cette direction
Moyenne quadratique sur le ciel (intégrale numérique)
Valeur exacte √(2/5)
Retard maximal entre les détecteurs, base/c
Avec deux détecteurs, GW150914 a été localisée à environ 600 degrés carrés. Avec Virgo en plus, GW170817 à 28 degrés carrés : cela a suffi pour trouver sa galaxie en 11 heures. Il y a aussi une direction aveugle, mais aucun point du ciel n'est masqué par la Terre : l'onde la traverse.

Le signal doit émerger du bruit. Aux basses fréquences, les vibrations du sol dominent, malgré une isolation extrême. Au milieu, le bruit thermique des miroirs. Aux hautes fréquences, le bruit quantique des photons. Le rapport signal sur bruit se calcule en intégrant le signal contre le bruit, fréquence par fréquence.

Courbe de bruit simplifiée (bleu) et signal d'une spirale (rouge), en strain par √Hz. La courbe de bruit est un modèle à trois termes, et non une courbe mesurée.

Rapport signal sur bruit ρ (Simpson)
Distance où ρ = 8, orientation idéale (dichotomie)
La même par la loi ρ ∝ 1/D
Volume surveillé, rapporté à celui du détecteur actuel
Doubler la sensibilité permet d'aller deux fois plus loin, donc d'observer huit fois plus d'événements. C'est pourquoi LIGO et Virgo ont passé de 3 candidats en 2015 à plus de 200 en dix ans. Les valeurs de ρ sont calculées pour un seul détecteur, une orientation idéale et un modèle simple de bruit : elles ne remplacent pas les résultats publiés. En moyenne sur les orientations, il faut diviser la distance par environ 2,26.
Planche 5 · signature

Les informations apportées par les observations

Un signal à peine plus gros que le bruit livre pourtant les masses, la distance, l'orientation et un test de la relativité générale. La clé est le filtrage adapté : comparer les données à des milliers de formes d'onde prévues.

Voici 2 secondes de données. Un chirp y est caché, invisible à l'œil. Essayez un gabarit, et le calcul le fait glisser sur les données en mesurant la corrélation. Quand gabarit et signal coïncident, un pic apparaît. La hauteur indique la force du signal, la position le temps et la masse de chirp celle qui donne le plus haut pic.

En haut, les données brutes (bruit blanc simulé), et le signal caché si vous le révélez. En bas à gauche, ρ(t) pour le gabarit choisi. En bas à droite, le pic selon la masse de chirp du gabarit.

Pic de ρ avec ce gabarit
Meilleure masse de chirp trouvée (balayage)
Vraie masse de chirp (une fois révélée)
Temps de fusion trouvé / vrai
Plus grand pic dû au seul bruit (mêmes gabarits, données sans signal)
Essayez ρ = 4 ou 5 : le pic se confond avec les fluctuations du bruit. C'est pourquoi les détections exigent d'ordinaire ρ ≈ 8, et deux détecteurs en coïncidence. Le véritable calcul se fait en fréquence, avec un bruit coloré et des gabarits complets, mais l'idée est bien celle-ci.

L'amplitude d'un chirp dépend de la masse de chirp, que l'on connaît par la fréquence, et de la distance. Le signal donne donc la distance directement, sans échelle des distances cosmiques : c'est une sirène standard. Si la galaxie est identifiée, on obtient v/D, la constante de Hubble.

Valeurs de H₀ et leurs barres d'erreur. Le point rouge est votre calcul. Les deux lignes vertes rappellent l'écart entre mesures précoces et tardives.

H₀ = v/D, avec v = 3 017 km/s (GW170817)
Amplitude h₊ relative pour cette inclinaison, (1 + cos²ι)/2
Erreur relative attendue avec N sirènes (∝ 1/√N)
Sirènes pour départager 67,4 et 73 à 2σ (dichotomie)
Publiées : GW170817 ; GWTC-470⁺¹²₋₈ ; 76,6⁺¹³,⁰₋₉,₅
L'inclinaison du plan et la distance sont liées : une orbite vue en biais paraît plus faible, donc plus lointaine. C'est la principale source d'erreur. Le catalogue GWTC-4 utilise aussi des événements sans contrepartie lumineuse, d'où un résultat encore peu contraignant : il n'a pas tranché le désaccord, mais il est indépendant.

En septembre 2015, GW150914 ouvrait une nouvelle astronomie. Le catalogue GWTC-4 compte 218 candidats. Il montre des trous noirs bien plus massifs que ceux connus par les rayons X, et qui semblent se regrouper autour de quelques masses, ce qui suggère plusieurs voies de formation.

Quelques événements remarquables, en masses (échelle logarithmique). Les valeurs sont les médianes publiées.

Masse de chirp
Rapport des masses
Énergie rayonnée (ajustement, sans spin)
Ce qu'il a apporté
Chaque événement est aussi un test. La vitesse, la forme de l'onde et les oscillations finales n'ont, à ce jour, montré aucun écart à la relativité générale. Aucun trou noir n'a été vu autrement que selon la solution de Kerr. La question de la formation (étoiles isolées, amas, fusions successives) reste en revanche ouverte.

Comme pour la lumière, chaque gamme de fréquence a ses sources et ses instruments. Le fond à nanohertz serait dû à des couples de trous noirs géants : plusieurs réseaux de pulsars en ont annoncé des indices en 2023. LISA, mission spatiale de l'ESA, doit écouter les millihertz. Le Télescope Einstein et d'autres projets visent à aller plus loin depuis le sol.

Les bandes de fréquence des instruments, en échelle logarithmique. La ligne rouge est la fréquence choisie.

Période
Longueur d'onde c/f
Instrument
Séparation d'un couple de 10⁶ M☉ à cette fréquence

Un siècle d'attente

L'onde

Une déformation transverse de l'espace-temps, à deux polarisations, qui se propage à la vitesse c et traverse la matière presque sans perte.

Les sources

Les couples compacts dessinent un chirp dont la pente donne la masse de chirp. La fin du signal révèle l'objet final, ou la matière d'une étoile à neutrons.

La mesure

Deux bras, un laser et le filtrage adapté extraient h ≈ 10⁻²¹. On y lit masses, distance, H₀, éléments lourds et test de la gravité.
Pour aller plus loin

Pourquoi deux polarisations ? La métrique a dix composantes, mais quatre sont liées par la jauge, quatre par les contraintes, et il en reste deux qui se propagent. En jauge dite transverse et sans trace (TT), seules h₊ et h× restent.

Redshift. Pour une source lointaine, la masse apparente est multipliée par (1 + z), et le signal ne permet pas seul de les distinguer. Les masses citées ici sont celles de la source, déjà corrigées.

Limites de cette page. La spirale est à l'ordre le plus bas (relativité restreinte du mouvement, ni spin ni ordres supérieurs). Le filtrage adapté utilise du bruit blanc et des gabarits simples.

Dans la même série. trous-noirs-v1.html pour les orbites et l'horizon ; equations-einstein-v1.html pour les équations ; tests-relativite-generale-v1.html pour les autres tests.

Sources principales. LIGO-Virgo, Physical Review Letters 116 et 119 ; Weisberg et Huang, 2016 ; LIGO-Virgo-KAGRA, GWTC-4.0, 2025 ; Aasi et coll., Advanced LIGO, 2015.