🔭 astrophysics

Revisiting GW170817 at milliarcsecond scale: high-precision constraints on jet geometry and $H_0$

Cet article présente une analyse bayésienne robuste de toutes les données VLBI disponibles pour GW170817, utilisant de nouveaux modèles hydrodynamiques pour contraindre précisément l'angle de vision du jet et déduire une constante de Hubble de $65.5\pm4.4$ km s $^{-1}$ Mpc $^{-1}$ , une valeur qui s'aligne plus étroitement sur les mesures Planck de l'Univers primordial que sur les résultats SH0ES de l'Univers tardif.

Auteurs originaux : Kelly Gourdji, Adam T. Deller, Chris Flynn, Taya Govreen-Segal, Cullan Howlett, Kunal P. Mooley, Ehud Nakar

Publié 2026-05-13

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Kelly Gourdji, Adam T. Deller, Chris Flynn, Taya Govreen-Segal, Cullan Howlett, Kunal P. Mooley, Ehud Nakar

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez l'univers comme un ballon géant en expansion. Depuis des décennies, les scientifiques débattent de la vitesse exacte à laquelle ce ballon se gonfle. Ce taux est appelé la constante de Hubble ( $H_0$ ). Le problème est que, lorsque vous la mesurez en utilisant l'« univers primordial » (comme le rayonnement fossile du Big Bang), vous obtenez un chiffre. Lorsque vous la mesurez en utilisant l'« univers tardif » (comme des étoiles explosives à proximité), vous obtenez un chiffre différent, légèrement plus élevé. Ce désaccord est connu sous le nom de « tension de Hubble ».

En 2017, les astronomes ont réalisé une percée. Ils ont détecté un crash cosmique : la collision de deux étoiles à neutrons. Cet événement, nommé GW170817, était spécial car il a émis à la fois des ondes gravitationnelles (des ondulations dans l'espace-temps) et de la lumière (un flash brillant). Cela en a fait une « sirène standard » — un phare cosmique capable d'aider à mesurer le taux d'expansion de l'univers.

Cependant, il y avait un piège. Pour mesurer la distance du phare, vous devez connaître exactement son inclinaison par rapport à nous. Si vous regardez une toupie en rotation de profil, elle a une apparence différente de celle que vous avez en la regardant de dessus. Si vous devinez mal l'inclinaison, vous devinez mal la distance, et votre calcul du taux d'expansion de l'univers est erroné.

Le problème du « point aveugle »
Pour GW170817, les détecteurs d'ondes gravitationnelles ne pouvaient pas nous indiquer l'inclinaison très précisément, car l'événement était observé presque de face (comme regarder directement dans le canon d'une arme à feu). Il était difficile de distinguer la rotation « gauche-droite » de la rotation « haut-bas ».

Pour résoudre ce problème, les scientifiques devaient observer le « rayonnement fossile » — la lumière résiduelle qui s'estompe après le crash. Alors que les débris du crash s'éloignent à grande vitesse, ils créent un jet de particules. Si nous pouvons observer ce jet se déplacer dans le ciel au fil du temps, nous pouvons déterminer l'angle.

La nouvelle approche de l'article : une carte haute précision
Cet article, rédigé par Kelly Gourdji et ses collègues, revient à passer d'une carte de rue floue à une image satellite haute définition.

L'ancienne méthode : Les études précédentes tentaient de deviner la forme du jet en observant l'évolution de la luminosité au fil du temps (la « courbe de lumière »). Mais c'est comme essayer de deviner la forme d'une voiture en écoutant seulement le bruit de son moteur. C'est délicat car le bruit du moteur dépend à la fois de la forme de la voiture et de sa vitesse.
La nouvelle méthode : Les auteurs ont utilisé une technique appelée VLBI (Interférométrie à très longue base). Imaginez relier des radiotélescopes à travers toute la Terre pour créer un seul télescope géant de la taille de la planète. Cela leur offre une vue si nette qu'ils peuvent discerner des détails de la taille d'un milliarcseconde (ce qui revient à voir une pièce de monnaie sur la Lune depuis la Terre).

Ils n'ont pas simplement observé le « centre » du blob lumineux ; ils ont analysé la forme et le mouvement de l'ensemble des ondes radio directement, en utilisant un modèle informatique sophistiqué qui simule le comportement d'un jet de gaz lorsqu'il percute l'espace. Ils ont combiné des données provenant de quatre moments différents :

Jour 0 : Le moment exact du crash (repéré par le télescope spatial Hubble).
Jours 75 et 230 : Observations radio provenant d'un réseau de télescopes.
Jour 207 : Une observation radio haute résolution qui a servi de « point de contrôle » crucial.

Les résultats : un angle plus précis
En ajustant leurs modèles 3D complexes à ces données ultra-nettes, ils ont déterminé l'angle de vue avec une précision bien supérieure à celle d'avant.

L'angle : Ils ont constaté que le jet est incliné d'environ 18 à 20 degrés par rapport à notre ligne de visée. (Les estimations précédentes étaient un peu plus floues, allant de 14 à 28 degrés).
La distance : Parce qu'ils connaissaient mieux l'angle, ils ont pu calculer la distance plus précisément. Ils ont déterminé que le crash s'est produit à environ 44 millions d'années-lumière.
Le taux d'expansion ( $H_0$ ) : En utilisant cette nouvelle distance précise, ils ont calculé la constante de Hubble à 65,5 km/s/Mpc.

Pourquoi cela compte
Ce résultat est un « arbitre » dans le débat cosmique.

L'« univers primordial » (satellite Planck) indique un taux de 67,4.
L'« univers tardif » (équipe SH0ES) indique un taux de 73,0.
Cette nouvelle mesure indique 65,5.

Les auteurs notent que leur résultat est très proche du chiffre de l'« univers primordial » (dans une marge d'erreur infime) mais significativement plus bas que celui de l'« univers tardif ». Cela suggère que peut-être les mesures de l'« univers tardif » sont légèrement erronées, ou que notre compréhension de l'expansion de l'univers doit être repensée.

Les mises en garde
Les auteurs notent prudemment que leur résultat dépend de plusieurs facteurs :

Le modèle du jet : Ils supposent que le jet se comporte exactement comme leurs simulations informatiques. Si le jet réel est étrange ou possède une source d'énergie cachée (comme un résidu du crash qui brillerait plus tard que prévu), le calcul de l'angle pourrait changer.
Le mouvement local : Les galaxies ne se déplacent pas seulement avec l'expansion de l'univers ; elles ont aussi leurs propres « vitesses particulières » (dérive due à la gravité locale). Les auteurs ont dû estimer et moyenner ces dérives locales, ce qui ajoute une certaine incertitude.

En résumé
Cet article revient à prendre une photo floue d'une toupie en rotation et à utiliser une super-lentille pour voir exactement comment elle est inclinée. Ce faisant, ils ont mesuré la distance du crash avec plus de précision et ont découvert que l'univers s'étend à un taux qui correspond davantage aux mesures du « Big Bang » qu'aux mesures des « étoiles proches ». Cela ne résout pas complètement le mystère de la tension de Hubble, mais il fournit une preuve très solide et indépendante qui penche vers un taux d'expansion « plus lent ».

Résumé technique : Réexamen de GW170817 à l'échelle de la milliarcseconde

Énoncé du problème
La détection de la fusion d'étoiles à neutrons binaires GW170817 a offert la première opportunité de mesurer la constante de Hubble ( $H_0$ ) grâce à une « sirène standard ». Cependant, la précision de cette mesure est fondamentalement limitée par la forte covariance entre la distance de luminosité ( $D_L$ ) et l'angle d'inclinaison ( $\iota$ ) du système binaire. Alors que les données d'ondes gravitationnelles (OG) seules ne peuvent pas contraindre précisément l'inclinaison pour des événements de face comme GW170817, le jet relativiste lancé par la fusion offre une prise électromagnétique sur cette géométrie. Les tentatives précédentes de contraindre l'angle de vue ( $\theta_v$ ) en utilisant la courbe de lumière du sursaut ont été limitées par des dégénérescences géométriques entre l'angle de vue et l'angle d'ouverture du jet ( $\theta_{cp}$ ). À l'inverse, les observations par interférométrie à très longue base (VLBI) du mouvement propre superluminal du jet fournissent des contraintes orthogonales, mais ont historiquement reposé sur l'ajustement de modèles de jets discrets et semi-analytiques ou de sous-ensembles limités de l'espace des paramètres, négligeant souvent l'information spatiale complète (telle que la taille de la source) ou les incertitudes systématiques de manière rigoureuse.

Méthodologie
Cette étude présente un cadre complet d'ajustement de modèle dans le plan des visibilités bayésiennes conçu pour extraire des contraintes robustes sur la géométrie de vue de GW170817 et, par conséquent, sur $H_0$ . La méthodologie intègre les composants suivants :

Intégration des données : L'analyse incorpore les quatre ensembles de données astrométriques disponibles à l'échelle de la milliarcseconde :
- Deux époques VLBI du High Sensitivity Array (HSA) (jours 75 et 230 post-fusion) de Mooley et al. (2018).
- Une époque Global-VLBI (jour 207) de Ghirlanda et al. (2019), que les auteurs ont réduite indépendamment pour assurer la cohérence avec l'étalonnage HSA.
- Une mesure du télescope spatial Hubble (HST) de la position de la kilonova (jour 8), servant d'ancrage pour la fusion (jour 0).
Modélisation hydrodynamique : Au lieu de s'appuyer sur des modèles précalculés discrets, les auteurs emploient un formalisme d'échelle appliqué à un continuum de nouveaux modèles hydrodynamiques de sursauts de jets basés sur des simulations relativistes 2D (Govreen-Segal & Nakar 2023). Ces modèles couvrent quatre distributions d'énergie de jet (top-hat et lois de puissance avec des indices $b=3, 4, 6$ ). L'échelle permet d'ajuster les modèles pour des angles de vue ( $\theta_v$ ), des angles d'ouverture ( $\theta_{cp}$ ) et des temps d'observateur arbitraires, en utilisant l'invariance des modèles sous des rapports de paramètres spécifiques.
Inférence bayésienne : Le cadre opère directement sur les visibilités observées par interféromètre radio en utilisant le logiciel accéléré par GPU GALARIO et l'échantillonneur imbriqué Dynesty (via Bilby). Cette approche évite la perte d'information associée à la déconvolution dans le plan de l'image. Le modèle inclut des paramètres de nuisance pour traiter rigoureusement les incertitudes systématiques, notamment les décalages astrométriques ( $\delta_{ra}, \delta_{dec}$ ) et les incertitudes d'échelle de flux ( $\delta_F$ ) pour chaque époque.
Contraintes cosmologiques : Pour mesurer $H_0$ , le cadre incorpore les vraisemblances des OG (via un estimateur de densité de noyau de l'inclinaison et de la distance) et les corrections de vitesse particulière. Les auteurs considèrent 28 combinaisons de catalogues de groupes de galaxies et de cartes de reconstruction de vitesse particulière, en employant la moyenne de modèles bayésiens (BMA) pour marginaliser sur ces incertitudes systématiques.

Contributions clés

Espace de paramètres continu : L'étude va au-delà de l'ajustement de modèles discrets pour explorer une plage continue de géométries de jets ( $\theta_v$ et $\theta_{cp}$ ), évitant ainsi les biais introduits par la sélection de structures de jets spécifiques.
Ajustement dans le plan des visibilités : En ajustant directement aux visibilités plutôt qu'aux centroïdes dans le plan de l'image, la méthode utilise l'information spatiale complète des données, y compris les contraintes de taille de la source, ce qui est crucial pour briser les dégénérescences.
Traitement des incertitudes systématiques : L'inclusion de paramètres de nuisance pour les incertitudes d'étalonnage astrométrique et de flux, couplée à l'approche BMA pour les corrections de vitesse particulière, fournit un budget d'erreur plus robuste que les études précédentes qui supposaient souvent des corrections fixes.
Réduction indépendante des données Global-VLBI : Les auteurs ont effectué une réduction indépendante des données Global-VLBI pour assurer la cohérence du référentiel avec les données HSA, en identifiant et en corrigeant les décalages astrométriques par rapport aux réductions précédentes.

Résultats

Contraintes sur l'angle de vue :
- En supposant une distance de luminosité fixe de $D_L = 40,7$ Mpc (basée sur des mesures de fluctuations de brillance de surface), l'analyse donne un angle de vue de $\theta_v = 18,3^\circ - 20,3^\circ$ (intervalle de crédibilité à 68 %) lors de l'incorporation de priors sur la géométrie du jet informés par la courbe de lumière.
- En marginalisant sur $D_L$ et les corrections de vitesse particulière, l'angle de vue préféré se déplace légèrement vers $\theta_v = 16,8^\circ - 19,2^\circ$ .
Distance de luminosité : L'ajustement direct pour la distance donne $D_L = 44,0 \pm 1,6$ Mpc. Cela est environ $1\sigma$ plus élevé que la valeur de $40,7 \pm 2,4$ Mpc dérivée de la méthode des fluctuations de brillance de surface, soulignant la sensibilité des contraintes d'angle de vue à la distance supposée.
Constante de Hubble ( $H_0$ ) : La mesure moyennée par modèles bayésiens sur 28 scénarios de correction de vitesse particulière est $H_0 = 65,5 \pm 4,4$ $H_{0} = 65, 5 \pm 4, 4$ km s $^{-1}$ $^{- 1}$ Mpc $^{-1}$ $^{- 1}$ .
- Cette valeur est à moins de $0,5\sigma$ de la valeur de la Collaboration Planck pour l'Univers primordial ( $67,4 \pm 0,5$ km s $^{-1}$ Mpc $^{-1}$ ).
- Elle est inférieure de $1,7\sigma$ à la mesure SH0ES pour l'Univers tardif ( $73,04 \pm 1,04$ km s $^{-1}$ Mpc $^{-1}$ ).

Signification et revendications
L'article revendique que cette mesure indépendante de $H_0$ dans l'Univers tardif est plus cohérente avec le résultat Planck de l'Univers primordial qu'avec le résultat SH0ES. Les auteurs suggèrent que, si l'on prend cette découverte au pied de la lettre, elle pourrait remettre en question la nécessité de modèles de cosmologie non standard pour résoudre la tension de Hubble. Cependant, les auteurs restent modestes quant à cette revendication, notant explicitement plusieurs mises en garde :

Mauvaise spécification du modèle : L'évolution tardive de la courbe de lumière du sursaut ( $\gtrsim 350$ jours) n'est pas bien ajustée par les modèles synchrotrons standards, et la source de l'émission excédentaire reste inconnue. Bien que les auteurs soutiennent que leur analyse (axée sur des époques antérieures) est dominée par la géométrie plutôt que par des imprécisions du modèle, ils reconnaissent que des physiques non modélisées pourraient biaiser les résultats.
Alignement du jet : L'analyse suppose que le jet est lancé le long de l'axe du moment angulaire orbital. Si un désalignement significatif existe, il pourrait biaiser de manière non linéaire l'estimation de $H_0$ .
Dominance de la vitesse particulière : Les auteurs identifient les corrections de vitesse particulière comme la source dominante d'incertitude dans la mesure de $H_0$ (contribuant à environ 5 % d'incertitude fractionnelle), plutôt que la contrainte sur l'angle de vue elle-même.

L'étude conclut que, bien que les observations VLBI de haute précision des sursauts de jets restent un outil unique et puissant pour contraindre les angles d'inclinaison à court terme, parvenir à une résolution définitive de la tension de Hubble nécessitera un plus grand échantillon de sirènes brillantes et une amélioration de la modélisation de la physique des jets et des vitesses particulières.

Articles similaires