ICM-SHOX III. The case of MACS J0018.5+1626, a radio relic that looks like a radio halo?

Resumen Técnico: ICM-SHOX III. El caso de MACS J0018.5+1626, ¿un reliquia de radio que parece un halo de radio?

Planteamiento del Problema
La emisión de radio difusa en cúmulos de galaxias, típicamente clasificada como halos de radio (morfología central, redondeada) o reliquias de radio (periférica, alargada, impulsada por choques), proporciona información crítica sobre el presupuesto de energía no térmica y la dinámica de la fusión del medio intracúmulo (ICM). Sin embargo, la interpretación de dicha emisión suele verse complicada por los efectos de proyección y la falta de restricciones multi-sonda sobre los parámetros de la fusión (masa, parámetro de impacto, ángulo de visión, etapa de la fusión). El cúmulo de galaxias MACS J0018.5+1626 presenta un desafío específico: las observaciones de LOFAR a 144 MHz revelan una fuente de radio difusa, ubicada centralmente, con una morfología típica de un halo de radio. Por el contrario, el flujo de trabajo ICM-SHOX (Improved Constraints on Mergers with SZ, Hydrodynamical simulations, Optical, and X-ray) sugiere que el cúmulo está experimentando una fusión binaria cerca del paso por el pericentro, vista casi a lo largo del eje de la fusión. Esta configuración produce teóricamente dos choques impulsados por la fusión que, si se observan de frente (face-on), podrían imitar un halo de radio central. El problema central abordado es si la aceleración de choque difusiva (DSA) estándar actuando sobre electrones térmicos en tal configuración de choque "de frente" puede reproducir la potencia y morfología de radio observadas en MACS J0018.5+1626, o si se requieren mecanismos alternativos (por ejemplo, re-aceleración turbulenta o electrones fósiles).

Metodología
Los autores emplean un marco numérico híbrido que combina simulaciones magnetohidrodinámicas (MHD) en 3D con un resolvedor de Fokker–Planck para modelar la aceleración de partículas y la emisión de sincrotrón.

• Simulaciones MHD: Utilizando el código de malla móvil AREPO, los autores simulan fusiones binarias idealizadas de cúmulos de galaxias. Las condiciones iniciales están restringidas por los resultados del flujo de trabajo ICM-SHOX para MACS J0018.5+1626, específicamente estableciendo el desplazamiento al rojo ($z \approx 0.5456$), las masas de los componentes y los perfiles de gas. Las simulaciones exploran variaciones en las velocidades relativas iniciales ($v_i = 2400$ y $3000$ km s$^{-1}$), parámetros de impacto ($b = 100$ y $250$ kpc) y el beta de plasma inicial ($\beta_p = 100$ y $200$).
• Seguimiento de Partículas: Se advectan partículas trazadoras con el flujo del fluido para rastrear la evolución lagrangiana del gas. Solo los trazadores que encuentran celdas con choque (número de Mach $M > 1.3$) son considerados para la aceleración de partículas.
• Resolvedor de Fokker–Planck: Un resolvedor de Fokker–Planck modela la evolución de la distribución de momento de los electrones de rayos cósmicos (CRe). El resolvedor contabiliza la aceleración de choque difusiva (DSA), la expansión adiabática y las pérdidas radiativas (sincrotrón, Compton inverso, Coulomb, Bremsstrahlung).
• Modelos de Aceleración de Choque: El estudio pone a prueba dos modelos de eficiencia de aceleración de choque: una eficiencia constante ($\eta_0$) y una eficiencia dependiente del número de Mach basada en Kang et al. (2007). La eficiencia de inyección de electrones se trata como un parámetro libre ($\eta_e$) escalado por la eficiencia de los protones, explorando valores de $10^{-3}$ a $10^{-1}$.
• Emisión de Sincrotrón: El resolvedor computa la emisividad de sincrotrón para los trazadores, generando mapas 3D proyectados a lo largo de varias líneas de visión (LOS). Estos mapas se conlucionan con los tamaños de haz de LOFAR para comparar directamente con las observaciones a 144 MHz.
• Comparación Observacional: Los mapas simulados de radio y rayos X se comparan con los datos de LOFAR (144 MHz) y Chandra (0.5–7 keV). Las métricas incluyen la potencia de radio total, perfiles de brillo superficial radial, índices espectrales y coeficientes de correlación punto a punto de radio-rayos X.

Resultados Clave

1. Configuración de la Fusión y Propiedades de los Choques: Las simulaciones confirman que una fusión vista casi a lo largo del eje (de frente) produce dos choques que se superponen significativamente en el centro del cúmulo. Los choques resultantes, impulsados por la fusión, tienen números de Mach promedio típicos de $M_s \sim 2–3$ (con $\sigma_M \sim 0.5–1.5$).
2. Reproducción de la Potencia de Radio: La DSA estándar puede reproducir la potencia de radio observada de LOFAR ($\sim 37$ mJy a 144 MHz) solo bajo condiciones específicas:
   • Eficiencia: Se requiere una eficiencia de aceleración de choque de electrones de $\eta_e \sim 10^{-3}–10^{-2}$. Esto corresponde al extremo optimista de las estimaciones teóricas, pero es consistente con algunos resultados de simulaciones de Celda de Partícula en Celda (PIC) para choques débiles.
   • Campo Magnético: Un beta de plasma inicial de $\beta_p = 100$ es suficiente para generar la amplificación de campo magnético necesaria ($\sim 0.8–1$ $\mu$G en las regiones de choque) para coincidir con las observaciones. Los modelos con $\beta_p = 200$ no logran alcanzar la potencia observada.
   • Época: El mejor ajuste ocurre poco después del paso por el pericentro (aproximadamente 40–200 Myr post-pericentro), coincidiendo con el pico de disipación de energía del choque.
3. Morfología y Ángulo de Visión: La proyección de frente de los dos choques superpuestos produce naturalmente una morfología de radio circular y centralmente ubicada que se asemeja estrechamente a un halo de radio. El estudio encuentra que los ángulos de visión $\lesssim 30^\circ$ (para $v_i=2400$ km s$^{-1}$) o $\lesssim 36^\circ$ (para $v_i=3000$ km s$^{-1}$) proporcionan el mejor ajuste a los perfiles radiales y la distribución de brillo superficial observados. Los ángulos de visión más grandes resultan en estructuras tipo reliquia no superpuestas, inconsistentes con las observaciones.
4. Propiedades Espectrales y de Correlación:
   • El índice espectral integrado ($\alpha_{int}$) oscila entre $\sim -1.3$ y $-1.8$, consistente con los valores típicos de halos y reliquias de radio.
   • La correlación punto a punto entre el brillo superficial de radio y rayos X se encuentra variable. Dependiendo de la etapa de la fusión, el ángulo de visión y la resolución, la correlación puede ser superlineal, sublineal o lineal. Esto desafía la noción de que las correlaciones sublineales son una firma única de los halos de radio (origen turbulento) frente a las reliquias (origen por choque).

Significación y Reivindicaciones
El artículo afirma que el aparente "halo de radio" en MACS J0018.5+1626 es probablemente el resultado de la superposición de dos reliquias de radio impulsadas por la fusión vistas casi de frente, en lugar de un halo turbulento ubicado centralmente. Esta interpretación está respaldada por las restricciones multi-longitud de onda del flujo de trabajo ICM-SHOX, que favorecen una geometría de fusión que se alinea con el escenario de choque de frente.

El estudio demuestra que la DSA estándar, cuando se aplica a una configuración de choque de frente con eficiencias de aceleración optimistas pero físicamente plausibles ($\eta_e \sim 10^{-3}–10^{-2}$) y fuerzas de campo magnético, puede reproducir la potencia y morfología de radio del cúmulo sin invocar re-aceleración turbulenta o poblaciones de electrones fósiles. Sin embargo, los autores señalan que las eficiencias requeridas son del lado optimista, lo que sugiere que mecanismos adicionales (como la re-aceleración turbulenta o la presencia de electrones fósiles) aún pueden jugar un papel.

El trabajo destaca la importancia crítica del ángulo de visión y los efectos de proyección en la clasificación de la emisión de radio difusa. Además, establece el flujo de trabajo ICM-SHOX como una herramienta poderosa para vincular observables térmicos (rayos X, SZ) y no térmicos (radio) para restringir el estado físico de los cúmulos en fusión. Los autores concluyen que, aunque su modelo explica con éxito las observaciones bajo la DSA estándar, se requieren observaciones de radio de alta frecuencia y datos de polarización para distinguir mejor entre los orígenes impulsados por choques y por turbulencia, y para restringir mejor la eficiencia de la aceleración de electrones en el ICM.