Proton Acceleration by Collisionless Shocks in Supermassive Black Hole Coronae: Implications for High-Energy Neutrinos

Resumen Técnico: Aceleración de Protones mediante Choques Sin Colisiones en Coronas de Agujeros Negros Supermasivos

Planteamiento del Problema
Observaciones recientes del Observatorio de Neutrinos IceCube han detectado un exceso significativo de neutrinos de alta energía provenientes de núcleos galácticos activos (AGN) cercanos y de baja luminosidad de radio, notablemente la galaxia Seyfert NGC 1068. Estas observaciones presentan una discrepancia: mientras que se detectan neutrinos de alta energía, el flujo correspondiente de rayos gamma de alta energía es significativamente menor de lo predicho por los modelos de interacción hadrónica estándar. Esto sugiere que los neutrinos se producen en una región opaca para los rayos gamma, probablemente en la corona caliente que rodea al agujero negro supermasivo (SMBH) central, donde los rayos gamma son atenuados mediante la producción de pares $\gamma\gamma$ con fotones de rayos X.

Para explicar el espectro de neutrinos observado, los protones parentales deben ser acelerados a energías de $\sim 100$ TeV dentro de esta corona. Sin embargo, el mecanismo de aceleración dominante en las coronas de SMBH sigue sin establecerse. Aunque la Aceleración de Choques Difusiva (DSA, por sus siglas en inglés) es un candidato principal, su viabilidad es incierta bajo las condiciones de la corona. Específicamente, se espera que los choques coronales tengan números de Mach sónico bajos ($M_s \sim 2-10$) y velocidades transrelativistas, regímenes que tradicionalmente se consideran ineficientes para la aceleración de partículas. Además, el comportamiento de los plasmas de dos temperaturas ($T_i \gg T_e$) y la eficiencia específica de la aceleración de protones frente a la de electrones en estos entornos carecen de restricciones basadas en primeros principios.

Metodología
Para abordar estas incertidumbres, los autores emplean simulaciones de partículas en celda (PIC) unidimensionales basadas en primeros principios utilizando el código SMILEI. El estudio se centra en choques casi paralelos ($\theta = 30^\circ$) ajustados a los parámetros físicos de las coronas de SMBH.

El muestreo de simulaciones explora sistemáticamente un amplio espacio de parámetros motivado por las restricciones observacionales para las galaxias Seyfert:

• Número de Mach Sónico ($M_s$): Desde transónico ($M_s \approx 2$) hasta valores más altos ($M_s \approx 30$).
• Número de Mach Alvénico ($M_A$): Desde $4$ hasta $46$.
• Velocidad del Choque ($v_{sh}$): Velocidades transrelativistas ($v_{sh}/c \approx 0.17 - 0.66$).
• Relación de Temperatura Ión-Electrón ($T_i/T_e$): Desde $0.01$ hasta $50$, reflejando los estados térmicos desacoplados esperados en las coronas sin colisiones.

Las simulaciones utilizan una relación de masa ión-electrón reducida ($m_i/m_e = 100$) para asegurar la viabilidad computacional preservando al mismo tiempo los parámetros adimensionales clave. Las ejecuciones se evolucionaron durante un tiempo suficiente ($t \gtrsim 400 \Omega_{ci}^{-1}$) para permitir la saturación de las inestabilidades electromagnéticas y el desarrollo de colas de potencia no térmicas.

Resultados Clave

1. Aceleración de Protones Robusta: El estudio encuentra que la DSA es un mecanismo robusto y eficiente para la aceleración de protones en las coronas de SMBH. A pesar de los bajos números de Mach sónico, los choques canalizan consistentemente aproximadamente el 10% de la energía cinética del choque hacia iones no térmicos. Esta eficiencia permanece estable en todo el muestreo de parámetros, incluyendo casos con $M_s \approx 2$, desafiando la noción de que los choques de bajo $M_s$ son aceleradores ineficaces.
2. Dinámica de Aceleración: La energía máxima de los protones acelerados aumenta linealmente con el tiempo después de que la estructura del choque alcanza un estado de cuasi-equilibrio. La tasa de aceleración sigue la escala predicha por la teoría de la DSA, $t_{acc} \propto \gamma \beta_{sh}^{-2} \Omega_{ci}^{-1}$. Los autores derivan un coeficiente de aceleración normalizado, $C^{-1} \equiv d(\gamma \beta_{sh}^{-2})/d(t\Omega_{ci})$, que exhibe un valor universal de $C^{-1} \approx 0.1$, variando por menos de un factor de dos en todos los parámetros probados.
3. Ineficiencia de la Aceleración de Electrones: En contraste con los protones, la aceleración de electrones es altamente ineficiente, con fracciones de energía electrónica no térmica que permanecen por debajo del 1% (típicamente $\sim 10^{-3}$ para relaciones $T_i/T_e$ similares a las coronales). La eficiencia de la aceleración de electrones disminuye a medida que la relación $T_i/T_e$ aumenta. En casos transrelativistas específicos con amplificación de campo magnético fuerte, la cola de potencia de los electrones no logra extenderse a energías más altas con el tiempo, lo que sugiere que la DSA sostenida se ve inhibida cuando los campos amplificados vuelven al choque efectivamente superlumínico para los electrones.
4. Mecanismos de Inestabilidad: La turbulencia magnética de aguas arriba requerida para la dispersión es impulsada principalmente por la inestabilidad no resonante (NRI, o inestabilidad de Bell) en la mayor parte del espacio de parámetros. La inestabilidad resonante (RI) se vuelve dominante solo a números de Mach alvénicos muy bajos ($M_A \lesssim 8$).

Implicaciones para NGC 1068 y Producción de Neutrinos
Aplicando estos resultados a NGC 1068, los autores demuestran que:

• Escalas de Tiempo: Los protones pueden ser acelerados a $\sim 100$ TeV en escalas de tiempo más cortas que las escalas de enfriamiento relevantes (Bethe-Heitler, fotomesón, $pp$) y las escalas dinámicas (caída libre).
• Presupuesto de Energía: La potencia de protones no térmicos necesaria para explicar la luminosidad de neutrinos observada se estima en $\sim 3\%$ de la potencia total del choque. Esto está bien dentro del $\sim 10\%$ de eficiencia de aceleración medido en las simulaciones.
• Déficit de Rayos Gamma: La baja eficiencia de la aceleración de electrones y la alta opacidad de la corona a los rayos gamma explican naturalmente el déficit observado de rayos gamma de alta energía mientras permiten que los neutrinos escapen.

Significado y Reivindicaciones
El artículo afirma proporcionar las primeras restricciones consistentes con los primeros principios sobre la DSA en el contexto específico de las coronas de SMBH. La significancia principal radica en establecer que los choques sin colisiones son un mecanismo viable y eficiente para producir los protones de alta energía necesarios para explicar la señal de neutrinos de IceCube de NGC 1068, incluso bajo las condiciones desafiantes de bajos números de Mach sónico y altas relaciones de temperatura ión-electrón.

Los autores declaran explícitamente que, si bien sus resultados apoyan fuertemente el modelo hadrónico para la producción de neutrinos en los AGN, la aplicación de los resultados de simulaciones cinéticas locales a un objeto astrofísico complejo como una corona de AGN requiere una consideración cuidadosa. El análisis del presupuesto de energía se presenta como una evaluación de orden de magnitud basada en la microfísica establecida en el documento, cerrando la brecha entre la física de plasma microscópica y los fenómenos macroscópicos de alta energía. El estudio no pretende resolver todas las incertidumbres respecto a la energía máxima de los protones o la transición detallada de las inestabilidades en 3D, señalando que estas requieren mayor investigación.