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

Technische Samenvatting: Zwart-gat-neutronenster-binairiteiten nabij de disruptielimiet van de neutronenster in het massaregieme van GW230529

Probleemstelling
De detectie van de zwaartekrachtgolfgebeurtenis GW230529 door de LIGO–Virgo–KAGRA (LVK) Collaboratie markeerde de eerste waarschijnlijke zwart-gat–neutronenster (BHNS) fusie met een aanzienlijk potentieel voor het produceren van elektromagnetische (EM) tegenhangers. Deze gebeurtenis betreft een primair zwart gat (BH) van 2,5–4,5 $M_\odot$ en een secundair compact object van 1,2–2,0 $M_\odot$, wat het in een massaregieme plaatst waar de neutronenster (NS) getidale disrupteerd kan worden in plaats van direct in het zwarte gat te vallen. Tidale disruptie is een vereiste voor het genereren van observeerbare EM-signalen, zoals kilonovae en compact-binair gamma-ray bursts (cbGRBs), die de vorming van een accretieschijf en de uitstoot van neutronenrijk materiaal vereisen.

Echter, het voorspellen of een specifieke BHNS-fusie zal resulteren in disruptie, is sterk afhankelijk van analytische passende formules, zoals het FNH18-model (Foucart et al. 2018), dat de massa van de post-merger remnant schijf schat. Deze modellen zijn primair gekalibreerd op simulaties met massaverhoudingen $q \gtrsim 3$ en spins die ver verwijderd zijn van de kritieke drempel die vereist is voor disruptie. Bijgevolg is er een gebrek aan hoogwaardige numerieke relativiteit (NR) data voor systemen nabij de disruptielimiet, met name voor de massa- en spinparameters die consistent zijn met GW230529. Deze kloof beperkt de mogelijkheid om de prevalentie van EM-signalen nauwkeurig te voorspellen en de neutronenster-toestandsvergelijking (EOS) vanuit EM-observaties te beperken.

Methodologie
De auteurs maken gebruik van de Spectral Einstein Code (SpEC) om hoog-resolutie numerieke simulaties van BHNS-fusies uit te voeren. De studie richt zich op een massaverhouding van $q=3$ en een neutronenster massa van $M_{NS} = 1.35 M_\odot$, parameters die consistent zijn met de GW230529-gebeurtenis. Twee verschillende toestandsvergelijkingen worden gebruikt om verschillende niveaus van neutronenster-compactheid te testen:

1. SFHo: Een zachtere EOS met een straal $R_{NS} \approx 11,9$ km.
2. DD2: Een stijvere EOS met een straal $R_{NS} \approx 13,2$ km.

De simulaties variëren de spin-magnitude van het zwarte gat ($\chi_{BH}$) om de regio direct rond de kritieke spin ($\chi_{crit}$) te verkennen die door het FNH18-model wordt voorspeld voor tidale disruptie. Voor SFHo is $\chi_{crit} \approx 0,037$; voor DD2 is $\chi_{crit} \approx -0,23$. De studie omvat simulaties met spins zowel iets onder als boven deze kritieke drempels.

Belangrijke methodologische kenmerken zijn:

• Monte Carlo Neutrino Transport: Dit markeert de eerste toepassing van Monte Carlo neutrino-radiatietransport in BHNS-fusiesimulaties. In tegen tegenstelling tot leakage-modellen, volgt deze methode neutrino-trajecten in positie- en momentumruimte om absorptie, emissie en verstrooiing nauwkeurig te modelleren. Dit is cruciaal voor het bepalen van de compositie (elektronenfractie, $Y_e$) van de post-merger remnant en outflows, wat de potentie voor r-proces nucleosynthese dicteert.
• Resolutiestudies: Convergentietests werden uitgevoerd voor de SFHo-geval bij drie resoluties om numerieke onzekerheden te kwantificeren.
• Initiële condities: De initiële data werden gegenereerd met de Spells solver, waarbij de orbitale eccentriciteit werd verminderd tot $e \lesssim 0,009$.

Belangrijkste Resultaten

1. Schijfmassa versus Analytische Voorspellingen:

   • SFHo EOS: De numerieke resultaten voor de totale remnant massa vallen over het algemeen binnen de verwachte modelleringsfouten van de FNH18-formule. Echter, de simulaties laten consistent iets hogere remnant massa's zien dan het model voorspelt.
   • DD2 EOS: Er wordt een significante discrepantie waargenomen. Voor spins nabij de disruptielimiet onderschat het FNH18-model de remnant schijfmassa met een factor die de opgegeven foutmargen van het model overschrijdt. Bijvoorbeeld, bij $\chi_{BH} = -0,030$ voorspelt het model een remnant massa van $\sim 0,012 M_\odot$, terwijl de simulatie $\sim 0,049 M_\odot$ oplevert.
   • Kritieke Spin Drempel: Ondanks de massa-discrepanties wordt de kritieke spinwaarde die nodig is om disruptie te triggeren (waar de remnant massa overgaat van nul naar niet-nul) gevonden als consistent tussen het analytische model en de simulaties. Het primaire verschil ligt in de helling van de remnant massa-curve ($dM_{rem}/d\chi_{BH}$) nabij de drempel; simulaties tonen een steilere stijging dan het model.

2. Opkomst van Laag-massa Schijven:

   • De simulaties onthullen de vorming van niet-nul accretieschijven zelfs voor spins die iets onder de voorspelde disruptielimiet liggen (bijv. de SFHo-geval met $\chi_{BH} = -0,0130$). Hoewel deze schijven een lage massa hebben ($\sim 10^{-3} M_\odot$), zijn ze potentieel voldoende om korte cbGRBs aan te drijven, zoals gesuggereerd door eerdere studies (bijv. Gottlieb et al.).

3. Massa-uitstoot en Compositie:

   • Dynamische Ejecta: De hoeveelheid massa die dynamisch wordt uitgestoten tijdens de fusie is zeer laag, doorgaans $\lesssim 0,004 M_\odot$. Voor de laagste spinconfiguraties ligt de ongebonden massa nabij de limieten van de numerieke resolutie, waarbij onzekerheden potentieel $\sim 100\%$ kunnen bereiken.
   • Schijfwinden: De studie sugggeert dat voor deze systemen de totale massa uitgestoten via latere-tijd schijfwinden (viskeus en magnetisch gedreven) waarschijnlijk zal domineren over de initiële dynamische ejecta.
   • Elektronenfractie ($Y_e$): De dynamische ejecta is zeer neutronenrijk ($Y_e < 0,1$). Echter, de outflows en het schijfmateriaal vertonen een breder bereik van $Y_e$ (tot $\sim 0,3$) door zwakke interacties (neutrino-verhitting). De DD2 EOS produceert meer neutronenrijke ejecta vergeleken met SFHo.

4. Observationele Implicaties:

   • Elk kilonova-signaal van deze specifieke GW230529-achtige fusies wordt voorspeld dim te zijn, waarschijnlijk gedomineerd door de optische emissie van schijfwinden in plaats van de infrarode emissie van dynamische ejecta.
   • De totale beschikbare ejecta-massa voor r-proces nucleosynthese is $\lesssim 0,02 M_\odot$.

Betekenis en Claims
Het artikel claimt dat hoewel het FNH18 analytische model breed genomen betrouwbaar blijft voor het voorspellen van de existentie van een disruptiedrempel (de kritieke spin), het systematisch de massa van de resulterende remnant schijf onderschat, met name voor stijvere toestandsvergelijkingen (DD2) nabij de disruptielimiet.

De auteurs benadrukken dat hun resultaten de grenzen van de huidige kalibratie belichten wanneer deze wordt toegepast op BHNS-fusies dicht bij de disruptiedrempel. Specifiek:

• De aanname van het model van een "nul" remnant massa onder de kritieke spin is niet strikt accuraat; schijven met een lage massa kunnen zelfs onder de voorspelde limiet worden gevormd.
• Deze laag-massa schijven zijn fysiek significant als potentiële producenten van korte gamma-ray bursts, zelfs als ze geen heldere kilonovae produceren.
• De bevindingen suggereren dat een bescheiden herkalibratie van de passende formule noodzakelijk kan zijn naarmate er meer simulaties in het laag-massaverhouding, nabij-kritieke-spin regime beschikbaar komen.

De studie claimt niet definitief de waarschijnlijkheidsbeoordeling van de LVK voor GW230529 als een disruptieve fusie te hebben gewijzigd, waarbij wordt opgemerkt dat de geobserveerde verschillen de resultaten van de LVK-analyse waarschijnlijk kwalitatief niet zullen veranderen. Het biedt echter een rigoureuzere numerieke basis voor het begrijpen van het gedrag van BHNS-binairiteiten in deze specifieke, observationeel cruciale parameterruimte.