Modeling the Cosmological Lyman-$α$ Forest at the Field Level

Resumo Técnico: Modelagem da Floresta Lyman-$\alpha$ Cosmológica ao Nível de Campo

Definição do Problema
A floresta Lyman-$\alpha$ (Ly-$\alpha$), observada como linhas de absorção nos espectros de quasares distantes, serve como uma sonda única do meio intergaláctico (IGM) e da estrutura cosmológica em altos redshifts ($2 \lesssim z \lesssim 4$) e em pequenas escalas. Embora os levantamentos atuais e futuros, como o Dark Energy Spectroscopic Instrument (DESI), estejam gerando volumes de dados sem precedentes, as capacidades de modelagem teórica existentes são insuficientes. Historicamente, as análises foram restritas à função de dois pontos (espectro de potência) ou à sua contraparte de Fourier. No entanto, comprimir um campo tridimensional em uma função de dois pontos é inerentemente uma compressão com perda, descartando correlações de ordem superior e informações de fase. Além disso, embora simulações hidrodinâmicas numéricas possam modelar a floresta Ly-$\alpha$, elas são computacionalmente proibitivas para gerar os grandes volumes cosmológicos necessários para a inferência ao nível de levantamento. Há uma necessidade crítica de um modelo direto analítico e perturbativo que opere ao nível de campo para preencher a lacuna entre a teoria analítica de grande escala e as simulações de pequena escala, permitindo a inferência ao nível de campo e o uso de priors baseados em simulações.

Metodologia
Os autores desenvolvem um modelo direto para as flutuações do fluxo da floresta Ly-$\alpha$ ($\delta_F$) baseado na Teoria de Campo Efetiva (EFT) da Estrutura de Grande Escala. O núcleo da abordagem envolve modelar o campo de decremento de fluxo como uma expansão determinística do campo de densidade de matéria escura subjacente e seus gradientes de velocidade, suplementado por um termo de ruído estocástico.

1. Estrutura EFT: O modelo utiliza o princípio da equivalência para estabelecer que as flutuações de fluxo rastreiam a superdensidade de matéria escura ($\delta$) e o gradiente de velocidade na linha de visada ($\eta$). Os autores estendem isso além da teoria linear, incluindo operadores dependentes da linha de visada ditados pelas simetrias da floresta Ly-$\alpha$.
2. Expansão de Viés (Bias Expansion): O modelo é construído como uma expansão de viés até a ordem cúbica. Para evitar redundância, o conjunto de operadores de viés é ortogonalizado usando o algoritmo de Gram-Schmidt. O modelo resultante expressa o campo de fluxo $\delta_F$ como uma soma de funções de transferência ($\beta$) atuando sobre operadores ortogonalizados (incluindo densidade linear, densidade quadrática, velocidade e termos de maré).
3. Ressumação IR: Para contabilizar as dependências não perturbativas de deslocamentos lineares relevantes para as Oscilações Acústicas de Bárions (BAO), os operadores são deslocados por deslocamentos de Zel'dovich.
4. Estocasticidade: O modelo inclui um campo estocástico ($\epsilon$) para capturar o erro irredutível ou "ruído de choque" (shot noise) decorrente da natureza não determinística do fluxo de Ly-$\alpha$ em pequenas escalas. O espectro de potência deste ruído é modelado como uma expansão dependente da escala ($n_0 + \alpha_1 k^2 + \alpha_2 \mu^2 k^2$).
5. Validação: O modelo é testado contra simulações hidrodinâmicas de alta fidelidade da suíte Sherwood (volume de $160^3 \, (h^{-1}\text{Mpc})^3$, $2048^3$ partículas) no redshift $z=2.8$. Os autores comparam os campos modelados diretamente contra os "verdadeiros" campos simulados, avaliando a concordância via espectros de potência, correlações cruzadas com halos e funções de densidade de probabilidade (PDF) de um ponto.

Principais Contribuições

• Modelo Direto ao Nível de Campo: Este trabalho apresenta o primeiro modelo direto analítico para a floresta Ly-$\alpha$ que opera ao nível de campo, em vez de apenas em estatísticas de resumo.
• Descrição Sistemática de EFT: Fornece uma descrição sistemática e de primeiros princípios para o campo de Ly-$\alpha$ em grandes escalas usando EFT, incorporando os operadores necessários dependentes da linha de visada.
• Quantificação de Estocasticidade: O artigo fornece as primeiras estimativas da estocasticidade (espectro de potência de ruído) das flutuações 3D de Ly-$\alpha$, caracterizando sua amplitude e dependência de escala.
• Modelagem de Correlação Cruzada: O arcabouço é estendido para modelar consistentemente a correlação cruzada entre a floresta Ly-$\alpha$ e halos de matéria escura (proxies para quasares e galáxias de alto redshift).

Resultados

• Precisão: O modelo direto de EFT reproduz os resultados da simulação Sherwood com alta fidelidade. Em escalas maiores que alguns Mpc, o modelo coincide com a distribuição de fluxo, correlações cruzadas Ly-$\alpha$–halo e estatísticas de contagem em células com precisão de nível percentual.
• Desempenho do Espectro de Potência: O modelo completo de viés cúbico da EFT combina com o espectro de potência da simulação até $k \approx 0.6 \, h \text{ Mpc}^{-1}$ com resíduos abaixo de 5%. Em contraste, a teoria linear suplementada com uma função de transferência falha significativamente, produzindo erros superiores a 5% mesmo em grandes escalas.
• Concordância ao Nível de Campo: O modelo reproduz com sucesso as amplitudes e fases de todos os modos de Fourier, um teste mais rigoroso do que comparar estatísticas de resumo. Comparações visuais de fatias de decremento de fluxo mostram excelente concordância qualitativa, com resíduos concentrados principalmente em regiões subdensas próximas a halos.
• Caracterização do Ruído: O espectro de potência do ruído é encontrado como branco em grandes escalas, consistente com as previsões de EFT. Os autores estimam a amplitude do ruído de choque $n_0 \approx 0.18 \, [h^{-1}\text{Mpc}]^3$. A dependência de escala do ruído sugere uma conexão com halos de matéria escura, onde o gás frio se acumula ao redor de halos massivos, mas é expelido de halos mais leves.
• Limitações: A inclinação do espectro de potência do ruído torna-se mais rasa em torno de $k \approx 0.6 \, h \text{ Mpc}^{-1}$, indicando uma quebra da expansão de gradiente para o componente estocástico nessas escalas.

Significância e Alegações
Os autores afirmam que este trabalho fornece ferramentas cruciais para unir a modelagem analítica em grandes escalas com simulações em pequenas escalas. Ao permitir a inferência ao nível de campo, o modelo permite a análise de dados da floresta Ly-$\alpha$ de maneiras anteriormente impossíveis, como utilizar todo o conteúdo de informação do campo 3D, em vez de apenas o espectro de potência.

O artigo afirma que esta abordagem é oportuna para maximizar o potencial científico das medições do DESI. Especificamente, o modelo pode:

• Unir mapas de Ly-$\alpha$ em grandes e pequenas escalas ao gerar analiticamente modelos de grande escala que correspondam às propriedades de simulações de alto volume e alta resolução (ex: Illustris TNG, CAMELS).
• Remover o gargalo computacional das simulações hidrodinâmicas de grande volume, permitindo uma inferência baseada em simulações eficiente.
• Gerar simulações semianalíticas para estimar matrizes de covariância necessárias para análises conjuntas de dados de floresta Ly-$\alpha$ e quasares.
• Facilitar medições de precisão dos parâmetros de EFT de Ly-$\alpha$ ao nível de campo sem variância de amostra.

Em última análise, os autores declaram que esta abordagem avança a compreensão teórica da conexão da floresta Ly-$\alpha$ com a distribuição de matéria subjacente, aproximando o campo da resolução dos enigmas da formação de estruturas no Universo de alto redshift.