main concepts scientists sections
SimplePopulaireAvancé

En rotation dans la brèche : l'effondrement direct comme source de GW231123 ✛

Original: "Spinning into the Gap: Direct-Horizon Collapse as the Origin of GW231123 from End-to-End GRMHD Simulations"
· Ore Gottlieb, Brian D. Metzger, Danat Issa, Sean E. Li, Mathieu Renzo, Maximiliano Isi
arXiv:2508.15887 · 2025-08-21 · CC BY 4.0 · 6 min
Phénomènes astrophysiques de haute énergie
Des étoiles massives avec un champ magnétique modéré produisent, lors de leur effondrement, des trous noirs dont les masses et les spins correspondent à GW231123.

Contexte

GW231123, détecté le 23 novembre 2023 par la collaboration LIGO-Virgo-KAGRA, est le système binaire de trous noirs le plus massif jamais observé : les masses des composants sont de 137 et 103 masses solaires, et les spins sans dimension sont de 0.9 et 0.8. Les deux trous noirs se situent dans la lacune de l'instabilité de paire (plage de masses de 69 à 139 masses solaires), où l'évolution stellaire standard prédit l'absence de trous noirs. Les fusions hiérarchiques dans les amas denses peuvent produire de telles masses, mais donnent rarement des spins élevés en raison de l'orientation aléatoire. L'évolution isolée des étoiles massives est confrontée à un problème : l'effondrement direct d'un cœur au-dessus de la lacune conduit généralement à un trou noir dont la masse est égale à celle du cœur, c'est-à-dire au-dessus de la lacune. Ce travail s'inscrit dans la lignée des études sur la formation des trous noirs massifs dans les systèmes binaires, mais contrairement aux travaux voisins qui considèrent des scénarios dynamiques ou une évolution sans champs magnétiques, il modélise ici pour la première fois tout le chemin, du cœur d'hélium au trou noir, en tenant compte de la rotation et des champs magnétiques.

Méthodes

Les auteurs ont commencé avec un modèle de cœur d'hélium de 250.75 masses solaires et une métallicité Z=0.001, dont l'évolution a été calculée avec le code MESA jusqu'au moment de l'effondrement, lorsque la masse du cœur a diminué à 149 masses solaires en raison des vents stellaires. L'effondrement jusqu'à la formation du trou noir a été modélisé avec le code hydrodynamique relativiste général unidimensionnel GR1D, avec l'équation d'état LS220 ; un trou noir de 40 masses solaires s'est formé 4.29 secondes après le début de l'effondrement. Ensuite, les profils de densité et de vitesse ont été transférés dans le code magnétohydrodynamique relativiste général tridimensionnel H-AMR, qui a fait évoluer le système après la formation du trou noir avec une masse fixe et un spin a=0.5. La rotation a été définie par un profil obtenu en moyennant 125 étoiles massives à évolution chimiquement homogène, et le champ magnétique a été initialisé sous la forme d'une configuration axisymétrique avec différentes intensités : de très faible (Bw) à très forte (Bvs), avec une intensité au rayon de 100 rayons gravitationnels allant de $10^{11}$.6 à $10^{15}$.5 G. La grille comportait jusqu'à trois niveaux de raffinement adaptatif, avec une résolution allant jusqu'à 6144×6144×8192 cellules pour le champ faible, ce qui permettait de résoudre l'instabilité magnéto-rotationnelle avec un paramètre de qualité Q>10.

Résultats

Les simulations ont montré que, immédiatement après la formation du trou noir, le taux d'accrétion atteint des valeurs extrêmes supérieures à 100 masses solaires par seconde, puis se stabilise à un niveau de 10 à 40 masses solaires par seconde jusqu'à ce que la masse principale de l'étoile soit épuisée, environ 1 seconde plus tard. Dans les modèles avec un champ magnétique fort et très fort (Bs et Bvs), le flux magnétique sans dimension atteint une saturation d'environ 50, ce qui indique un état de freinage magnétique et une extraction efficace de l'énergie de rotation du trou noir, générant des jets avec une puissance allant jusqu'à $10^{54}$ erg/s. Dans le modèle avec un champ faible (Bw), aucun jet ne se forme, et le trou noir accrète presque toute l'étoile, atteignant une masse d'environ 140 masses solaires et un spin maximal a≈1. Dans le modèle avec un champ modéré (Bm), les jets apparaissent tardivement, lorsque le taux d'accrétion diminue, et le trou noir se stabilise à une masse d'environ 100 masses solaires et un spin a≈0.8. Dans les modèles avec des champs forts, le spin diminue jusqu'à des valeurs d'équilibre, allant jusqu'à un spin rétrograde a≈-0.1 dans le modèle Bvs en raison d'une accrétion stochastique avec un moment cinétique opposé. La corrélation finale masse-spin montre que les champs faibles produisent des trous noirs massifs et en rotation rapide, les champs forts produisent des trous noirs moins massifs et en rotation lente, et les champs modérés se situent dans la lacune de masse avec un spin élevé, ce qui est cohérent avec les paramètres de GW231123 : le trou noir secondaire correspond au modèle Bm, le primaire au modèle Bw.

Implications

Ce travail démontre pour la première fois que l'effondrement direct d'une étoile massive en rotation et magnétisée peut naturellement produire des trous noirs dans la lacune de l'instabilité de paire avec des spins élevés, ce qui lève la principale objection contre ce canal de formation pour GW231123. Contrairement aux travaux voisins qui considèrent soit des scénarios dynamiques avec de faibles spins, soit des modèles semi-analytiques sans champs magnétiques, une approche cohérente de bout en bout, de l'étoile au trou noir, avec une magnétohydrodynamique tridimensionnelle, est appliquée ici. Cela établit une nouvelle corrélation entre la masse et le spin du trou noir en fonction de l'intensité du champ magnétique du progéniteur, ce qui peut servir de diagnostic pour les études de population. De plus, le modèle prédit des transitoires électromagnétiques puissants, comparables aux sursauts gamma les plus énergétiques, ce qui ouvre la possibilité d'observations multi-messagers.

Développements futurs

Les recherches futures devraient inclure l'évolution dynamique de la masse et du spin du trou noir dans les simulations, ainsi que le transfert radiatif et le refroidissement par neutrinos pour déterminer plus précisément les propriétés des jets et des éjections. Il est nécessaire d'étudier une gamme plus large de masses de progéniteurs et de profils de rotation afin de construire une image complète de la formation des trous noirs massifs. Il est également important de calculer la dynamique détaillée des kicks, qui peuvent conduire à une désorientation des spins et à une modification des orbites dans les systèmes binaires. Les recherches observationnelles d'analogues électromagnétiques de tels événements dans l'Univers primordial aideront à vérifier les prédictions du modèle.

Impact

Les résultats influenceront l'interprétation des données des observatoires d'ondes gravitationnelles, en particulier pour les événements massifs dans la lacune de masse, et les modèles de synthèse de population des trous noirs binaires. Ils stimuleront également la recherche de transitoires associés à l'effondrement d'étoiles massives dans les domaines optique et gamma.

Prochaines étapes

Il est prévu de réaliser des simulations avec une évolution auto-cohérente du spin du trou noir et d'inclure le calcul des kicks pour déterminer le destin du système binaire. L'influence de différentes équations d'état et configurations magnétiques initiales sera également étudiée.

Problèmes ouverts clés

Ce travail aborde directement le problème de l'origine des trous noirs dans la lacune de l'instabilité de paire, qui reste l'une des principales énigmes de l'astrophysique des objets compacts. Il relie également l'incertitude sur les profils de rotation et les champs magnétiques des étoiles massives aux propriétés finales des trous noirs, soulignant la nécessité d'améliorer les modèles d'évolution stellaire.

Le regard de la machine du savoir

Cette section s'adresse à l'auteur du travail et à ses collègues. Les travaux voisins ont été trouvés par une recherche sémantique dans notre archive, et non par correspondance de mots. Nous n'évaluons pas : nous montrons le champ autour.

Ce que fait ce travail

Les auteurs ont modélisé l'évolution d'un cœur d'hélium de 250,75 masses solaires, de l'effondrement à la formation d'un trou noir, à l'aide d'une magnétohydrodynamique relativiste générale tridimensionnelle, en incluant pour la première fois la rotation et les champs magnétiques. Ils ont montré que les vents de disque et les jets, régulés par le champ magnétique, permettent la formation de trous noirs dans la fenêtre d'instabilité de paire avec des spins élevés, reproduisant les masses et les spins de GW231123 et en accord avec GW190521.

Où est sa force

La modélisation de bout en bout, de l'évolution du cœur à l'effondrement et à l'accrétion, en tenant compte des champs magnétiques, constitue une avancée dans la compréhension de la formation des trous noirs massifs dans les systèmes isolés, offrant une explication naturelle des spins élevés.

Vers où avancer

À essayer

  • Vérifier si les spins élevés sont préservés en faisant varier le profil initial de rotation et le flux magnétique.
  • Estimer la fréquence d'un tel scénario dans la population d'étoiles massives en fonction de la métallicité.
  • Comparer les taux d'événements prédits avec ceux observés par LIGO-Virgo-KAGRA pour les masses dans la fenêtre.

Travaux voisins dans notre archive

Le travail le plus proche de l'archive se situe à 0.64 par le sens — le niveau habituel de notre corpus.

Le graphe entier →

Dans le modèle avec un champ magnétique modéré, le trou noir perd environ 40 masses solaires par les jets et les vents, ce qui est comparable à la masse d'une étoile entière. Les jets dans le modèle le plus puissant atteignent une puissance isotrope de $10^{56}$ erg/s, ce qui est des ordres de grandeur plus brillant que les sursauts gamma les plus énergétiques connus.

Le scénario d'effondrement direct d'une étoile massive avec formation d'un trou noir en rotation rapide et de jets puissants fait écho à la description de la formation des trous noirs dans le roman « Lumière des étoiles lointaines » de Greg Egan, où l'effondrement des étoiles de première génération engendre des objets extrêmes.

\ell(r,\theta) = \begin{cases} 0.55 r^2 \sin^2\theta \, \text{s}^{-1} & r < r_c \\ 0.55 r_c^2 \sin^2\theta \, \text{s}^{-1} & r_c < r < R_{\text{He}} \end{cases}
Profil de rotation : la partie interne tourne comme un corps solide, la partie externe conserve un moment cinétique constant.
a = \frac{c}{G M_{\text{BH}}^2} \int_0^{M_{\text{BH}}} \ell(m) \, dm
Le spin est calculé comme l'intégrale du moment cinétique spécifique sur la masse accrétée par le trou noir.
\dot{a}(t) = \dot{a}_{\text{acc}}(t) + \dot{a}_{\text{out}}(t)
La variation totale du spin est la somme de la contribution de l'accrétion (augmentation) et de l'extraction magnétique d'énergie (freinage).

Chiffres clés

  • masse du cœur d'hélium: 250.75 masses solaires
  • masse du trou noir lors de la formation: 40 masses solaires
  • taux d'accrétion juste après la formation: >100 masses solaires par seconde
  • spin du trou noir dans le modèle Bm: 0.8
  • masse du trou noir dans le modèle Bm: environ 100 masses solaires
Scientifiques
John Archibald WheelerJacob BekensteinFritz ZwickyKarl SchwarzschildStephen HawkingSubrahmanyan Chandrasekhar
Notions
trous noirs à effondrement direct trou noir binaire mécanisme de Hills GW190814 fusion de sursaut gamma combustion de l'hélium fusion de binaires d'étoiles à neutrons ringdown effondrement collision de bulles inflation de masse déficit de masse de réaction nucléaire inspiral étoile à hélium spin effectif événement de rupture de marée instabilité gravitationnelle horizon interne bulbe galactique m87