Angular bispectrum of matter number counts in cosmic structures

Resumo Técnico: Bispectro Angular das Contagens de Número de Matéria em Estruturas Cósmicas

Enunciado do Problema
A função de correlação de dois pontos (espectro de potência) continua sendo a principal ferramenta para restringir o modelo padrão $\Lambda$CDM. No entanto, à medida que a estrutura de grande escala (LSS) evolui para o regime fracamente não linear, estatísticas de ordem superior, especificamente a função de correlação de três pontos e sua contraparte no espaço de Fourier, o bispectro, tornam-se essenciais para sondar as correlações de fase de Fourier e os acoplamentos de modos induzidos pela evolução gravitacional não linear e não-gaussianidade primordial (PNG).

As análises atuais de bispectros de galáxias dependem fortemente da aproximação de céu plano local, da aproximação de Limber e da dinâmica Newtoniana. Essas aproximações falham para levantamentos de próxima geração (ex: Euclid, DESI, SPHEREx) que cobrem áreas amplas e utilizam bins de redshift finitos (tomografia). Além disso, embora os efeitos relativísticos (projeção, Relatividade Geral dinâmica e radiação) sejam teoricamente conhecidos, seu impacto preciso no bispectro angular em uma configuração de bin de redshift finito sem a aproximação de Limber ainda não foi totalmente computado. Adicionalmente, extrair esses sinais relativísticos sutis de simulações é desafiador devido a artefatos numéricos e à dominância dos termos Newtonianos.

Metodologia
Os autores apresentam um cálculo de céu total do bispectro angular para contagens de número de matéria usando teoria de perturbação de segunda ordem. A metodologia é caracterizada pelas seguintes escolhas técnicas:

1. Estrutura (Framework): O cálculo utiliza "harmônicos esféricos tomográficos" (TSH), empregando funções de janela para definir bins de redshift finitos ($1.75 \le z \le 2.25$ e $0.55 \le z \le 0.65$). Esta abordagem evita a aproximação de Limber e a suposição de céu plano.
2. Expansão de Perturbação: A perturbação da contagem de número $\Delta$ é expandida até segunda ordem. Os autores organizam os termos de acordo com o parâmetro de expansão de campo fraco $H/k$:
   • Termos Newtonianos ($O((H/k)^0)$): Inclui densidade linear, distorções no espaço de Redshift (RSD) e termos quadráticos (acoplamentos densidade-densidade, densidade-velocidade).
   • Efeitos de Projeção Relativística ($O(H/k)$): Efeitos de projeção de ordem líder (Doppler, Sachs-Wolfe, etc.) até primeira ordem na expansão de campo fraco.
   • Efeitos de GR Dinâmica e Radiação ($O((H/k)^2)$ e superiores): Inclui correções de GR dinâmica não linear e correções radiativas (seguindo Tram et al. 2016).
   • Termos de Acoplamento: Inclui acoplamentos entre efeitos de projeção e dinâmicos até $O((H/k)^3)$.
3. Exclusões: Para simplificação e para permitir verificações de consistência direta, os autores negligenciam o viés de traçadores (tracer bias), efeitos integrados ao longo da linha de visada (lente gravitacional, efeito Sachs-Wolfe integrado) e termos da posição do observador. Eles também omitem efeitos de projeção de próxima ordem devido à complexidade e falta de consenso na literatura.
4. Implementação Numérica: Os autores estendem o framework matemático eficiente de Assassi et al. (2017), utilizando decomposições FFTLog e avaliações de funções hipergeométricas para reduzir o problema a integrais unidimensionais. Eles fornecem um código público, ang_bispec, para computar as 14 contribuições distintas de árvore (tree-level) do bispectro.
5. Validação: As previsões teóricas são comparadas com simulações de N-corpos relativísticas (Montandon et al. 2023). Para garantir uma comparação justa, os autores removem efeitos de projeção integrados dos dados da simulação e aplicam binagem e suavização Gaussiana às previsões teóricas para corresponder ao estimador da simulação (Bucher et al. 2016).

Principais Contribuições

• Primeiro Cálculo de Céu Total: Este trabalho fornece o primeiro cálculo de céu total do bispectro angular em teoria de perturbação de segunda ordem sem a aproximação de Limber, formulado para bins de redshift finitos.
• Novo Cálculo Newtoniano: Mesmo a parte Newtoniana do bispectro dentro deste setup específico (bins finitos, sem Limber) é reportada como novel.
• Inclusão Relativística Abrangente: Os autores incluem, até segunda ordem, todos os efeitos de GR dinâmica, correções de radiação e efeitos de projeção relativística de ordem líder dentro de um framework unificado.
• Liberação de Código: O código ang_bispec é disponibilizado publicamente, estendendo computações anteriores para todos os termos Newtonianos dominantes e efeitos relativísticos.

Resultados
Os autores avaliaram o bispectro para dois bins de redshift e compararam as contribuições de vários termos:

1. Dominância de Termos Newtonianos: Como esperado, as contribuições Newtonianas (densidade, RSD, quadrática) são tipicamente uma ou mais ordens de magnitude maiores que os sinais relativísticos em todo o espectro para ambos os redshifts.
2. Hierarquia Relativística em $z \approx 2$:
   • Efeitos de projeção e de GR dinâmica têm amplitudes comparáveis em grandes escalas.
   • Achado Inesperado: No limite comprimido (squeezed limit), os efeitos de radiação dominam o sinal relativístico, excedendo outros termos relativísticos por uma ordem de magnitude.
3. Hierarquia Relativística em $z \approx 0.6$: A hierarquia esperada é recuperada, onde os efeitos de projeção dominam. No entanto, as correções dinâmicas (da evolução não linear) permanecem não desprezíveis, sendo apenas um fator de 2 a 3 menores que os efeitos de projeção.
4. Comparação com Simulações:
   • Bispectro Total: Os resultados teóricos concordam razoavelmente bem com as medições de simulação para o bispectro total, exceto em pequenas escalas onde a resolução da simulação limita o poder.
   • Sinal Relativístico: Ao comparar o sinal relativístico "puro" (soma de termos de GR e radiação) com as medições de simulação (derivadas pela subtração de Newtoniano de simulações Relativísticas), os autores encontram uma discrepância. A soma teórica é menor por um fator de $\sim 5$ do que a medição da simulação em alto redshift.
   • Interpretação da Discrepância: Os autores atribuem essa lacuna a artefatos numéricos nas simulações, especificamente efeitos de gauge residuais (descompasso entre gauges de Poisson e N-body) e efeitos de backscaling (como a radiação é incluída em execuções Newtonianas). Eles concluem que esses erros numéricos são comparáveis em amplitude aos efeitos físicos de GR dinâmica.
   • Efeitos de Suavização: A análise destaca que a suavização (necessária para detectar sinais relativísticos minúsculos) pode introduzir vieses sistemáticos ao vazar erros numéricos de pequena escala para grandes escalas, potencialmente explicando desvios observados em comparações de baixo redshift.

Significância e Alegações
O artigo afirma estabelecer um pipeline teórico robusto para analisar o bispectro angular das contagens de número de matéria em levantamentos fotométricos, preenchendo a lacuna entre observáveis de lente gravitacional e agrupamento de galáxias (galaxy clustering).

• Validação de Amplitudes: O trabalho confirma as estimativas de amplitude dos efeitos relativísticos derivadas de simulações por Montandon et al. (2023), fornecendo simultaneamente uma decomposição termo a termo que as simulações não conseguem isolar facilmente.
• Relevância de Efeitos Dinâmicos: Contrariando a suposição de que os efeitos de projeção são a única contribuição relativística significativa, os autores demonstram que os efeitos de GR dinâmica e radiação são da mesma ordem de magnitude que os efeitos de projeção em alto redshift e permanecem significativos (dentro de um fator de 2) em redshifts mais baixos.
• Cautela sobre Extração de Simulações: Um achado significativo é a identificação de artefatos numéricos (gauge e backscaling) nas simulações relativísticas atuais que são tão grandes quanto os próprios sinais físicos buscados. Isso sugere que extrair sinais relativísticos "puros" dos métodos atuais de subtração de simulação exige extrema cautela.
• Utilidade Futura: Os autores postulam que seu framework é essencial para futuros levantamentos de próxima geração (como o Euclid) para distinguir a não-gaussianidade primordial da não-gaussianidade intrínseca e para testar teorias da gravidade no setor não linear, desde que simplificações como viés de galáxia e integrais ao longo da linha de visada sejam eventualmente incorporadas.