Constraining Galaxy Cluster Triaxiality via Weak Lensing -- I. Preparation for the Rubin Data Beyond Leading Order

Resumo Técnico: Constraindo a Triaxialidade de Aglomerados de Galáxias via Lente Fraca — I. Preparação para os Dados do Rubin Além da Ordem Dominante

Declaração do Problema
Os aglomerados de galáxias, as maiores estruturas gravitacionalmente colapsadas no universo, possuem distribuições de massa 3D triaxiais em vez de esféricas. Essa geometria introduz uma dispersão significativa (aproximadamente 20%) nas medições de massa por lente fraca (LF) de aglomerados individuais se um modelo de densidade esférico for assumido. Além disso, a triaxialidade acoplada a efeitos de projeção pode induzir vieses de detecção e seleção, onde aglomerados alinhados ao longo da linha de visão parecem mais ricos e mais fáceis de detectar, enquanto aqueles perpendiculares parecem mais arredondados e mais difíceis de detectar. Constringir com precisão as formas projetadas (elipticidades) de conjuntos de aglomerados é crucial para reduzir a dispersão nas medições de massa, corrigir vieses de seleção e sondar parâmetros cosmológicos e física da matéria escura (por exemplo, restrições sobre matéria escura auto-interagente). Embora estudos anteriores tenham medido as formas de aglomerados, uma análise abrangente de conjunto dos dados do Ano 3 (Y3) do Dark Energy Survey (DES) usando modelagem avançada de multipolos não havia sido realizada.

Metodologia
Os autores utilizam o software Cluster Lensing Mass Modeling (CLMM), desenvolvido pela Colaboração de Ciência da Energia Escura (DESC), para analisar a liberação pública de dados DES Y3. O pipeline de análise envolve:

1. Seleção de Dados:

   • Aglomerados: 9.962 aglomerados do catálogo redMaPPer (riqueza $20 < \lambda < 222$, desvio para o vermelho $0.2 < z < 0.65$) com alta probabilidade de centralização ($p_{cen} > 0.9$) para garantir sistemas relaxados e centralização precisa.
   • Galáxias de Fundo: Formas do catálogo de metacalibração DES Y3 e desvios para o vermelho fotométricos (photo-z) do catálogo Directional Neighborhood Fitting (DNF).
   • Validação: Catálogos simulados simulando halos Navaro–Frenk–White (NFW) elípticos e um subconjunto de aglomerados redMaPPer do DES Y1 para validação cruzada.

2. Determinação da Orientação:

   • A orientação do eixo maior de cada aglomerado é determinada pela distribuição de suas galáxias satélites membros (usando momentos de segunda ordem ponderados), em vez de apenas a galáxia central, para minimizar efeitos de má centralização.
   • As estimativas de cisalhamento de lente fraca são rotacionadas para alinhar com este eixo maior antes da empilhamento.

3. Modelagem Teórica (Além da Ordem Dominante):

   • A densidade de massa superficial é expandida usando uma decomposição em multipolos até a segunda ordem (quadrupolo).
   • Monopolo ($\Sigma_0$): Representa o componente circular. Os autores incluem um termo de correção de segunda ordem ($\epsilon^2$) ao monopolo, que accounts para a redistribuição de massa devido à elipticidade.
   • Quadrupolo ($\Sigma_2$): Representa a deformação elíptica.
   • A análise ajusta simultaneamente os sinais de monopolo e quadrupolo para restringir a massa, a concentração e a elipticidade projetada ($\epsilon$).

4. Correções Sistemáticas:

   • Fator de Diluição: Corrige o desalinhamento entre a distribuição de galáxias satélites (usada para definir a orientação) e a verdadeira orientação do halo de matéria escura. Isso é estimado via amostragem de Monte Carlo baseada na riqueza do aglomerado.
   • Fator de Boost: Corrige a contaminação de galáxias de fundo por membros do aglomerado não lensados (erros de desvio para o vermelho fotométrico), estimado usando a distribuição de desvio para o vermelho de galáxias em bins radiais.
   • Cortes Radiais: Fontes dentro de 300 kpc são excluídas para mitigar efeitos de mistura e de lente forte.

Principais Contribuições

• Primeira Medição de Forma de Conjunto DES Y3: Este trabalho apresenta a primeira medição da elipticidade projetada efetiva para um grande conjunto de aglomerados do catálogo redMaPPer do DES Y3 usando lente fraca.
• Análise Além da Ordem Dominante: A análise estende a expansão em multipolos até a segunda ordem da elipticidade para o termo de monopolo. Isso fornece uma correção de nível percentual ao perfil de massa, melhorando a precisão das medições de massa em comparação com aproximações de ordem dominante.
• Tratamento Robusto de Sistemáticos: O estudo implementa e valida um pipeline que corrige explicitamente os fatores de diluição (desalinhamento) e boost (contaminação), demonstrando a prontidão da ferramenta CLMM para futuras pesquisas em grande escala.
• Validação: A metodologia é validada contra catálogos simulados com formas de entrada conhecidas e contra dados do DES Y1, mostrando consistência e robustez contra variações nos métodos de medição (resposta por objeto vs. por HEALPix).

Resultados

• Eltipticidade Projetada: Para a amostra de referência, a elipticidade projetada do halo medida é $\epsilon = 0.310^{+0.017}_{-0.016}$, correspondendo a uma razão de eixos de $q = 0.527^{+0.018}_{-0.019}$. Isso representa uma detecção de $\sim 19\sigma$ de elipticidade não nula.
• Dependência de Propriedades: A análise não encontra nenhuma dependência estatisticamente significativa da elipticidade projetada na massa do aglomerado (riqueza) ou no desvio para o vermelho dentro do tamanho de amostra e precisão estatística atuais. As diferenças entre bins de alta/baixa riqueza e alta/baixa desvio para o vermelho são de aproximadamente $0.4\sigma$.
• Viés de Massa: Ao comparar o ajuste do modelo triaxial com um modelo esférico, o modelo triaxial produz uma estimativa de massa aproximadamente 5% menor para o mesmo sinal de monopolo, consistente com expectativas teóricas para perfis NFW elípticos.
• Validação com Simulações: O algoritmo de medição recupera as razões de eixos de entrada com um viés inferior a 5% para razões de eixos verdadeiras $q > 0.35$. A análise é encontrada insensível a variações na concentração do halo (testado para $c=2, 4, 6$).

Significância e Alegações
O artigo posiciona este trabalho como um passo preparatório crítico para o Levantamento de Legado de Espaço e Tempo (LSST) do Observatório Vera C. Rubin. Os autores afirmam que:

1. A ferramenta CLMM, agora equipada com um módulo de triaxialidade, é capaz de analisar de forma eficiente e consistente dados de lente fraca de aglomerados em grandes volumes.
2. A metodologia demonstrada é diretamente aplicável a futuras pesquisas de grande área (LSST, Euclid, Telescópio Espacial Roman).
3. Conjuntos de dados futuros com tamanhos de amostra maiores e observações mais profundas restringirão significativamente as elipticidades dos aglomerados (estimando uma redução de 82% nas barras de erro para o LSST Y10 em comparação com o DES Y3), permitindo testes mais precisos da dependência da forma do halo em relação à massa e ao desvio para o vermelho, e um melhor controle de sistemáticos na cosmologia de aglomerados.
4. Este artigo serve como o primeiro de uma série, com trabalhos futuros planejados para utilizar simulações cosmológicas (SkySim5000) para quantificar ainda mais vieses de seleção, má centralização e a relação entre a triaxialidade do halo e a estrutura em grande escala circundante.