Formation and X-ray emission from hot bubbles in planetary nebulae - III. The impact of [Wolf-Rayet]-type winds

Resumen Técnico: Formación y Emisión de Rayos X de Burbujas Calientes en Nebulosas Planetarias - III. El Impacto de los Vientos de tipo [Wolf-Rayet]

Planteamiento del Problema
Las nebulosas planetarias (PN) exhiben burbujas calientes que emiten rayos X, formadas por la interacción de vientos rápidos de la fase post-Asintótica de la Rama Gigante (post-AGB) con eyecciones de la fase AGB, que son más lentas y densas. Mientras que los modelos estándar asumen estrellas centrales ricas en hidrógeno (H-ricas), una fracción significativa de las PN que emiten rayos X alberga estrellas centrales deficientes en hidrógeno, de tipo [Wolf-Rayet] ([WR]). Estudios previos sugirieron que la alta metalicidad (rica en C) de los vientos [WR] potencia el enfriamiento radiativo, lo que potencialmente retrasa la formación de la burbuja y altera las estructuras térmicas. Sin embargo, los modelos existentes a menudo dependieron de simulaciones 1D con aumentos artificiales del viento o asumieron la conducción térmica como el mecanismo principal para enfriar el plasma caliente a las temperaturas de rayos X suaves observadas ($T_X \approx 1-3 \times 10^6$ K). El impacto específico de los parámetros de viento [WR] consistentes con uno mismo, derivados de prescripciones empíricas recientes de pérdida de masa, sobre la evolución hidrodinámica y la emisión de rayos X sintética de las burbujas calientes aún debe ser plenamente explorado en un contexto multidimensional.

Metodología
Los autores emplean simulaciones de radiación-hidrodinámica utilizando el código PLUTO para modelar el medio circunestelar de las PN. La metodología involucra un enfoque de múltiples etapas:

1. Evolución Estelar: Se generan modelos de evolución para progenitores de $1, 2$y$3 M_\odot$ utilizando MESA. Se computan dos trayectorias evolutivas para cada masa: una trayectoria post-AGB estándar rica en H (utilizando la prescripción de pérdida de masa de Miller Bertolami 2016) y una trayectoria [WR] utilizando una prescripción reciente de tasa de pérdida de masa para estrellas deficientes en H derivada de Toalá, Todt & Sander (2024).
2. Configuración Hidrodinámica:
   • Simulaciones 1D: Realizadas en mallas radiales para estudiar la formación y evolución global de la burbuja, incorporando curvas de enfriamiento específicas para abundancias ricas en H y deficientes en H (ricas en C) generadas mediante Cloudy.
   • Simulaciones 2D: Realizadas en mallas esféricas para capturar inestabilidades hidrodinámicas (capa delgada de Vishniac y Rayleigh-Taylor) en la interfaz de interacción viento-viento. Estas simulaciones rastrean la mezcla de material AGB y del viento post-AGB utilizando un parámetro escalar advectado.
3. Síntesis de Rayos X: Los espectros sintéticos de rayos X (0.3–2.0 keV) se calculan mediante post-procesamiento usando la base de datos atómica CHIANTI. La Medida de Emisión Diferencial (DEM) se calcula por separado para las fases de gas ricas en H y pobres en H basándose en el escalar advectado, lo que permite el cálculo de las temperaturas de plasma ($T_X$) y luminosidades ($L_X$) ponderadas por emisividad que contabilizan la composición mixta de la burbuja.

Contribuciones Clave

• Evolución [WR] Consistente con uno Mismo: A diferencia de trabajos anteriores que aumentaban artificialmente los parámetros del viento para una única masa estelar, este estudio utiliza trayectorias MESA con tasas de pérdida de masa [WR] consistentes con uno mismo a través de un rango de masas progenitoras ($1-3 M_\odot$).
• Mezcla Multidimensional: El estudio evita explícitamente la conducción térmica, demostrando en su lugar que las inestabilidades hidrodinámicas en simulaciones 2D generan naturalmente capas de mezcla turbulenta que regulan las temperaturas de las burbujas.
• Análisis de Rayos X de Doble Componente: Los autores introducen un método riguroso para separar y recombinar las contribuciones de rayos X de las fases de gas ricas en H y pobres en H, reconociendo que las burbujas [WR] contienen una mezcla de materiales químicamente distintos.

Resultos

• Retraso en la Formación de la Burbuja: El aumento del enfriamiento radiativo en los vientos [WR] ricos en C retrasa significativamente la formación de burbujas calientes en comparación con sus contrapartes ricas en H. Para un modelo de $1 M_\odot$, la formación de la burbuja se retrasa aproximadamente 6,000 años tras el inicio de la fase post-AGB, mientras que los modelos de mayor masa muestran retrasos menos dramáticos. Este retraso se atribuye al requisito de que la velocidad del viento debe superar un umbral crítico ($v_{crit}$) para vencer el rápido enfriamiento radiativo.
• Luminosidad de Rayos X Mejorada: Los modelos que incorporan vientos [WR] exhiben luminosidades de rayos X ($L_X$) significativamente mayores que los modelos ricos en H. Esto es impulsado por la mayor luminosidad mecánica del viento ($L_{wind}$) resultante de las elevadas tasas de pérdida de masa de las estrellas [WR].
• Convergencia de la Temperatura del Plasma: A pesar de la mayor potencia del viento y las diferentes composiciones químicas, la temperatura del plasma ponderada por emisividad ($T_X$) del gas que emite rayos X converge a valores de $1-3 \times 10^6$ K tanto para los modelos [WR] como para los ricos en H una vez que se establece la mezcla.
• Papel de la Mezcla vs. Conducción Térmica: Las simulaciones 2D muestran que la mezcla turbulenta transporta eficientemente material nebular más frío hacia la burbuja caliente, regulando la temperatura a los niveles observados sin requerir conducción térmica. Las comparaciones con modelos 1D revelan que las simulaciones 1D sobreestiman sistemáticamente $T_X$ (por factores de 1.2–1.6, hasta ~5 en modelos masivos) debido a la ausencia de este mecanismo de mezcla.
• Brillo Superficial: El brillo superficial de rayos X alcanza su máximo temprano y declina rápidamente (por $\sim 4$ órdenes de magnitud) a medida que la burbuja se expande, lo que sugiere que la emisión difusa detectable de rayos X está restringida a una ventana evolutiva específica donde la burbuja se ha expandido lo suficiente pero antes de que se enfríe y se diluya demasiado.

Significado y Reivindicaciones
El artículo afirma reforzar la hipótesis de que la mezcla hidrodinámica es un mecanismo clave para generar la emisión de rayos X suaves en las PN, incluso para sistemas que albergan estrellas centrales [WR]. Los autores aseveran que sus resultados:

1. Explican la ausencia de emisión de rayos X en las PN [WR] jóvenes y compactas como una consecuencia del retraso en la formación de la burbuja debido al fuerte enfriamiento radiativo.
2. Demuestran que las temperaturas de rayos X suaves observadas ($T_X \approx 1-3 \times 10^6$ K) son un resultado natural de la mezcla multidimensional, independientemente de si la estrella central sigue una trayectoria evolutiva rica en H o [WR].
3. Proporcionan una alternativa robusta a los modelos de conducción térmica, sugiriendo que la emisión de rayos X suaves puede surgir incluso en presencia de campos magnéticos que podrían suprimir la conducción térmica electrónica.
4. Validan el uso de prescripciones empíricas de pérdida de masa [WR] en el modelado hidrodinámico, mostrando que reproducen las propiedades del viento estelar observadas ($L, T_{eff}, \dot{M}, v_\infty$) de las CSPNe que emiten rayos X.

Los autores concluyen que, aunque su aproximación de viento es práctica, el trabajo futuro debería modelar explícitamente el escenario de "renacimiento" (born-again) para capturar plenamente la transición a vientos deficientes en H. Sin embargo, el enfoque actual proporciona un marco físicamente motivado para comprender las propiedades de rayos X de las PN [WR].