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

Resumen Técnico: Simulaciones Hidrodinámicas de la Nova Recurrente T Coronae Borealis: Predicciones de Nucleosíntesis

Problema y Contexto
Las novas recurrentes (RN) son sistemas binarios caracterizados por desatadas térmicas nucleares (TNR) que ocurren en la superficie de enanas blancas (WD) con tiempos de recurrencia notablemente cortos, típicamente entre 1 y 100 años. Estos intervalos cortos requieren tasas de acreción de masa muy altas, masas de la enana blanca que se aproximan al límite de Chandrasekhar y luminosidades iniciales de la enana blanca elevadas. T Coronae Borealis (T CrB) es una prominente RN Galáctica que estalló en 1866 y 1946, con un periodo de recurrencia predicho de aproximadamente 80 años. Dada la inminencia de su próximo estallido, existe una necesidad crítica de restringir los parámetros físicos del sistema —específicamente la masa de la enana blanca, la luminosidad inicial, la tasa de acreción y la metalicidad— para comprender las condiciones que impulsan tales erupciones. Además, las RN son potenciales progenitoras de supernovas de Tipo Ia, por lo que la caracterización de su nucleosíntesis y acumulación de masa es vital.

Metodología
Los autores emplearon el código hidrodinámico lagrangiano de diferencias finitas SHIVA para simular la explosión de T CrB. Las simulaciones utilizaron una ecuación de estado general que incluye gas de electrones degenerado, plasma de iones multicomponente y radiación, con opacidades de Rosseland OPAL. La generación de energía nuclear se calculó mediante una red de reacción de 120 especies (desde $^1$H hasta $^{48}$Ti) vinculadas por 630 procesos, incorporando tasas actualizadas de la base de datos STARLIB.

El estudio presenta 11 nuevos modelos hidrodinámicos diseñados para reproducir un periodo de recurrencia de 80 años ($\tau_{rec}$) bajo diversas condiciones:

• Masas de la Enana Blanca ($M_{WD}$): 1.30, 1.35 y 1.38 $M_{\odot}$.
• Luminosidades Iniciales de la Enana Blanca ($L_{WD}$): 0.01, 0.1 y 1.0 $L_{\odot}$.
• Metalicidad Acrecionada ($Z_{acc}$): Solar ($1.0 Z_{\odot}$), subsolar ($0.1 Z_{\odot}$) y supersolar ($10 Z_{\odot}$).
• Mezcla: Se computaron 8 modelos adicionales para enanas blancas de 1.35 y 1.38 $M_{\odot}$ para evaluar el impacto de la mezcla (10% y 25%) entre la envoltura acrecionada y las capas subyacentes de la enana blanca (composiciones de Carbono-Oxígeno y Oxígeno-Neón).
• Enanas Blancas de Menor Masa: Se computaron 3 modelos para una enana blanca de 1.20 $M_{\odot}$ para dar cuenta de las RN con estimaciones de masa más bajas.

Para cada configuración, la tasa de acreción de masa ($\dot{M}_{acc}$) se ajustó para desencadenar una explosión tras exactamente 80 años de acreción.

Resultados Clave

• Tasas de Acreción y Condiciones Físicas: Para lograr una recurrencia de 80 años con material de composición solar, se requieren tasas de acreción de masa de $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$ para enanas blancas con masas de $1.30 - 1.38 M_{\odot}$ y luminosidades de $0.01 - 1.0 L_{\odot}$. Menores luminosidades iniciales de la enana blanca, menores masas de la enana blanca o una metalicidad reducida en el material acrecionado requieren tasas de acreción más altas para alcanzar la presión de ignición crítica dentro de la ventana de 80 años.
• Dinámica de la Explosión: Los modelos con menor $L_{WD}$ o menor $Z_{acc}$ resultan en masas de envoltura acumuladas más altas y presiones de ignición más altas. En consecuencia, estos modelos producen estallidos más violentos caracterizados por temperaturas pico ($T_{max}$) más altas, masas eyectadas ($\Delta M_{eje}$) más grandes y energías cinéticas mayores. Por ejemplo, una enana blanca de 1.38 $M_{\odot}$ con $L_{WD}=0.01 L_{\odot}$ y metalicidad solar alcanza $T_{max} \approx 2.92 \times 10^8$ K, mientras que una enana blanca de 1.30 $M_{\odot}$ con $L_{WD}=1.0 L_{\odot}$ alcanza solo $2.01 \times 10^8$ K.
• Nucleosíntesis y Abundancias Elementales:
  • Razones Isotópicas: Los eyectos están dominados por isótopos del ciclo CNO ($^{12,13}$C, $^{14,15}$N). Se observan variaciones significativas en las razones isotópicas (por ejemplo, $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si) a través de diferentes modelos.
  • Dependencia de la Masa de la Enana Blanca: Los modelos con diferentes masas de enana blanca pero con la misma luminosidad inicial exhiben diferencias significativas en las abundancias de Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca y Sc. Específicamente, los modelos de 1.38 $M_{\odot}$ muestran un enriquecimiento significativo en S, Cl, K, Ca y Sc en relación con las abundancias solares, mientras que los modelos de 1.20 $M_{\odot}$ y 1.30 $M_{\odot}$ generalmente muestran abundancias solares o ligeramente subsolares para estos elementos (con la excepción del modelo 130C). El silicio está enriquecido en los modelos de 1.20 y 1.30 $M_{\odot}$, pero se encuentra depletado o cerca de los valores solares en los modelos de 1.38 $M_{\odot}$.
  • Efectos de la Metalicidad: En los modelos de 1.38 $M_{\odot}$, variar la metalicidad altera drásticamente la composición. Una metalicidad baja ($0.1 Z_{\odot}$) conduce a una producción significativa de $^{40}$Ca, mientras que una metalicidad alta ($10 Z_{\odot}$) desplaza la dominancia hacia N, C, Ne y Si.
• Escenarios de Mezcla: Los modelos que incorporan la mezcla entre la envoltura y el sustrato de la enana blanca (CO o ONe) muestran fases de acreción acortadas y masas eyectadas totales reducidas en comparación con los modelos sin mezcla. Crucialmente, la mezcla aumenta significativamente la producción de emisores de rayos gamma ($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al). Mientras que los modelos sin mezcla predicen una emisión de rayos gamma insignificante, los modelos con mezcla podrían producir firmas detectables, lo que permitiría distinguir entre enanas blancas de tipo CO y ONe.
• Emisión de Rayos Gamma: Para los modelos estándar sin mezcla, los autores predicen que las líneas de rayos gamma (511 keV y 1275 keV) y el continuo asociado estarán órdenes de magnitud por debajo de los límites de detección de los instrumentos actuales debido a los bajos rendimientos de los isótopos radiactivos.

Significancia y Afirmaciones
El artículo afirma que las diferencias composicionales en los eyectos, particularmente las abundancias elementales de especies tales como Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca y Sc, sirven como una herramienta diagnóstica potencial. Estas variaciones pueden restringir el espacio de parámetros y discriminar entre los diversos modelos de T CrB propuestos en el estudio. Específicamente, los autores afirman que las mediciones espectroscópicas tempranas del próximo estallido de T CrB podrían determinar la masa de la enana blanca y el grado de mezcla, aclarando así las condiciones físicas que gobiernan las erupciones del sistema. El estudio concluye que, si bien la detección de rayos gamma es poco probable sin una mezcla significativa, los patrones de abundancia elemental ofrecen un método robusto para identificar la naturaleza de la enana blanca subyacente y los procesos de mezcla en juego.