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

Technische Zusammenfassung: Ableitung der nuklearen Eigenschaften von Neutronensternen aus Gravitationswellen- und Gammastrahlenausbrüsten

Problemstellung
Jüngste Beobachtungen von langgestreckten Gammastrahlenausbrüsten (lGRBs), die mit Kilonova-Emissionen einhergehen (z. B. GRB 211211A, GRB 230307A), stellen das traditionelle Paradigma infrage, wonach lGRBs ausschließlich aus dem Kollaps massereicher Sterne resultieren. Falls diese Ereignisse aus Verschmelzungen von binären Neutronensternen (BNS) oder Neutronenstern–Schwarzen-Löchern (NSBH) stammen, stellt sich die kritische Frage: Welche physikalische Eigenschaft bestimmt die Dauer des Ausbruchs? Theoretische Modelle legen nahe, dass die Gesamtmasse der Verschmelzung die Lebensdauer des Verschmelzungsremnants bestimmt. Leichtere Verschmelzungen können langlebige Neutronenstern-Remnants (NS) bilden, die kurze Gammastrahlenausbrüche (sGRBs) antreiben, während schwerere Verschmelzungen kurzlebige Remnants erzeugen, die in Schwarze Löcher (BHs) kollabieren, was potenziell lGRBs über massive Akkretionsscheiben antreibt. Die spezifischen Massenschwellenwerte, die diese Ergebnisse voneinander trennen, sowie deren Beziehung zu grundlegenden nuklearen Eigenschaften, insbesondere der Tolman-Oppenheimer-Volkoff-Masse ($M_{\text{TOV}}$), bleiben jedoch beobachtungstechnisch unbeschränkt.

Methodik
Die Autoren entwickeln ein Bayessches Framework, um charakteristische Massenskalen abzuleiten, indem sie die Gravitationswellen-Verschmelzungsraten der LIGO-Virgo-Kagra-Kollaboration (LVK) mit den beobachteten Raten von sGRBs und kilonova-assoziierten lGRBs vergleichen.

1. Definition der charakteristischen Massen: Die Studie parametrisiert drei kritische Gesamtmassen-Schwellenwerte relativ zur TOV-Masse ($M_{\text{TOV}}$):

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$: Die Grenze zwischen sehr langlebigen und langlebigen Remnants.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$: Die Grenze zwischen langlebigen und kurzlebigen Remnants. Verschmelzungen mit $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$ werden hypothetisch sGRBs erzeugen.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$: Die Grenze zwischen kurzlebigen Remnants und promptem Kollaps. Verschmelzungen mit $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$ werden hypothetisch lGRBs erzeugen. Verschmelzungen mit $M_{\text{tot}} > M_{\text{PC}}$ unterliegen einem prompten Kollaps; lGRBs werden nur produziert, wenn das Massenverhältnis ($q$) eine ausreichende Scheibenbildung ermöglicht ($1,2 \le q \le 3$ für BNS, $q \le 3$ für NSBH).

2. Datenquellen:

   • GW-Raten: Der LVK GWTC-4 Katalog (Abac et al. 2025a) und das FullPop-4.0 Populationsmodell (Abac et al. 2025b) liefern lokale Verschmelzungsraten als Funktion der Komponenteneigenschaften.
   • GRB-Raten: Die intrinsischen Raten für sGRBs ($R_{\text{sGRB}}$) und das Ratenverhältnis $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ werden als Beobachtungsbeschränkungen mit signifikanten Unsicherheiten behandelt, die über Normalverteilungen modelliert werden.

3. Inferenz: Unter Verwendung des Bayes-Theorems inferieren die Autoren die Posterior-Verteilungen für den Parametersatz $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$. Die Likelihood-Funktion vergleicht die vorhergesagten GW-induzierten GRB-Raten (abgeleitet aus den Massenschwellenwerten) mit den beobachteten GRB-Raten. Priors für $M_{\text{TOV}}$ und den NS-Radius bei der TOV-Masse ($R_{\text{TOV}}$) werden aus Legred et al. (2021) gezogen, unter Berücksichtigung von Pulsar- und GW-Beobachtungen.

Wichtigste Ergebnisse

• Charakteristisches Massenverhältnis ($k_2$): Die Analyse liefert einen Medianwert von $k_2 = 1,36^{+0,08}_{-0,09}$ (68 % Konfidenzintervall), wobei $k_2 > 1,24$ bei 90 % Konfidenz gilt. Dies deutet darauf hin, dass der Übergang von langlebigen zu kurzlebigen Remnants bei etwa $1,36 M_{\text{TOV}}$ stattfindet.
• Implikation für die Lebensdauer der Remnants: Der hohe Wert von $k_2$ legt nahe, dass massive Neutronensterne nach der Verschmelzung über längere Zeiträume (Hunderte von Millisekunden) überleben können, bevor sie kollabieren, anstatt sofort zu kollabieren.
• Korrelationen: Die Studie identifiziert eine starke Korrelation zwischen $k_2$ und $k_3$ sowie eine Anti-Korrelation zwischen $k_2$ und $M_{\text{TOV}}$. Diese Anti-Korrelation ermöglicht es, die Beschränkungen auf das charakteristische Massenverhältnis $k_2$ direkt auf die Obergrenzen der TOV-Masse abzubilden.
• Robustheit: Die Inferenz von $k_2$ bleibt robust über einen weiten Bereich an angenommenen sGRB-Raten ($30$ bis $2000 \text{ Gpc}^{-3}\text{yr}^{-1}$) und Ratenverhältnissen ($\eta$ von $0,01$ bis $3$). Variationen in den Massenverhältnis-Schwellenwerten ($q_{\text{low}}$ und $q_{\text{high}}$) und die Einbeziehung oder der Ausschluss der $M_{\text{VL}}$-Grenze haben minimale Auswirkungen auf das inferierte $k_2$.

Bedeutung und Behauptungen
Die Arbeit etabliert eine neue, unabhängige Methode zur Beschränkung der Zustandsgleichung (EoS) von Neutronensternen, indem sie die GW-Verschmelzungspopulationen mit der Dauer von GRBs verknüpft. Durch die Identifizierung einer charakteristischen Massenskala ($1,36 M_{\text{TOV}}$), die sGRB- und lGRB-Progenitoren trennt, liefern die Autoren eine direkte beobachtbare Beschränkung auf die maximale Masse von Neutronensternen.

Die Autoren merken an, dass ihre Ergebnisse konsistent mit numerischen Relativitätssimulationen (z. B. Perna et al. 2025) sind, die einen Übergang zwischen $1,3$ und $1,4 M_{\text{TOV}}$ vorhersagen. Sie betonen, dass, während die aktuellen Ergebnisse gegenüber Beobachtungsunsicherheiten robust sind, zukünftige Verbesserungen der GW- und GRB-Ratenmessungen, insbesondere gemeinsame Detektionen von GRBs mit GW-Ereignissen, die Korrelation zwischen charakteristischen Massen und $M_{\text{TOV}}$ verschärfen werden, was potenziell die Unsicherheit auf der TOV-Masse auf $\sim 0,1 M_{\odot}$ reduzieren wird. Dieser Ansatz ergänzt bestehende Techniken und bietet einen Weg, engere Grenzen für die NS-Zustandsgleichung festzulegen.