Analytical method for computing the covariance matrix of cosmic shear two-point correlation function

Resumo Técnico: Método Analítico para Computação da Matriz de Covariância da Função de Correlação de Dois Pontos do Cisalhamento Cósmico

Enunciação do Problema
A estimativa precisa da matriz de covariância para estatísticas de cisalhamento cósmico é um pré-requisito para análises cosmológicas utilizando levantamentos de lentes gravitacionais fracas de grande área atuais e futuros. Embora a covariância de amostra possa ser medida a partir de numerosas realizações de simulação, essa abordagem torna-se computacionalmente proibitiva para levantamentos do Estágio IV com áreas de céu superiores a 10.000 deg², onde poucas realizações independentes podem ser extraídas de simulações de céu completo. Consequentemente, a modelagem analítica da covariância é necessária. No entanto, máscaras de levantamento introduzem efeitos de acoplamento de modos que complicam o cálculo do termo de covariância Gaussiana, particularmente para funções de correlação de dois pontos (2PCFs) no espaço real. Métodos analíticos existentes frequentemente lutam para reproduzir com precisão a estrutura de covariância na presença de geometrias de levantamento complexas.

Metodologia
Os autores investigam métodos analíticos para computar a matriz de covariância Gaussiana da 2PCF do cisalhamento cósmico, $\xi(\theta)$, levando em conta a geometria finita do levantamento (efeitos da janela do levantamento). O estudo utiliza 1.000 realizações Gaussianas de campos de cisalhamento cósmico geradas com máscaras do levantamento HSC Ano 3 (HSC-Y3) como referência para testar previsões analíticas.

O quadro analítico procede da seguinte forma:

1. Mapeamento do Harmônico para o Espaço Real: O estudo estabelece a relação entre o espectro de potência convolvido pela janela, $\tilde{C}_\ell$, e a função de correlação medida, $\xi_i$. A função de correlação observada é tratada como um estimador não viesado corrigido por um fator dependente do número de pares de galáxias na região do levantamento ($N_{pair}$) e da densidade numérica média ($n$), derivado do espectro de potência afetado pela janela $\tilde{C}_\ell$.
2. Projeção da Covariância: A covariância da 2PCF no espaço real, $\text{Cov}(\xi_i, \xi_j)$, é formalmente expressa como uma transformada dupla de Legendre da covariância no espaço harmônico, $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$.
3. Estimadores Analíticos: Três abordagens específicas para estimar a covariância no espaço harmônico são avaliadas:
   • Aproximação de Núcleo Estreito Melhorada (iNKA): Assume que a matriz de acoplamento de modos é fortemente localizada em torno de $\ell \simeq \ell'$, permitindo o cálculo rápido de $\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$.
   • Aproximação $f_{sky}$: Uma aproximação padrão assumindo que a covariância é diagonal e escalada pela fração do céu.
   • Método de Contagens Quarticas Ponderadas: Calcula a covariância da 2PCF diretamente no espaço de configuração usando contagens quarticas de pares de galáxias ponderadas por $\xi \times \xi$, seguindo o formalismo de Schneider et al. (2002).

Principais Resultados
O estudo compara as previsões analíticas contra a covariância de amostra medida diretamente a partir das 1.000 simulações Gaussianas:

• Validação do Quadro: A relação entre o espectro de potência convolvido pela janela e a função de correlação (Eq. 2.6) foi validada, confirmando que o quadro pseudo-$C_\ell$ descreve corretamente os efeitos da janela do levantamento e fornece um mapeamento preciso entre o espaço harmônico e o espaço real.
• Limitações da iNKA: Embora a iNKA reproduza com precisão os componentes diagonais da covariância no espaço harmônico ($\text{Cov}(\tilde{C}_\ell, \tilde{C}_{\ell'})$) dentro de 10%, ela falha em modelar com precisão os componentes fora da diagonal. Quando essas covariâncias no espaço harmônico são projetadas no espaço real via transformação de Legendre, as imprecisões residuais nos termos fora da diagonal propagam-se significativamente. Consequentemente, a covariância derivada da iNKA para $\xi(\theta)$ não reproduz a covariância medida a partir das simulações, mesmo para componentes diagonais da matriz de covariância da 2PCF.
• Superioridade das Contagens Quarticas: Em contraste, o método de contagens quarticas ponderadas demonstra concordância significativamente melhor com os resultados das simulações. Este método, que modela explicitamente correlações baseadas em pares no espaço de configuração, reproduz com sucesso a estrutura de covariância sem depender da suposição de dominância diagonal da iNKA.
• Aproximação $f_{sky}$: A aproximação $f_{sky}$ também mostra discrepâncias em comparação com as simulações, embora o método de contagens quarticas seja destacado como o mais preciso entre os estimadores analíticos testados.

Significado e Alegações
O artigo conclui que a modelagem precisa da estrutura fora da diagonal da covariância no espaço harmônico é crucial para obter estimativas de covariância confiáveis de estatísticas de lentes gravitacionais fracas no espaço real quando efeitos de janela do levantamento estão presentes. Os autores demonstram que simplesmente estender a iNKA para estatísticas no espaço real é insuficiente devido à amplificação de erros fora da diagonal durante a transformação de Legendre.

A contribuição primária deste trabalho é a identificação do modo de falha específico da iNKA no contexto da covariância da 2PCF e a validação do método de contagens quarticas ponderadas como uma alternativa analítica mais robusta para levantamentos do Estágio IV. Os resultados implicam que futuros levantamentos de grande área que exigem inferência precisa de parâmetros cosmológicos devem empregar métodos analíticos que capturem o acoplamento de modos além de aproximações diagonais simples no espaço harmônico.