Self-Consistent Modelling of Neutrino Production in Turbulent Black Hole Coronae

Título: Modelado Autoconsistente de la Producción de Neutrinos en Coronas de Agujeros Negros Turbulentos

1. El Problema

En 2022, la colaboración IceCube reportó un exceso estadístico significativo de neutrinos de alta energía (1-10 TeV) asociados a la galaxia Seyfert 2 NGC 1068. Este hallazgo sugiere que los núcleos galácticos activos (AGN) no jeteados son fuentes importantes de neutrinos extragalácticos. Sin embargo, existe una discrepancia crítica: la emisión de rayos gamma observada es mucho menor que la inferida para los neutrinos, lo que implica que los neutrinos se producen en una región compacta y opaca a fotones de alta energía (debido a la absorción $\gamma\gamma$).

El desafío principal para modelar este fenómeno radica en la autoconsistencia dinámica. Los modelos existentes a menudo utilizan aproximaciones estáticas o separan el proceso de aceleración de las pérdidas radiativas y la retroalimentación sobre el medio. Específicamente, se requiere un tratamiento que considere simultáneamente:

• La aceleración estocástica de partículas en plasmas turbulentos.
• La retroalimentación de las partículas aceleradas sobre la cascada de turbulencia (amortiguamiento de la turbulencia).
• Las pérdidas radiativas y hadrónicas, incluyendo el reprocesamiento de la radiación electromagnética en entornos densos.
• La evolución temporal de las distribuciones de partículas, ya que los protones evolucionan en escalas de tiempo macroscópicas, a diferencia de los electrones.

2. Metodología: El Código Híbrido Turb-AM3

Los autores presentan Turb-AM3, un código numérico híbrido diseñado para abordar estas complejidades. Su arquitectura combina dos componentes principales:

• Solucionador Radiativo Lepto-Hadrónico (AM3): Un código de estado del arte que resuelve las ecuaciones de transporte para fotones, neutrinos, electrones, positrones y partículas secundarias (piones, muones). Calcula pérdidas radiativas (Compton inverso, sincrotrón), hadrónicas ($pp$, $p\gamma$) y la producción de cascadas electromagnéticas.
• Módulo de Aceleración Estocástica: Acoplado dinámicamente a AM3, este módulo modela la aceleración de protones en plasmas turbulentos magnetizados.
  • Ecuación de Transporte: Utiliza una ecuación de transporte en el espacio de momentos que incluye operadores de difusión (Fokker-Planck) y un operador generalizado de Fermi (LGF) para regímenes de turbulencia fuerte, basado en simulaciones recientes de PIC (Partícula en Celda).
  • Acoplamiento No Lineal: Implementa un mecanismo de amortiguamiento de la turbulencia. La energía extraída por los protones acelerados de la cascada turbulenta reduce la potencia turbulenta disponible, lo que a su vez disminuye la tasa de aceleración. Esto crea un mecanismo de autorregulación.
  • Escalas de Tiempo: El modelo integra la evolución temporal hasta el tiempo de advección ($t_{adv}$), considerando la inyección de protones, la difusión espacial y la pérdida por advección hacia el agujero negro.

3. Contribuciones Clave

• Desarrollo de Turb-AM3: La creación de un marco numérico capaz de tratar la evolución temporal acoplada de la aceleración de protones y la disipación de turbulencia, algo que las simulaciones PIC no pueden alcanzar debido a la enorme gama dinámica de energías requerida (5 órdenes de magnitud).
• Modelado Autoconsistente de Retroalimentación: La inclusión explícita de cómo la densidad de energía de los protones no térmicos afecta la magnetización y la tasa de aceleración, evitando soluciones no físicas donde la energía de las partículas excedería la disponible en el campo magnético.
• Predicción de Formas Espectrales: Generación de espectros de neutrinos "plantilla" que no son meras soluciones de estado estacionario, sino promedios temporales que reflejan la dinámica real de la corona.

4. Resultados Principales

• Reproducción de NGC 1068: Con un conjunto estándar de parámetros físicos ($R_{cor} \sim 15 r_g$, $v_A \sim 0.25c$, $\ell_c \sim 4 r_g$, $\tau_T \sim 0.5$), el modelo reproduce satisfactoriamente el flujo de neutrinos observado por IceCube en el rango de 1-10 TeV sin violar las restricciones de rayos gamma.
• Mecanismo de Autorregulación: Se encuentra que la densidad de energía de los protones no térmicos se estabiliza cerca de la densidad de energía turbulenta debido al amortiguamiento. Esto explica naturalmente la normalización del flujo de neutrinos sin necesidad de ajustar parámetros arbitrarios.
• Estructura del Espectro de Neutrinos:
  • El espectro predicho muestra un pico alrededor de 1-10 TeV.
  • Presenta una extensión plana a bajas energías (dominada por interacciones $pp$), que es una firma distintiva de la autorregulación de la turbulencia.
  • Un corte a altas energías (>10 TeV) impuesto por las pérdidas foto-hadrónicas ($p\gamma$).
• Aplicación a Otros AGN: El modelo se extiende exitosamente a otras galaxias Seyfert (NGC 4151 y NGC 7469). Al ajustar únicamente la luminosidad del disco (y por tanto la densidad de fotones objetivo), el modelo reproduce los niveles de emisión de neutrinos y la tendencia observada de que fuentes con neutrinos de mayor energía tienden a tener flujos más bajos.
• Dependencia de Parámetros: El estudio de sensibilidad revela que la forma espectral es robusta frente a variaciones en la fracción de inyección ($\xi_p$) siempre que sea suficientemente alta ($\gtrsim 10^{-2}$), pero es muy sensible a la relación entre los tiempos de escape y aceleración ($t_{esc}/t_{acc}$).

5. Significado e Impacto

• Validación del Escenario de Corona: Los resultados apoyan fuertemente la hipótesis de que las coronas de agujeros negros turbulentos son los sitios de producción de neutrinos en AGN no jeteados, resolviendo la discrepancia entre neutrinos y rayos gamma mediante la opacidad de la región de producción.
• Herramienta para Futuras Observaciones: El modelo proporciona espectros "plantilla" detallados que mejorarán la sensibilidad de los futuros telescopios de neutrinos y permitirán distinguir entre diferentes mecanismos de aceleración.
• Nueva Física de Plasmas: Demuestra que la interacción no lineal entre partículas aceleradas y turbulencia es crucial para entender la emisión multi-mensajero, ofreciendo una nueva vía para inferir las condiciones del plasma (como el factor de carga de pares) a partir de observaciones de neutrinos.
• Escalabilidad: El marco Turb-AM3 es aplicable a otros entornos astrofísicos compactos, como binarias de rayos X, eventos de disrupción de marea y chorros de AGN, consolidándose como una herramienta estándar para la astrofísica de multi-mensajeros.

En resumen, este trabajo establece un nuevo estándar en el modelado de fuentes de neutrinos de alta energía, demostrando que la aceleración estocástica autoconsistente en coronas turbulentas puede explicar las observaciones de IceCube de manera física y robusta.