Euclid preparation. Full-shape modelling of 2-point and 3-point correlation functions in real space

Resumo Técnico: Modelagem de Forma Completa (Full-Shape) de Funções de Correlação de 2 e 3 Pontos no Espaço Real para a Missão Euclid

1. Problema e Contexto

A missão Euclid da Agência Espacial Europeia visa mapear a estrutura em grande escala do universo com precisão sem precedentes, utilizando aglomerados de galáxias e lentes gravitacionais fracas. Para extrair informações cosmológicas robustas, é necessário analisar estatisticamente as flutuações na distribuição de galáxias.

• O Desafio: Tradicionalmente, análises focam na função de correlação de 2 pontos (2PCF) ou no espectro de potência. A inclusão de estatísticas de ordem superior, como a função de correlação de 3 pontos (3PCF) no espaço de configuração (ou o espectro cruzado no espaço de Fourier), é crucial para capturar propriedades não gaussianas e quebrar degenerescências entre parâmetros cosmológicos e de viés (bias).
• Obstáculo Computacional: A análise de "forma completa" (full-shape) da 3PCF no espaço de configuração é computacionalmente proibitiva devido ao custo de avaliar modelos teóricos (perturbação cosmológica) e à complexidade de estimadores diretos em grandes volumes de dados. Além disso, a validade de modelos perturbativos em pequenas escalas (regime não linear) precisa ser rigorosamente testada antes da aplicação aos dados reais do Euclid.

2. Metodologia

Os autores desenvolveram uma estrutura de análise integrada para validar modelos teóricos e realizar uma análise conjunta de 2PCF e 3PCF no espaço real.

• Dados de Simulação:

  • Utilizaram quatro "snapshots" (instantâneos) da simulação Flagship I (N-body), com caixas de $L = 3780 \, h^{-1}$ Mpc, em redshifts $z = \{0.9, 1.2, 1.5, 1.8\}$.
  • As simulações foram preenchidas com galáxias sintéticas usando uma distribuição de ocupação de halos (HOD) que reproduz a população esperada de galáxias emissoras de $H\alpha$ do Euclid (fluxo limite $f_{H\alpha} > 2 \times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$).
  • O estudo considera o caso ideal, sem efeitos observacionais (como pureza ou completude da amostra).

• Modelos Teóricos:

  • 2PCF: Modelada no espaço real até a ordem de um laço (one-loop) dentro da Teoria de Perturbação Cosmológica (PT) e da Teoria de Campo Efetivo da Estrutura em Grande Escala (EFTofLSS). Inclui expansão de viés até termos não locais e derivativos de ordem superior, além de contratermos.
  • 3PCF: Modelada no nível de árvore (tree-level) da PT.
  • Re-sumação Infravermelha (IR): Implementada para tratar o alargamento das oscilações acústicas de bárions (BAO) devido à evolução não linear.

• Estimadores e Emuladores:

  • Para a 3PCF, utilizaram o estimador baseado em decomposição em harmônicos esféricos (SHD) de Slepian & Eisenstein (2015), que reduz a complexidade computacional de $O(N_g^3)$ para $O(N_g^2)$.
  • Emulador: Desenvolveram um emulador baseado em PyTorch capaz de gerar previsões rápidas e robustas para 2PCF e 3PCF em função de três parâmetros cosmológicos ($A_s$, $\omega_{cdm}$, $h$) e parâmetros de viés. Isso permitiu a amostragem eficiente da função de verossimilhança conjunta.

• Análise Estatística:

  • Realizaram uma análise de verossimilhança conjunta (joint likelihood) combinando 2PCF e 3PCF.
  • Avaliaram a validade do modelo variando o corte de escala mínima ($r_{min}$) e excluindo configurações de triângulos quase isósceles através de um parâmetro $\eta_{min}$.
  • Métricas de desempenho: Goodness-of-Fit ($\chi^2$), Figure of Bias (FoB) e Figure of Merit (FoM).

3. Principais Contribuições

1. Primeira Análise Conjunta Full-Shape no Espaço Real: Este trabalho apresenta a primeira análise conjunta de 2PCF e 3PCF no espaço real para restringir parâmetros cosmológicos, superando a barreira computacional anterior através do uso de emuladores.
2. Validação de Modelos Perturbativos: Determinaram o limite de validade dos modelos perturbativos para a 3PCF no espaço real, identificando escalas onde a teoria falha devido a efeitos não lineares e limitações na expansão em polinômios de Legendre.
3. Desenvolvimento de Emulador: Criaram e validaram um emulador eficiente para a 3PCF, essencial para futuras análises de dados do Euclid que exigem amostragem de espaços de parâmetros de alta dimensão.
4. Estudo de Viés e Degenerescências: Demonstraram como a inclusão da 3PCF quebra a degenerescência entre a amplitude das perturbações escalares ($A_s$) e o viés linear ($b_1$), um problema persistente nas análises apenas de 2PCF.

4. Resultados Chave

• Validade do Modelo (Goodness-of-Fit):

  • O modelo de 2PCF (um laço) ajusta bem os dados em todas as escalas consideradas.
  • O modelo de 3PCF (nível de árvore) falha para separações abaixo de $20 \, h^{-1}$ Mpc em $z=0.9$.
  • Para a 3PCF, a exclusão de triângulos quase isósceles (via $\eta_{min}$) é crítica. Um corte de $r_{min} = 30 \, h^{-1}$ Mpc com $\eta_{min} = 3$ (excluindo triângulos isósceles) fornece os resultados mais robustos e sem viés sistemático.
  • Uma configuração mais agressiva ($r_{min} = 30 \, h^{-1}$ Mpc, $\eta_{min} = 2$) é aceitável, mas mostra tensões de ~68% no parâmetro $A_s$ em baixos redshifts.

• Restrições Cosmológicas:

  • A análise conjunta (2PCF + 3PCF) fornece restrições significativamente mais apertadas (maior FoM) em comparação com a 2PCF isolada.
  • A degenerescência entre $A_s$ e $b_1$ é efetivamente quebrada sem a necessidade de priores informativos externos sobre o viés.
  • Os parâmetros cosmológicos ($A_s$, $\omega_{cdm}$, $h$) são recuperados com precisão, com viés (FoB) dentro dos limites de confiança esperados para cortes de escala conservadores.

• Relações de Viés:

  • A imposição de relações teóricas entre parâmetros de viés (ex: $b_{G2}$ como função de $b_1$) reduz o espaço de parâmetros e melhora a precisão das restrições cosmológicas, especialmente para $A_s$, sem introduzir viés significativo nos parâmetros principais.

5. Significado e Impacto

Este trabalho é um passo fundamental na preparação para a análise de dados da missão Euclid.

• Viabilidade Técnica: Demonstra que a análise conjunta de 2PCF e 3PCF no espaço real é viável e computacionalmente tratável graças aos emuladores desenvolvidos.
• Guia para Dados Reais: Estabelece cortes de escala recomendados ($r_{min} \approx 30-40 \, h^{-1}$ Mpc para 3PCF) e estratégias para lidar com configurações de triângulos problemáticas, garantindo que as futuras medições do Euclid não sejam contaminadas por erros sistemáticos do modelo teórico.
• Melhoria Cosmológica: Confirma que a adição da 3PCF ao conjunto de dados melhora a precisão na determinação de parâmetros cosmológicos fundamentais, tornando a missão Euclid ainda mais poderosa para investigar a natureza da energia escura e a gravidade em grande escala.
• Base para Futuras Análises: A metodologia e o emulador desenvolvidos servirão de base para análises futuras no espaço de redshift (que incluem distorções de espaço de redshift) e para a combinação com lentes gravitacionais fracas.

Em suma, o artigo valida a infraestrutura teórica e computacional necessária para explorar o potencial máximo da estatística de 3 pontos na cosmologia de precisão do futuro.