Gamma-Ray Constraints on Heavy Axion-Like-Particle Decays from Fermi-LAT and H.E.S.S. Blazar Spectra

Resumen Técnico: Restricciones de Rayos Gamma sobre la Desintegración de Axiones-Similares Pesados a partir de los Espectros de Blázares de Fermi-LAT y H.E.S.S.

Planteamiento del Problema
Los rayos gamma de muy alta energía (VHE; $E_\gamma \ge 100$ GeV) provenientes de fuentes extragalácticas, tales como los blázares, se atenúan durante su propagación a través del espacio intergaláctico debido a la producción de pares ($\gamma\gamma \to e^+e^-$) con los fotones del Fondo Cósmico Extragaláctico (EBL, por sus siglas en inglés). Si bien los modelos estándar del EBL, basados en la evolución de las galaxias y sondeos profundos, describen con éxito esta atenuación, las posibles anomalías en la transparencia del Universo han motivado la exploración de la física más allá del Modelo Estándar (BSM). Específicamente, este estudio investiga un escenario que involucra Axiones-Similares (ALPs) pesados con masas $m_a \sim 2.5\text{--}20$ eV. A diferencia de los ALPs más ligeros donde las oscilaciones fotón-ALP podrían aumentar la transparencia, se hipotetiza que los ALPs pesados en este rango de masa se desintegran en dos fotones ($a \to \gamma\gamma$) a lo largo de escalas de tiempo cosmológicas. Si los ALPs constituyeran una fracción significativa de la materia oscura, su desintegración inyectaría fotones adicionales en el fondo cósmico infrarrojo/óptico, aumentando así la intensidad del EBL y elevando la profundidad óptica de los rayos gamma. Este artículo busca restringir la masa ($m_a$) y el acoplamiento fotón-ALP ($g_{a\gamma}$) de tales ALPs pesados, probando si los espectros observados de los blázares son consistentes con un modelo de EBL aumentado por las contribuciones de la desintegración de los ALPs.

Metodología
El análisis utiliza un conjunto de datos conjunto que comprende observaciones del Telescopio de Gran Área Fermi (Fermi-LAT) y del Sistema de Telescopios H.E.S.S. (High Energy Stereoscopic System).

• Selección de Datos: La muestra consiste en 11 blázares seleccionados basándose en su detección a energías VHE y en desplazamientos al rojo donde la absorción del EBL es significativa ($\tau > 2$). El conjunto de datos incluye 14 años de datos de Fermi-LAT (2008–2022) y observaciones de H.E.S.S.
• Segmentación Temporal: Para dar cuenta de la variabilidad de las fuentes, los autores emplearon el algoritmo de bloques de Bayes para segmentar las curvas de luz. El análisis espectral se realizó en estos segmentos para asegurar que la evolución espectral intrínseca no sesmara los resultados. Para las fuentes variables, se emparejaron los datos contemporáneos de Fermi-LAT y H.E.S.S.; para los segmentos previos a Fermi, solo se utilizaron datos de H.E.S.S.
• Modelado Espectral: Se realizó un análisis de máxima verosimilitud conjunta utilizando el paquete Gammapy. Los espectros intrínsecos de las fuentes se modelaron utilizando funciones de ley de potencia o de log-parábola. El análisis incluyó un parámetro de molestia para dar cuenta de un desplazamiento coherente de la escala de energía en los datos de H.E.S.S. ($\sim 10\text{--}15\%$ de incertidumbre).
• Formalismo de Desintegración de ALPs: La contribución de la desintegración de los ALPs al EBL se modeló utilizando el formalismo de [54–56]. La intensidad específica del EBL inducido por ALPs depende de la masa del ALP, el acoplamiento y la densidad ($\Omega_a$). La profundidad óptica total se calculó como la suma de la contribución estándar del EBL (usando el modelo de Domínguez et al. [9] como base, con comparaciones con otros modelos) y la atenuación adicional inducida por los fotones de la desintegración de los ALPs.
• Marco Estadístico: Los autores construyeron una razón de verosimilitud de perfil, $\lambda(m_a, g_{a\gamma})$, sumando los estadísticos C a través de todas las fuentes y segmentos de tiempo. Los límites de exclusión se derivaron al nivel de confianza (C.L.) del 95% bajo la suposición de que los ALPs constituyen la totalidad de la densidad de materia oscura ($\Omega_a = \Omega_{DM}$).

Contribuciones Clave

• Primer Análisis Conjunto: Este trabajo representa el primer análisis de las firmas de desintegración de ALPs pesados utilizando un marco de verosimilitud unificado que combina datos de Fermi-LAT y H.E.S.S., maximizando la robustez estadística a través de un amplio rango de energía (MeV a TeV).
• Cobertura del Espacio de Parámetros: El estudio escanea el rango de masa de los ALPs de $m_a \approx 2.5\text{--}20$ eV. Este rango está motivado físicamente porque los fotones de desintegración en esta ventana de masa caen en la banda óptica a infrarrojo cercano, lo cual afecta más fuertemente la atenuación de los rayos gamma TeV mediante la producción de pares.
• Manejo de Sistemas: El análisis incorpora explícitamente las incertidumbres de la escala de energía de H.E.S.S. y prueba la robustez de los resultados frente a cinco modelos diferentes de EBL (Domínguez, Franceschini, Finke, Saldana-Lopez, y Porras-Bedmar & Meyer).

Resultados

• Límites de Exclusión: El análisis combinado produce límites de exclusión de 95% C.L. sobre el acoplamiento fotón-ALP $g_{a\gamma}$. Las restricciones más fuertes se encuentran en las masas más altas donde la vida media de desintegración es más corta, lo que conduce a una mayor emisión de fotones.
  • En $m_a \approx 15$ eV, el límite alcanza $g_{a\gamma} \sim 7 \times 10^{-12}$ GeV$^{-1}$.
  • En $m_a \approx 20$ eV, el límite es aproximadamente $g_{a\gamma} \sim 5.4 \times 10^{-12}$ GeV$^{-1}$.
• Consistencia con la Hipótesis Nula: Los datos no muestran una preferencia estadística significativa por un componente inducido por ALPs sobre el modelo estándar de EBL. La gran mayoría de los segmentos de fuentes individuales y el conjunto de datos combinado son consistentes con la hipótesis nula (sin desintegración de ALPs).
• Robustez del Modelo: Los límites derivados son robustos frente a las incertidumbres en el modelado del EBL, con variaciones entre los diferentes plantillas de EBL que se mantienen dentro de un factor de aproximadamente 1.5.
• Comparación con Reclamaciones Previas: Los resultados desfavorecen el espacio de parámetros sugerido en trabajos previos (por ejemplo, [55, 67]) que proponían que las desintegraciones de ALPs podrían explicar excesos en el fondo óptico cósmico o en los espectros del EBL. Específicamente, la región sugerida por [55] ($m_a \sim 8\text{--}20$ eV, $g_{a\gamma} \sim (3\text{--}6) \times 10^{-11}$ GeV$^{-1}$) queda excluida bajo la suposición de $\Omega_a = \Omega_{DM}$.

Significancia y Reclamaciones
El artículo afirma proporcionar las restricciones más fuertes hasta la fecha sobre ALPs pesados en el rango $m_a \sim 15\text{--}20$ eV basadas en observaciones de rayos gamma. Al cerrar una región de parámetros previamente no restringida, estos resultados ofrecen una prueba complementaria a las búsquedas espectroscópicas de líneas de desintegración estrechas y experimentos de laboratorio. Los autores enfatizan que su método depende de la contribución integrada de la desintegración al EBL y la resultante modificación de la profundidad óptica de los rayos gamma, lo que lo hace independiente de las suposiciones respecto a factores D específicos de la fuente o fondos de líneas estrechas. Aunque el análisis no encuentra evidencia de la desintegración de ALPs pesados, establece límites superiores estrictos que restringen la viabilidad de tales partículas como un componente dominante de la materia oscura capaz de alterar significativamente la transparencia de los rayos gamma del Universo. El trabajo demuestra la eficacia de los análisis de rayos gamma multi-instrumento para sondear la física BSM y establece una línea base para futuras búsquedas más sensibles con instrumentos de próxima generación como el Cherenkov Telescope Array (CTAO).