On gamma rays as predictors of UHECR flux in AGNs

Resumen Técnico: Sobre los rayos $\gamma$ como predictores del flujo de UHECR en AGNs

Planteamiento del Problema
Los Núcleos Galácticos Activos (AGNs) son candidatos primarios para las fuentes de Rayos Cósmicos de Ultra Alta Energía (UHECRs). Sin embargo, modelos fenomenológicos recientes que utilizan el flujo observado de rayos $\gamma$ de muy alta energía (VHE) como un proxy directo de la luminosidad de los UHECR ($L_{CR}$) no logran reproducir los datos experimentales del Observatorio Pierre Auger. Específicamente, estos modelos tienen dificultades para ajustar simultáneamente el espectro de energía, la composición de masa y el mapa de dirección de llegada. Un punto crítico de falla es la sobreestimación de la contribución del blázar Markarian 421 (Mkn 421), el cual crea un "hotspot" (punto caliente) y un dipolo de dirección inconsistentes con las observaciones. Además, estos modelos a menudo producen un dipolo de dirección "espurio" y un ajuste deficiente a la forma del espectro, particularmente cuando se mezclan diferentes clases de AGNs (por ejemplo, galaxias de radio y blázares) en un mismo análisis.

Metodología
Los autores proponen que la discrepancia surge de una suposición implícita en trabajos previos: que la emisión de UHECR está irradiada (beamed) de la misma manera que la emisión de rayos $\gamma$ de los chorros relativistas. El artículo argumenta que, mientras los rayos $\gamma$ están sujetos a un fuerte dopaje Doppler ($D^4$) debido al movimiento relativista en "blobs" dentro del chorro, los UHECR probablemente experimentan isotropización (decollimación) antes de escapar del entorno de la fuente.

El estudio procede a través de los siguientes pasos:

1. Revisión Teórica: Los autores revisan los mecanismos de aceleración en los chorros de AGNs (choques, reconexión, aceleración de cizalla) y la propagación de partículas a través de los lóbulos de los chorros. Argumentan que los campos magnéticos en los lóbulos y el capullo (cocoon) de la fuente dispersan los UHECR, isotropizando efectivamente su emisión incluso si se originan en flujos relativistas.
2. Corrección del Proxy: En lugar de usar la luminosidad de rayos $\gamma$ observada ($L^{obs}_{\gamma}$), los autores proponen usar la luminosidad intrínseca de rayos $\gamma$ ($L^{int}_{\gamma}$) o la potencia radiativa del chorro como el proxy para $L_{CR}$. Esto implica corregir el flujo observado por el factor Doppler ($D$) elevado a la potencia $q$, donde $L^{int}_{\gamma} = D^{-q}L^{obs}_{\gamma}$. Se prueban dos escenarios: $q=4$ (asumiendo que $L_{CR} \propto L^{int}_{\gamma}$) y $q=2$ (asumiendo que $L_{CR} \propto$ potencia radiativa del chorro).
3. Simulación y Ajuste: Utilizando el marco de trabajo CRPropa3, los autores realizan simulaciones de propagación 1-D para cinco núcleos primarios (H, He, N, Si, Fe) a través de una cuadrícula de energías y distancias. Comparan el flujo simulado contra los datos de Pierre Auger (espectro y momentos de $X_{max}$) y mapas de dirección de llegada.
4. Análisis Comparativo: Se comparan tres casos:
   • Caso (a): Usando $L^{obs}_{\gamma}$ como el proxy (enfoque estándar).
   • Caso (b): Usando $D^{-2}L^{obs}_{\gamma}$ como el proxy.
   • Caso (c): Usando $D^{-4}L^{obs}_{\gamma}$ como el proxy.

Contribuciones Clave y Resultados

• Supresión de la Dominancia de los Blázares: La corrección reduce significativamente el peso de los blázares altamente irradiados (como Mkn 421) en relación con las galaxias de radio (como Cen A, M87 y Fornax A). En el modelo estándar, Mkn 421 domina el flujo; en los modelos corregidos, su contribución se suprime de 2 a 5 órdenes de magnitud, volviéndolo insignificante.
• Mejora en el Ajuste de Espectro y Composición: El ajuste combinado de los datos del espectro de energía y de composición mejora significativamente cuando se utiliza la luminosidad intrínseca como proxy. El chi-cuadrado reducido ($\chi^2/NDF$) disminuye de 4.6 (usando $L^{obs}_{\gamma}$) a 3.1 (usando $D^{-2}L^{obs}_{\gamma}$). Esta mejora se atribuye a una contribución relativa más realista de las fuentes cercanas frente a las fuentes distantes, evitando las distorsiones espectrales causadas por la sobrecontribución de fuentes distantes y altamente irradiadas.
• Resolución de la Tensión del Dipolo: La tensión entre las direcciones del dipolo modeladas y observadas se reduce. Para el campo de visión de Auger a energías $E > 8$ EeV, la discrepancia cae de 5.9$\sigma$ a 3.5$\sigma$. Para $E > 32$ EeV, la tensión se reduce de 2.1$\sigma$ a 1.1$\sigma$.
• Realineación de Hotspots: Los mapas de dirección de llegada predichos cambian de estar dominados por Mkn 421 a mostrar hotspots alrededor de Cen A, Fornax A, M87 y ciertos blázares distantes (Mkn 180, 1ES2344+514). Estos nuevos hotspots se alinean más estrechamente con los dos hotspots reportados por el Observatorio Pierre Auger y el uno reportado por el Telescope Array.

Significancia y Reivindicaciones
El artículo afirma ofrecer una solución a las contradicciones encontradas en trabajos previos de modelos de fuentes de UHECR basados en AGNs al desacoplar los supuestos de irradiación (beaming) de los rayos $\gamma$ y los UHECR. Los autores sostienen que:

1. La asociación entre los rayos $\gamma$ y los UHECR es débil pero no nula; aunque los rayos $\gamma$ trazan las regiones de aceleración de partículas, los UHECR son probablemente isotropizados por los campos magnéticos en los lóbulos de la fuente antes de escapar.
2. Usar el flujo de rayos $\gamma$ observado como un proxy asume implícitamente que los UHECR están irradiados, un supuesto que los autores argumentan es físicamente improbable para los UHECR que escapan de entornos magnéticos complejos.
3. Corregir por el dopaje Doppler ($L^{int}_{\gamma}$ o potencia radiativa del chorro) proporciona un proxy más robusto que reconcilia los datos de espectro, composición y dirección de llegada del Observatorio Pierre Auger.

Los autores concluyen que, si bien su modelo es un ajuste combinado simplificado (que no tiene en cuenta las distribuciones completas de $X_{max}$ ni las deflexiones detalladas de los campos magnéticos), los resultados fortalecen la hipótesis de los AGNs como fuentes de UHECR y proporcionan un marco más consistente para interpretar los datos de mensajeros múltiples de diferentes clases de AGNs. Señalan que el trabajo futuro debe abordar la decollimación dependiente de la energía y los efectos específicos de los campos magnéticos galácticos y extragalácticos regulares.