Forward analytical model for the optical selection bias on galaxy cluster lensing profiles

Título: Um modelo analítico forward para o viés de seleção óptica em perfis de lentes de aglomerados de galáxias

Autores: M. Costanzi et al.

---

1. O Problema

Os catálogos de aglomerados de galáxias selecionados por propriedades ópticas (baseados em contagem de galáxias ou "riqueza", $\lambda$) são suscetíveis a viéses de seleção, principalmente decorrentes de efeitos de projeção ao longo da linha de visada.

• Mecanismo do Viés: Estruturas não resolvidas (outros aglomerados ou halos de matéria escura) que se sobrepõem à linha de visada do aglomerado alvo podem ser erroneamente contadas como membros do aglomerado. Isso infla a riqueza observada ($\lambda_{ob}$) em relação à riqueza intrínseca ($\lambda_{tr}$).
• Consequência Cosmológica: Esse efeito não apenas altera a contagem de aglomerados, mas também distorce o perfil de densidade projetada medido via lentes gravitacionais fracas. Aglomerados selecionados opticamente tendem a residir em ambientes mais densos ao longo da linha de visada do que uma amostra aleatória de mesma massa, levando a uma superestimação do sinal de lente (bias de 10-20% em escalas relevantes para calibração de massa).
• Limitação Atual: Abordagens anteriores utilizavam simulações baseadas ou assumiam templates analíticos com parâmetros livres marginalizados. Falta um modelo analítico preditivo, sem parâmetros livres, que ligue diretamente os viéses de seleção às especificações do levantamento e ao modelo cosmológico.

2. Metodologia

Os autores desenvolveram um modelo analítico forward (de frente para trás) que quantifica o impacto dos efeitos de projeção nos perfis de densidade. O modelo baseia-se na premissa de que as mesmas estruturas que contribuem para o aumento da riqueza observada também contribuem para o perfil de densidade projetado.

Componentes Principais do Modelo:

1. Definição de Riqueza: A riqueza observada é modelada como $\lambda_{ob} = \lambda_{tr} + \Delta_{prj}$, onde $\Delta_{prj}$ é a contaminação devido a halos projetados.
2. Viés de Aglomerado Óptico ($b_{\lambda_{ob}}$):
   • O modelo introduz uma parametrização dependente da escala para o viés dos aglomerados selecionados opticamente.
   • Utiliza uma função sigmoide para descrever a transição do viés em pequenas escalas (dentro do raio do aglomerado) para grandes escalas (ambiente externo).
   • O viés em pequena escala ($b^0_{\lambda_{ob}}$) depende do "boost" de riqueza ($\Delta_{prj}$) em relação ao fundo esperado.
   • O viés em grande escala ($b^\infty_{\lambda_{ob}}$) é modelado como uma função linear do excesso de contaminação de riqueza ($\delta_{prj}$) em relação a uma amostra aleatória.
3. Perfil de Densidade Projetada ($\langle \Sigma \rangle$):
   • O perfil total é decomposto em dois termos:
     • Termo de 1-halo: Contribuição do próprio halo do aglomerado.
     • Termo de 2-halos ($\langle \Sigma_{prj} \rangle$): Contribuição integrada de halos correlacionados ao longo da linha de visada, ponderados pelo viés $b_{\lambda_{ob}}$ e pela geometria de sobreposição.
4. Validação: O modelo foi testado contra um catálogo sintético (mock) baseado na simulação N-body de DeRose et al., configurado para imitar o catálogo de aglomerados do Dark Energy Survey (DES) Ano 3 (Y3) usando o algoritmo redMaPPer.

3. Contribuições Chave

• Modelo Preditivo Sem Parâmetros Livres: Diferente de abordagens anteriores que marginalizavam sobre parâmetros de viés, este modelo deriva o viés diretamente do modelo de ocupação de halos (HOD), cosmologia e especificações do algoritmo de detecção (redMaPPer).
• Parametrização do Viés de Escala: A introdução de uma forma funcional sigmoide para $b_{\lambda_{ob}}(\theta)$, onde a transição é governada pelo tamanho angular do aglomerado e pela magnitude dos efeitos de projeção.
• Conexão Direta entre Riqueza e Ambiente: O modelo demonstra explicitamente como o excesso de riqueza ($\Delta_{prj}$) se correlaciona com um ambiente mais denso ao longo da linha de visada, afetando o perfil de lentes.
• Ferramenta para Inferência Cosmológica: O trabalho fornece uma metodologia consistente para corrigir viéses sistemáticos em análises cosmológicas baseadas em aglomerados ópticos, ligando as especificações do levantamento à cosmologia subjacente.

4. Resultados

A validação contra os dados sintéticos do DES Y3 mostrou que o modelo:

• Reproduz a Dependência de Escala: Captura com precisão a dependência do termo de 2-halos em relação aos aumentos de riqueza induzidos por projeção.
• Recupera o Perfil de Densidade: O modelo prevê corretamente o perfil de densidade projetado total, incluindo o viés relativo a uma amostra aleatória de mesma distribuição de massa. O viés apresenta uma forma "em sino" característica, atingindo o pico em torno de $R \sim 1 h^{-1}$ Mpc (escala de transição entre 1-halo e 2-halos).
• Evolução com Redshift e Riqueza: O modelo reproduz a evolução do viés com o redshift e a riqueza. Aglomerados mais ricos (e, portanto, mais massivos) residem em ambientes mais densos, resultando em um viés maior.
• Aproximação Analítica: Os autores fornecem uma função de ajuste analítica (dupla lei de potência com transição suave) para o viés no perfil de densidade, facilitando sua implementação em pipelines de inferência cosmológica onde o cálculo integral completo seria computacionalmente proibitivo.

5. Significado e Conclusões

Este trabalho estabelece um framework analítico robusto para interpretar e prever viéses de seleção em catálogos de aglomerados fotométricos.

• Impacto na Cosmologia: Ao permitir a modelagem consistente dos efeitos de projeção, o modelo reduz as incertezas sistemáticas na calibração de massa de aglomerados, essencial para restringir parâmetros cosmológicos (como $\sigma_8$ e $\Omega_m$) a partir de levantamentos ópticos.
• Generalização: Embora focado no redMaPPer e no DES, o framework é generalizável para outros proxies de massa óptica (como massa estelar ponderada) e outros algoritmos de detecção, desde que o kernel de projeção e a distribuição de probabilidade sejam calibrados.
• Limitações e Futuro: O modelo atual negligencia efeitos de desalinhamento do centro (miscentering) e viéses de orientação/triaxialidade dos halos. No entanto, ele fornece a base física necessária para incorporar essas complexidades em futuras análises de precisão.

Em resumo, o artigo oferece uma solução analítica elegante e validada para um dos principais desafios sistemáticos na cosmologia de aglomerados de galáxias: a correção dos efeitos de projeção que distorcem tanto a contagem quanto a distribuição de massa observada.