Assessing the sensitivity to Axion-Like-Particle Dark Matter with very-high-energy gamma-ray observations of selected AGN and galaxy cluster pairs

Resumen Técnico: Evaluación de la Sensibilidad a la Materia Oscura de Partículas Similares a Axiones con Observaciones de Gamma de Muy Alta Energía

Planteamiento del Problema
Las Partículas Similares a Axiones (ALPs, por sus siglas en inglés) son candidatas pseudoescalares para la materia oscura ligera que pueden acoplarse a los fotones, lo que conduce a oscilaciones entre ALPs y fotones en campos magnéticos externos. Si bien estas oscilaciones pueden dejar irregularidades espectrales distintivas en los rayos gamma de muy alta energía (VHE) de fuentes astrofísicas, su detección es un desafío. Específicamente, cuando los rayos gamma VHE de Núcleos Activos de Galaxias (AGNs) atraviesan los campos magnéticos de cúmulos de galaxias intervinientes, las características de absorción resultantes dependen fuertemente de la realización específica del campo magnético (orientación y longitud de coherencia). Para pares de fuente-cúmulo individuales, las propiedades del campo magnético suelen estar mal restringidas en escalas pequeñas, lo que dificulta la predicción de la impronta espectral y conduce a restricciones estadísticas débiles al analizar fuentes únicas.

Metodología
Para superar las limitaciones del análisis de una sola fuente, los autores proponen un análisis de apilamiento (stacking analysis) de múltiples pares de AGN-cúmulo de galaxias (AGN-GC). La metodología implica los siguientes pasos:

1. Selección de Fuentes: El estudio cruza la información del catálogo 4LAC-DR3-h de Fermi-LAT de AGNs brillantes en rayos gamma con catálogos de Sunyaev-Zeldovich de cúmulos de galaxias. Los criterios de selección son:

   • El desplazamiento al rojo (redshift) del AGN debe ser mayor que el del cúmulo.
   • La línea de visión debe pasar a menos de 500 kpc (comóvil) del centro del cúmulo.
   • El espectro del AGN debe ajustarse bien a una ley de potencia.
   • El par debe ser visible para los Telescopios de Imágenes Cherenkov actuales (IACTs: H.E.S.S., MAGIC, VERITAS) con un ángulo cenital inferior a 60°.
   • Este proceso genera 11 pares para H.E.S.S. y 15 pares para MAGIC/VERITAS, sumando un total de 41 conjuntos de datos simulados (mock datasets) independientes.

2. Modelado Físico:

   • Campos Magnéticos: El estudio se centra en cúmulos de núcleo frío donde el campo magnético escala con la densidad electrónica ($B \propto n_e^{0.67}$). La densidad electrónica sigue un modelo $\beta$.
   • Probabilidad de Oscilación: La probabilidad de supervivencia del fotón ($P_{\gamma\gamma}$) se calcula basándose en la masa del ALP ($m_a$) y el acoplamiento ($g_{a\gamma\gamma}$). Mientras que las realizaciones individuales de los dominios del campo magnético producen patrones de oscilación erráticos, los autores demuestran que el apilamiento de múltiples pares suaviza estas fluctuaciones. La probabilidad de supervivencia promediada se parametriza como una función de paso: $\langle P_{\gamma\gamma} \rangle = 1 - p_0 / (1 + (E/E_c)^k)$.
   • Simulación de Flujo: El flujo observado esperado se modela combinando el espectro intrínseco del AGN, la probabilidad de supervivencia inducida por ALPs y la absorción por la Radiación de Fondo Extragaláctica (EBL). Se generan datos simulados para 50 horas de observación por par en un rango de energía de 100 GeV a 60 TeV, incorporando las áreas efectivas y la resolución energética del detector.

3. Análisis Estadístico:

   • Se emplea una prueba de razón de log-verosimilitud de Poisson por intervalos para comparar dos hipótesis: la hipótesis nula (sin oscilación de ALP) y la hipótesis de ALP (espectros modulados por oscilaciones).
   • El análisis margina sobre los parámetros de molestia (nuisance parameters) para derivar límites superiores al Nivel de Confianza (C.L.) del 95% sobre el acoplamiento $g_{a\gamma\gamma}$.
   • Las incertidumbres sistemáticas, particularmente respecto a la selección del modelo de EBL (Dominguez, Franceschini, Finke), se evalúan mediante la prueba de escenarios de modelado erróneo.

Resultados Clave

• Mejora de la Sensibilidad: La combinación de conjuntos de datos agudiza significativamente la sensibilidad. Para una masa de ALP de $m_a = 3 \times 10^{-8}$ eV, el límite de acoplamiento alcanza $g_{a\gamma\gamma} = 6 \times 10^{-13}$ GeV$^{-1}$. Esto representa una mejora de un factor de tres al aumentar el conjunto de datos de 11 pares (solo H.E.S.S.) a 41 pares (combinando H.E.S.S., MAGIC y VERITAS).
• Rango de Masa: El método es más sensible en el rango de masa de 10 a 100 neV. La sensibilidad se degrada cerca de $10^{-8}$ eV debido al umbral de baja energía de los IACT y en masas superiores ($>10^{-7}$ eV) debido a la fuerte absorción de la EBL y la reducción de la estadística de fotones.
• Efectos Sistemáticos:
  • Modelos de EBL: La elección del modelo de EBL afecta la sensibilidad hasta en un 3% para $m_a < 7 \times 10^{-8}$ eV y hasta un 6% para $m_a \le 2 \times 10^{-7}$ eV.
  • Modelado Erróneo: Utilizar un modelo de EBL de Finke para analizar datos generados con un modelo de Franceschini puede inducir una señal falsa (TS negativo) en conjuntos de datos pequeños (11 pares) debido a las diferentes dependencias de la energía por debajo de 10 TeV para fuentes de bajo desplazamiento al rojo. Este sesgo disminuye con el conjunto de datos más grande de 41 pares, ya que las fuentes con mayor desplazamiento al rojo reducen el impacto.
  • Homogeneidad: Los autores señalan que un conjunto de datos estadísticamente homogéneo, que evite la sobrerrepresentación de objetos brillantes y cercanos, minimiza los efectos del modelado erróneo de la EBL.

Significancia y Reivindicaciones
El artículo afirma que esta estrategia de análisis de apilamiento ofrece un método robusto para sondear las interacciones entre ALPs y fotones en un espacio de parámetros previamente inexplorado. Específicamente, las observaciones simuladas combinadas de 41 pares de AGN-GC permiten la exploración de masas de ALP entre 10 y 100 neV con acoplamientos por encima de $5 \times 10^{-13}$ GeV$^{-1}$. Esta región específica es significativa porque corresponde al espacio de parámetros donde las ALPs podrían constituir toda la materia oscura del Universo.

Si bien el estudio reconoce que el enfoque requiere tiempos de observación sustanciales (550 horas para H.E.S.S. y 750 horas para MAGIC/VERITAS), afirma que tales campañas son factibles a lo largo de varios años. Los autores concluyen que, aunque las incertidumbres estadísticas dominan actualmente, los efectos sistemáticos pueden mitigarse mediante futuras observaciones multi-longitud de onda y mejores simulaciones. Posicionan este método como una vía prometedora para la búsqueda de materia oscura de tipo ALP, especialmente con la llegada de futuros observatorios como el Observatorio del Telescopio Cherenkov Array (CTAO).