The Very Late Time Afterglow of GW170817 Favors a Wobbling Jet

Résumé technique : La rémanence très tardive de GW170817 favorise un jet oscillant

Énoncé du problème
La fusion d'étoiles à neutrons binaires (BNS) GW170817 demeure le seul événement ayant fait l'objet d'une détection simultanée d'ondes gravitationnelles, d'un sursaut gamma (GRB 170817A) et d'une kilonova. Bien qu'un modèle de jet structuré collimaté ait réussi à expliquer la phase de montée de la rémanence et le mouvement apparent superluminal, il fait face à une incohérence critique dans la rémanence très tardive (observée pendant près d'une décennie). Les modèles standards prédisent qu'après la rupture du jet, un jet collimaté s'étale latéralement, entraînant un déclin temporel de $F_\nu \propto t^{-p}$ (où $p > 2$ est l'indice de puissance des électrons). Cependant, les observations de GW170817 révèlent un taux de déclin nettement plus faible ($\alpha \lesssim 1.8$). Résoudre cette divergence sans invoquer de composantes supplémentaires non observées (telles qu'une rémanence de kilonova ou des chocs rafraîchis) ou des hypothèses physiques extrêmes reste un défi.

Méthodologie
Motivés par des simulations magnétohydrodynamiques relativistes générales (GRMHD) de collapsars et de fusions étoile à neutrons-trou noir, qui montrent que les disques d'accrétion peuvent être inclinés en raison d'interactions avec du gaz en chute ou d'une éjection asymétrique de moment angulaire, les auteurs proposent que le jet de GW170817 est « oscillant ».

• Modèle théorique : Les auteurs approximent le jet oscillant comme une structure en forme d'anneau traînant une bande continue sur le ciel. Ils dérivent analytiquement l'évolution temporelle de la rémanence, en considérant trois étapes : une décélération initiale, une phase intermédiaire où l'arc visible s'allonge, et une phase finale où l'anneau complet se fond dans un jet large.
• Dérivation analytique : Au stade intermédiaire, la région émettrice visible augmente dans une dimension (longueur de l'arc) tandis que le cône de faisceau s'étend dans deux dimensions. Cela résulte en un indice de déclin unique de $F_\nu \propto t^{-3(2p-1)/7}$. Pour des valeurs typiques de $p$ (2 à 2,5), cela produit des indices de déclin compris entre 1,3 et 1,7, correspondant à la pente faible observée.
• Simulation numérique et analyse bayésienne : Les auteurs utilisent le code hydrodynamique jetsimpy pour modéliser la dynamique de l'onde de choc, ce qui permet de prendre en compte plus précisément l'étalement complexe dans une structure en anneau que les modèles semi-analytiques. Ils réalisent une analyse bayésienne sur l'ensemble des données multimessagers, incluant les courbes de lumière radio et X, les mesures de mouvement apparent superluminal, et les échantillons postérieurs des ondes gravitationnelles pour l'inclinaison et la distance de luminosité. Ils comparent deux modèles : un jet structuré à loi de puissance standard et le jet en forme d'anneau (oscillant) proposé.

Contributions et résultats clés

• Validation du modèle : Le modèle de jet en forme d'anneau reproduit avec succès le déclin tardif faible ($\alpha \approx 1,3-1,7$) tout en restant cohérent avec les contraintes de mouvement apparent superluminal, qui nécessitent une structure « d'aile » collimatée.
• Signification statistique : La comparaison bayésienne des modèles produit un facteur de Bayes logarithmique de $\log_{10} \text{BF} = 6,16$, correspondant à un niveau de signification de $4,8\sigma$. Cela favorise fortement le modèle de jet en forme d'anneau par rapport au jet structuré à loi de puissance standard.
• Estimation des paramètres : L'analyse infère un angle d'oscillation de $\theta_{wo} \approx 27^\circ$. Les paramètres d'ajustement optimal pour le jet en forme d'anneau incluent un demi-angle d'ouverture $\theta_c \approx 0,03$ rad et un angle d'observateur $\theta_{obs} \approx 0,73$ rad.
• Implications physiques : Le grand angle d'oscillation inféré suggère une inclinaison significative du disque d'accrétion, potentiellement résultant d'interactions disque-gaz en chute ou d'une éjection asymétrique de moment angulaire lors de la fusion BNS.

Signification et affirmations
L'article affirme que le modèle de jet oscillant fournit une explication auto-cohérente aux données multimessagers de GW170817 sans invoquer de composants supplémentaires ni de processus physiques non standards, adhérant au rasoir d'Occam comme extension minimale du modèle de jet structuré.

• Résolution du « manque de rupture de jet » : Les auteurs suggèrent que ce modèle pourrait résoudre le problème du « manque de rupture de jet » dans d'autres sursauts gamma présentant des pentes de déclin de 1,4 à 1,8, impliquant que les jets oscillants pourraient être courants.
• Impact cosmologique : Le modèle produit un angle de vue plus large comparé aux modèles de jet collimaté, résultant en une contrainte réduite sur la distance de luminosité ($d_L \approx 36,4$ Mpc). Cela conduit à une estimation de la constante de Hubble de $H_0 \approx 80,8$ km/s/Mpc, qui favorise légèrement les inférences des supernovae de type Ia, bien que les auteurs notent que cette valeur puisse souffrir d'erreurs systématiques si l'axe d'oscillation est mal aligné avec l'axe orbital.
• Inférence du taux de fusion : Une structure en forme d'anneau implique un angle solide plus grand ($4\pi\theta_c\theta_{wo}$) comparé à un jet collimaté ($\pi\theta_c^2$). Les auteurs soutiennent que cela pourrait résoudre la divergence entre les taux de fusion BNS inférés des ondes gravitationnelles et ceux inférés des sursauts gamma courts, suggérant que ces derniers pourraient être surestimés.
• Limites : Les auteurs déclarent explicitement que le jet en forme d'anneau est une approximation de « modèle jouet ». Ils reconnaissent que les structures d'écoulement réelles révélées par les simulations sont hautement asymétriques et complexes, et ils appellent au développement d'outils de modélisation de rémanence tridimensionnelle efficaces pour investiguer davantage ces signatures.