CoLoRe-2LPT: Lyman-$\alpha$ mock catalogues for the validation of DESI cosmological analyses

Declaração do Problema
A floresta de Lyman-$\alpha$ (Ly$\alpha$) é uma sonda crítica para o estudo da estrutura em grande escala do universo em altos redshifts ($z > 2$), oferecendo restrições sobre as Oscilações Acústicas de Bárions (BAO) e o agrupamento de forma completa (FS). No entanto, à medida que o Dark Energy Spectroscopic Instrument (DESI) entra em uma era de precisão sem precedentes com seu Data Release 2 (DR2), a sensibilidade às suposições de modelagem, aos sistemas de erro astrofísicos e aos efeitos observacionais aumentou. Catálogos de mocks realistas são essenciais para validar pipelines de análise, quantificar vieses sistemáticos e verificar a recuperação de parâmetros. Métodos anteriores de geração de mocks para a floresta de Ly$\alpha$ baseavam-se em realizações log-normais do campo de densidade. Embora eficazes para estudos de BAO em grande escala, essas abordagens são não confiáveis para escalas levemente não lineares ($k \sim 0,1 - 0,5 \, h/\text{Mpc}$) e falham em reproduzir com precisão o alargamento não linear do pico de BAO ou o agrupamento em pequena escala necessário para análises robustas de FS. Além disso, os mocks anteriores frequentemente careciam de uma modelagem realista de quasares (QSOs) e contaminantes astrofísicos como sistemas de Alta Coluna de Densidade (HCDs) e linhas metálicas.

Metodologia
Os autores apresentam o "CoLoRe-2LPT", um novo conjunto de catálogos de mocks de Ly$\alpha$ de custo computacional rápido, projetados para abordar as limitações das abordagens log-normais. A metodologia envolve um pipeline de múltiplos estágios:

1. Campos de Densidade e Velocidade 2LPT: A inovação central é o uso da Teoria de Perturbação Lagrangiana de Segunda Ordem (2LPT) para evoluir o campo de densidade de matéria, substituindo a transformação log-normal. O pipeline utiliza o código CoLoRe para gerar cones de luz cosmológicos. Ele começa com um campo aleatório Gaussiano em $z=0$, evolui partículas de teste usando 2LPT para obter deslocamentos e velocidades, e as reescala para gerar cones de luz cobrindo a faixa de redshift de QSOs do DESI ($0 < z < 3,8$).
2. Calibração de Agrupamento de QSO: Para povoar os cones de luz com QSOs, os autores empregam um modelo de viés de limiar dependente da densidade. Os parâmetros deste modelo são calibrados empiricamente contra simulações de N-corpos Abacus de alta resolução. Isso garante que o viés linear do QSO e as distorções no espaço de redshift (RSD) correspondam aos dados observacionais. O modelo inclui a adição de efeitos de "Dedos de Deus" (Fingers-of-God, FoG) via velocidades aleatórias para reproduzir o amortecimento visto em simulações de N-corpos.
3. Geração de Transmissão de Ly$\alpha$: O módulo LyaCoLoRe transforma os espetros (skewers) de densidade e velocidade 2LPT em espetros de transmissão de Ly$\alpha$. Este processo envolve:
   • Adicionar potência em pequena escala aos campos de densidade 2LPT usando um espectro de potência 1D para corresponder ao agrupamento característico da floresta até $\sim 100 \, \text{kpc}/h$.
   • Converter densidade em profundidade óptica usando a Aproximação de Gunn-Peterson Fluctuante (FGPA).
   • Aplicar RSDs para mapear a profundidade óptica no espaço real para o espaço de redshift.
   • Computar o fluxo transmitido.
   • Os parâmetros para potência em pequena escala, profundidade óptica e escala de velocidade são ajustados para reproduzir as medições do DESI DR1 do espectro de potência 1D ($P_{1D}$), fluxo médio, viés de Ly$\alpha$ ($b_{\text{Ly}\alpha}$) e o parâmetro de RSD ($\beta_{\text{Ly}\alpha}$).
4. Contaminantes Astrofísicos e Ruído: Para criar mocks semelhantes ao DESI, o pipeline incorpora:
   • HCDs: Modelados ao povoar células acima de um limiar de densidade, com parâmetros ajustados para corresponder ao agrupamento e taxas de incidência de HCD.
   • Linhas Metálicas: Adicionadas após as RSDs usando o script quickquasars, modelando transições (ex: Si II, Si III) que criam picos de correlação espúrios.
   • Contínuo e Ruído: Os contínuos de QSO são gerados usando SIMQSO, e o ruído instrumental é adicionado via specsim para corresponder à pegada do DESI DR2 e características de exposição.

Principais Contribuições

• Primeira Aplicação de 2LPT a Mocks de Cone de Luz Completo de Ly$\alpha$: Embora métodos LPT tenham sido usados para levantamentos de galáxias e mocks de Ly$\alpha$ baseados em caixas, este trabalho constitui a primeira aplicação de 2LPT a mocks de cone de luz completo para um levantamento de floresta de Ly$\alpha$.
• Modelagem Não Linear Melhorada: A abordagem 2LPT captura naturalmente não linearidades leves e o alargamento não linear do pico de BAO, que anteriormente eram modelados de forma ad-hoc em mocks log-normais.
• Validação Abrangente: Os mocks são validados contra múltiplas estatísticas do DESI, incluindo o espectro de potência 1D, fluxo médio, vieses de Ly$\alpha$ e QSO, e o parâmetro $\beta$.
• Modelagem Realista de Contaminantes: A inclusão de metais e HCDs permite a exploração de sistemáticas induzidas por esses contaminantes astrofísicos, o que é crucial para análises de FS.

Resultados

• Estatísticas de Agrupamento: Os novos mocks reproduzem as estatísticas de Ly$\alpha$ dentro de 10% das medições mais recentes do DESI. Especificamente, o viés de Ly$\alpha$ e o parâmetro $\beta$, o fluxo transmitido médio e o espectro de potência 1D são bem reproduzidos.
• Agrupamento de QSO: O viés linear de QSO é recuperado dentro de 5% dos valores observacionais. A inclusão de velocidades 2LPT e efeitos de FoG melhora significativamente a recuperação de distorções no espaço de redshift em comparação com abordagens anteriores.
• Validação de BAO e FS: Quando ajustados com o modelo de BAO completo do DESI DR2, as correlações dos mocks mostram excelente concordância com os dados observacionais. Os parâmetros cosmológicos derivados dos mocks desviam-se da unidade em no máximo 0,003, um desvio negligenciável comparado às incertezas estatísticas. Os mocks reproduzem com sucesso o amortecimento anisotrópico do pico de BAO naturalmente, sem alargamento ad-hoc.
• Validação de Contaminantes: A inclusão de metais e HCDs resulta em correlações de mocks que coincidem qualitativa e quantitativamente com os picos espúrios e distorções observados nos dados do DESI DR2. Os vieses de metais concordam com os parâmetros do DESI DR2 dentro de $\sim 20\%$.

Significância
O artigo afirma que este novo conjunto de mocks é uma ferramenta "chave" para aumentar a modelagem e validação da análise cosmológica de forma completa (full-shape) de Ly$\alpha$ do DESI DR2. Ao fornecer uma ferramenta computacionalmente eficiente e motivada fisicamente, que captura com precisão tanto o agrupamento em grande escala quanto o levemente não linear, os mocks CoLoRe-2LPT permitem a validação robusta de análises de FS e BAO para levantamentos de Ly$\alpha$ atuais e de próxima geração. O trabalho demonstra que mocks baseados em 2LPT podem substituir abordagens log-normais para fornecer características de pequena escala e alargamento de BAO mais realistas, o que é essencial conforme as análises cosmológicas avançam além de simples medições de BAO para restrições de forma completa (full-shape).