Forecast for the detectability of patchy hydrogen reionization in WEAVE-QSO measurements of the Lyman-$\alpha$ forest power spectrum at redshift $z \geq 4$

Resumo Técnico: Previsão para a Detectabilidade da Reionização de Hidrogênio Fragmentada em Medições de WEAVE-QSO

Definição do Problema
O meio intergaláctico (IGM) retém "relíquias térmicas" da época da reionização do hidrogênio. Durante a reionização, o gás fotoaquecido em bolhas ionizadas resfriou-se adiabaticamente ao longo de longas escalas de tempo, imprimindo inhomogeneidades espaciais na distribuição de temperatura do IGM. Essas inhomogeneidades persistem até $z \sim 4$ e são previstas para aumentar o espectro de potência de fluxo do Lyman-$\alpha$ ($P_{Ly\alpha}$) em grandes escalas (baixo número de onda, $k \lesssim 10^{-3} \, \mathrm{s\,km^{-1}}$). Embora estudos anteriores utilizando dados de alta resolução (ex: KODIAK, SQUAD, XQ-100) tenham encontrado evidências tentativas ($\simeq 2.7\sigma$) para esse aumento, a precisão estatística limitada e o comprimento de caminho de redshift impediram uma detecção de alta significância. Os autores visam prever se o próximo levantamento WEAVE-QSO, com seu massivo tamanho de amostra, poderá detectar e caracterizar definitivamente este rastro de relíquia em grande escala.

Metodologia
Os autores desenvolveram um pipeline abrangente para gerar espectros simulados da floresta $Ly\alpha$ adaptados à configuração do levantamento WEAVE-QSO e para prever as incertezas associadas.

1. Simulações Hidrodinâmicas: A estrutura subjacente do IGM e sua história térmica foram modeladas usando a suíte de simulações Sherwood-Relics. Estas simulações hidrodinâmicas cosmológicas resolvem especificamente o IGM durante e após a reionização, incorporando efeitos de reionização fragmentada que geram o aumento previsto de potência em grande escala.
2. Geração de Espectros Simulados (WEAVEify): Um pacote de software customizado, WEAVEify, foi desenvolvido para modelar o output do Sherwood-Relics para frente em espectros realistas de WEAVE-QSO. Este processo inclui:
   • Mapeamento de skews (linhas de visão) de profundidade óptica para os referenciais de repouso de QSO.
   • Convolução com um kernel Gaussiano dependente do comprimento de onda para simular a função de espalhamento de linha instrumental (assumindo $R \sim 5000$).
   • Posicionamento de contínuo usando um método de Análise de Componentes Principais (PCA) derivado de SDSS QSOs.
   • Adição de ruído realista limitado por fótons baseado nas características do detector, brilho do céu e tempos de exposição ($3060$ s), contabilizando estatística de Poisson, corrente escura e ruído de leitura.
3. Amostragem do Levantamento: Os autores geraram $10.000$ espectros de QSO simulados para quatro bins de redshift centrados em $z = 4.0, 4.2, 4.4,$ e $4.6$. A amostra foi construída usando um método de Monte Carlo baseado na função de luminosidade de QSO, aplicando um corte de magnitude de banda $r$ de $m_r < 21.9$. Isso resulta em tamanhos de amostra esperados de $2.942$, $1.822$, $1.385$, e$1.155$ QSOs por bin para uma área de levantamento de $10.000 \, \mathrm{deg^2}$.
4. Estimativa do Espectro de Potência e Modelagem de Incerteza: O espectro de potência 1D foi estimado usando um estimador de Transformada Rápida de Fourier (FFT). Os autores construíram uma matriz de covariância completa contabilizando:
   • Incertezas Estatísticas: Derivadas via bootstrapping do tamanho da amostra simulada.
   • Incertezas Sistemáticas: Modeladas separadamente para correções de resolução espectral, subtração de ruído, contaminação não correlacionada de linhas metálicas (C IV, Si IV, Mg II), incompletude de absorvedores de alta densidade de coluna (HCD) e erros de posicionamento de contínuo.
   • A análise foi restrita a números de onda $k < 0.05 \, \mathrm{s\,km^{-1}}$, onde o fator de correção de resolução permanece acima de $0.2$.
5. Inferência de Parâmetros: Um framework de Cadeia de Markov Monte Carlo (MCMC) Bayesiano foi empregado para comparar os dados simulados contra modelos teóricos. O modelo incluiu parâmetros padrão do IGM (profundidade óptica efetiva $\tau_{\mathrm{eff}}$, temperatura $T_0$, inclinação $\gamma$, energia térmica $u_0$) e um parâmetro de amplitude adicional, $A_p$, representando a correção de reionização fragmentada. Dois cenários foram testados: um ajuste "fiducial" permitindo que $A_p$ varie, e um ajuste "sem $A_p$" onde a correção de fragmentação é desativada.

Principais Contribuições

• Primeira Previsão Detalhada para WEAVE-QSO: Este trabalho fornece a primeira previsão específica para a precisão das medições de $P_{Ly\alpha}$ pelo levantamento WEAVE-QSO em $z \ge 4$.
• Modelagem Abrangente de Sistemas: O artigo introduz um tratamento rigoroso de incertezas sistemáticas específicas para levantamentos de baixa resolução e grande amostra, incluindo o impacto das asas de amortecimento de HCD e o erro de posicionamento de contínuo no poder de grande escala.
• Validação do Potencial de Detecção: Ao integrar a física do Sherwood-Relics com modelagem instrumental realista, o estudo demonstra que a assinatura específica da reionização fragmentada é distinta das variações térmicas padrão do IGM.

Resultos

• Orçamento de Incerteza: A incerteza relativa total prevista para $P_{Ly\alpha}$ é de aproximadamente $5,2\%$ quando média sobre todas as escalas e bins de redshift abaixo de $k_{\mathrm{max}}$. Em grandes escalas ($k \lesssim 0.01 \, \mathrm{s\,km^{-1}}$), as incertezas permanecem abaixo de $10\%$, dominadas por erros estatísticos no intervalo intermediário e erros sistemáticos (especificamente HCDs e posicionamento de contínuo) nas escalas mais amplas.
• Significância de Detecção: Quando a análise MCMC é aplicada aos dados simulados contendo o sinal de reionização fragmentada:
  • O modelo fiducial (com $A_p$) fornece um ajuste significativamente melhor do que o modelo sem reionização fragmentada ($\chi^2/\mathrm{dof}$ melhora de $65.28/32$ para $36.07/28$).
  • A melhoria corresponde a $\Delta\chi^2 = 29.21$ para quatro parâmetros adicionais, resultando em uma significância equivalente a uma Gaussiana de dois lados de $\simeq 4.5\sigma$.
  • As distribuições posteriores para $A_p$ recuperam com sucesso as amplitudes de entrada em todos os bins de redshift.
• Sensibilidade à Escala: A significância de detecção é fortemente dependente das maiores escalas. Quando o bin de menor $k$ ($k \sim 10^{-3} \, \mathrm{s\,km^{-1}}$) é removido, a significância cai para $\simeq 1.2\sigma$, indicando que a detecção depende dos modos mais afetados pelo rastro térmico da fragmentação.

Significância e Alegações
O artigo alega que o levantamento WEAVE-QSO, com seu tamanho de amostra previsto e propriedades de covariância, deve ser capaz de detectar diretamente o rastro de relíquia de grande escala da reionização de hidrogênio fragmentada no espectro de potência da floresta $Ly\alpha$ em $z \ge 4$ com alta significância estatística ($\simeq 4.5\sigma$). Isso representa uma melhoria substancial em relação às detecções tentativas anteriores ($\simeq 2.7\sigma$) de conjuntos de dados menores e de alta resolução. Os autores concluem que o framework desenvolvido permite a avaliação do tempo de reionização através de medições do espectro de potência de grande escala e que a configuração do WEAVE-QSO é suficiente para transformar evidências tentativas em uma detecção robusta, desde que as incertezas sistemáticas sejam controladas nos níveis previstos. O estudo não pretende medir o redshift exato da reionização, mas sim a amplitude do rastro de fragmentação em relação a um modelo de fundo UV uniforme.