High-energy Multi-messenger Emission from Galaxy Clusters in the Local Universe

Resumen Técnico: Emisión Multi-mensajero de Alta Energía de Cúmulos de Galaxias en el Universo Local

Planteamiento del Problema
El origen de los fondos difusos de neutrinos y rayos $\gamma$ de alta energía sigue siendo un problema mayor sin resolver en la astrofísica. Aunque se sabe que los núcleos galácticos activos (NGA) y las galaxias con brote estelar son contribuyentes, el papel específico de los cúmulos de galaxias como fuentes de estos mensajeros cósmicos no está completamente restringido. Los cúmulos de galaxias, caracterizados por un gas medio intra-cumular (MIC) caliente ($\sim 10^8$ K) y campos magnéticos fuertes ($\sim 1$ $\mu$G), son teóricamente capaces de acelerar rayos cósmicos (RC) mediante choques de fusión, choques de acreción y turbulencia. Estos protones acelerados pueden interactuar con el MIC para producir rayos $\gamma$ y neutrinos secundarios mediante procesos hadrónicos. Sin embargo, estudios anteriores han dependido a menudo de modelos simplificados y esféricamente simétricos o han asumido la inyección de RC únicamente desde los núcleos centrales de los cúmulos, lo que podría distorsionar la distribución espacial y el flujo total de la emisión difusa. Además, distinguir la emisión difusa de los cúmulos de la emisión de fuentes puntuales (por ejemplo, de NGA centrales como NGC 1275 o M87) sigue siendo un desafío debido a la resolución angular limitada en las observaciones actuales.

Metodología
Este trabajo emplea un marco numérico riguroso que combina simulaciones magnetohidrodinámicas (MHD) de alta resolución con métodos de transporte Monte Carlo para modelar la emisión multi-mensajero de cúmulos de galaxias en el Universo local ($z \lesssim 0.1$).

1. Simulaciones MHD (SLOW): El estudio utiliza la simulación "Simulating the LOcal Web" (SLOW), que incorpora condiciones iniciales restringidas observacionalmente basadas en el catálogo CosmicFLow. A diferencia de modelos anteriores, SLOW proporciona una descripción 3D autoconsistente y de alta resolución de las distribuciones de densidad de gas, temperatura y campo magnético dentro de los cúmulos (específicamente Virgo, Perseo y Coma), capturando morfologías no esféricas y subestructuras derivadas de fusiones.
2. Inyección de Rayos Cósmicos: Se adopta un modelo de inyección motivado físicamente donde la inyección de RC sigue la distribución espacial de la densidad de gas del cúmulo. Este enfoque considera la aceleración de RC en regiones de alta densidad asociadas a turbulencia y choques en todo el volumen del cúmulo, en lugar de restringir la inyección al núcleo del cúmulo. El presupuesto total de energía de RC se parametriza mediante la relación $X_{CR} = E_{CR}/E_{cluster} \simeq 0.01$, consistente con los límites superiores observacionales.
3. Transporte de RC y Producción Secundaria: La propagación de RC (asumidos como protones) dentro del MIC y el medio intergaláctico (MIG) se simula utilizando el código Monte Carlo CRPropa. El código resuelve la ecuación de transporte de Parker, teniendo en cuenta la advección, la difusión espacial, la difusión de momento y las pérdidas de energía. Los canales de interacción clave incluyen colisiones protón-protón ($pp$) (dominantes por debajo de $10^{17}$ eV), producción de fotopiones y cascadas electromagnéticas (dispersión Compton inversa sobre el CMB y el EBL, pérdidas por sincrotrón). Los rayos $\gamma$ y neutrinos secundarios se generan de manera autoconsistente a partir de estas interacciones.

Contribuciones Clave

• Modelado Ambiental Realista: La integración de simulaciones MHD SLOW permite modelar la propagación de RC en entornos de cúmulos realistas y no esféricos, superando las suposiciones simplificadas de simetría esférica de estudios analíticos o semi-analíticos anteriores.
• Esquema de Inyección Extendido: La adopción de un modelo de inyección de RC espacialmente extendido, que rastrea el perfil de densidad de gas, representa una mejora metodológica significativa sobre trabajos anteriores que asumían inyección exclusiva en el centro. Esto reduce el tiempo de confinamiento efectivo de los RC y las tasas de interacción en el núcleo, lo que conduce a un perfil de emisión más difuso.
• Estimación Integral de Flujo: El estudio calcula tanto los flujos individuales de cúmulos (Virgo, Perseo, Coma) como el flujo difuso acumulativo de toda la población de cúmulos masivos ($M \gtrsim 5 \times 10^{13} M_\odot$) dentro de una distancia comóvil de 500 Mpc.

Resultados

• Cúmulos Individuales: Los flujos de rayos $\gamma$ difusos predichos para los cúmulos de Virgo, Perseo y Coma están muy por debajo de los límites superiores actuales reportados por la colaboración LHAASO. Los flujos predichos también son inferiores a la emisión de rayos $\gamma$ observada por MAGIC desde la fuente central NGC 1275 en Perseo, lo cual es consistente con la expectativa de que MAGIC detecta una fuente puntual mientras que el modelo predice emisión difusa.
• Flujo Acumulativo: El flujo de rayos $\gamma$ integrado de todos los cúmulos dentro de 500 Mpc es significativamente menor que el fondo difuso de rayos $\gamma$ medido por Fermi-LAT. El estudio señala que, a energías por debajo de $\sim 100$ GeV, otras fuentes (blázares, NGA desalineados, galaxias con brote estelar) dominan el fondo.
• Emisión de Neutrinos: El flujo de neutrinos de todos los sabores predicho para la población local de cúmulos está muy por debajo del flujo difuso observado por IceCube, pero es comparable a los límites superiores existentes de IceCube específicamente para cúmulos de galaxias ($M > 10^{14} M_\odot$). Se encuentra que el flujo de neutrinos predicho para el cúmulo de Virgo es comparable a la sensibilidad proyectada del detector KM3NeT/ARCA.
• Características Espectrales: Los índices espectrales de los flujos predichos ($\alpha \approx 2.0 - 2.5$) son consistentes con las expectativas para la aceleración por choques de fusión y acreción. Los resultados indican que los cúmulos de menor masa contribuyen eficientemente a la emisión de rayos $\gamma$ de $\sim 10$ GeV, mientras que los sistemas masivos dominan el flujo de neutrinos de alta energía ($\gtrsim 10$ TeV) debido a sus superiores capacidades de confinamiento de RC.

Significado y Afirmaciones
El artículo afirma que sus resultados son consistentes con las restricciones observacionales actuales, particularmente la no detección de emisión difusa significativa por parte de LHAASO y los límites superiores establecidos por IceCube. Los autores postulan que la tensión entre las predicciones de modelos anteriores (que a menudo asumían flujos más altos) y los límites superiores de IceCube se reduce sustancialmente mediante la adopción del modelo de inyección espacialmente extendido y el uso de simulaciones SLOW de alta resolución.

El estudio concluye que, aunque es poco probable que los cúmulos de galaxias dominen el fondo total de neutrinos de alta energía difuso, podrían representar una fracción no despreciable del mismo. El trabajo proporciona un marco teórico que motiva búsquedas dedicadas de neutrinos difusos desde el cúmulo de Virgo utilizando KM3NeT e IceCube-Gen2, y de rayos $\gamma$ difusos utilizando CTA y LHAASO. Los autores enfatizan que sus hallazgos se alinean con las perspectivas de los futuros observatorios multi-mensajero, ofreciendo una prueba oportuna de los procesos de aceleración hadrónica de RC en estructuras a gran escala.