Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

Résumé technique : Binaires trou noir-étoile à neutrons proches de la limite de rupture de l'étoile à neutrons dans le régime de masse de l'événement GW230529

Énoncé du problème
La détection de l'événement d'ondes gravitationnelles GW230529 par la collaboration LIGO–Virgo–KAGRA (LVK) a marqué le premier cas probable de fusion trou noir–étoile à neutrons (BHNS) avec un potentiel significatif de produire des contreparties électromagnétiques (EM). Cet événement implique un trou noir primaire de 2,5–4,5 $M_\odot$ et un objet compact secondaire de 1,2–2,0 $M_\odot$, ce qui le place dans un régime de masse où l'étoile à neutrons (NS) peut être soumise à une rupture de marée plutôt que de plonger directement dans le trou noir (BH). La rupture de marée est une condition préalable à la génération de signaux EM observables, tels que les kilonovae et les sursauts gamma de binaires compactes (cbGRB), qui nécessitent la formation d'un disque d'accrétion et l'éjection de matière riche en neutrons.

Cependant, prédire si une fusion BHNS spécifique entraînera une rupture repose fortement sur des formules d'ajustement analytiques, telles que le modèle FNH18 (Foucart et al. 2018), qui estime la masse du disque résiduel post-fusion. Ces modèles ont été principalement calibrés sur des simulations avec des rapports de masse $q \gtrsim 3$ et des spins éloignés du seuil critique requis pour la rupture. Par conséquent, il existe un manque de données de haute fidélité issues de la relativité numérique (NR) pour les systèmes proches de la limite de rupture, particulièrement pour les paramètres de masse et de spin cohérents avec GW230529. Cette lacune limite la capacité de prédire avec précision la prévalence des signaux EM et de contraindre l'équation d'état (EOS) de l'étoile à neutrons à partir des observations EM.

Méthodologie
Les auteurs utilisent le code Spectral Einstein (SpEC) pour effectuer des simulations numériques à haute résolution de fusions BHNS. L'étude se concentre sur un rapport de masse de $q=3$ et une masse d'étoile à neutrons de $M_{NS} = 1,35 M_\odot$, des paramètres cohérents avec l'événement GW230529. Deux équations d'état distinctes sont employées pour tester différents niveaux de compacité de l'étoile à neutrons :

1. SFHo : Une EOS plus molle avec un rayon $R_{NS} \approx 11,9$ km.
2. DD2 : Une EOS plus raide avec un rayon $R_{NS} \approx 13,2$ km.

Les simulations font varier l'amplitude du spin du trou noir ($\chi_{BH}$) pour sonder la région immédiatement entourant le spin critique ($\chi_{crit}$) prédit par le modèle FNH18 pour la rupture de marée. Pour SFHo, $\chi_{crit} \approx 0,037$ ; pour DD2, $\chi_{crit} \approx -0,23$. L'étude inclut des simulations avec des spins légèrement en dessous et au-dessus de ces seuils critiques.

Caractéristiques méthodologiques clés :

• Transport de neutrinos de Monte Carlo : Cela marque la première application du transport de rayonnement de neutrinos de Monte Carlo dans les simulations de fusion BHNS. Contra o aux modèles de fuite (leakage models), cette méthode suit les trajectoires des neutrinos en position et en quantité de mouvement pour modéliser avec précision l'absorption, l'émission et la diffusion. Ceci est crucial pour déterminer la composition (fraction électronique, $Y_e$) du résidu post-fusion et des éjectas, ce qui dicte le potentiel pour la nucléosynthèse r.
• Études de résolution : Des tests de convergence ont été effectués pour le cas SFHo à trois résolutions afin de quantifier les incertitudes numériques.
• Conditions initiales : Les données initiales ont été générées à l'aide du solveur Spells, avec des excentricités orbitales réduites à $e \lesssim 0,009$.

Résultats clés

1. Masse du disque vs prédictions analytiques :

   • EOS SFHo : Les résultats numériques pour la masse totale du résidu tombent généralement dans les erreurs de modélisation attendues de la formule FNH18. Cependant, les simulations montrent systématiquement des masses résiduelles légèrement supérieures à ce que prédit le modèle.
   • EOS DD2 : Une divergence significative est observée. Pour les spins proches de la limite de rupture, le modèle FNH18 sous-estime la masse du disque résiduel par un facteur qui dépasse les marges d'erreur déclarées du modèle. Par exemple, à $\chi_{BH} = -0,030$, le modèle prédit une masse résiduelle de $\sim 0,012 M_\odot$, tandis que la simulation produit $\sim 0,049 M_\odot$.
   • Seuil de spin critique : Malgré les divergences de masse, la valeur du spin critique requis pour déclencher la rupture (où la masse résiduelle passe de zéro à non-zéro) est jugée cohérente entre le modèle analytique et les simulations. La principale différence réside dans la pente de la courbe de masse résiduelle ($dM_{rem}/d\chi_{BH}$) près du seuil ; les simulations montrent une montée plus abrupte que le modèle.

2. Émergence de disques de faible masse :

   • Les simulations révèlent la formation de disques d'accrétion non nuls même pour des spins légèrement en dessous de la limite de rupture prédite (par exemple, le cas SFHo avec $\chi_{BH} = -0,0130$). Bien que ces disques soient de faible masse ($\sim 10^{-3} M_\odot$), ils sont potentiellement suffisants pour alimenter de courts cbGRB, comme suggéré par des études antérieures (ex. Gottlieb et al.).

3. Éjection de masse et composition :

   • Éjectas dynamiques : La quantité de masse éjectée dynamiquement lors de la fusion est très faible, typiquement $\lesssim 0,004 M_\odot$. Pour les configurations de spin les plus basses, la masse non liée est proche des limites de la résolution numérique, avec des incertitudes pouvant atteindre $\sim 100\%$.
   • Vents de disque : L'étude suggère que pour ces systèmes, la masse totale éjectée via les vents de disque tardifs (visqueux et magnétiquement pilotés) sera probablement dominante par rapport aux éjectas dynamiques initiaux.
   • Fraction électronique ($Y_e$) : Les éjectas dynamiques sont hautement riches en neutrons ($Y_e < 0,1$). Cependant, les flux sortants et le matériau du disque présentent une plage de $Y_e$ plus large (jusqu'à $\sim 0,3$) en raison des interactions faibles (chauffage par neutrinos). L'EOS DD2 produit des éjectas plus riches en neutrons que l'EOS SFHo.

4. Implications observationnelles :

   • Tout signal de kilonova provenant de ces fusions spécifiques de type GW230529 est prédit comme étant faible, probablement dominé par l'émission optique des vents de disque plutôt que par l'émission infrarouge des éjectas dynamiques.
   • La masse totale d'éjectas disponible pour la nucléosynthèse r est $\lesssim 0,02 M_\odot$.

Signification et affirmations
L'article affirme que bien que le modèle analytique FNH18 reste globalement fiable pour prédire l'existence d'un seuil de rupture (le spin critique), il sous-estime systématiquement la masse du disque résiduel résultant, particulièrement pour les équations d'état plus rigides (DD2) près de la limite de rupture.

Les auteurs soulignent que leurs résultats mettent en évidence les limites de la calibration actuelle lorsqu'elle est appliquée aux fusions BHNS proches du seuil de rupture. Plus précisément :

• L'hypothèse du modèle d'une masse résiduelle « zéro » en dessous du spin critique n'est pas strictement exacte ; des disques de faible masse peuvent se former même en dessous de la limite prédite.
• Ces disques de faible masse sont physiquement significatifs en tant que progéniteurs potentiels pour les sursauts gamma courts, même s'ils ne produisent pas de kilonovae brillants.
• Les conclusions suggèrent qu'un recalibrage modeste de la formule de l'ajustement pourrait être nécessaire à mesure que de nouvelles simulations dans le régime de faible rapport de masse et de spin proche du seuil deviennent disponibles.

L'étude ne prétend pas modifier de manière définitive l'évaluation de la probabilité par la LVK concernant une fusion par rupture pour GW230529, notant que les différences observées sont peu susceptibles de changer qualitativement les résultats de l'analyse de la LVK. Elle fournit toutefois une base numérique plus rigoureuse pour comprendre le comportement des binaires BHNS dans cet espace de paramètres spécifiquement critique pour l'observation.