The Light Curve of Wind-Reprocessed Tidal Disruption Events

Resumen Técnico: La curva de luz de los eventos de disrupción de marea con reprocesamiento por viento

Planteamiento del problema
El origen de la ubicua emisión óptica y ultravioleta (UV) en los eventos de disrupción de marea (TDE, por sus siglas en inglés) sigue siendo una cuestión central en el campo. Aunque se predice que los TDE serán brillantes en rayos X debido a la acreción de agujeros negros supermasivos (SMBH), la mayoría se descubren mediante sondeos ópticos, exhibiendo a menudo una alta relación óptico-rayos X. Trabajos previos han sugerido que la emisión de alta energía cerca del agujero negro es reprocesada por el gas circundante y los flujos de salida (outflows), pero el modelado detallado de este proceso ha sido limitado. Las simulaciones existentes se centran o bien en la hidrodinámica sin un transporte de radiación dependiente del tiempo detallado, o bien asumen condiciones de estado estacionario que no logran capturar la evolución dramática del flujo de gas y los estados de ionización durante un TDE. En consecuencia, la conexión entre la física de la acreción subyacente y las curvas de luz observadas, particularmente en cuanto a escalas de tiempo y presupuestos energéticos, no se comprende plenamente.

Metodología
Los autores emplean simulaciones de radiación hidrodinámica 1D dependientes del tiempo utilizando el código de transferencia de radiación de Monte Carlo Sedona, extendido para incluir hidrodinámica de malla móvil (moving-mesh). El marco de la simulación modela un TDE fiduciario que involucra un agujero negro de $3 \times 10^6 M_\odot$ que disrumpe una estrella de $1 M_\odot$.

• Fuente: Se asume una fuente caliente, brillante en rayos X/EUV, en el radio de circularización ($R_{source} \approx 10^{13}$ cm), con una luminosidad que evoluciona según la tasa de retorno de masa (escalada a $\sim 2.5 \times$ Eddington en su pico).
• Flujo de salida (Outflow): Un viento es lanzado desde el radio de la fuente a $v_{wind} = 0.01c$, transportando el 75% de la masa que retorna.
• Física: El código resuelve las ecuaciones de transporte de radiación asumiendo el Equilibrio Termodinámico Local (LTE) para las poblaciones de niveles y los estados de ionización, mientras resuelve el no-equilibrio de la energía del gas. Las opacidades incluyen dispersión de electrones, procesos bound-free, free-free y un tratamiento simplificado de las transiciones bound-bound utilizando un parámetro de eficiencia de absorción de líneas ($\epsilon = 0.1$).
• Alcance: La simulación rastrea el ascenso, el pico y el declive inicial del evento durante varios meses, cubriendo la transición de un entorno transparente a uno ópticamente grueso.

Resultos clave

1. Eficiencia de reprocesamiento y retrasos en la curva de luz: Las simulaciones demuestran que un viento de flujo de salida puede reprocesar con éxito la emisión de rayos X y EUV en las bandas ópticas/UV. Un hallazgo crítico es un retraso temporal significativo entre el pico bolométrico y el pico óptico/UV. La curva de luz óptica/UV presenta un retraso respecto a la curva bolométrica de aproximadamente 3 a 4 semanas (específicamente $\sim 20$ días en el modelo). Este retraso corresponde al tiempo requerido para que el material se acumule lo suficiente y construya una capa de reprocesamiento ópticamente gruesa.
2. Escalas de tiempo dependientes de la longitud de onda: Las bandas de mayor energía alcanzan su pico más temprano que las bandas de menor energía. El modelo predice un estallido (flare) de rayos X temprano y de corta duración (que alcanza su pico en $\sim 7$ días) que es absorbido rápidamente a medida que la densidad del viento aumenta. Por el contrario, la luminosidad óptica/UV aumenta a medida que esta radiación de alta energía es reprocesada. Dentro del rango óptico/UV, las bandas más azules (por ejemplo, UVW2) alcanzan su pico ligeramente antes que las bandas más rojas (por ejemplo, la banda $i$), lo cual es consistente con las observaciones de eventos como AT2020ocn.
3. Evolución de la opacidad: La opacidad del viento es altamente dependiente del tiempo y del radio. Al inicio de la simulación, el gas está altamente ionizado, lo que permite que los rayos X escapen. A medida que el viento se expande y se enfría, aparecen estados de ionización más bajos de oxígeno (por ejemplo, OVI, OV) en radios grandes, creando altas opacidades bound-free que absorben los rayos X blandos. Los autores encuentran que las transiciones bound-bound (líneas) se vuelven cada vez más importantes después del pico, aumentando el continuo óptico/UV hasta en un factor de 5 en comparación con los modelos que carecen de opacidades de línea.
4. Presupuesto de energía y pérdidas adiabáticas: La distribución espectral de energía (SED) de la simulación alcanza su máximo en el EUV (ultravioleta extremo), no en el óptico. Por consiguiente, ajustar cuerpos negros a los datos ópticos/UV subestima la energía total emitida por un orden de magnitud o más. El estudio cuantifica las pérdidas adiabáticas debidas a la expansión del viento, encontrando que estas reducen la luminosidad de la fuente de entrada por un factor de $\sim 6$ en el pico. Esta pérdida es mayor que las estimaciones analíticas simples basadas en el radio de atrapamiento, lo que sugiere que una pérdida de energía significativa ocurre incluso fuera del radio de atrapamiento nominal.
5. Irradiación de rayos X: Incluso cuando el continuo de rayos X está completamente oscurecido y es inobservable, las capas externas del viento permanecen irradiadas por fotones de rayos X. Esto implica que las líneas espectrales excitadas por rayos X podrían ser observables incluso en ausencia de un continuo de rayos X detectado.

Significancia y afirmaciones
El artículo afirma que el reprocesamiento por viento proporciona un mecanismo viable para explicar las curvas de luz ópticas/UV observadas en los TDE, incluyendo sus formas y la evolución de sus colores. La principal significancia de este trabajo radica en demostrar que la emisión de alta energía en etapas tempranas puede perderse en los sondeos ópticos porque el pico óptico tiene un retraso respecto al pico bolométrico y a la disrupción inicial.

Los autores argumentan que:

• El "problema de la energía faltante" en los TDE (donde la luminosidad óptica observada es menor de lo esperado para la tasa de acreción) se explica naturalmente porque la SED alcanza su máximo en el EUV no observado y por las pérdidas adiabáticas.
• El retraso temporal entre los picos bolométrico y óptico implica que la disrupción inicial y el retorno de masa al agujero negro ocurren antes de lo que sugieren las curvas de luz ópticas.
• La detección de estallidos tempranos de rayos X (como se ve en AT2022dsb) es consistente con un escenario en el que la capa de reprocesamiento aún no se ha construido hasta alcanzar la densidad óptica.
• Aunque la simulación 1D representa un ángulo de visión intermedio (análogo a vistas "de lado" en 3D), los resultados ofrecen una base para comprender cómo el reprocesamiento dependiente del tiempo moldea los observables, motivando observaciones UV de alta cadencia para restringir las propiedades de la fuente y los presupuestos de energía.

Los autores señalan que sus resultados son una prueba de concepto y que trabajos futuros explorarán los espacios de parámetros (por ejemplo, variando las fracciones de eyección de masa, las luminosidades de la fuente y las geometrías 3D) para refinar estas predicciones.