Black Hole Spin-down in Collapsars in 3D Neutrino Transport GRMHD Simulations

Probleemstelling

Langdurige gammaflitsen (LGRBs) worden geacht te worden aangedreven door "collapsars": massieve sterren waarvan de kern instort tot een zwart gat (BH). De energie van de daaruit voortkomende relativistische jets hangt sterk af van de spin van het zwarte gat ($a$) en de magnetische flux ($\Phi$).
Eerdere studies, gebaseerd op niet-stralende (adiabatische) schijven, toonden aan dat magnetisch gearresteerde schijven (MADs) het zwarte gat snel tot een zeer lage evenwichtsspin kunnen afremmen ($a_{eq} \approx 0.035 - 0.07$). Dit resulteert in jets met zo weinig energie dat ze moeite hebben om uit de ster te ontsnappen.
Echter, de werkelijke omstandigheden in een collapsar zijn extreem dicht en heet, waardoor de accretieschijf koelt via neutrino-emissie. Deze koeling verandert de thermodynamica en de dikte van de schijf, wat de hoekmomentumtransport en de spin-evolutie van het zwarte gat kan beïnvloeden. Tot nu toe ontbrak het aan simulaties die neutrino-transport en GRMHD (General Relativistic Magnetohydrodynamics) gecombineerd gebruiken om dit effect te bestuderen.

Methodologie

De auteurs hebben een reeks van 3D GRMHD-simulaties uitgevoerd met de GPU-versnelde code $\nu$H-AMR.

• Fysica: De simulaties omvatten twee-moment neutrino-transport (voor elektron-neutrino's, antineutrino's en zware lepton-neutrino's) gekoppeld aan ideale GRMHD. De vergelijkingen van Einstein en de magnetohydrodynamica worden opgelost in Kerr-Schild-coördinaten.
• Modelopzet: Er zijn twee soorten progenitor-sterrenmodellen gebruikt met verschillende dichtheidsprofielen:
  1. p#-modellen: Een typisch power-law profiel ($\rho \propto r^{-1.5}$), resulterend in accretiesnelheden van $\dot{m} \sim 0.1 - 1 M_\odot/s$.
  2. c#-modellen: Een model met een constante dichtheidskern, resulterend in veel hogere initiële accretiesnelheden ($\dot{m} \gtrsim 10 M_\odot/s$).
• Variatie: Voor beide modellen zijn simulaties uitgevoerd met verschillende initiële spins ($a$ van 0.0 tot 0.8) en met/zonder neutrino-transport (voor vergelijking).
• Analyse: De auteurs meten de spin-up parameter ($s$), de magnetische flux ($\phi$), en de efficiënties van elektromagnetische en neutrino-uitstroom. Ze passen een semi-analytisch model toe om de spin-evolutie te extrapoleren naar de tijdschaal van een volledige GRB ($\sim 30-40$ seconden).

Belangrijkste Bijdragen

1. Eerste 3D Simulaties: Dit is de eerste studie die succesvolle collapsar-jets simuleert met volledige neutrino-transport in 3D GRMHD.
2. Spin-Equilibrium Model: De auteurs hebben een nieuw evenwichtspunt voor de spin van het zwarte gat bepaald voor neutrino-gekoelde schijven, wat afwijkt van eerdere modellen voor niet-stralende schijven.
3. Koppeling met Observaties: Ze verbinden de theoretische spin-evolutie direct met waarnemingen van LIGO/Virgo/KAGRA (zwaartekrachtgolven) en de luminositeit van GRBs.

Resultaten

• Evenwichtsspin ($a_{eq}$): In tegenstelling tot niet-stralende MADs ($a_{eq} \approx 0.035-0.07$), remmen neutrino-gekoelde MADs het zwarte gat af tot een hogere evenwichtsspin van $a_{eq} \approx 0.13$.
  • Dit komt doordat neutrino-koeling de schijf dunner maakt, wat de spin-down efficiëntie vermindert.
  • Deze waarde is 2 tot 4 keer hoger dan bij niet-stralende schijven en resulteert in jets die 4 tot 16 keer krachtiger zijn.
• MAD-toestand: De systemen bereiken snel de MAD-toestand (binnen $0.3 - 0.4$ seconden), waarbij de magnetische flux saturatie bereikt ($\phi \approx 50$). Zolang de flux aanwezig blijft, blijft het systeem in deze toestand.
• Invloed van Neutrino's: Hoewel de directe bijdrage van neutrino's aan het koppel (torque) op de spin verwaarloosbaar klein is, zijn ze cruciaal omdat ze de schijfstructuur (dikte) bepalen, wat op zijn beurt de elektromagnetische en hydrodynamische koppels beïnvloedt.
• Stabiliteit: Bij zeer lage spins ($a \approx 0.1$) en ondieke dichtheidsprofielen (p#-modellen) vertonen de jets variabiliteit en instabiliteit (kink-instabiliteit), wat kan leiden tot disk-misalignement en tijdelijk verlies van de MAD-toestand.
• Tijdschaal: Voor typische GRB-duraties ($t \gtrsim 30$ s) hebben de zwarte gaten voldoende tijd om massa te accretteren en de evenwichtsspin van $0.13$ te bereiken, ongeacht de initiële spin.

Betekenis en Conclusie

De resultaten van dit onderzoek hebben belangrijke implicaties voor de astrofysica:

1. GRB Luminositeit: De hogere evenwichtsspin ($a \approx 0.13$) verklaart waarom collapsars krachtige jets kunnen produceren die de ster kunnen doorboren, terwijl eerdere modellen met $a \approx 0.04$ te zwakke jets voorspelden.
2. Zwaartekrachtgolven (GW): De voorspelde spinverdeling ($a \approx 0.1 - 0.2$) komt uitstekend overeen met de waarnemingen van LIGO/Virgo/KAGRA, die aangeven dat zwarte gaten in binaire systemen vaak lage spins hebben. Dit ondersteunt het idee dat collapsars de progenitors kunnen zijn van sommige binaire zwarte gat-samensmeltingen.
3. Beperkingen: De studie suggereert dat de maximale luminositeit van GRBs wordt beperkt door de spin-dynamiek van het zwarte gat; als de spin te snel afremt voordat de jet de ster verlaat, wordt de burst minder krachtig.

Samenvattend biedt dit werk een zelfconsistentere benadering van collapsar-modellering, waarbij neutrino-koeling een sleutelrol speelt in het bepalen van de uiteindelijke spin van het zwarte gat en de daaruit voortkomende energie van gammaflitsen.