On the impact of the carbon fusion rate over the properties of superbursts -- Numerical simulations of superbursts with MESA

Technische Samenvatting: Over de impact van de koolstoffusierate op de eigenschappen van superbursts

Probleemstelling
Superbursts zijn energetische thermonucleaire explosies die plaatsvinden in de korsten van neutronensterren binnen Low-Mass X-ray Binaries (LMXBs), getriggerd door instabiele koolstofverbranding bij temperaturen $T \le 10^9$ K. Hoewel eerdere multi-zone simulaties (bijv. Keek & Heger 2011, 2012) succesvol de recurrentieperioden en lichtcurves van superbursts hebben gemodelleerd, vertrouwden zij op de klassieke koolstoffusierate van Caughlan & Fowler (1988, CF88). Recente nucleaire natuurkundige onderzoeken hebben deze rate uitgedaagd, wat suggereert dat de werkelijke rate bij relevante temperaturen ($T \le 10^9$ K) aanzienlijk anders zou kunnen zijn. Specifiek suggereert experimentele data van de LUNA-collaboratie (Tumino et al. 2018) dat de rate met een factor $\sim 10^3$ verhoogd kan worden, terwijl hinder-modellen (Monpribat et al. 2022) suggereren dat de rate met een vergelijkbare factor verlaagd kan worden. De impact van deze substantiële variaties op superburst-observabelen—zoals recessietijden, vervaltijden, ontbrandingsdieptes en nucleosynthese-opbrengsten—blijft onverkend in de context van moderne stellaire evolutiecodes.

Methodologie
De auteurs gebruikten de publieke stellaire evolutiecode MESA v.24.08.1 om de evolutie van neutronenster-enveloppen te simuleren tijdens massatoevoeging (accretie). De simulaties maakten gebruik van een gereduceerd nucleair reactienetwerk, Approx_A56 (dekking tot nucliden met $A=56$), met selecte runs met Approx_A64 om effecten van afkap (truncation) te verifiëren.

Belangrijke parameters en configuraties omvatten:

• Accretiesnelheden: Twee fiduciaire snelheden werden getest: een "gemiddelde" rate ($\dot{M}_{moderate} \approx 5.88 \times 10^{-9} M_\odot \text{yr}^{-1}$) en een "hoge" rate ($\dot{M}_{high} \approx 2 \times 10^{-8} M_\odot \text{yr}^{-1}$).
• Samenstelling: Het geaccreteerde materiaal werd gemodelleerd als 20% $^{12}\text{C}$ en 80% $^{56}\text{Fe}$, consistent met de literatuurvereisten voor het triggeren van superbursts. De initiële envelop bestond uit zuiver $^{56}\text{Fe}$.
• Basale Verwarming: Een vaste basisluminositeit ($L_b$) werd toegepast aan de onderzijde van de envelop om instabiele koolstofontbranding te waarborgen, gekalibreerd om overeen te komen met eerdere KEPLER-simulaties.
• Koolstoffusierates: Vier verschillende reactierates werden vergeleken:
  1. CF88: De klassieke rate (fiduciaal model).
  2. LUNA: Gebaseerd op LUNA-metingen, wat een rate voorspelt die $\sim 10^3$ keer hoger is dan CF88 bij $T \sim 2\text{--}3 \times 10^8$ K.
  3. HIN: Gebaseerd op het hinder-model, wat een rate voorspelt die $\sim 10^3$ keer lager is dan CF88.
  4. HIN-RES: Het HIN-model inclusief een laagliggende resonantie, wat resultaten oplevert die tussen HIN en CF88 in liggen.
• Fysieke Behandelingen: De simulaties bevatten radiatieve opaciteiten (Wolf-voorschrift), Mixing Length Theory voor convectie, en vaste oppervlakte-drukrandvoorwaarden. Algemene relativistische correcties werden weggelaten.

Belangrijkste Resultaten
De studie kwantificeert hoe variaties in de koolstoffusierate de eigenschappen van superbursts veranderen:

1. Tijdschalen en Ontbrandingsdiepte:

   • Verhoogde Rate (LUNA): Leidt tot kortere recessietijden en vervaltijden, en triggert ontbranding bij een ondiepere kolomdiepte ($y_{ignition}$). Voor $\dot{M}_{moderate}$ daalt de recessietijd van $\sim 241$ dagen (CF88) naar $\sim 100$ dagen.
   • Verlaagde Rate (HIN): Resulteert in langere recessie- en vervaltijden, waarbij ontbranding plaatsvindt bij diepere kolomdieptes. Voor $\dot{M}_{moderate}$ neemt de recessietijd toe naar $\sim 364$ dagen.
   • Intermediaire Rate (HIN-RES): Produceert resultaten die bijna ononderscheidbaar zijn van het CF88-fiduciale model.

2. Temperatuur en Lichtcurves:

   • De maximale temperatuur die bereikt wordt tijdens de explosie is gevoelig voor de rate. De HIN-rate levert de hoogste piektemperaturen op ($\sim 5 \times 10^9$ K), terwijl LUNA lagere pieken levert ($\sim 3 \times 10^9$ K).
   • De vervalprofielen van de luminositeit verschuiven dienovereenkomstig. LUNA-simulaties vertonen een sneller verval, terwijl HIN-simulaties een langzamer verval vertonen.
   • Bij $\dot{M}_{high}$ is de piek-luminositeit grotendeels ongevoelig voor de keuze van de rate, maar de vervaltijden en recessie-intervallen variëren nog steeds significant.

3. Nucleosynthese en As-samenstelling:

   • De distributie van de as vertoont een karakteristieke dubbele piekstructuur (rond $^{28}\text{Si}$ en $^{56}\text{Fe}$) in alle modellen.
   • Rate-gevoeligheid: Voor $A \le 44$ is de abundantie van $\alpha$-nucliden zeer gevoelig voor de rate. De LUNA-rate (hogere temperatuur/efficiëntie) produceert aanzienlijk meer $\alpha$-nucliden vergeleken met CF88/HIN-RES, terwijl de HIN-rate er minder produceert.
   • Effecten van Hoge Accratiesnelheid: Bij $\dot{M}_{high}$ is er een merkbare versterking in de synthese van kernen met $28 \le A \le 50$ en een specifieke accumulatie van $\sim 0.1\%$ $^{16}\text{O}$, een kenmerk dat minder prominent is bij gematigde rates.
   • Fotodisintegratie: Simulaties waarbij de inerte $^{56}\text{Fe}$ werd vervangen door $^{72}\text{Ge}$, bevestigden dat de ijzer-piek as-distributie grotendeels wordt gedreven door $^{56}\text{Fe}$ fotodisintegratie in plaats van enkel door koolstoffusie. Het verwijderen van fotodisintegratiekanalen stopte de nucleosynthese voorbij $^{40}\text{Ca}$.

4. Degeneratie met Basale Verwarming:

   • De studie toont aan dat het veranderen van de koolstoffusierate een effect heeft op de superburst-tijdschalen en ontbrandingsdiepte dat vergelijkbaar is met het veranderen van de basisluminositeit ($L_b$) van de envelop. Bijvoorbeeld, het verlagen van de rate (HIN) en het verhogen van $L_b$ produceert vergelijkbare recessietijden. Dit suggereert een degeneratie waarbij het fitten van de lichtcurve alleen niet uniek zowel de fusierate als de basale verhitting kan vaststellen.

Betekenis en Conclusies
Het artikel concludeert dat de koolstoffusierate bij $T \le 10^9$ K een kritische parameter is die de fysica van superbursts bepaalt. De auteurs stellen dat:

• Onzekerheden in de koolstoffusierate variaties in geobserveerde recessie- en vervaltijden kunnen verklaren zonder noodzakelijkerwijs een verandering in de massatoevoersnelheid of basale verhitting aan te nemen.
• De "canonieke" as-samenstelling afgeleid van eerdere KEPLER-simulaties (gebruikmakend van CF88) robuust is wat betreft de dubbele piekstructuur, maar gevoelig is voor de specifieke reactierate wat betreft de abundanties van $\alpha$-nucliden.
• De degeneratie tussen de koolstoffusierate en de basale verhitting impliceert dat het fitten van een enkele superburst-lichtcurve de onderliggende nucleaire fysica of de thermische staat van de neutronenster-korst niet ondubbelzinnig kan bepalen.
• De resultaten valideren het gebruik van MESA voor superburst-modellering, waarbij consistentie wordt getoond met eerdere KEPLER-resultaten, terwijl een publieke code-optie wordt geboden voor het verkennen van nucleaire fysica-onzekerheden.

De auteurs merken bescheiden op dat hoewel hun bevindingen suggereren dat de koolstoffusierate de noodzaak voor hoge basale verhitting in modellen zou kunnen verminderen (indien de LUNA-rate correct is), de onzekerheid in massatoevoersnelheden voor specifieke bronnen (zoals 4U 1636-53) definitieve conclusies over welke rate fysiek gerealiseerd wordt, verhindert. Verder verduidelijken zij dat hun behandeling van koolstofverbranding als puur thermonucleair gerechtvaardigd is door de dichtheden en temperaturen die bereikt worden, verschillend van pycnonucleaire regimes.