Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

Technische Samenvatting: Hydrodynamische Simulaties van de Recurrente Nova T Coronae Borealis: Voorspellingen voor Nucleosynthese

Probleem en Context
Recurrente nova's (RN's) zijn binaire systemen die worden gekenmerkt door thermonucleaire runaways (TNR's) die plaatsvinden op het oppervlak van witte dwergen (WD's) met opmerkelijk korte herhalingsintervallen, doorgaans tussen 1 en 100 jaar. Deze korte intervallen vereisen zeer hoge massa-accretiecijfers, witte dwerfmassa's die de Chandrasekhar-limiet naderen, en hoge initiële WD-luminositeiten. T Coronae Borealis (T CrB) is een prominente Galactische RN die uitbarstte in 1866 en 1946, met een voorspelde herhalingsperiode van ongeveer 80 jaar. Gezien de nabijheid van de volgende uitbarsting is er een kritische behoefte om de fysische parameters van het systeem te beperken—met name de WD-massa, de initiële luminositeit, het accretiecijfer en de metalliciteit—om de condities te begrijpen die dergelijke erupties aandrijven. Bovendien zijn RN's potentiële voorlopers van Type Ia supernova's, waardoor de karakterisering van hun nucleosynthese en massaccumulatie essentieel is.

Methodologie
De auteurs gebruikten de hydrodynamische, Lagrangiaanse, eindige-verschil code SHIVA om de explosie van T CrB te simuleren. De simulaties maakten gebruik van een algemene toestandsvergelijking inclusief gedegenereerd elektrengas, multicomponent ionenplasma en straling, met OPAL Rosseland gemiddelde opaciteiten. De nucleaire energiegeneratie werd berekend via een reactienetwerk van 120 soorten (van $^1$H tot $^{48}$Ti) verbonden door 630 processen, waarbij gebruik is gemaakt van bijgewerkte snelheden uit de STARLIB-database.

De studie presenteert 11 nieuwe hydrodynamische modellen die ontworpen zijn om een herhalingsperiode van 80 jaar ($\tau_{rec}$) te reproduceren onder variërende condities:

• WD-massa's ($M_{WD}$): 1.30, 1.35, en 1.38 $M_{\odot}$.
• Initiële WD-luminositeiten ($L_{WD}$): 0.01, 0.1, en 1.0 $L_{\odot}$.
• Geaccreteerde Metalliciteit ($Z_{acc}$): Solair ($1.0 Z_{\odot}$), sub-solair ($0.1 Z_{\odot}$), en super-solair ($10 Z_{\odot}$).
• Mixing (Mengeling): 8 aanvullende modellen werden berekend voor 1.35 en 1.38 $M_{\odot}$ WD's om de impact van mixing (10% en 25%) tussen de geaccreteerde envelop en de onderliggende WD-lagen (Carbon-Oxygen en Oxygen-Neon composities) te evalueren.
• Lagere Massa WD's: 3 modellen werden berekend voor een 1.20 $M_{\odot}$ WD om RN's met lagere massa-inschattingen te accounten.

Voor elke configuratie werd de massa-accretiesnelheid ($\dot{M}_{acc}$) afgestemd om een explosie na exact 80 jaar accretie te triggeren.

Belangrijkste Resultaten

• Accretiecijfers en Fysische Condities: Om een 80-jarige herhaling te bereiken met materiaal van solaire compositie, zijn massa-accretiecijfers van $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$ vereist voor WD's met massa's van $1.30 - 1.38 M_{\odot}$ en luminositeiten van $0.01 - 1.0 L_{\odot}$. Lagere initiële WD-luminositeiten, lagere WD-massa's of een verminderde metalliciteit in het geaccreteerde materiaal vereisen hogere accretiecijfers om de kritische ontbrandingsdruk binnen het 80-jarige venster te bereiken.
• Explosiedynamiek: Modellen met een lagere $L_{WD}$ of een lagere $Z_{acc}$ resulteren in hogere geaccumuleerde envelopmassa's en hogere ontbrandingsdrukken. Bijgevolg produceren deze modellen meer gewelddadige uitbarstingen gekenmerkt door hogere piektemperaturen ($T_{max}$), grotere geëjecteerde massa's ($\Delta M_{eje}$) en grotere kinetische energieën. Bijvoorbeeld, een 1.38 $M_{\odot}$ WD met $L_{WD}=0.01 L_{\odot}$ en solaire metalliciteit bereikt $T_{max} \approx 2.92 \times 10^8$ K, terwijl een 1.30 $M_{\odot}$ WD met $L_{WD}=1.0 L_{\odot}$ slechts $2.01 \times 10^8$ K bereikt.
• Nucleosynthese en Elementaire Abundanties:
  • Isotopische Ratio's: Het ejecta wordt gedomineerd door CNO-cyclus isotopen ($^{12,13}$C, $^{14,15}$N). Significante variaties in isotopische ratio's (bijv. $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si) worden waargenomen over verschillende modellen.
  • Afhankelijkheid van WD-massa: Modellen met verschillende WD-massa's maar identieke initiële luminositeiten vertonen significante verschillen in de abundanties van Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, en Sc. Specifiek vertonen 1.38 $M_{\odot}$ modellen een significante verrijking in S, Cl, K, Ca, en Sc ten opzichte van solaire abundanties, terwijl 1.20 $M_{\odot}$ en 1.30 $M_{\odot}$ modellen over het algemeen solaire of licht sub-solaire abundanties vertonen voor deze elementen (met uitzondering van model 130C). Silicon is verhoogd in de 1.20 en 1.30 $M_{\odot}$ modellen, maar gedepleteerd of nabij-solair in de 1.38 $M_{\odot}$ modellen.
  • Effecten van Metalliciteit: In 1.38 $M_{\odot}$ modellen verandert variërende metalliciteit de compositie drastisch. Lage metalliciteit ($0.1 Z_{\odot}$) leidt tot significante $^{40}$Ca productie, terwijl hoge metalliciteit ($10 Z_{\odot}$) de dominantie verschuift naar N, C, Ne, en Si.
• Mixing Scenario's: Modellen die mixing bevatten tussen de envelop en de WD-substraat (CO of ONe) vertonen kortere accretiefases en lagere totale geëjecteerde massa's vergeleken met non-mixing modellen. Cruciaal is dat mixing de productie van gamma-straaleremissies ($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al) aanzienlijk verhoogt. Waar non-mixing modellen verwaarloosbare gamma-emissie voorspellen, zouden mixing-modellen detecteerbare signaturen kunnen produceren, wat potentieel kan helpen bij het onderscheiden tussen CO en ONe WD's.
• Gamma-Straling: Voor de standaard non-mixing modellen voorspellen de auteurs dat gamma-lijnen (511 keV en 1275 keV) en de bijbehorende continuüm orden van grootte onder de detectielimieten van huidige instrumenten liggen vanwege de lage opbrengsten van radioactieve isotopen.

Betekenis en Claims
Het artikel stelt dat de compositionele verschillen in het ejecta, met name de elementaire abundanties van soorten zoals Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca, en Sc, dienen als een potentieel diagnostisch instrument. Deze variaties kunnen de parameterruimte beperken en de verschillende T CrB-modellen die in de studie worden voorgesteld, onderscheiden. Specifiek claimen de auteurs dat vroege spectroscopische metingen van de komende T CrB uitbarsting de WD-massa en de mate van mixing kunnen bepalen, waardoor de fysische condities die de erupties van het systeem beheersen, kunnen worden verduidelijkt. De studie concludeert dat hoewel gamma-detectie onwaarschijnlijk is zonder significante mixing, de elementaire abundantiepatronen een robuuste methode bieden voor het identificeren van de aard van de onderliggende WD en de mixing-processen die aan het werk zijn.