Black hole-neutron star binaries with high spins and large mass asymmetries: III. Properties of the ejected material and its electromagnetic signatures

Technische Samenvatting: Zwart gat–neutronenster-binair met hoge spins en grote massa-asymmetrieën: III. Eigenschappen van het uitgestoten materiaal en de elektromagnetische signaturen

Probleem en Context
Zwart gat–neutronenster (BHNS) binaire systemen zijn primaire kandidaten voor multi-messenger astronomie, omdat ze signalen bieden in zowel gravitatiegolf- (GW) als elektromagnetische (EM) kanalen. Hoewel de GW170817-gebeurtenis ongekende beperkingen oplegde aan binaire neutronenster-fusies (BNS), blijven bevestigde EM-tegenhangers voor BHNS-fusies onbevestigd, ondanks verschillende GW-detecties (bijv. GW200105, GW200115). Het begrijpen van de eigenschappen van dynamische ejecta in BHNS-systemen is cruciaal voor het vaststellen van de toestandsvergelijking (EOS) van dicht nucleair materiaal, het verifiëren van relativistische jet-lanceermechanismen en het evalueren van de bijdrage van BHNS-fusies aan de kosmische abundantie van zware elementen via het $r$-proces.

De bestaande literatuur leunt vaak op geïdealiseerde ejecta-distributies afgeleid van analytische modellen of vereenvoudigde numerieke resultaten. Dit artikel, het derde in een reeks, pakt dit aan door een ab-initio benadering te hanteren. Het onderzoekt de geometrische, thermodynamische en kinematische eigenschappen van puur dynamische ejecta uit algemene-relativistische magnetohydrodynamische (GRMHD) simulaties van BHNS-systemen die gekenmerkt worden door hoge zwarte gat (BH) spins en grote massaverhoudingen ($Q = M_{BH}/M_{NS}$).

Methodologie
De studie maakt gebruik van een uitgebreide parameterruimte die zes BHNS-configuraties beslaat:

1. Vaste massaverhouding $Q=4$ met variërende BH spins $\chi_{BH} \in \{0.4, 0.6, 0.8\}$.
2. Vaste BH spin $\chi_{BH}=0.8$ met variërende massaverhoudingen $Q \in \{4, 5, 6, 7\}$.

De simulaties maken gebruik van de DD2 EOS voor de neutronenster (NS) met een baryonische massa van $1.53 M_\odot$. De numerieke infrastructuur omvat de EinsteinToolkit met de Antelope code voor metriek-evolutie (gebruikmakend van CCZ4 of Z4c formuleringen) en de Frankfurt-IllinoisGRMHD (FIL) code voor het gemagnetiseerde fluïdum. Initiële data wordt gegenereerd met de FUKA solver.

De analyse verloopt in drie fasen:

1. Karakterisering van Dynamische Ejecta: De auteurs analyseren de energie- en impulsoverdracht naar het gedisrumpeerde materiaal, waarbij ze specifiek kijken naar de fase-ruimte distributie van specifieke energie ($\tilde{E}$) en genormaliseerd specifiek impulsmoment ($\tilde{J}$). Ze definiëren dynamische ejecta met behulp van zowel geodesische ($u_t < -1$) als Bernoulli ($h u_t < -1$) ongebondenheidscriteria.
2. Nucleosynthese: Kunstmatige tracer-deeltjes worden geëxtraheerd uit de simulatiedata op een detectorstraal van $400 M_\odot$. Deze tracers, gewogen naar massa, dienen als input voor het SkyNet nucleaire reactienetwerk om de elementen-abundanties te berekenen die voortkomen uit het $r$-proces. Het netwerk draait gedurende $\sim 32$ jaar om de evolutie van zware kernen te vangen.
3. Elektromagnetische Signaturen: De resulterende ejecta-eigenschappen (massa, snelheid, entropie, geometrie) worden ingevoerd in de radiatieve-overdrachtscode POSSIS om kilonova (KN) lichtcurves te generen. Deze worden vergeleken met de AT2017gfo data en de bovengrenzen van S190814bv.

Belangrijkste Bijdragen en Resultaten

• Mechanisme van Massa-ejectie: De studie beschrijft het mechanisme van dynamische massa-ejectie voor hoge-massaverhouding, hoge-spin binaire systemen. In tegen afwezigheid van lagere-massaverhouding gevallen, vertonen hoge-$Q$ binaire systemen met hoge spins een vertraagde vorming van een accretieschijf. De getijdenstaart voert meer rotaties rond het BH uit voordat de disruptie voltooid is, wat de overdracht van impulsmoment faciliteert en de accumulatie van gebonden materiaal vertraagt. In het $Q=7, \chi_{BH}=0.8$ geval wordt de ISCO naar buiten geduwd, en wordt materie dichter dan $\approx 120$ km binnen 5 ms na de fusie uitgeput.

• Bulk-eigenschappen van Ejecta:

  • Massa: Sterke getijden-disruptie (STD) gevallen (bijv. $Q=4, \chi_{BH}=0.8$) ejecteren een substantiële massa ($\sim 0.04 - 0.06 M_\odot$), die afneemt naarmate de massaverhouding groter wordt. De $Q=4, \chi_{BH}=0.8$ en $Q=5, \chi_{BH}=0.8$ configuraties vertegenwoordigen een plateau van de maximale ejecta-massa voor de gegeven spin.
  • Thermodynamica: Sterkere getijden-disrupties leiden tot een hogere gemiddelde entropie door frequentere botsingen in de getijdenstaart. Entropiewaarden kunnen de $50 k_B$/baryon overschrijden in de staarten van de distributie.
  • Kinematica: Een toenemende massa-asymmetrie verhoogt over het algemeen de gemiddelde snelheid van de ejecta. De kinetische energie wordt versterkt door het impulsmoment van het systeem. De gemiddelde snelheid $\langle v \rangle$ varieert van $\sim 0.18c$ tot $0.29c$, met een duidelijke trend van toenemende snelheid bij een hogere $Q$.
  • Compositie: De elektronfractie ($Y_e$) blijft laag ($\approx 0.04 - 0.06$) over alle simulaties, wat consistent is met de afwezigheid van neutrino-stralings-effecten in de simulaties. Deze lage $Y_e$ drijft de productie van zware elementen aan.

• Ruimtelijke en Angulaire Distributie:

  • De ejecta vertoont een significante anisotropie. Sterke getijden-disrupties resulteren in meer gekollimeerde ejecta met een polaire halve openingshoek van $\sim 5^\circ$, terwijl zwakke disrupties breder zijn ($\sim 10^\circ$).
  • De azimutale omvang varieert. Sterke disrupties beslaan $\sim 120^\circ$, terwijl zwakke disrupties en hoge-$Q$ gevallen de $360^\circ$ kunnen benaderen door het wikkelen van de getijdenstaart.
  • Het 90%-confinement gebied van de ejecta beslaat tussen 4% en 8% van de hemel, afhankelijk van de sterkte van de disruptie.

• Nucleosynthese-opbrengsten:

  • De $r$-proces opbrengsten worden gedomineerd door zeer neutronenrijk materiaal ($Y_e \approx 0.05$).
  • De resulterende abundantiepatronen tonen een klassieke signatuur van neutronenrijke dynamische ejecta: de tweede en derde $r$-proces pieken worden met een orde van grootte meer geproduceerd vergeleken met de zonnere $r$-proces data, terwijl de eerste piek ernstig is uitgeput.
  • Variaties in de abundantiepatronen tussen simulaties zijn het meest uitgesproken tussen de tweede en derde pieken ($140 < A < 190$) en in de actiniden-regio, wat correleert met verschillen in entropie en expansietijdsschalen.

• Elektromagnetische Signaturen:

  • Kilonova lichtcurves berekend met POSSIS laten zien dat de bolometrische en optische/nabij-infrarode evolutie afhankelijk is van de kijkhoek en de binaire parameters ($Q, \chi_{BH}$).
  • De modellen zijn over het algemeen consistent met de AT2017gfo data en de bovengrenzen voor S190814bv, wat suggereert dat BHNS-fusies met deze parameters waarneembare transiënten kunnen produceren, hoewel ze mogelijk minder helder en langduriger zijn dan hun BNS-tegenhangers.

Significantie
Het artikel claimt significantie door een gedetailleerde, ab-initio link te leggen tussen de onderliggende binaire parameters (massaverhouding en BH spin) en de observeerbare eigenschappen van de ejecta en de kilonova. Door verder te gaan dan geïdealiseerde aannames, demonstreren de auteurs dat:

1. Hoge BH spins en grote massaverhoudingen de dynamica van getijden-disruptie aanzienlijk veranderen, de vorming van de accretieschijf vertragen en de angulaire distributie van de ejecta wijzigen.
2. De geometrische en thermodynamische eigenschappen van de ejecta direct verbonden zijn met de binaire configuratie, wat de resulterende $r$-proces opbrengsten en kilonova lichtcurves beïnvloedt.
3. De gedetailleerde beschrijving van de angulaire omvang (polaire en azimutale) die hier wordt gegeven, noodzakelijke input biedt voor het verfijnen van synthetische kilonova-modellen, die vaak steunen op vaste geometrieën.
4. De bevindingen ondersteunen de levensvatbaarheid van BHNS-fusies als bronnen van zware elementen en multi-messenger signalen, zelfs in afwezigheid van een bevestigde EM-tegenhanger tot nu toe.

Het werk concludeert dat hoewel de dynamische ejecta alleen niet volledig het zonnere abundantiepatroon kunnen reproduceren (door het ontbreken van de eerste piek en de bijdragen van seculiere uitstromingen), het wel een robuuste baseline biedt voor het begrijpen van de zware-element productie in de meest extreme BHNS-fusies.