A MINOT-based Study of Gamma-ray emission from SPT-CL J2012-5649/Abell 3667

Resumen Técnico: Un Estudio Basado en MINOT de la Emisión de Rayos Gamma de SPT-CL J2012-5649/Abell 3667

Planteamiento del Problema
Los cúmulos de galaxias son las estructuras virializadas más grandes del universo, actuando como reservorios de rayos cósmicos (RC) de alta energía acelerados por choques impulsados por fusiones. Si bien las búsquedas iniciales con Fermi-LAT de emisión de rayos gamma procedente de cúmulos arrojaron resultados nulos, estudios recientes han reportado detecciones tentativas, incluyendo una señal de $\sim6.1\sigma$ del cúmulo en fusión SPT-CL J2012-5649 (Abell 3667, en adelante A3667). Sin embargo, el origen de esta emisión permanece incierto debido a la presencia de galaxias de radio cercanas dentro de la función de dispersión puntual de Fermi-LAT y a una tensión significativa entre el índice espectral observado ($\Gamma \approx -3.6$) y las predicciones hadrónicas teóricas ($\Gamma \approx -2.4$). Además, las contribuciones relativas de los procesos hadrónicos (interacciones protón-protón) frente a los procesos leptónicos (dispersión Compton inversa de electrones de rayos cósmicos) requieren una cuantificación rigurosa para interpretar la física no térmica del medio intracúmulo (MIC).

Metodología
Los autores emplean el marco MINOT (Modelado e Interpretación de la emisión No Térmica de cúmulos de galaxias del Tipo de la Vía Láctea) para construir un modelo físico autoconsistente de A3667. El estudio utiliza las siguientes entradas y procedimientos:

• Datos Termodinámicos: Los perfiles de densidad y temperatura de electrones se derivan de observaciones de rayos X de archivo de Chandra (catálogo ACCEPT), abarcando desde el núcleo del cúmulo hasta $0.22 R_{500}$.
• Extrapolación de Perfiles: Los perfiles termodinámicos observados se extrapolan a un radio de truncamiento ($R_{trunc} = 5000$ kpc $\approx 3.7 R_{500}$) utilizando un modelo $\beta$ para la densidad de electrones y un modelo de ley de potencia para la temperatura.
• Configuración de Física No Térmica:
  • Campo Magnético: Modelado como escalando con la densidad de electrones ($B \propto n_e^{\eta_B}$) con $\eta_B = 0.5$ y un campo de referencia de $5 \mu$G a 100 kpc.
  • Rayos Cósmicos: Las distribuciones espaciales de protones ($CR_p$) y electrones ($CRe$) se escalan en relación con la densidad térmica ($n_{CR} \propto n_e^{\eta_{CR}}$). La línea base asume un escalado isodens ($\eta_{CR}=1$).
  • Espectros: Los $CR_p$ siguen una ley de potencia ($dN/dp \propto p^{-2.4}$). Los $CRe$ primarios siguen un espectro de inyección continua con una ruptura de enfriamiento a 5 GeV. Los $CRe$ secundarios se calculan de manera autoconsistente mediante interacciones $pp$ utilizando el modelo Pythia8.
• Cálculo de Emisión: El marco calcula los flujos de rayos gamma desde dos canales:
  1. Hadrónico: Desintegración de piones neutros ($\pi^0 \to \gamma\gamma$) a partir de colisiones $pp$.
  2. Leptónico: Dispersión Compton inversa (CI) de $CRe$ sobre el Fondo Cósmico de Microondas (CMB).
• Integración: Los flujos se integran sobre volúmenes esféricos desde 1 kpc hasta $R_{trunc}$ y se comparan con las observaciones de Fermi-LAT (1–300 GeV) y los límites superiores de DAMPE, COMPTEL e INTEGRAL.

Resultados Clave

1. Predicción de Flujo Hadrónico: El flujo de rayos gamma hadrónico predicho en la banda de 1–300 GeV es de $2.82 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$ dentro de $R_{500}$. Cuando se integra hasta el radio de truncamiento ($3.7 R_{500}$), el flujo aumenta a $1.15 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$. Este valor concuerda en orden de magnitud con el flujo reportado por Fermi-LAT de $1.3 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$.
2. Distribución Espacial: Aproximadamente el 76% del flujo hadrónico total predicho proviene de regiones más allá de $R_{500}$, consecuencia del perfil de densidad exterior poco profundo ($\beta \approx 0.33$) que mantiene una columna de gas significativa a grandes radios.
3. Subdominancia Leptónica: La contribución de CI tanto de electrones primarios como secundarios es subdominante frente al componente hadrónico por un factor de $\sim20$ dentro de $R_{500}$ en la banda de 1–300 GeV. Incluso al integrar hasta $R_{trunc}$, el flujo de CI ($1.49 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$) permanece un orden de magnitud por debajo del flujo observado por Fermi-LAT y de la predicción hadrónica.
4. Tensión Espectral: Si bien la normalización del flujo en $R_{trunc}$ coincide con las observaciones, el índice espectral presenta una discrepancia significativa. El espectro observado es empinado ($\Gamma = -3.61 \pm 0.33$), mientras que el modelo hadrónico predice un espectro más duro ($\Gamma \approx -2.4$ a $-2.6$ para $\alpha_p = 2.4$). Incluso suavizar el espectro de protones a $\alpha_p = 2.8$ no logra reproducir la pendiente observada.

Significado y Conclusiones
El artículo concluye que A3667 es un objetivo prometedor pero complejo para estudios de emisión no térmica. El significado principal del trabajo radica en demostrar que:

• El nivel de flujo de rayos gamma observado puede reproducirse mediante un modelo hadrónico estándar si la integración se extiende a las afueras del cúmulo ($R_{trunc}$), destacando la importancia de modelar la emisión más allá de $R_{500}$.
• El canal de CI es despreciable para explicar la señal de Fermi-LAT en A3667, confirmando que cualquier señal detectable es de origen hadrónico bajo suposiciones estándar de RC.
• Existe una tensión entre la forma espectral observada y las predicciones hadrónicas. Los autores afirman que un modelo de MIC puramente hadrónico con protones de RC isobáricos estándar no puede explicar completamente la pendiente espectral observada. Sugieren que esta discrepancia podría indicar contaminación por fuentes de fondo o que el espectro de protones de RC ha sido modificado por procesos de reaceleración y transporte no capturados por un modelo simple de ley de potencia.

Los autores enfatizan que sus resultados son estimaciones de orden de magnitud sujetas a incertidumbres en la distribución de RC, la estructura del campo magnético y la asunción de simetría esférica, la cual no captura la morfología asimétrica de los choques de fusión y frentes fríos de A3667. Se sugiere trabajo futuro para incorporar datos de radio resueltos espacialmente y modelado 3D para refinar las predicciones de brillo superficial.