Interaction between the ejecta, the accretion disk, and the secondary star in the recurrent nova system U Sco

Resumen Técnico: Interacción entre los eyectos, el disco de acreción y la estrella secundaria en el sistema de nova recurrente U Sco

Planteamiento del Problema
Si bien se ha dedicado un esfuerzo significativo al modelado de las fases iniciales de la fuga termonuclear descontrolada y de la eyección en las novas recurrentes (RNe), la interacción hidrodinámica posterior entre los eyectos de la nova, el disco de acreción que orbita la enana blanca y la compañera estelar secundaria permanece insuficientemente explorada. Esta brecha es crítica para comprender la evolución a largo plazo de las RNe, específicamente respecto a la integridad del disco de acreción (que dicta los periodos de recurrencia), el presupuesto de masa del sistema binario (escape vs. retención) y la potencial contaminación química de la estrella secundaria. El artículo se centra en U Sco, una nova recurrente prototípica con una enana blanca masiva, una compañera subgigante y un periodo orbital corto, para abordar estas dinámicas.

Metodología
El estudio emplea un enfoque de modelado híbrido que combina hidrodinámica 1D y 3D:

1. Condiciones Iniciales: La explosión y la eyección inicial se modelan utilizando el código lagrangiano de simetría esférica SHIVA. La estructura resultante de los eyectos (masa, radio, densidad, temperatura y perfiles de velocidad) se mapea en un dominio 3D.
2. Simulación 3D: La evolución subsiguiente se simula utilizando el código de hidrodinámica de partículas suavizadas (SPH) GADGET-2. El sistema incluye:
   • Una enana blanca puntual ($1.38 M_\odot$).
   • Una secundaria subgigante ($0.88 M_\odot$), modelada solo con sus capas más externas ($0.0054 M_\odot$) resueltas mediante SPH para gestionar la carga computacional, con el núcleo tratado como una masa puntual.
   • Un disco de acreción que orbita la enana blanca, probado en diversas geometrías (ensanchado/flared vs. en forma de V) y perfiles de densidad (derivados de Shakura-Sunyaev vs. constantes).
   • Eyectos de nova con masas que oscilan entre $1.1 \times 10^{-6}$ y $2.1 \times 10^{-6} M_\odot$ y velocidades de hasta $10,000$ km s$^{-1}$.
3. Diseño Experimental: Se llevó a cabo una serie de diez simulaciones (Modelos A–H) para aislar los efectos de parámetros específicos: masa de los eyectos, velocidad de los eyectos, perfil de densidad de los eyectos, masa del disco y geometría del disco. Crucialmente, estas simulaciones incorporan la rotación orbital alrededor del centro de masa del sistema, una característica que suele omitirse en estudios previos de binarias más anchas.

Resultados Clave
Las simulaciones revelan resultados distintos basados en la relación entre la masa de los eyectos y la masa del disco ($M_{ejecta}/M_{disk}$) y la geometría del disco:

• Disrupción del Disco: Para los modelos que presentan discos ensanchados donde $M_{ejecta}/M_{disk} \geq 1$ (Modelos A, B, C, D, G, H), la explosión de la nova resulta en la disrupción completa y el barrido del disco de acreción. La energía cinética de los eyectos excede significativamente la energía de ligadura gravitatoria del disco.
• Supervivencia del Disco: En contraste, los modelos con un disco masivo (M_{ejecta}/M_{disk} = 0.48}, Modelo E) o discos en forma de V con alta densidad inicial constante (Modelo F) muestran una supervivencia parcial del material del disco tras el impacto.
• Pérdida de Masa y Retención:
  • Las altas velocidades de los eyectos (características de las RNe, ej. $10,000$ km s$^{-1}$) impulsan a una fracción significativa de la masa combinada del disco y los eyectos a exceder la velocidad de escape del sistema. En el Modelo A fiduciario, aproximadamente el $97\%$ del material abandona el sistema binario, mientras que solo un $\sim3\%$ permanece gravitacionalmente ligado a la enana blanca.
  • Esto contrasta con hallazgos previos para novas clásicas (Artículo I), donde velocidades menores resultaban en una mayor retención.
  • Aumentar la masa o la velocidad de los eyectos incrementa la masa total perdida del sistema y disminuye la masa acrecentada por la enana blanca.
  • Por el contrario, eyectos más densos (densidad constante vs. perfil) y discos más masivos aumentan la fracción de material retenido por la enana blanca debido a una desaceleración mejorada durante la colisión.
• Interacción con la Estrella Secundaria:
  • El impacto de los eyectos y el material del disco en la compañera subgigante es mínimo. Solo una fracción ínfima de partículas ($\sim0.07\%$ de la masa total de eyectos/disco en el Modelo A) se incorpora en las capas exteriores de la secundaria, lo que sugiere una contaminación química insignificante.
  • Se observa una eyección menor de masa de las capas exteriores de la subgigante durante la colisión en etapa tardía. Parte de esta masa secundaria es expulsada del sistema, mientras que una pequeña fracción es acrecentada por la enana blanca.
• Papel de la Rotación: En la órbita ancha de U Sco, la rotación juega un papel marginal en la pérdida o retención de masa total de la mezcla de eyectos/disco ($<1\%$ de diferencia). Sin embargo, la rotación reduce la cantidad de masa perdida por la subgigante que es posteriormente acrecentada por la enana blanca en aproximadamente un $20\%$.

Significancia y Reivindicaciones
Los autores posicionan este trabajo como un paso cualitativo necesario para llenar el vacío en la comprensión de la evolución post-estallido de las novas recurrentes. El estudio afirma proporcionar el primer análisis sistemático en 3D mediante SPH de la interacción entre los eyectos de la nova, los discos de acreción y las compañeras en un sistema binario en rotación con parámetros específicos para U Sco.

El artículo concluye modestamente que sus hallazgos son cualitativos debido a las limitaciones de resolución (específicamente el bajo conteo de partículas para los eyectos de baja densidad y el disco). No obstante, afirma que estos resultados ofrecen una guía valiosa para investigaciones futuras, particularmente en:

1. Identificar las condiciones (razones de masa, geometría, densidad) que conducen a la supervivencia del disco frente a la disrupción total.
2. Destacar la sensibilidad de la retención de masa a la velocidad de los eyectos, un diferenciador clave entre novas recurrentes y clásicas.
3. Enfatizar la necesidad de modelos detallados de la geometría del disco (ensanchados vs. en forma de V) para simular con precisión la evolución a largo plazo de estos sistemas.

Los autores señalan que, aunque sus simulaciones sugieren la disrupción completa del disco en muchos escenarios estándar, las afirmaciones observacionales de supervivencia del disco (ej. en el estallido de U Sco de 2022) pueden implicar la presencia de discos extendidos o configuraciones físicas diferentes no capturadas totalmente en el espacio de parámetros actual.