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

Technische Zusammenfassung: Schwarze-Loch-Neutronenstern-Binärsysteme nahe der Grenze der Neutronenstern-Disruption im Massenregime von GW230529

Problemstellung
Die Detektion des Gravitationswellenereignisses GW230529 durch die LIGO–Virgo–KAGRA (LVK) Kollaboration markierte das erste wahrscheinliche Schwarzes-Loch–Neutronenstern-Verschmelzung (BHNS) Ereignis mit einem signifikanten Potenzial für die Erzeugung elektromagnetischer (EM) Gegenstücke. Dieses Ereignis beinhaltet ein primäres Schwarzes Loch (BH) von 2,5–4,5 $M_\odot$ und ein sekundäres kompaktes Objekt von 1,2–2,0 $M_\odot$, was es in ein Massenregime platziert, in dem der Neutronenstern (NS) eher durch Gezeitenkräfte disruptiert wird, anstatt direkt in das Schwarze Loch zu stürzen. Die Gezeiten-Disruption ist eine Voraussetzung für die Erzeugung beobachtbarer EM-Signale, wie etwa Kilonovae und kompakte-Binär-Gamma-Ray-Bursts (cbGRBs), welche die Bildung einer Akkretionsscheibe und den Ausstoß neutronenreicher Materie erfordern.

Die Vorhersage, ob eine spezifische BHNS-Verschmelzung jedoch zu einer Disruption führt, hängt stark von analytischen Anpassungsformeln ab, wie etwa dem FNH18-Modell (Foucart et al. 2018), welches die Masse der post-merger Remnant-Scheibe schätzt. Diese Modelle wurden primär auf Simulationen mit Massenverhältnissen $q \gtrsim 3$ und Spins kalibriert, die weit von der kritischen Schwelle für eine Disruption entfernt liegen. Infolgedessen mangelt es an hochpräzisen numerischen Relativität (NR) Daten für Systeme nahe der Disruptionsgrenze, insbesondere für die Spin- und Massenparameter, die konsistent mit GW230529 sind. Diese Lücke schränkt die Fähigkeit ein, die Häufigkeit von EM-Signalen genau vorherzusagen und die Zustandsgleichung (EOS) des Neutronensterns aus EM-Beobachtungen einzugrenzen.

Methodik
Die Autoren nutzen den Spectral Einstein Code (SpEC), um hochauflösende numerische Simulationen von BHNS-Verschmelzungen durchzuführen. Die Studie konzentriert sich auf ein Massenverhältnis von $q=3$ und eine Neutronensternmasse von $M_{NS} = 1,35 M_\odot$, Parameter, die konsistent mit dem GW230529-Ereignis sind. Es werden zwei unterschiedliche Zustandsgleichungen verwendet, um verschiedene Grade der Neutronenstern-Kompaktheit zu testen:

1. SFHo: Eine weichere EOS mit einem Radius $R_{NS} \approx 11,9$ km.
2. DD2: Eine steifere EOS mit einem Radius $R_{NS} \approx 13,2$ km.

Die Simulationen variieren die Spin-Magnitude des Schwarzen Lochs ($\chi_{BH}$), um die Region unmittelbar um den kritischen Spin ($\chi_{crit}$) zu untersuchen, der vom FNH18-Modell für die Gezeiten-Disruption vorhergesagt wird. Für SFHo ist $\chi_{crit} \approx 0,037$; für DD2 ist $\chi_{crit} \approx -0,23$. Die Studie umfasst Simulationen mit Spins sowohl leicht unterhalb als auch oberhalb dieser kritischen Schwellenwerte.

Wesentliche methodische Merkmale sind:

• Monte Carlo Neutrino-Transport: Dies stellt die erste Anwendung von Monte-Carlo-Neutrino-Strahlungstransport in BHNS-Verschmelungssimulationen dar. Im Gegensatz zu Leakage-Modellen verfolgt diese Methode Neutrino-Trajektorien in Positions- und Impulsraum, um Absorption, Emission und Streuung genau zu modellieren. Dies ist entscheidend für die Bestimmung der Zusammensetzung (Elektronenfraktion, $Y_e$) des Post-Merger-Remnants und der Ausflüsse, was die potenzielle r-Prozess-Nukleosynthese bestimmt.
• Auflösungsstudien: Konvergenztests wurden für den SFHo-Fall bei drei Auflösungen durchgeführt, um numerische Unsicherheiten zu quantifizieren.
• Anfangsbedingungen: Die Anfangsdaten wurden mit dem Spells-Solver generiert, wobei die orbitale Exzentrizität auf $e \lesssim 0,009$ reduziert wurde.

Hauptergebnisse

1. Scheibenmasse vs. analytische Vorhersagen:

   • SFHo EOS: Die numerischen Ergebnisse für die gesamte Remnant-Masse fallen im Allgemeinen innerhalb der erwarteten Modellierungsfehler der FNH18-Formel. Dennoch zeigen die Simulationen konsistent etwas höhere Remnant-Massen, als das Modell vorhersagt.
   • DD2 EOS: Es wird eine signifikante Diskrepanz beobachtet. Bei Spins nahe der Disruptionsgrenze unterschätzt das FNH18-Modell die Remnant-Scheibenmasse um einen Faktor, der die im Modell angegebenen Fehlermargen übersteigt. Beispielsweise sagt das Modell bei $\chi_{BH} = -0,030$ eine Remnant-Masse von $\sim 0,012 M_\odot$ voraus, während die Simulation $\sim 0,049 M_\odot$ liefert.
   • Kritischer Spin-Schwellenwert: Trotz der Massendiskrepanzen ist der kritische Spin-Wert, der die Disruption auslöst (wo die Remnant-Masse von Null auf Nicht-Null übergeht), zwischen dem analytischen Modell und den Simulationen konsistent. Der primäre Unterschied liegt in der Steigung der Remnant-Masse-Kurve ($dM_{rem}/d\chi_{BH}$) nahe der Schwelle; Simulationen zeigen einen steileren Anstieg als das Modell.

2. Entstehung von Niedrigmasse-Scheiben:

   • Die Simulationen zeigen die Bildung von nicht-null Akkretionsscheiben selbst für Spins, die leicht unter der vorhergesagten Disruptionsgrenze liegen (z. B. der SFHo-Fall mit $\chi_{BH} = -0,0130$). Obwohl diese Scheiben eine geringe Masse aufweisen ($\sim 10^{-3} M_\odot$), sind sie potenziell ausreichend, um kurze cbGRBs anzutreiben, wie bereits in früheren Studien (z. B. Gottlieb et al.) suggeriert wurde.

3. Masse-Ausstoß und Zusammensetzung:

   • Dynamischer Ejecta: Die Menge der während der Verschmelzung dynamisch ausgestoßenen Masse ist sehr gering, typischerweise $\lesssim 0,004 M_\odot$. Für die niedrigsten Spin-Konfigurationen liegt die ungebundene Masse nahe an den Grenzen der numerischen Auflösung, wobei die Unsicherheiten potenziell $\sim 100\%$ erreichen können.
   • Scheibenwinde (Disk Winds): Die Studie legt nahe, dass für diese Systeme die durch spätere Zeiträume (viskose und magnetisch getriebene) Scheibenwinde ausgestoßene Gesamtmasse wahrscheinlich den anfänglichen dynamischen Ejecta dominieren wird.
   • Elektronenfraktion ($Y_e$): Der dynamische Ejecta ist hochgradig neutronenreich ($Y_e < 0,1$). Die Ausflüsse und das Scheibenmaterial zeigen jedoch einen breiteren Bereich von $Y_e$ (bis zu $\sim 0,3$) aufgrund schwacher Wechselwirkungen (Neutrino-Heating). Die DD2 EOS produziert einen neutronenreicheren Ejecta im Vergleich zu SFHo.

4. Beobachtungsimplikationen:

   • Jegliches Kilonova-Signal von diesen spezifischen GW230529-ähnlichen Verschmelzungen wird als schwach vorhergesagt und wird wahrscheinlich von der optischen Emission der Scheibenwinde dominiert, statt von der Infrarot-Emission des dynamischen Ejecta.
   • Die zur r-Prozess-Nukleosynthese verfügbare Gesamt-Ejecta-Masse beträgt $\lesssim 0,02 M_\odot$.

Bedeutung und Behauptungen
Das Paper behauptet, dass das analytische FNH18-Modell zwar im Großen und Ganzen zuverlässig bleibt, um die Existenz einer Disruption-Schwelle (den kritischen Spin) vorherzusagen, aber die resultierende Remnant-Scheibenmasse systematisch unterschätzt, insbesondere bei steiferen Zustandsgleichungen (DD2) nahe der Disruptionsgrenze.

Die Autoren betonen, dass ihre Ergebnisse die Grenzen der aktuellen Kalibrierung verdeutlichen, wenn diese auf BHNS-Verschmelzungen nahe der Disruptionsgrenze angewendet wird. Speziell:

• Die Annahme des Modells einer "Null"-Remnant-Masse unterhalb des kritischen Spins ist nicht streng korrekt; auch unter dem vorhergesagten Limit können Niedrigmasse-Scheiben entstehen.
• Diese Niedrigmasse-Scheiben sind physikalisch signifikant als potenzielle Progenitoren für kurze Gamma-Ray-Bursts, selbst wenn sie keine hellen Kilonovae erzeugen.
• Die Ergebnisse legen nahe, dass eine moderate Rekalibrierung der Anpassungsformel notwendig sein könnte, sobald mehr Simulationen im Bereich niedriger Massenverhältnisse und nahe der kritischen Spin-Schwelle verfügbar sind.

Die Studie beansprucht nicht, die Wahrscheinlichkeitseinschätzung der LVK bezüglich der Disruption eines GW230529-ähnlichen Mergers definitiv zu ändern, und stellt fest, dass die beobachteten Unterschiede die Ergebnisse der LVK-Analyse qualitativ unwahrscheinlich verändern werden. Sie bietet jedoch eine fundiertere numerische Basis für das Verständnis des Verhaltens von BHNS-Binärsystemen in diesem spezifischen, beobachtungstechnisch kritischen Parameterraum.