Impacts of radiative cooling on the images of a black hole shadow and extended jets in two-temperature GRMHD simulations

Resumen Técnico: Impactos del Enfriamiento Radiativo en las Imágenes de la Sombra de un Agujero Negro y Chorros Extendidos en Simulaciones GRMHD de Dos Temperaturas

Planteamiento del Problema
Las observaciones recientes del Telescopio del Horizonte de Sucesos (EHT) de la sombra del agujero negro M87* han revelado una eficiencia radiativa en el flujo de acreción caliente que sugiere que el enfriamiento radiativo no es despreciable. Si bien las simulaciones de Magnetohidrodinámica Relativista General (GRMHD) son herramientas estándar para modelar la física de acreción, muchos modelos existentes omiten el enfriamiento radiativo o se basan en aproximaciones de un solo fluido con prescripciones paramétricas de temperatura electrónica (por ejemplo, el modelo "R-β"). Además, estudios anteriores han asumido frecuentemente funciones de distribución electrónica térmicas (eDF), ignorando potencialmente procesos no térmicos impulsados por la reconexión magnética y la turbulencia. Este estudio aborda la brecha en la comprensión de cómo el enfriamiento radiativo, cuando se acopla a un modelo de dos temperaturas impulsado físicamente y a eDFs no térmicas, altera la termodinámica electrónica, la estructura del disco y las imágenes sintéticas resultantes de las sombras de agujeros negros y los chorros extendidos.

Metodología
Los autores realizaron simulaciones de Magnetohidrodinámica Relativista General (GRMHD) de flujos de acreción magnetizados hacia un agujero negro en rotación ($a = 0.9375$) utilizando el código BHAC. Las simulaciones emplearon un marco de dos temperaturas donde la termodinámica electrónica evoluciona por separado de la de los iones, incorporando:

• Mecanismos de Calentamiento: Se probaron dos prescripciones distintas: calentamiento turbulento (Kawazura et al. 2019) y calentamiento por reconexión magnética (Rowan et al. 2017).
• Mecanismos de Enfriamiento: Se incluyeron el enfriamiento radiativo (bremsstrahlung, ciclotrón-sincrotrón y dispersión Compton inversa) y las interacciones de Coulomb entre protones y electrones como términos fuente/sumidero.
• Tasas de Acreción: Las simulaciones se ejecutaron en un rango de tasas de acreción de masa normalizadas a la tasa de Eddington: $\dot{m} = (1.0 - 10) \times 10^{-6} \dot{M}_{Edd}$.
• Transferencia Radiativa: Se realizaron cálculos de Transferencia Radiativa Relativista General (GRRT) utilizando el código BHOSS para generar imágenes a 230 GHz con una inclinación de $163^\circ$.
• Distribución Electrónica: El estudio utilizó una función de distribución electrónica híbrida que combina una distribución Maxwell-Jüttner térmica y una distribución no térmica $\kappa$ variable. El parámetro $\kappa$ se definió paramétricamente en función de la magnetización ($\sigma$) y el beta del plasma ($\beta$) utilizando modelos derivados de simulaciones de Partículas en Celda (PIC) (Ball et al. 2018; Meringolo et al. 2023).

Contribuciones y Resultados Clave
Al comparar modelos con y sin enfriamiento radiativo a través de diferentes tasas de acreción y prescripciones de calentamiento, el estudio arroja los siguientes hallazgos:

1. Estructura de Temperatura y Densidad Electrónica: El enfriamiento radiativo disminuye drásticamente la temperatura electrónica en el disco interno ($r \lesssim 10 r_g$) alrededor del plano ecuatorial, donde la densidad es más alta. En contraste, la temperatura en la vaina del chorro disminuye solo ligeramente. La inclusión del enfriamiento conduce a la formación de una estructura de disco más delgada a tasas de acreción más altas debido al aumento de la eficiencia de enfriamiento.
2. Morfología de la Imagen y Flujo: El enfriamiento radiativo resulta en un anillo de acreción (sombra) más tenue y una reducción en el flujo total emitido. Sin embargo, produce estructuras de chorro más extendidas y relativamente más brillantes. La contribución del chorro del lado cercano a la imagen total aumenta en los modelos con enfriamiento en comparación con los modelos sin enfriamiento.
3. Distribuciones de Energía Espectral (SEDs): La inclusión del enfriamiento desplaza la frecuencia de inversión (pico de la SED) hacia frecuencias más bajas para una tasa de acreción dada, lo que indica que la radiación sincrotrón se vuelve menos autoabsorbida. A 230 GHz, los modelos sin enfriamiento a tasas de acreción más altas se vuelven ópticamente gruesos, mientras que los modelos con enfriamiento permanecen en el régimen ópticamente delgado.
4. Variabilidad Temporal: El estudio encuentra que la variabilidad temporal en las curvas de luz a 230 GHz se origina principalmente en la región del plano medio tanto para los casos con enfriamiento como sin enfriamiento. A medida que aumenta la tasa de acreción de masa, la variabilidad temporal (índice de modulación) disminuye para ambos escenarios.
5. Prescripciones de Calentamiento: Bajo la misma tasa de acreción, los modelos de calentamiento por reconexión producen generalmente flujos totales más altos que los modelos de calentamiento turbulento, tanto con como sin enfriamiento. Sin embargo, el impacto cualitativo del enfriamiento radiativo (disco más tenue, chorros más brillantes) permanece consistente a través de las prescripciones de calentamiento.

Significado y Afirmaciones
Los autores afirman que incorporar el enfriamiento radiativo en simulaciones GRMHD de dos temperaturas con distribuciones electrónicas no térmicas es crucial para modelar con precisión la sombra del agujero negro y la estructura del chorro a gran escala de M87*. El estudio demuestra que omitir el enfriamiento puede conducir a modelos incompletos o predicciones inconsistentes con respecto a la temperatura electrónica y la variabilidad del flujo.

El artículo concluye que, si bien las características específicas identificadas (como la visibilidad mejorada del chorro extendido y reducciones específicas de flujo) están actualmente por debajo del rango dinámico de las observaciones actuales del EHT, representan predicciones comprobables para arrays de próxima generación como el ngEHT. Los resultados se presentan como generalmente aplicables a otros núcleos galácticos activos de baja luminosidad (LLAGNs) en el estado de Disco Arrestado Magnéticamente (MAD), sugiriendo que el enfriamiento radiativo conduce universalmente a discos más fríos, chorros más brillantes, flujo total reducido y curvas de luz más suaves, independientemente de la masa del agujero negro y la distancia. Los autores mantienen una postura modesta, señalando que sus simulaciones aún no están ajustadas para ajustar perfectamente los datos de M87* y que el trabajo futuro debe integrar de manera más autoconsistente los términos de enfriamiento para electrones no térmicos para reproducir completamente los espectros de banda ancha.