Control variates from Eulerian and Lagrangian perturbation theory: Application to the bispectrum

Resumo Técnico: Variáveis de Controle a partir da Teoria de Perturbação Euleriana e Lagrangiana: Aplicação ao Bispectro

Declaração do Problema
Simulações cosmológicas de N-corpos são essenciais para modelar a estrutura não gaussiana da teia cósmica, particularmente para estatísticas de ordem superior como o bispectro. No entanto, essas simulações sofrem de "variância cósmica" devido à aleatoriedade de suas condições iniciais. Para alcançar a alta precisão exigida por levantamentos de galáxias de Fase IV, tipicamente é necessário fazer uma média sobre muitas realizações de simulação, o que é computacionalmente proibitivo. Embora existam técnicas de redução de variância como "emparelhamento e fixação", elas são menos eficazes para bispectros do que para espectros de potência. O método de variáveis de controle oferece uma alternativa promissora, explorando simulações mais baratas e correlacionadas (substitutas) para reduzir a variância amostral de observáveis. Um desafio crítico é identificar modelos substitutos que sejam tanto tratáveis analiticamente (para determinar sua média sem ensembles caros) quanto altamente correlacionados com o bispectro completo de N-corpos em uma ampla gama de escalas.

Metodologia
Os autores investigam variáveis de controle perturbativas derivadas tanto da Teoria de Perturbação Euleriana (EPT) quanto da Teoria de Perturbação Lagrangiana (LPT), utilizando o bispectro da matéria em $z=0.5$ como estudo de caso. O estudo utiliza a suíte Quijote de simulações de N-corpos ($N=1000$ realizações) como verdade fundamental.

1. Análise de Modelo de Brinquedo: Os autores analisam primeiro um modelo simplificado onde o campo "verdadeiro" é gerado pela Aproximação de Zeldovich (ZA). Eles demonstram que variáveis de controle "lineares" (gaussianas), embora tenham média zero, exibem um desencontro exponencial com o campo não linear proporcional a $\exp[-\frac{1}{2}\Sigma^2(k_1^2 + k_2^2 + k_3^2)]$, onde $\Sigma$ é a dispersão de deslocamento.
2. Variáveis de Controle Eulerianas: Os autores constroem variáveis de controle usando campos EPT até a quinta ordem na densidade ($O(\delta^5)$). Eles empregam um parâmetro de contagem $\epsilon$ para isolar contribuições de nível de árvore ($O(\epsilon)$) e diagramas de ordem superior de média zero ($O(\epsilon^2)$) da realização completa da rede. Eles abordam o desafio de ajustar previsões analíticas a medições de rede desenvolvendo esquemas de "triângulo efetivo" e média discreta para corrigir efeitos de agrupamento (binning).
3. Variáveis de Controle Lagrangianas: Os autores avaliam os bispectros da ZA e da Teoria de Perturbação Lagrangiana de Segunda Ordem (2LPT). Eles mostram que, embora a 2LPT capture a estrutura correta de nível de árvore e os fluxos de massa em grande escala, seu bispectro completo não é calculável analiticamente, exigindo estimativa empírica da média. Por outro lado, o bispectro da ZA é calculável analiticamente, mas possui uma estrutura de nível de árvore incorreta em comparação com a EPT, levando a uma baixa correlação em baixos-$k$.
4. Variáveis de Controle Deslocadas: Para combinar as vantagens de ambas as abordagens, os autores introduzem uma "variável de controle deslocada". Esta é construída usando a base de "operador deslocado", onde o campo é definido como um traçador enviesado lagrangiano com apenas o deslocamento linear (Zeldovich) exponenciado. Esta construção garante que o campo retenha o bispectro de nível de árvore correto (correspondendo à EPT) enquanto resume deslocamentos lineares (mantendo as propriedades de correlação da LPT).
5. Cálculo Analítico: Para o bispectro da ZA, os autores desenvolvem um método para calcular previsões analíticas médias por bin usando pesos derivados de cálculos de nível de árvore e amostragem Monte-Carlo para bins largos, alcançando concordância sub-1% com medições de rede.

Principais Contribuições e Resultados

• Limites de Desencontro: O estudo confirma que variáveis de controle baseadas em EPT sofrem de desencontro exponencial em escalas menores que a dispersão de deslocamento ($k \gtrsim \Sigma^{-1}$), limitando sua utilidade a grandes escalas.
• ZA vs. 2LPT: A aproximação de Zeldovich, embora tratável analiticamente, é uma variável de controle subótima para o bispectro em comparação com o espectro de potência, porque seu bispectro de nível de árvore difere da previsão da EPT. A 2LPT oferece melhor correlação, mas carece de tratabilidade analítica para a média.
• Desempenho da Variável de Controle Deslocada: A "variável de controle deslocada" proposta supera todos os outros métodos perturbativos. Ela combina a estrutura correta de nível de árvore da EPT com a resomação de deslocamentos lineares da LPT.
  • Redução de Variância: Aplicar a variável de controle deslocada a uma única simulação de N-corpos com $V = 1 (\text{Gpc}/h)^3$ reduz a variância do bispectro da matéria por um fator equivalente a fazer uma média sobre $10^4$ simulações para os triângulos de menor-$k$.
  • Precisão: Isso permite medições de precisão sub-2% do bispectro da matéria para todas as configurações de triângulo na faixa de $k_{\text{min}} = 0.04$ a $k_{\text{max}} = 0.47 \, h\text{Mpc}^{-1}$ usando uma única simulação.
• Tratabilidade Analítica: Os autores demonstram que a variável de controle deslocada é, em princípio, calculável analiticamente em todas as ordens, embora a implementação numérica eficiente para o bispectro completo permaneça uma tarefa para trabalhos futuros. Atualmente, eles estimam seu desempenho numericamente.

Significância e Afirmações
O artigo afirma que o desenvolvimento de emuladores precisos de bispectro é agora viável sem o custo proibitivo de fazer uma média de milhares de simulações de N-corpos. Ao utilizar a variável de controle deslocada, a variância inerente à modelagem baseada em simulação pode ser suprimida a um nível onde a precisão é limitada apenas pela medição de configurações de triângulo em pequena escala, e não pela variância amostral.

Os autores posicionam este trabalho como estabelecendo a base teórica para estender variáveis de controle a toda a hierarquia $n$-pontos. Eles afirmam explicitamente que a metodologia é suficiente para se estender a traçadores enviesados (galáxias) e espaço de redshift, embora essas aplicações sejam deixadas para trabalhos futuros. O artigo não afirma ter resolvido o cálculo analítico do bispectro deslocado completo, mas estabelece sua viabilidade teórica e desempenho superior em testes numéricos. A significância primária reside em permitir modelagem cosmológica baseada em simulação de alta precisão, adequada para futuros levantamentos de estrutura em grande escala.