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

Resumo Técnico: CHEX-MATE: CLUMP-3D II. Análise Combinada de Gás e Matéria Escura a partir de Raios-X, SZE e WL

Problema
Sob o modelo padrão de formação de estruturas hierárquicas, os aglomerados de galáxias são intrinsecamente triaxiais em vez de esféricos. Aplicar modelos com simetria esférica a essas estruturas introduz vieses significativos nas propriedades derivadas, particularmente nas estimativas de massa via lentes gravitacionais fracas (WL), que são altamente sensíveis à orientação e ao alongamento do aglomerado ao longo da linha de visada (LOS). Embora modelos triaxiais possam mitigar esses vieses, eles requerem dados de múltiplas sondas profundas e modelos fisicamente motivados para restringir a geometria tridimensional. Um trabalho anterior (Paper I, Kim et al. 2024, dorante K24) estabeleceu um formalismo de ajuste triaxial utilizando apenas dados de Raios-X e do Efeito Sunyaev-Zel'dovich (SZ) para restringir a geometria do gás. No entanto, essa abordagem focada apenas no gás não podia medir diretamente a massa total do aglomerado ou o suporte de pressão não térmica no meio intrac de aglomerado (ICM) sem assumir relações sobre a relação entre as distribuições de gás e de massa total.

Metodologia
Este artigo apresenta uma extensão do pipeline de ajuste triaxial do K24, integrando dados de lentes gravitacionais fracas (WL) à análise existente de Raios-X e SZ. A metodologia envolve:

1. Geometria Triaxial: O modelo assume um elipsoide triaxial para a distribuição do gás do ICM e para o potencial gravitacional. A geometria intrínseca é definida por razões axiais ($q_1, q_2$) e ângulos de Euler ($\theta, \phi, \psi$) em relação ao observador.
2. Modelagem do ICM: A densidade eletrônica segue um perfil $\beta$ modificado (Eq. 3), e a pressão eletrônica segue um perfil Navarro-Frenk-White generalizado (gNFW) (Eq. 4). Esses perfis são projetados ao longo da LOS para gerar mapas de modelo para brilho de superfície (SB) de Raios-X, parâmetro Compton-$y$ de SZ e temperatura espectroscópica.
3. Reconstrução de Massa: Diferente de abordagens anteriores que assumiam simetria elíptica para o perfil de densidade total, este trabalho assume que o potencial gravitacional segue a geometria triaxial do ICM ($q_{pot} = q_{ICM}$). Isso é motivado por simulações que sugerem que o potencial é mais suave do que a densidade total. O perfil de massa total é restringido pela forma NFW em uma base esférica, que é então invertida numericamente via equação de Poisson na grade triaxial para gerar o potencial elíptico. Este potencial é projetado para produzir observáveis de WL (mapas de cisalhamento reduzido).
4. Estrutura de Ajuste: É empregado um método de Monte Carlo via Cadeias de Markov (MCMC) usando emcee para ajustar o conjunto de dados combinado. A verossimilhança é definida por uma estatística $\chi^2$ total combinando mapas de SZ, SB de Raios-X (usando uma abordagem dual 2D/1D para mitigar o viés de agrupamento/clumping) e perfis de temperatura de Raios-X, além de mapas de cisalhamento 2D de WL.
5. Pressão Não Térmica: Uma vez que os perfis de densidade de massa total e de densidade de gás são restringidos, a fração de pressão não térmica ($P_{nt}/P_{tot}$) é calculada assumindo um estado de equilíbrio hidrostático generalizado. A pressão total necessária para sustentar o aglomerado contra a gravidade é derivada do potencial gravitacional, e o componente não térmico é a diferença entre esta pressão total e a pressão térmica derivada da temperatura do gás.

Principais Contribuições

• Integração de Pipeline: O artigo integra com sucesso os dados de WL no pipeline de ajuste triaxial de apenas gás do K24, permitindo a medição simultânea da geometria do aglomerado, massa total e concentração sem depender de suposições de simetria esférica.
• Modelagem Baseada no Potencial: Os autores introduzem uma abordagem inovadora onde a triaxialidade é aplicada ao potencial gravitacional em vez ao perfil de densidade total, permitindo um potencial mais suave que se alinha com as superfícies de isodensidade do ICM.
• Validação: A metodologia é validada usando observações simuladas (mock) de aglomerados com parâmetros de entrada conhecidos ($M_{200c}$, $c_{200c}$ e geometria). Os testes demonstram que o pipeline pode recuperar com precisão as massas e geometrias de entrada, com apenas vieses menores na concentração quando o valor de entrada desvia significativamente da distribuição a priori.

Resultados: Aplicação ao Abell 1689
O pipeline foi aplicado ao aglomerado CHEX-MATE PSZ2 G313.33+61.13 (Abell 1689), utilizando dados de Raios-X do XMM-Newton, dados de SZ de ACT/Planck e dados de WL do Subaru.

• Geometria: O aglomerado é encontrado como alongado ao longo da linha de visada em relação ao plano do céu, com um fator de alongamento de $R_{LP} = 1,27 \pm 0,02$.
• Massa: Devido a este alongamento, o ajuste triaxial resulta em uma massa total de $M_{200c} = (13,69^{+1,56}_{-1,41}) \times 10^{14} M_\odot$. Este valor é significativamente menor do que o de $(17,77^{+2,00}_{-1,75}) \times 10^{14} M_\odot$ obtido de um ajuste esférico simétrico utilizando a mesma metodologia, destacando o viés introduzido por suposições esféricas.
• Concentração: O ajuste triaxial resulta em uma concentração de $c_{200c} = 8,55^{+2,20}_{-1,61}$. Isto é consistente com o valor esférico simétrico de $9,99^{+2,26}_{-1,78}$, sugerindo que a alta concentração anteriormente observada em Abell 1689 não é um artefato de triaxialidade ou orientação.
• Pressão Não Térmica: A fração de pressão não térmica foi medida na faixa radial de 0,18–1,37 Mpc. A fração exibe um mínimo de aproximadamente 20% em raios intermediários, aumentando para perto de 30% tanto nos raios menores quanto nos maiores, com uma precisão de medição típica de $\pm 5\%$.

Significância
Este artigo demonstra que combinar dados de Raios-X, SZ e WL dentro de um arcabouço triaxial é essencial para obter medições não enviesadas de massa e concentração para aglomerados de galáxias individuais. Ao modelar explicitamente a geometria 3D, o estudo resolve a degenerescência entre massa e alongamento que assola as análises esféricas. Além disso, a capacidade de reconstruir o perfil de massa total em uma base triaxial permite uma medição robusta do suporte de pressão não térmica nas periferias do aglomerado, fornecendo restrições observacionais sobre a eficiência de termalização do material de acreção que podem ser comparadas com simulações hidrodinâmicas. Os resultados para Abell 1689 especificamente indicam que sua alta concentração é uma propriedade intrínseca e não um efeito de projeção, e que sua massa é menor do que previamente estimado quando a simetria esférica é assumida.