FolpsD: combining EFT and phenomenological approaches for joint power spectrum and bispectrum analyses

1. Problema e Motivação

A extração de informação cosmológica a partir de dados de agrupamento de galáxias (como os do instrumento DESI - Dark Energy Spectroscopic Instrument) enfrenta um desafio central: como estender a faixa de escalas analisadas para além do regime linear, sem introduzir vieses sistemáticos significativos.

• Limitações Atuais: As análises de "forma completa" (full-shape) baseadas na Teoria de Perturbação (PT) e na Teoria de Campo Efetivo (EFT) são rigorosas, mas geralmente restritas a escalas de grande comprimento de onda (baixo $k$) para evitar efeitos não lineares e distorções no espaço de redshift (RSD) não modelados.
• O Dilema: Modelos fenomenológicos de amortecimento (damping) são frequentemente usados para estender a análise a escalas menores (maior $k$), capturando o efeito "Dedo de Deus" (FoG), mas podem violar princípios de invariância de Galileu ou introduzir degenerescências indesejadas com parâmetros cosmológicos.
• Objetivo: Desenvolver uma abordagem que combine a precisão da EFT com a flexibilidade de modelos fenomenológicos, utilizando estatísticas de ordem superior (bispectro) para quebrar degenerescências e melhorar as restrições cosmológicas.

2. Metodologia

Os autores desenvolveram e testaram o código FolpsD, uma extensão do código existente Folps (Fast Optimized Large-scale structure Perturbation Solver).

A. Modelo Teórico (FolpsD)

O modelo combina três componentes principais:

1. Espectro de Potência EFT com Amortecimento (LoS Damping):

   • Utiliza o espectro de potência de galáxias em um loop da EFT.
   • Introduz um fator de amortecimento fenomenológico de Lorentziana ao longo da linha de visada (LoS) para modelar a supressão de potência em pequenas escalas devido à dispersão de velocidades (efeito FoG).
   • A função de amortecimento é $D(x^2) = 1/(1+x^2)$, onde $x \propto k_\parallel \sigma_v$. Isso permite modelar o espectro até $k \sim 0.3 \, h\text{Mpc}^{-1}$, além do limite típico da EFT padrão ($\sim 0.2 \, h\text{Mpc}^{-1}$).
   • Diferencia-se de modelos puramente fenomenológicos (como TNS) ao manter a estrutura completa dos correlatores de um-loop da EFT, garantindo consistência com a conservação de momento e invariância de Galileu em grandes escalas.

2. Bispectro de Galáxias (Nível Árvore):

   • Implementa o bispectro de galáxias no nível de árvore (tree-level), projetado na base de harmônicos esféricos tripolares (base de Sugiyama).
   • Inclui termos de contramedida de ordem superior (NLO) e fatores de amortecimento projetados em cada vetor de onda do triângulo de configuração.
   • O uso da base de Sugiyama permite uma análise eficiente dos multipolos do bispectro, focando principalmente no monopolo ($B_{000}$) e no quadrupolo ($B_{202}$).

3. Efeitos Alcock-Paczynski e Viés de Galáxia:

   • O modelo incorpora corretamente o efeito Alcock-Paczynski (AP) para corrigir distorções geométricas entre a cosmologia fiducial e a verdadeira.
   • Utiliza expansões de viés de galáxia de segunda ordem (parâmetros $b_1, b_2, b_s$) e inclui tratamentos para diferentes esquemas de viés (McDonald, class-pt).

B. Validação e Dados

• Simulações: O modelo foi testado usando 25 realizações independentes de catálogos mocks de segunda geração da simulação AbacusSummit, configurados para replicar as amostras de galáxias Luminosas Vermelhas (LRG) e Quasares (QSO) do DESI.
• Dados Reais: A componente de espectro de potência foi aplicada aos dados reais do DESI Data Release 1 (DR1), cobrindo um volume efetivo de 17.5 Gpc$^3$.
• Análise Estatística: Utilização de cadeias de Markov (MCMC) via Cobaya e desilike, com marginais analíticas sobre parâmetros de ruído e contramedidas. Foram calculados o Figure of Merit (FoM) e o Figure of Bias (FoB) para avaliar precisão e viés.

3. Principais Contribuições

1. Desenvolvimento do FolpsD: Criação de um código unificado que integra EFT de um-loop, amortecimento fenomenológico LoS e bispectro de árvore na base de Sugiyama.
2. Análise Conjunta Pk+Bk: Demonstração de que a combinação de espectro de potência e bispectro reduz significativamente as degenerescências entre parâmetros de viés e cosmológicos, superando a análise de apenas espectro de potência.
3. Avaliação de Limites do Amortecimento: Identificação crítica de que, embora o amortecimento melhore as restrições para amostras de alto sinal-ruído (LRGs), ele pode introduzir viés em amostras dominadas por ruído (QSOs) ou em cenários onde a forma do espectro degenera com o amortecimento (ex: massa de neutrinos).
4. Extensão da Faixa de Escalas: Validação de que o bispectro pode ser incluído até escalas maiores ($k_{max} \approx 0.24 \, h\text{Mpc}^{-1}$ para LRGs) sem viés significativo, especialmente quando combinado com o espectro de potência amortecido.

4. Resultados Chave

Em Simulações (Mocks Abacus)

• LRGs (Alto Sinal-Ruído):
  • A análise conjunta FolpsD-Pk+Bk (com amortecimento) fornece restrições ~30% mais apertadas nos parâmetros $A_s$ (amplitude das flutuações) e $\omega_{cdm}$ (densidade de matéria escura fria) em comparação com a EFT padrão sem bispectro.
  • O modelo é robusto até $k \approx 0.3 \, h\text{Mpc}^{-1}$ no espectro de potência e $k \approx 0.24 \, h\text{Mpc}^{-1}$ no bispectro, mantendo o viés abaixo de $1\sigma$.
  • O amortecimento ajuda a restringir fortemente os parâmetros de viés de segunda ordem ($b_2, b_s$), o que, por sua vez, melhora a inferência cosmológica.
• QSOs (Baixo Sinal-Ruído):
  • Para quasares, o modelo com amortecimento introduz viés significativo nos parâmetros cosmológicos (desvios $>1\sigma$), mesmo com cortes de escala conservadores.
  • O parâmetro de amortecimento absorve flutuações de ruído, distorcendo a inferência. Recomenda-se cautela ou uso de modelos EFT padrão para amostras de baixo sinal-ruído.
• Massa de Neutrinos:
  • O amortecimento fenomenológico degrada as restrições sobre a soma das massas de neutrinos ($\sum m_\nu$) devido à degenerescência entre o amortecimento e a supressão do espectro de potência causada pelos neutrinos. A inclusão do bispectro ajuda a mitigar, mas não eliminar totalmente, esse efeito.

Em Dados Reais (DESI DR1)

• Cosmologia $\Lambda$CDM:
  • Aplicando o FolpsD ao espectro de potência do DESI DR1, obtiveram-se melhorias de até 20% na precisão de $\sigma_8$ e $\Omega_m$ em comparação com a análise EFT padrão.
• Cosmologia $w_0w_a$CDM (Energia Escura Dinâmica):
  • Utilizando o método ShapeFit (que comprime a informação de forma em parâmetros físicos), a análise com FolpsD resultou em restrições 15% mais apertadas em $w_0$ e 21% mais apertadas em $w_a$.
  • O Figure of Merit (FoM) para energia escura aumentou de 68 (EFT) para 161 (FolpsD).
  • Descoberta Importante: A análise de forma completa com FolpsD indicou um desvio da energia escura constante ($w = -1$) em um nível ligeiramente superior a $1\sigma$, usando apenas dados de estrutura em grande escala (LSS).

5. Significância e Conclusão

O trabalho demonstra que a combinação de estatísticas de dois pontos (espectro de potência) e três pontos (bispectro), aliada a modelos fenomenológicos controlados de amortecimento, é uma ferramenta poderosa para a cosmologia de próxima geração.

• Eficiência: Permite extrair mais informação cosmológica dos mesmos dados, estendendo a faixa de escalas válidas.
• Robustez: Identifica claramente onde os modelos fenomenológicos falham (amostras ruidosas, neutrinos), fornecendo diretrizes claras para a análise de dados futuros do DESI (DR2 e além).
• Impacto Imediato: Os resultados sugerem que a energia escura dinâmica pode ser detectada com maior significância estatística ao se utilizar abordagens de "forma completa" avançadas, como o FolpsD, antes mesmo da conclusão completa da missão DESI.

O código FolpsD é disponibilizado publicamente, facilitando a adoção dessas técnicas pela comunidade científica.