One with HI: Modelling HI Intensity Mapping one-point statistics including systematics

Resumo Técnico: One with HI: Modellando a Estatística de Um Ponto do Mapeamento de Intensidade de HI Incluindo Sistemáticas

Definição do Problema
O mapeamento de intensidade (IM) de hidrogênio neutro (HI) oferece uma janela única para o Universo pós-reionização ($z \leq 6$) e será um dos principais motores científicos para levantamentos futuros, como o Square Kilometre Array Observatory (SKAO). Embora o espectro de potência ($P(k)$) seja a ferramenta padrão para extrair informações cosmológicas, ele captura apenas a informação Gaussiana. Além das estatísticas de dois pontos, como a função de densidade de probabilidade (PDF) de um ponto, prometem extrair informações não-Gaussianas adicionais. No entanto, modelar essas estatísticas com precisão é desafiador devido às sistemáticas observacionais inerentes aos levantamentos de rádio, especificamente os efeitos do feixe do telescópio, a remoção de foregrounds e o ruído térmico. Essas sistemáticas amortecem modos específicos e suavizam o sinal, podendo obscurecer as assinaturas não-Gaussianas necessárias para quebrar degenerescências entre parâmetros cosmológicos (ex: $\Omega_M$ e $\sigma_8$) e parâmetros de viés ($b_1$).

Metodologia
Os autores desenvolvem um arcabouço teórico para modelar a PDF de HI em células esféricas, incorporando sistemáticas observacionais e validando-a contra simulações de alta resolução.

1. Dados de Simulação: O estudo utiliza a suíte de simulações UNIT (N-corpos de matéria escura apenas) com um volume de $1 \, h^{-3}\text{Gpc}^3$. As propriedades das galáxias e as massas de gás frio são atribuídas aos halos usando o modelo de Evolução de Galáxias Semi-Analítico (SAGE). As massas de HI são derivadas das massas de gás frio e mapeadas para mapas de intensidade usando um esquema de Vizinho Mais Próximo (NGP) em voxels cúbicos ($1 \, \text{Mpc}/h$).
2. Modelagem de Sistemáticas:
   • Remoção de Foreground: Modelada como uma supressão exponencial de grandes modos radiais ($k_\parallel$) para emular a remoção da emissão de sincrotrão galáctico suave.
   • Feixe do Telescópio: Modelado como uma convolução Gaussiana no plano transversal, suprimindo pequenas escalas angulares.
   • Ruído Térmico: Adicionado como um campo aleatório Gaussiano com variância derivada das especificações do instrumento (SKA1-MID).
   • Esses efeitos são aplicados no espaço de Fourier ao campo de intensidade de HI.
3. Arcabouço Teórico:
   • PDF de Matéria: Modelada usando Teoria de Grandes Desvios (LDT) e a aproximação de Colapso Esférico (SC) para prever a PDF de densidade de matéria não-linear.
   • Viés de HI e Estocasticidade: A distribuição condicional de HI dado a densidade de matéria escura, $P(\delta_{HI}|\delta_m)$, é modelada usando uma expansão de viés Euleriano quadrático e uma parametrização de estocasticidade não-Poissoniana. Os autores utilizam o código CosMomentum para computar a PDF do traçador por meio da convolução da PDF de matéria com esta distribuição condicional.
   • Espectro de Potência: Modelado usando hmcode-2020 para o espectro de potência de matéria não-linear, combinado com um modelo de viés linear, ruído de shot e os kernels de sistemáticas (feixe e foreground) aplicados como uma função de modulação $S_{eff}(k)$.
4. Validação e Previsão: Os modelos teóricos são validados contra as estatísticas baseadas em simulações usando uma técnica de reamostragem jackknife para estimar covariâncias. Uma análise de matriz de informação de Fisher é realizada para prever as restrições sobre parâmetros cosmológicos ($\Omega_M, \sigma_8$) e o viés linear ($b_1$) para uma configuração de levantamento SKA1-MID.

Principais Contribuições

• Modelo de PDF Integrado com Sistemáticas: O artigo apresenta um modelo teórico para a PDF de um ponto de HI que incorpora explicitamente os efeitos dos feixes de telescópios, remoção de foregrounds e ruído térmico. Isso é alcançado modificando a função de taxa no arcabouço LDT para incluir kernels de sistemáticas e ajustando os parâmetros de viés/estocasticidade condicional.
• Caracterização de Viés e Estocasticidade: Os autores quantificam o viés dependente de escala e a estocasticidade dos traçadores de HI na presença de sistemáticas, demonstrando que um modelo de viés quadrático e uma expansão de segunda ordem para estocasticidade fornecem um bom ajuste aos dados da simulação, embora termos de ordem superior possam ser necessários em escalas menores.
• Quebra de Degenerescência: O trabalho demonstra como combinar a PDF de HI com o espectro de potência pode quebrar a degenerescência entre o viés linear ($b_1$) e a amplitude das flutuações ($\sigma_8$), uma limitação inerente às estatísticas de dois pontos isoladas.

Resultos

• Precisão do Modelo: As previsões teóricas da PDF concordam com as medições baseadas em simulações com precisão de nível percentual (tipicamente $<4\%$ de desvio na menor escala de suavização $R=15.5 \, \text{Mpc}/h$) ao usar momentos condicionais simulados. Ao usar modelos de viés paramétricos, as discrepâncias são ligeiramente maiores, mas permanecem dentro de limites aceitáveis para previsões de Fisher.
• Impacto das Sistemáticas: O feixe do telescópio é identificado como a sistemática dominante, reduzindo significativamente a variância da PDF. A remoção de foreground modifica levemente a forma, enquanto o ruído térmico é amplamente suprimido pelo suavizamento esférico.
• Restrições Cosmológicas:
  • Plano $\Omega_M - \sigma_8$: Combinar o espectro de potência com a PDF em uma grande escala de suavização ($R=31.0 \, \text{Mpc}/h$) reduz significativamente a degenerescência entre $\Omega_M$ e $\sigma_8$, estreitando as incertezas em aproximadamente 64% para $\Omega_M$ e 49% para $\sigma_8$ comparado ao espectro de potência sozinho. Isso é atribuído ao desalinhamento das direções de degenerescência entre as duas estatísticas.
  • Degenerescência $\sigma_8 - b_1$: Para quebrar a degenerescência entre $\sigma_8$ e o viés $b_1$, a PDF em uma escala de suavização menor ($R=15.5 \, \text{Mpc}/h$) prova ser mais eficaz. Isso ocorre porque a redução da assimetria (skewness) da PDF é mais sensível ao parâmetro de viés em escalas menores e mais não-Gaussianas.
• Viés de Parâmetros: Os autores descobrem que o viés nos parâmetros recuperados devido às diferenças entre o modelo teórico e as simulações é geralmente inferior a $1\sigma$ para o volume simulado. No entanto, eles observam que, para os volumes muito maiores esperados do SKA, esses vieses podem se tornar significativos se não forem corrigidos (ex: via correções de Edgeworth para a assimetria).

Significância
O artigo estabelece que a PDF de um ponto de intensidade de HI, mesmo quando degradada por sistemáticas observacionais realistas, retém informação não-Gaussiana suficiente para complementar o espectro de potência. Ao fornecer um modelo teórico validado que leva em conta essas sistemáticas, o trabalho permite a extração de restrições cosmológicas mais rigorosas e a quebra de degenerescências críticas de parâmetros (especificamente $b_1$ vs. $\sigma_8$) em futuros levantamentos do SKAO. Os autores enfatizam que, embora seu modelo atual seja robusto para os volumes simulados, trabalhos futuros exigirão o refinamento da modelagem do redimensionamento da assimetria e volumes de simulação maiores para explorar plenamente as capacidades dos dados do SKA.