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

Technischer Bericht: Hydrodynamische Simulationen der rekurrenten Nova T Coronae Borealis: Vorhersagen zur Nukleosynthese

Problem und Kontext
Rekurrente Novae (RNe) sind Doppelsternsysteme, die durch thermonukleare Durchbrüche (Thermonuclear Runaways, TNRs) auf der Oberfläche von Weißen Zweren (WDs) gekennzeichnet sind, welche sich durch bemerkenswert kurze Wiederkehrzeiten auszeichnen, typischerweise zwischen 1 und 100 Jahren. Diese kurzen Intervalle erfordern sehr hohe Massen-Akkretionsraten, Weißer-Zwer-Massen nahe der Chandrasekhar-Grenze und hohe initiale Leuchtkräfte der Weißen Zwerge. T Coronae Borealis (T CrB) ist eine prominente galaktische RN, die in den Jahren 1866 und 1946 ausgebrochen ist und eine vorhergesagte Rekurrenzperiode von etwa 80 Jahren aufweist. Angesichts des unmittelbar bevorstehenden nächsten Ausbruchs besteht ein kritischer Bedarf daran, die physikalischen Parameter des Systems einzugrenzen – insbesondere die Masse des Weißen Zwergs, die initiale Leuchtkraft, die Akkretionsrate und die Metallizität –, um die Bedingungen zu verstehen, die solche Eruptionen antreiben. Darüber hinaus sind RNe potenzielle Vorläufer von Typ-Ia-Supernovae, weshalb die Charakterisierung ihrer Nukleosynthese und Massenakkumulation von entscheidender Bedeutung ist.

Methodik
Die Autoren verwendeten den hydrodynamischen, lagrange-basierten Finite-Differenzen-Code SHIVA, um die Explosion von T CrB zu simulieren. Die Simulationen nutzten eine allgemeine Zustandsgleichung, die ein entartetes Elektronengas, ein mehrkomponentiges Ionenplasma und Strahlung sowie OPAL-Rosseland-Mittelwerte der Opazität umfasst. Die nukleare Energieerzeugung wurde über ein Reaktionsnetzwerk von 120 Spezies (von $^1$H bis $^{48}$Ti) berechnet, das durch 630 Prozesse verknüpft ist und aktualisierte Raten aus der STARLIB-Datenbank verwendet.

Die Studie präsentiert 11 neue hydrodynamische Modelle, die darauf ausgelegt sind, eine 80-jährige Rekurrenzperiode ($\tau_{rec}$) unter variierenden Bedingungen zu reproduzieren:

• Weiße-Zwer-Massen ($M_{WD}$): 1,30, 1,35 und 1,38 $M_{\odot}$.
• Initiale Weiße-Zwer-Leuchtkräfte ($L_{WD}$): 0,01, 0,1 und 1,0 $L_{\odot}$.
• Akkretierte Metallizität ($Z_{acc}$): Solar ($1,0 Z_{\odot}$), sub-solar ($0,1 Z_{\odot}$) und super-solar ($10 Z_{\odot}$).
• Mischung: 8 zusätzliche Modelle wurden für Weiße Zwerge mit 1,35 und 1,38 $M_{\odot}$ berechnet, um den Einfluss der Mischung (10 % und 25 %) zwischen der akkretierte Hülle und den darunter liegenden Schichten des Weißen Zwergs (Kohlenstoff-Sauerstoff- und Sauerstoff-Neon-Zusammensetzung) zu bewerten.
• Geringere Massen der Weißen Zwerge: 3 Modelle wurden für einen Weißen Zwerg mit 1,20 $M_{\odot}$ berechnet, um RNe mit niedrigeren Massenschätzungen zu berücksichtigen.

Für jede Konfiguration wurde die Massen-Akkretionsrate ($\dot{M}_{acc}$) so abgestimmt, dass nach exakt 80 Jahren der Akkretion eine Explosion ausgelöst wird.

Wichtigste Ergebnisse

• Akkretionsraten und physikalische Bedingungen: Um eine 80-jährige Rekurrenz mit Material mit solarer Zusammensetzung zu erreichen, sind Massen-Akkretionsraten von $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$ für Weiße Zwerge mit Massen von $1,30 - 1,38 M_{\odot}$ und Leuchtkräften von $0,01 - 1,0 L_{\odot}$ erforderlich. Niedrigere initiale Leuchtkräfte der Weißen Zwerge, geringere Massen der Weißen Zwerge oder eine reduzierte Metallizität im akkretierte Material erfordern höhere Akkretionsraten, um innerhalb des 80-Jahres-Fensters den kritischen Zünddruck zu erreichen.
• Explosionsdynamik: Modelle mit niedrigerem $L_{WD}$ oder niedrigerem $Z_{acc}$ führen zu höheren akkumulierten Hüllmassen und höheren Zünddrücken. Folglich produzieren diese Modelle heftigere Ausbrüche, die durch höhere Spitzentemperaturen ($T_{max}$), größere ausgestoßene Massen ($\Delta M_{eje}$) und größere kinetische Energien gekennzeichnet sind. Beispielsweise erreicht ein 1,38 $M_{\odot}$ schwerer Weißer Zwerg mit $L_{WD}=0,01 L_{\odot}$ und solarer Metallizität ein $T_{max} \approx 2,92 \times 10^8$ K, während ein 1,30 $M_{\odot}$ schwerer Weißer Zwerg mit $L_{WD}=1,0 L_{\odot}$ nur $2,01 \times 10^8$ K erreicht.
• Nukleosynthese und Elementhäufigkeiten:
  • Isotopenverhältnisse: Das Ejekta wird von CNO-Zyklus-Isotopen ($^{12,13}$C, $^{14,15}$N) dominiert. Es werden signifikante Variationen in den Isotopenverhältnissen (z. B. $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si) über verschiedene Modelle hinweg beobachtet.
  • Abhängigkeit von der Weißer-Zwer-Masse: Modelle mit unterschiedlichen Massen der Weißen Zwerge, aber identischen initialen Leuchtkräften, weisen signifikante Unterschiede in den Häufigkeiten von Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca und Sc auf. Insbesondere zeigen 1,38 $M_{\odot}$ Modelle eine signifikante Anreicherung in S, Cl, K, Ca und Sc im Vergleich zu solaren Häufigkeiten, während 1,20 $M_{\odot}$ und 1,30 $M_{\odot}$ Modelle im Allgemeinen solare oder leicht sub-solare Häufigkeiten für diese Elemente aufweisen (mit Ausnahme des Modells 130C). Silizium ist in den 1,20 und 1,30 $M_{\odot}$ Modellen angereichert, aber in den 1,38 $M_{\odot}$ Modellen verarmt oder nahe bei solarer Häufigkeit.
  • Effekte der Metallizität: In den 1,38 $M_{\odot}$ Modellen verändert die variierende Metallizität die Zusammensetzung drastisch. Eine niedrige Metallizität ($0,1 Z_{\odot}$) führt zu einer signifikanten $^{40}$Ca-Produktion, während eine hohe Metallizität ($10 Z_{\odot}$) die Dominanz hin zu N, C, Ne und Si verschiebt.
• Mischungsszenarien: Modelle, die eine Mischung zwischen der Hülle und dem Substrat des Weißen Zwergs (CO oder ONe) beinhalten, zeigen verkürzte Akkretionsphasen und reduzierte Gesamtejektamassen im Vergleich zu Nicht-Mischungs-Modellen. Entscheidend ist, dass die Mischung die Produktion von Gamma-Strahler-Emittenten ($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al) signifikant erhöht. Während Nicht-Mischungs-Modelle eine vernachlässigbare Gamma-Strahlung vorhersagen, könnten Mischungsmodelle detektierbare Signaturen erzeugen, was potenziell die Unterscheidung zwischen CO- und ONe-Weißen Zwergen ermöglicht.
• Gamma-Strahlung: Für die Standard-Nicht-Mischungs-Modelle sagen die Autoren voraus, dass die Gamma-Linien (511 keV und 1275 keV) und das damit verbundene Kontinuum aufgrund der geringen Ausbeuten radioaktiver Isotope um Größenordnungen unter den Detektionslimits aktueller Instrumente liegen werden.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass die kompositorischen Unterschiede im Ejekta, insbesondere die Elementhäufigkeiten von Spezies wie Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca und Sc, als potenzielles diagnostisches Werkzeug dienen. Diese Variationen können den Parameterraum einschränken und die verschiedenen in der Studie vorgeschlagenen T CrB-Modelle differenzieren. Speziell behaupten die Autoren, dass frühe spektroskopische Messungen des bevorstehenden T CrB-Ausbruchs die Masse des Weißen Zwergs und den Grad der Mischung bestimmen könnten, wodurch die physikalischen Bedingungen, die die Eruptionen des Systems steuern, geklärt werden können. Die Studie kommt zu dem Schluss, dass, obwohl eine Gamma-Detektion ohne signifikante Mischung unwahrscheinlich ist, die Elementhäufigkeitsmuster eine robuste Methode zur Identifizierung der Natur des zugrunde liegenden Weißen Zwergs und der stattfindenden Mischungsprozesse bieten.