Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

Résumé Technique : Inférence des propriétés nucléaires des étoiles à neutrons à partir des observations d'ondes gravitationnelles et de sursauts gamma

Énoncé du problème
Des observations récentes de sursauts gamma à longue durée (lGRB) accompagnés d'une émission de kilonova (ex. GRB 211211A, GRB 230307A) remettent en question le paradigme traditionnel selon lequel les lGRB proviennent uniquement de l'effondrement d'étoiles massives. Si ces événements sont originaires de fusions d'étoiles à neutrons binaires (BNS) ou d'étoiles à neutrons-trous noirs (NSBH), une question critique se pose : quelle propriété physique régit la durée du sursaut ? Les modèles théoriques suggèrent que la masse totale de la fusion détermine la durée de vie du rémanent de la fusion. Les fusions plus légères peuvent former des rémanents d'étoiles à neutrons (NS) à longue durée de vie qui alimentent des sursauts gamma courts (sGRB), tandis que les fusions plus massives entraînent des rémanents à courte durée de vie qui s'effondrent en trous noirs (BH), pouvant potentiellement alimenter des lGRB via des disques d'accrétion massifs. Cependant, les seuils de masse spécifiques séparant ces résultats et leur relation avec les propriétés nucléaires fondamentales, spécifiquement la masse de Tolman-Oppenheimer-Volkoff ($M_{\text{TOV}}$), restent sans contraintes observationnelles.

Méthodologie
Les auteurs développent un cadre bayésien pour inférer des échelles de masse caractéristiques en comparant les taux de fusion par ondes gravitationnelles (GW) de la collaboration LIGO-Virgo-Kagra (LVK) avec les taux observés de sGRB et de lGRB associés à des kilonovae.

1. Définitions des masses caractéristiques : L'étude paramètre trois seuils de masse critiques relatifs à la masse TOV ($M_{\text{TOV}}$) :

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$ : la limite entre les rémanents très longue durée et longue durée.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$ : la limite entre les rémanents longue durée et courte durée. Les fusions dont $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$ sont supposées produire des sGRB.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$ : la limite entre les rémanents courte durée et l'effondrement immédiat (prompt collapse). Les fusions avec $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$ sont supposées produire des lGRB. Les lGRB ne sont produits que si le rapport de masse ($q$) permet une formation de disque suffisante ($1,2 \le q \le 3$ pour BNS, $q \le 3$ pour NSBH) ; sinon, les fusions avec $M_{\text{tot}} > M_{\text{PC}}$ subissent un effondrement immédiat.

2. Sources de données :

   • Taux de GW : Le catalogue GWTC-4 (Abac et al. 2025a) et le modèle de population FullPop-4.0 (Abac et al. 2025b) fournissent les taux de fusion locaux en fonction des masses des composantes.
   • Taux de GRB : Les taux intrinsèques pour les sGRB ($R_{\text{sGRB}}$) et le rapport de taux $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ sont traités comme des contraintes observationnelles présentant des incertitudes significatives, modélisées via des distributions normales.

3. Inférence : En utilisant le théorème de Bayes, les auteurs infèrent les distributions postérieures pour l'ensemble de paramètres $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$. La fonction de vraisemblance compare les taux de GRB induits par les GW (dérivés des seuils de masse) aux taux de GRB observés. Les priors pour $M_{\text{TOV}}$ et le rayon de l'étoile à neutrons à la masse TOV ($R_{\text{TOV}}$) sont tirés de Legred et al. (2021), conditionnés par les observations de pulsars et d'ondes gravitationnelles.

Résultats clés

• Rapport de masse caractéristique ($k_2$) : L'analyse produit une valeur médiane de $k_2 = 1,36^{+0,08}_{-0,09}$ (intervalle de confiance de 68 %), avec $k_2 > 1,24$ à un niveau de confiance de 90 %. Cela indique que la transition entre les rémanents de longue durée et de courte durée se produit à environ $1,36 M_{\text{TOV}}$.
• Implication pour la durée de vie des rémanents : La valeur élevée de $k_2$ suggère que les étoiles à neutrons massives peuvent survivre pendant des périodes prolongées (plusieurs centaines de millisecondes) après la fusion avant de s'effondrer, plutôt que de s'effondrer immédiatement.
• Corrélations : L'étude identifie une forte corrélation entre $k_2$ et $k_3$, ainsi qu'une anti-corrélation entre $k_2$ et $M_{\text{TOV}}$. Cette anti-corrélation permet de mapper directement les contraintes sur le rapport de masse caractéristique $k_2$ vers des limites supérieures sur la masse TOV.
• Robustesse : L'inférence de $k_2$ reste robuste à travers une large gamme de taux de sGRB supposés ($30$à$2000 \text{ Gpc}^{-3}\text{an}^{-1}$) et de rapports de taux $\eta$ (de $0,01$à$3$). Les variations des seuils de rapport de masse ($q_{\text{low}}$ et $q_{\text{high}}$) et l'inclusion ou l'exclusion de la limite $M_{\text{VL}}$ ont un impact minimal sur le $k_2$ inféré.

Signification et affirmations
Cet article établit une nouvelle méthode indépendante pour contraindre l'équation d'état (EoS) des étoiles à neutrons en liant les populations de fusions d'ondes gravitationnelles aux durées des sursauts gamma. En identifiant une échelle de masse caractéristique ($1,36 M_{\text{TOV}}$) qui sépare les progéniteurs de sGRB et de lGRB, les auteurs fournissent une contrainte observationnelle directe sur la masse maximale des étoiles à neutrons.

Les auteurs notent que leurs résultats sont cohérents avec les simulations de relativité numérique (ex. Perna et al. 2025) qui prédisent une transition entre $1,3$ et $1,4 M_{\text{TOV}}$. Ils soulignent que, bien que les résultats actuels soient robustes face aux incertitudes observationnelles, les améliorations futures des mesures de taux de GW et de GRB, particulièrement les détections conjointes de GRB avec des événements GW, permettront de resserrer la corrélation entre les masses caractéristiques et $M_{\text{TOV}}$, réduisant potentiellement l'incertitude sur la masse TOV à $\sim 0,1 M_{\odot}$. Cette approche complète les techniques existantes et offre une voie pour placer des limites plus strictes sur l'EoS des étoiles à neutrons.