CHEX-MATE: Cluster Multi-Probes in Three Dimensions (CLUMP-3D) II. Combined Gas and Dark Matter Analysis from X-ray, SZE, and WL

Resumen Técnico: CHEX-MATE: CLUMP-3D II. Análisis Combinado de Gas y Materia Oscura a partir de Rayos X, SZE y WL

Planteamiento del Problema
Bajo el modelo estándar de formación jerárquica de estructuras, los cúmulos de galaxias son intrínsecamente triaxiales en lugar de esféricos. La aplicación de modelos con simetría esférica a estas estructuras introduce sesgos significativos en las propiedades derivadas, particularmente en las estimaciones de masa mediante lentes gravitacionales débiles (WL), las cuales son altamente sensibles a la orientación y elongación del cúmulo a lo largo de la línea de visión (LOS). Si bien los modelos triaxiales pueden mitigar estos sesgos, requieren datos de múltiples sondas profundos y modelos motivados físicamente para restringir la geometría tridimensional. El trabajo previo (Paper I, Kim et al. 2024, en adelante K24) estableció un formalismo de ajuste triaxial utilizando únicamente datos de rayos X y de efecto Sunyaev-Zel'dovich (SZ) para restringir la geometría del gas. Sin embargo, este enfoque basado solo en el gas no podía medir directamente la masa total del cúmulo ni el soporte de presión no térmica en el medio intraclúster (ICM) sin asumir relaciones sobre la relación entre la distribución del gas y la de la masa total.

Metodología
Este artículo presenta una extensión del flujo de trabajo de ajuste triaxial de K24 integrando datos de lentes gravitacionales débiles (WL) con el análisis existente de rayos X y SZ. La metodología consiste en:

1. Geometría Triaxial: El modelo asume un elipsoide triaxial para la distribución del gas del ICM y el potencial gravitatorio. La geometría intrínseca se define mediante razones axiales ($q_1, q_2$) y ángulos de Euler ($\theta, \phi, \psi$) respecto al observador.
2. Modelado del ICM: La densidad electrónica sigue un perfil $\beta$ modificado (Ec. 3), y la presión electrónica sigue un perfil Navarro-Frenk-White generalizado (gNFW) (Ec. 4). Estos perfiles se proyectan a lo largo de la LOS para generar mapas de modelos para la brillantez superficial de rayos X (SB), el parámetro de Compton-$y$ de SZ y la temperatura espectroscópica.
3. Reconstrucción de Masa: A diferencia de enfoques previos que asumían simetría elíptica para el perfil de densidad total, este trabajo asume que el potencial gravitatorio sigue la geometría triaxial del ICM ($q_{pot} = q_{ICM}$). Esto está motivado por simulaciones que sugieren que el potencial es más suave que la densidad total. El perfil de masa total se restringe mediante la forma NFW en una base esférica, la cual es luego invertida numéricamente mediante la ecuación de Poisson sobre la rejilla triaxial para generar el potencial elíptico. Este potencial se proyecta para producir observables de WL (mapas de cizalladura reducida).
4. Marco de Ajuste: Se emplea un enfoque de Monte Carlo por Cadenas de Markov (MCMC) utilizando emcee para ajustar el conjunto de datos combinado. La verosimilitud se define mediante una estadística $\chi^2$ total que combina mapas de SZ, SB de rayos X (usando un enfoque dual 2D/1D para mitigar el sesgo de agrupamiento o clumping) y perfiles de temperatura de rayos X, así como mapas de cizalladura 2D de WL.
5. Presión No Térmica: Una vez que los perfiles de densidad de masa total y densidad de gas están restringidos, se calcula la fracción de presión no térmica ($P_{nt}/P_{tot}$) asumiendo un estado de equilibrio hidrostático generalizado. La presión total requerida para soportar el cúmulo contra la gravedad se deriva del potencial gravitatorio, y el componente no térmico es la diferencia entre esta presión total y la presión térmica derivada de la temperatura del gas.

Contribuciones Clave

• Integración del Flujo de Trabajo: El artículo integra con éxito los datos de WL en el flujo de trabajo de ajuste triaxial de solo gas de K24, permitiendo la medición simultánea de la geometría, la masa total y la concentración del cúmulo sin depender de supuestos de simetría esférica.
• Modelado Basado en el Potencial: Los autores introducen un enfoque novedoso donde la triaxialidad se aplica al potencial gravitatorio en lugar al perfil de densidad total, lo que permite un potencial más suave que se alinea con las superficies de isodensidad del ICM.
• Validación: La metodología se valida utilizando observaciones simuladas (mock) de cúmulos con parámetros de entrada conocidos ($M_{200c}$, $c_{200c}$ y geometría). Las pruebas demuestran que el flujo de trabajo puede recuperar con precisión las masas y geometrías de entrada, con solo sesgos menores en la concentración cuando el valor de entrada se desvía significativamente de la distribución previa.

Resultados: Aplicación a Abell 1689
El flujo de trabajo se aplicó al cúmulo CHEX-MATE PSZ2 G313.33+61.13 (Abell 1689), utilizando datos de rayos X de XMM-Newton, datos de SZ de ACT/Planck y datos de WL de Subaru.

• Geometría: Se encuentra que el cúmulo está elongado a lo largo de la línea de visión respecto al plano del cielo, con un factor de elongación de $R_{LP} = 1.27 \pm 0.02$.
• Masa: Debido a esta elongación, el ajuste triaxial arroja una masa total de $M_{200c} = (13.69^{+1.56}_{-1.41}) \times 10^{14} M_\odot$. Este valor es significativamente menor que el de $(17.77^{+2.00}_{-1.75}) \times 10^{14} M_\odot$ obtenido de un ajuste esféricamente simétrico utilizando la misma metodología, resaltando el sesgo introducido por los supuestos de simetría esférica.
• Concentración: El ajuste triaxial arroja una concentración de $c_{200c} = 8.55^{+2.20}_{-1.61}$. Esto es consistente con el valor esféricamente simétrico de $9.99^{+2.26}_{-1.78}$, lo que sugiere que la alta concentración previamente observada en Abell 1689 no es un artefacto de la triaxialidad o la orientación.
• Presión No Térmica: La fracción de presión no térmica fue medida en el rango radial de 0.18–1.37 Mpc. La fracción exhibe un mínimo de aproximadamente el 20% en radios intermedios, aumentando a cerca del 30% tanto en los radios más pequeños como en los más grandes, con una precisión de medición típica de $\pm 5\%$.

Significancia
Este artículo demuestra que combinar datos de rayos X, SZ y WL dentro de un marco triaxial es esencial para obtener mediciones imparciales de masa y concentración para cúmulos de galaxias individuales. Al modelar explícitamente la geometría 3D, el estudio resuelve la degeneración entre masa y elongación que afecta a los análisis esféricos. Además, la capacidad de reconstruir el perfil de masa total en una base triaxial permite una medición robusta del soporte de presión no térmica en las periferias del cúmulo, proporcionando restricciones observacionales sobre la eficiencia de termalización del material que se acrece, las cuales pueden compararse con simulaciones hidrodinámicas. Los resultados para Abell 1689 indican específicamente que su alta concentración es una propiedad intrínseca y no un efecto de proyección, y que su masa es menor de lo estimado previamente cuando se asume simetría esférica.