Euclid preparation: The flat-sky approximation for the clustering of Euclid's photometric galaxies

Título: Aproximação de Céu Plano para o Agrupamento das Galáxias Fotométricas do Euclid

Autores: Euclid Collaboration (W. L. Matthewson et al.)
Data: Outubro de 2025

---

1. Problema e Contexto

O satélite Euclid gerará um catálogo massivo de galáxias com redshifts fotométricos (aproximadamente $1,5 \times 10^9$ galáxias). Para extrair informações cosmológicas ótimas, é necessário calcular espectros de potência angulares teóricos que sejam comparáveis aos dados observados.

• Desafio Computacional: O cálculo completo (full calculation) dos espectros de potência é computacionalmente proibitivo para análises estatísticas extensas, como inferência de parâmetros via MCMC (Cadeia de Markov Monte Carlo), que exigem milhares de avaliações do modelo.
• Soluções Atuais: A Aproximação de Limber é amplamente utilizada por ser extremamente rápida (cerca de 100 vezes mais rápida que o cálculo completo), mas é conhecida por falhar em bins de redshift estreitos, subestimando ou superestimando significativamente os espectros em grandes escalas (baixos $\ell$).
• Objetivo: Avaliar se a Aproximação de Limber é suficientemente precisa para as futuras liberações de dados do Euclid (DR1 e DR3) ou se a Aproximação de Céu Plano (Flat-Sky) é necessária para evitar viéses nos resultados cosmológicos.

2. Metodologia

Os autores compararam três métodos de cálculo dos espectros de potência angular ($C_\ell$) para o número de contagens de galáxias:

1. Cálculo Completo: Utilizando o código CLASS (incluindo todos os termos relativísticos, distorções de espaço-vermelho e efeitos de lente).
2. Aproximação de Limber: Substitui a integral sobre o número de onda $k$ por uma função delta de Dirac, assumindo que apenas modos transversais contribuem.
3. Aproximação de Céu Plano (Flat-Sky): Uma extensão refinada (baseada em Gao et al., 2024) que trata o céu como um plano, substituindo harmônicos esféricos por transformadas de Fourier 2D. Esta abordagem mantém a componente paralela de $k$ ($k_\parallel$), corrigindo os erros da Limber em grandes escalas, mas ainda oferece ganhos de velocidade em relação ao cálculo completo.

Configurações de Dados Simuladas:

• DR1 (Primeira Liberação): 6 bins de redshift equi-populados, corte de magnitude $I_E < 23.5$. Representa um cenário otimista inicial.
• DR3 (Terceira e Última Liberação): 13 bins de redshift equi-populados, corte de magnitude mais profundo $I_E < 24.5$. Representa o cenários final da missão com bins mais estreitos.

Análise Estatística:

• Foram geradas 1000 realizações de ruído gaussiano baseadas na matriz de covariância teórica.
• Calculou-se o $\chi^2$ para cada aproximação em comparação com o "verdadeiro" espectro (cálculo completo).
• Avaliou-se a distribuição de $\Delta\chi^2$ (diferença entre o $\chi^2$ da aproximação e o do modelo completo) para determinar se as aproximações introduzem viéses estatisticamente significativos.

3. Principais Resultados

Cenário DR1 (6 bins)

• Desempenho da Limber: Aproximação de Limber é suficientemente precisa. As discrepâncias entre a Limber e o cálculo completo caem dentro da faixa de erro de $1\sigma$ esperada para a configuração do survey.
• Desempenho do Céu Plano: Aproximação de Céu Plano é altamente precisa (< 5% de erro), mas não oferece vantagem prática sobre a Limber neste caso, pois a precisão extra não é detectável dada a incerteza observacional e o menor número de bins.
• Conclusão DR1: A Limber pode ser usada sem viés significativo para a primeira liberação de dados.

Cenário DR3 (13 bins)

• Desempenho da Limber: Aproximação de Limber falha. Para bins estreitos (especialmente a partir do 4º bin e em correlações de alto redshift), a Limber não consegue recuperar o espectro dentro da margem de $1\sigma$.
  • O valor de $\Delta\chi^2$ para a Limber atinge picos de ~410, rejeitando a hipótese nula com $p < 0.001$ (confiança de 99%).
  • O erro é mais pronunciado em grandes escalas ($\ell \lesssim 100$) e em correlações entre bins de redshift diferentes.
  • A falha é atribuída principalmente à incapacidade da Limber de modelar corretamente as Distorções de Espaço-Verme (RSD) e termos de densidade em bins estreitos, onde a componente $k_\parallel$ não é desprezível.
• Desempenho do Céu Plano: Aproximação de Céu Plano mantém alta precisão (< 5% de erro) em todos os casos.
  • O $\Delta\chi^2$ médio é de apenas ~5, com um valor-p de 0.131, indicando que a distribuição é consistente com zero (sem viés significativo).
• Conclusão DR3: A Limber é inadequada para a análise final do Euclid; a Aproximação de Céu Plano é necessária.

4. Contribuições Técnicas e Significância

1. Validação da Aproximação de Céu Plano: O artigo demonstra que a versão "recalibrada" da aproximação de céu plano é robusta o suficiente para substituir o cálculo completo em análises de dados fotométricos de alta precisão, mantendo a precisão física necessária para a cosmologia de precisão.
2. Definição de Limites de Uso da Limber: Estabelece claramente que a Aproximação de Limber é válida para bins largos (DR1), mas torna-se um viés sistemático inaceitável para bins estreitos (DR3), orientando a escolha de métodos para futuras análises do Euclid.
3. Eficiência Computacional:
   • A Limber é ~200x mais rápida que o cálculo completo.
   • A Aproximação de Céu Plano (na implementação atual) é apenas ~2x mais rápida que o cálculo completo, o que ainda é lento para MCMC.
   • Solução Proposta: Os autores sugerem que, como as maiores discrepâncias ocorrem em baixos $\ell$, uma abordagem híbrida (usando Limber em altas escalas e Céu Plano em baixas escalas) ou otimizações futuras (como pré-cálculo de integrais, similar ao código BLAST) podem tornar a Aproximação de Céu Plano viável para inferência de parâmetros em tempo hábil.
4. Impacto na Cosmologia: A precisão na modelagem dos espectros de potência é crucial para evitar viéses na inferência de parâmetros cosmológicos (como $S_8$ e $H_0$). O uso incorreto da Limber em dados DR3 poderia levar a conclusões cosmológicas erradas.

Conclusão Final

O estudo conclui que, embora a Aproximação de Limber seja adequada para a primeira fase de dados do Euclid (DR1), ela é insuficiente para a análise final (DR3) devido à necessidade de bins de redshift mais estreitos. A Aproximação de Céu Plano é identificada como a alternativa precisa e necessária para a análise de dados do Euclid, embora ainda requeira otimizações de código para atingir a velocidade computacional necessária para análises de inferência de parâmetros em larga escala.