A Three-Dimensional Exploration of Magnetic Fields, Rotation, and Shock Revival in a $39 M_\odot$ Core-Collapse Supernova Progenitor

Technische Samenvatting: Een driedimensionale verkenning van magnetische velden, rotatie en shockherleving in een 39 M⊙-progenitor van een kerninstortings-supernova

Probleem en Motivatie
De herleving van de gestagneerde shock bij kerninstortings-supernova's (CCSNe) blijft een complex probleem dat de wisselwerking omvat tussen neutrino-verwarming, rotatie, magnetische velden en hydrodynamica. Hoewel het neutrino-mechanisme robuust is in veel progenitors, met name die met hoge compactheid, biedt het magnetorotatie-mechanisme (MR) een alternatief pad voor shockherleving en het aandrijven van energieke transiënten (bijv. hypernova's, Type Ic-BL-supernova's). Dit mechanisme vereist snelle rotatie en sterke magnetische velden om rotatievrije energie om te zetten in magnetische energie via winding en instabiliteiten, waardoor bipolaire uitstromen worden aangedreven. Eerdere studies hebben echter vaak een beroep gedaan op axiale symmetrie (2D), wat jets kunstmatig kan stabiliseren en niet-axiale symmetrische instabiliteiten zoals de kink-instabiliteit kan onderdrukken. Bovendien vereist de specifieke impact van magnetische velden op tijdschalen voor shockherleving, de morfologie van uitstromen en de stabiliteit van de proto-neutronenster (PNS) tegenover zwarte-gatvorming in omgevingen met hoge compactheid en snelle rotatie in 3D verder onderzoek.

Methodologie
De auteurs voerden driedimensionale hydrodynamische (HD) en magnetohydrodynamische (MHD) simulaties uit van een 39 M⊙-progenitor aan het begin van de hoofdreeks (ZAMS) met hoge compactheid ($\xi_{2.5} = 0.75$) en snelle rotatie ($\Omega_c \approx 0.58$ rad s$^{-1}$). De progenitor werd geëvolueerd met MESA met quasi-chemisch homogene evolutie. Drie distincte modellen werden gesimuleerd met de FLASH-code met de Spark-oplosser, waarbij adaptief meshverfijning (AMR) en een meerdimensionaal, energie-afhankelijk, twee-moment neutrino-transportschema met een M1-sluiting werden toegepast:

1. B1R1: Een roterend, gemagnetiseerd model. Een poloidaal zaadveld ($B_0 \approx 5 \times 10^{11}$ G) werd geïnitieerd, uitgelijnd met de rotatieas.
2. B0R1: Een roterend, niet-gemagnetiseerd controlemodel (alleen HD).
3. B0R0: Een niet-roterend, niet-gemagnetiseerd basismodel (alleen HD).

De simulaties evolueerden van vóór de instorting via de kernbouncing tot in de vroege post-bouncing-fase (tot $\sim$1 seconde), met gebruikmaking van een Cartesiaans rooster met een fijnste resolutie van $\sim$0,61 km nabij de kern. Algemene relativistische effecten werden benaderd met een effectief potentieel (Geval A van Marek et al. 2006).

Belangrijkste Resultaten
Alle drie de modellen bereikten een runaway-shockuitbreiding binnen het gesimuleerde interval, maar vertoonden duidelijk verschillende morfologieën, tijdschalen en lotgevallen van de restanten:

• Shockherleving en Morfologie:

  • B1R1 (Gemagnetiseerd/Roterend): Vertoonde de vroegste shockherleving ($t_{rev} \approx 150$ ms na bouncing) en ontwikkelde de duidelijkste bipolaire uitstroom. De uitstroom was aanvankelijk jet-achtig, maar werd verbreed en vervormd door niet-axiale symmetrische instabiliteiten.
  • B0R0 (Niet-roterend): Reviveerde later ($t_{rev} \approx 250$ ms) via sterke door neutrino's aangedreven turbulentie. De shock breidde zich uit in grote, ongeordende pluimen zonder een voorkeurs-polaire oriëntatie. Cruciaal instortte dit model ongeveer 1 seconde na bouncing tot een zwart gat vanwege het ontbreken van rotatieondersteuning.
  • B0R1 (Roterend/Niet-gemagnetiseerd): Ervond de laatste shockherleving ($t_{rev} \approx 370$ ms). Rotatie onderdrukte convectie en verlaagde de neutrino-luminositeiten, wat leidde tot een compacter evolutie vergeleken met B0R0. Het ontwikkelde geen duidelijke bipolaire uitstroom.

• Hoekmoment en Energietransport:

  • In B1R1 verdeelden Maxwell-spanningen het hoekmoment efficiënt naar buiten vanuit het binnenste PNS, waardoor de spin van de binnenste kern werd verminderd en rotatievrije energie werd onttrokken. Dit proces putte de reservoir van differentieel rotatie uit om de polaire uitstroom aan te drijven.
  • De begroting van het hoekmoment toonde aan dat, terwijl advektie hoekmoment aan de binnenste kern aanleverde, magnetische koppel momenten de dominante put vormden, wat de spin-up compenseerde. Neutrino-emissie fungeerde als een secundaire, kleinere put.
  • De diagnostische energie ($E_{diag}$) in B1R1 groeide het vroegst en snelst, in overeenstemming met de extractie van rotatie-energie. De begroting van de magnetische energie gaf aan dat een aanzienlijk deel van de onttrokken energie door een effectieve dissipatieterm ging ($\dot{Q}_{diss}$), waarschijnlijk geassocieerd met niet-opgeloste reconnectie en sub-roosteroverdracht, in plaats van dat het werd behouden als opgeloste magnetische energie.

• Instabiliteiten en 3D-dynamica:

  • De gemagnetiseerde uitstroom in B1R1 was onderhevig aan een niet-axiale symmetrische kink-instabiliteit ($m=1$-modus). Diagnostiek toonde een vroege exponentiële groei van de verplaatsing van het zwaartepunt van de magnetische energie (e-vouwingstijd $\sim$1,5 ms), in overeenstemming met theoretische kink-tijdschalen.
  • In tegenstelling tot axiaal-symmetrische simulaties die een stabiele, smalle jet zouden kunnen voorspellen, toonde de 3D-simulatie dat de uitstroomtrechter in de loop van de tijd kantelde (enkele graden) en verbreedde. De uitstroom overleefde en propageerde naar grote stralen, maar verloor zijn smalle collimatie, wat suggereert dat 3D-effecten cruciaal zijn voor het nauwkeurig modelleren van jet-stabiliteit.

• Nucleosynthese:

  • Naschakeling van het uitgeworpen materiaal toonde duidelijke samenstellingsverschillen. B0R0 produceerde proton-rijk uitgeworpen materiaal ($Y_e > 0.5$) dat typerend is voor door neutrino's aangedreven explosies, waarbij tot ijzergroep-elementen werd gesynthetiseerd ($A \sim 56$).
  • B1R1 produceerde een aanzienlijk neutron-rijk component ($Y_e$ tot $\sim$0,33) vanwege de snelle ontwikkeling van de uitstroom, waardoor een zwak r-proces mogelijk was tot $A \sim 90$. Echter, het initiële zaadveld en de door kink veroorzaakte verbreding voorkwamen de productie van het extreem neutron-rijke uitgeworpen materiaal ($Y_e < 0.32$) dat wordt gezien in modellen met meer extreme magnetische velden.

Betekenis en Claims
Het artikel beweert dat in een extreme, hoge-compactheid-progenitor, snelle rotatie in combinatie met sterke magnetische velden een vroege, magnetisch ondersteunde weg naar shockherleving kan openen die intrinsiek driedimensionaal is. Belangrijkste conclusies zijn:

1. Magnetische velden als drijvende kracht: Magnetische spanningen vormen niet alleen de uitstroom; ze onttrekken actief rotatievrije energie, vertragen de binnenste kern en faciliteren shockherleving eerder dan in niet-gemagnetiseerde roterende gevallen.
2. Noodzaak van 3D: De evolutie van de door magnetische velden aangedreven uitstroom is intrinsiek niet-axiaal symmetrisch. De kink-instabiliteit verbreedt en kantelt de uitstroom, waardoor de vorming van een aanhoudende, smalle jet wordt voorkomen. Axiaal-symmetrische simulaties kunnen zowel de stabiliteit als de collimatie van dergelijke uitstromen overschatten.
3. Stabiliteit van het restant: Rotatie biedt aanzienlijke ondersteuning tegen directe zwarte-gatvorming in progenitors met hoge compactheid, hoewel de langetermijnstabiliteit van het restant buiten het gesimuleerde interval onzeker blijft.
4. Energiepad: De energietransfer van differentieel rotatie naar de uitstroom behelst aanzienlijke niet-opgeloste dissipatie, wat het belang van sub-roosterfysica in MHD-CCSN-simulaties benadrukt.

De auteurs wijzen op beperkingen, waaronder het gebruik van een sterk zaadveld om niet-opgeloste kleinschalige dynamo-versterking te compenseren, resolutie-gevoeligheid in dissipatiediagnostiek, en de noodzaak van langere simulaties om het uiteindelijke lot van de roterende modellen te bepalen. Ze benadrukken dat, hoewel de simulaties een beeld van magnetorotatie-gestuurde shockherleving ondersteunen, de specifieke uitkomsten sterk afhankelijk zijn van de volledige 3D-dynamica van de uitstroom.