Dust Formation in Common Envelope Binary Interactions -- III. Lightcurves

Resumen Técnico: Formación de Polvo en Interacciones de Envoltura Común de Binarias — III. Curvas de Luz

Planteamiento del Problema
Las Novas Rojas Luminosas (LRNe) son eventos transitorios resultantes de interacciones de envoltura común (CE) en sistemas binarios, pero predecir sus curvas de luz observables a partir de simulaciones hidrodinámicas sigue siendo un desafío. Los intentos previos de modelar estos eventos a menudo trataban la formación de polvo como un proceso instantáneo o dependían de tablas de opacidad que no resolvían explícitamente la escala de tiempo de nucleación. Además, cerrar la brecha entre las simulaciones hidrodinámicas 3D y los datos observacionales requiere métodos de transferencia radiativa de post-procesamiento robustos que den cuenta de la naturaleza compleja, asimétrica y evolutiva de los eyectos de la CE, particularmente la formación de una cáscara de polvo ópticamente gruesa.

Metodología
Los autores realizan cálculos de curvas de luz de post-procesamiento en dos simulaciones 3D de hidrodinámica de partículas suavizadas (SPH) de interacciones de CE, descritas originalmente en Bermúdez-Bustamante et al. (2024). Las simulaciones modelan la interacción entre una compañera compacta ($0.6 M_\odot$) y estrellas de la rama asintótica de las gigantes (AGB) de dos masas distintas: $1.7 M_\odot$ y $3.7 M_\odot$. Crucialmente, estas simulaciones incluyen la nucleación explícita de polvo basado en carbono.

La metodología para derivar las curvas de luz implica:

1. Transferencia Radiativa: Aplicación de la ecuación de transferencia radiativa a lo largo de los rayos proyectados desde el dominio de la simulación hacia un observador. La intensidad específica se calcula utilizando una aproximación de opacidad gris asumiendo equilibrio termodinámico local (LTE).
2. Tratamiento de la Opacidad: La opacidad total es la suma de la opacidad del gas (derivada de las tablas de MESA para $T > 1450$ K) y la opacidad del polvo (calculada mediante la teoría de Mie para el polvo de carbono). Para temperaturas por debajo de 1450 K, donde se forma el polvo, la opacidad del gas se establece en un valor constante bajo, permitiendo que la opacidad del polvo domine.
3. Cálculo de la Luminosidad: La luminosidad bolométrica se deriva integrando la intensidad específica sobre una rejilla de $256 \times 256$ píxeles. La función fuente se trata como un cuerpo negro gris.
4. Resolución y Análisis de Incertidumbre: Los autores cuantifican explícitamente las incertidumbres derivadas de la baja resolución superficial (donde pocas partículas SPH definen la fotosfera) y la suposición adiabática (falta de enfriamiento radiativo). Introducen una métrica de resolución ($N_{res}$) basada en el número de partículas que contribuyen a la integración de la luminosidad por rayo para definir las barras de error.

Resultos Clave
Las simulaciones, ejecutadas durante 44 años bajo condiciones adiabáticas, arrojan las siguientes características de la curva de luz:

• Fase Temprana (Ascenso al Pico): Ambos modelos exhiben un pico brillante y caliente ($L \sim 2-4 \times 10^6 L_\odot$, $T_{eff} \sim 6000-7000$ K) que ocurre entre 1 y 5 años después del inicio del desbordamiento del lóbulo de Roche. Este ascenso es impulsado principalmente por la rápida expansión de la fotosfera debido a la espiral hacia adentro (inspiral) e inestabilidades superficiales, más que por cambios significativos de temperatura. El tiempo y la forma del pico varían según el ángulo de visión y la masa de la progenitora; el modelo más masivo ($3.7 M_\odot$) muestra un pico más ancho e irregular debido a una razón de masa menor ($q=0.16$) y una formación de polvo retrasada.
• Formación de Polvo y Declive: La nucleación de polvo comienza aproximadamente 1–3 años después de que la interacción inicia. Poco después, la luminosidad bolométrica experimenta un declive brusco a medida que la cáscara de polvo se vuelve ópticamente gruesa, oscureciendo la fotosfera de gas caliente. Esto es seguido por una recuperación parcial y una fase de meseta con una temperatura efectiva de $\sim 400$ K.
• Evolución de la Fotosfera: La fotosfera "geométrica" (definida por $\tau = 2/3$) se expande linealmente, alcanzando $\sim 250$ au al final de la simulación. Sin embargo, la fotosfera del "observador" (derivada del ajuste de SED) revela una evolución compleja donde un componente más frío, dominado por el polvo, eventualmente supera en brillo al componente de gas más caliente.
• Dispersión a Largo Plazo: La extrapolación de la expansión libre de la cáscara de polvo sugiere que el polvo se volverá ópticamente delgado en longitudes de onda visibles entre 100 y 200 años post-interacción, revelando nuevamente la fotosfera interna más cálida.
• Comparación con Observaciones: Los autores comparan sus resultados con OGLE-2002-BLG-360 y AT 2025abao. El modelo de $3.7 M_\odot$ proporciona una coincidencia más cercana con AT 2025abao en términos de masa de la progenitora y luminosidad pico, aunque el tiempo de ascenso observado para AT 2025abao es significativamente más rápido de lo que predice la simulación. El modelo de $1.7 M_\odot$ comparte algunas características con OGLE-2002, particularmente respecto a la masa de polvo y el radio fotosférico en tiempos tardíos.

Limitaciones e Incertidumbres
El artículo identifica dos fuentes primarias de incertidumbre que limitan la precisión cuantitativa de las curvas de luz:

1. Resolución Superficial: Los primeros 1–2 años de las curvas de luz están fuertemente afectados por la resolución insuficiente de SPH en la superficie estelar. La resolución actual (1.3 millones de partículas) resulta en solo 1–2 partículas definiendo la fotosfera en $t=0$, lo que genera grandes incertidumbres en la luminosidad inicial y el tiempo de ascenso. Los autores señalan que lograr una resolución convergente requeriría cientos de millones o miles de millones de partículas.
2. Suposición Adiabática: Las simulaciones son adiabáticas, careciendo de enfriamiento radiativo. Por consiguiente, la luminosidad de la meseta a largo plazo ($\sim 10^5 L_\odot$) es probablemente sobreestimada por un orden de magnitud respecto a las expectativas físicas, ya que el sistema no puede irradiar el exceso de energía interna.

Significado
Este trabajo representa un primer paso para cerrar la brecha entre las simulaciones hidrodinámicas de CE en 3D con nucleación de polvo explícita y las curvas de luz observables. Aunque los autores reconocen que sus modelos aún no están completamente convergidos o físicamente completos (debido a la suposición adiabática), demuestran con éxito que:

• La formación de polvo en las interacciones de CE es un proceso no instantáneo que altera fundamentalmente la morfología de la curva de luz, creando una fase distintiva de declive y meseta.
• La asimetría de la interacción conduce a dependencias significativas del ángulo de visión en la curva de luz observada.
• Es posible rastrear el tiempo de formación de polvo y la subsiguiente transición de una fotosfera dominada por gas a una dominada por polvo, proporcionando un marco para interpretar futuras observaciones de LRN.

El artículo concluye que, si bien las limitaciones actuales impiden un emparejamiento cuantitativo preciso con transitorios específicos, las características cualitativas de las curvas de luz (pico, declive, meseta y escalas de tiempo de dispersión de polvo) ofrecen una plantilla valiosa para comprender la evolución de los eventos de CE y la formación de nebulosas planetarias. Las mejoras futuras requerirán un refinamiento de partículas adaptativo para resolver la superficie y la implementación de enfriamiento radiativo para corregir la energética de tiempos tardíos.