Coupled Time-Dependent Proton Acceleration and Leptonic-Hadronic Radiation in Turbulent Supermassive Black Hole Coronae

Resumen Técnico: Aceleración de Protones Dependiente del Tiempo y Radiación Leptónica-Hadrónica Acoplada en Coronas de Agujeros Negros Supermasivos Turbulentas

Planteamiento del Problema
Las coronas turbulentas que rodean a los agujeros negros supermasivos (SMBH) son sitios potenciales para la aceleración de partículas no térmicas a altas energías, generando radiación multimensajero observable (neutrinos y ondas electromagnéticas). Sin embargo, modelar estos procesos es un desafío debido a las escalas de tiempo comparables de la aceleración de partículas, el enfriamiento radiativo y el desarrollo de cascadas electromagnéticas (EM). Se requiere un tratamiento consistente que incorpore la retroalimentación de las cascadas EM sobre el enfriamiento y la aceleración de protones, particularmente en regiones coronales compactas donde las densidades de fotones son altas. Aunque existen modelos de estado estacionario, hay una falta de marcos sistemáticos que acoplen la aceleración de protones dependiente del tiempo con la retroalimentación radiativa a través de múltiples escalas espaciales y temporales, especialmente para eventos transitorios como los eventos de disrupción de marea (TDE).

Metodología
Los autores presentan un marco numérico que acopla de manera consistente la aceleración de protones dependiente del tiempo con el modelado de radiación leptónica-hadrónica.

• Aceleración de Protones: La evolución de la función de distribución de protones $f(p, t)$ está gobernada por la ecuación de Fokker-Planck (FP) dependiente del tiempo. Esta ecuación incorpora términos de difusión de momento (aceleración), escape y enfriamiento. Los autores emplean el método de Chang-Cooper para la discretización espacial y el esquema de Crank-Nicolson para la discretización temporal para asegurar la estabilidad y la positividad, incluso en regímenes dominados por el enfriamiento.
• Modelado de Radiación: El resolvedor FP se acopla con AM3, un software de código abierto para el modelado de radiación leptónica-hadrónica dependiente del tiempo. AM3 calcula la evolución de las partículas secundarias (fotones, neutrinos, $e^\pm$, $\mu^\pm$, $\pi^\pm$, $\pi^0$) resultantes de interacciones $p\gamma$, $pp$, $\gamma\gamma$, sincrotrón de protones y Bethe-Heitler.
• Mecanismo de Acoplamiento: Para evitar la doble cuenta del enfriamiento, la ecuación FP determina exclusivamente la evolución del espectro de protones. En cada paso de tiempo, AM3 calcula las tasas de enfriamiento ($t_{cool}^{-1} = t_{p\gamma}^{-1} + t_{pp}^{-1} + t_{BH}^{-1} + t_{p-sy}^{-1}$) basadas en la distribución actual de protones y los campos de fotones. Estas tasas se retroalimentan al resolvedor FP para actualizar la evolución de los protones.
• Validación: El código se valida contra soluciones de estado estacionario (ecuaciones diferenciales ordinarias reducidas) y se prueba con funciones de inyección y perfiles de aceleración/enfriamiento específicos para asegurar la convergencia e independencia de los detalles de la inyección en el régimen de alta energía.

Aplicaciones Clave y Resultados
El marco se aplica a dos escenarios astrofísicos distintos:

1. Corona Estacionaria en NGC 1068:

   • Configuración: El modelo simula una corona magnetizada turbulenta alrededor del SMBH en la galaxia Seyfert NGC 1068, una fuente conocida de neutrinos de alta energía.
   • Resultados: La simulación dependiente del tiempo reproduce el espectro de neutrinos de IceCube (con un pico entre 1–100 TeV) mientras satisface las restricciones de rayos gamma impuestas por Fermi-LAT. El modelo demuestra que el sistema converge a un estado estacionario dentro de $\sim 300$ escalas de tiempo de caída libre ($t_{fs}$).
   • Hallazgos: En este entorno de alta luminosidad de rayos X, el enfriamiento de protones está dominado por las interacciones con los fotones de rayos X de la corona. Se encontró que la retroalimentación de las cascadas EM sobre la distribución de protones es insignificante, y la solución dependiente del tiempo se alinea estrechamente con los resultados previos de estado estacionario.

2. Corona Transitoria en el TDE AT 2019dsg:

   • Configuración: Se desarrolla un modelo fenomenológico para una corona transitoria impulsada por un TDE, caracterizada por una tasa de caída de masa $\propto t^{-5/3}$. El modelo vincula la tasa de acreción dependiente del tiempo con las propiedades de la corona (densidad, opacidad).
   • Resultados: La simulación predice emisiones multilongitud de onda y de neutrinos retrasadas. Específicamente, la retroalimentación de la cascada EM en etapas tempranas impacta significativamente el enfriamiento de protones en coronas más débiles (donde la luminosidad de rayos X es baja).
   • Hallazgos: La acumulación de cascadas EM, combinada con la acumulación de protones durante la fase super-Eddington, impulsa las curvas de luz óptica/UV, de rayos X y de neutrinos a alcanzar su máximo con un retraso de aproximadamente $O(100)$ días (específicamente $\sim 140$ días en el modelo). Esto ofrece una explicación potencial para la emisión de neutrinos con retraso temporal observada en los TDE.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera descripción sistemática que acopla la aceleración de protones dependiente del tiempo (vía Fokker-Planck) con la retroalimentación radiativa de manera totalmente dependiente del tiempo.

• Versatilidad: El marco se presenta como una herramienta flexible capaz de modelar tanto fuentes estables (ej. NGC 1068, otras galaxias Seyfert emisoras de neutrinos) como fenómenos transitorios (ej. TDE, brotes de rayos gamma, flares de AGN).
• Agnosticismo de Mecanismo: Aunque la aplicación principal se centra en la aceleración turbulenta, los autores señalan que el marco puede acomodar otros mecanismos como la reconexión magnética, flujos de cizalladura y aceleración de choque mediante la modificación de los términos de difusión y deriva.
• Perspectiva Multimensajero: Este trabajo une la dinámica de plasma (microfísica de aceleración) y el modelado de radiación, ofreciendo una herramienta para interpretar fenómenos astrofísicos de alta energía en el contexto de observatorios próximos como IceCube-Gen2 y KM3NeT.
• Limitaciones Moderadas: Los autores reconocen que el modelo actual de aceleración de protones depende de parametrizaciones existentes (ej. coeficientes de difusión derivados de simulaciones PIC) y aún no modela de forma consistente la retroalimentación de los protones acelerados sobre la cascada de turbulencia misma, ni el origen detallado de los rayos X de la corona. Sugieren estos puntos como direcciones para futuras refinaciones.