X-rays from shock-heated gas in recurrent-nova remnants: Nested nova shells in a structured circumstellar medium

Resumen Técnico: Rayos X de Gas Calentado por Choque en Remanentes de Nova Recurrente

Planteamiento del Problema
Las novae simbióticas recurrentes, como RS Ophiuchi (RS Oph), exhiben interacciones complejas entre los eyectos de la nova y el medio circunestelar (CSM). Mientras que los estudios iniciales post-erupción se centran en choques inmediatos en escalas de tiempo de días a meses, la evolución a largo plazo del CSM a lo largo de múltiples ciclos de erupción sigue siendo menos comprendida. Específicamente, el mecanismo que gobierna la emisión de rayos X extendida y difusa observada entre erupciones —distinta de la fuente de energía supersuave en la superficie de la enana blanca— no está completamente caracterizado. Observaciones previas sugieren que las novae recurrentes pueden haber experimentado muchas más erupciones de las registradas históricamente, creando una estructura de "capas anidadas" que influye en la evolución secular del choque y en el destino final del sistema (por ejemplo, como progenitor de una supernova de Tipo Ia).

Metodología
Los autores realizan simulaciones hidrodinámicas tridimensionales utilizando el código de refinamiento de malla adaptativa (AMR) FLASH 4.6 para modelar una secuencia fidedigna de nueve erupciones durante 130 años. La configuración adapta una prescripción de capas anidadas diseñada originalmente para remanentes de supernova al contexto de las novae recurrentes.

• Condiciones Iniciales: El CSM se inicializa como un viento de gigante roja ($\dot{M} = 10^{-7} M_\odot \text{ yr}^{-1}$, $v_\infty = 20 \text{ km s}^{-1}$) con una temperatura de $10^4$ K.
• Modelo de Erupción: Cada una de las nueve erupciones se modela como una capa bipolar supersónica que se expande desde un radio inicial de 0.009 pc. La masa de los eyectos se establece en $2 \times 10^{-6} M_\odot$ con una energía cinética de $4.02 \times 10^{44}$ erg, produciendo una velocidad de expansión de $\sim 4600 \text{ km s}^{-1}$. La morfología depende de la latitud, creando una estructura bipolar.
• Física: El enfriamiento radiativo incluye el enfriamiento de gas en equilibrio de ionización por colisión (CIE) ópticamente delgado y el enfriamiento inducido por polvo on-the-fly (relación polvo-gas $10^{-2}$). El dominio de la simulación abarca $(0.2 \text{ pc})^3$ con un refinamiento AMR de hasta $\sim 3.9 \times 10^{-4}$ pc.
• Síntesis de Rayos X: Los resultados hidrodinámicos se procesan posteriormente utilizando XSPEC v.12.12. Se aplica el modelo de choque plano-paralelo de ionización no equilibrio (NEI) vvpshock para calcular las emisividades en las bandas suaves (0.5–2.0 keV) y duras (2.0–10.0 keV). Los tiempos de ionización se asignan basados en una máscara de presión para distinguir los eyectos recién chocados del material más antiguo y mezclado.

Resultos Clave

1. Evolución Morfológica: La secuencia de erupciones excava una cavidad bipolar delimitada por capas anidadas. La estructura de densidad del gas muestra que los rayos X suaves trazan interfaces densas y comprimidas entre las capas, mientras que los rayos X duros se originan en el gas chocado más caliente, incluyendo los bordes de las capas y el interior del remanente excavado. A medida que el remanente se expande, la morfología de rayos X extendida se vuelve progresivamente más suave y de mayor llenado de volumen a medida que los contrastes de densidad entre sucesivas capas disminuyen.
2. Luminosidad de Rayos X Suaves: La emisión de la banda suave (0.5–2.0 keV) muestra un aumento inicial rápido seguido de un declive secular gradual asociado con la disminución de la medida de emisión a medida que el remanente se expande. Aunque las erupciones subsecuentes causan aumentos de corta duración, la tendencia general es un declive de aproximadamente 2.18 dex durante los 130 años de la simulación. Las luminosidades de rayos X suaves difusas predichas (que oscilan entre $\sim 10^{30}$ y unos pocos $\times 10^{31}$ erg s$^{-1}$ tras el pico inicial) son ampliamente consistentes con el orden de magnitud inferido para el componente extendido de RS Oph.
3. Luminosidad de Rayos X Duros: En contraste con la banda suave, la emisión de la banda dura (2.0–10.0 keV) evoluciona de forma episódica. Las erupciones individuales imprimen llamaradas distintas, con picos de luminosidad que declinan desde $\sim 10^{30}$ erg s$^{-1}$ después de la primera erupción hasta $\sim 10^{27}$ erg s$^{-1}$ en ciclos posteriores. Los valores entre erupciones caen a $\sim 10^{26}$ erg s$^{-1}$. Estas llamaradas corresponden a interacciones transitorias entre los eyectos rápidos y el sistema de capas acumuladas.
4. Comparación con Observaciones: Las tasas de declive de rayos X suaves simuladas y las escalas de luminosidad para los ciclos más brillantes se alinean estrechamente con las observaciones de Chandra de RS Oph (Montez et al. 2022) tomadas 3.4 y 5.3 años después de la erupción de 2006.

Significancia y Reivindicaciones
El artículo concluye que la emisión de choque de capas extendidas es una consecuencia natural de la evolución de capas anidadas en las novae recurrentes. Los autores identifican dos mecanismos reguladores distintos para la variabilidad de rayos X:

• La expansión del remanente regula el desvanecimiento lento del componente de rayos X suaves.
• Los nuevos encuentros entre eyectos y capas impulsan la variabilidad de la banda dura en escalas de tiempo de décadas.

El estudio postula que el gas calentado por choque puede contribuir apreciablemente al entorno observado de rayos X en las novae recurrentes entre erupciones. Además, el sistema de capas acumuladas crea un entorno circunestelar físicamente evolucionado que daría forma a la propagación del choque de cualquier evento explosivo futuro, incluyendo una potencial supernova de Tipo Ia. Los autores señalan que su modelo no incluye campos magnéticos ni conducción térmica anisotrópica, y que la simulación cubre solo nueve ciclos, mientras que sistemas como RS Oph pueden haber experimentado muchos más.