ECHO21: a tool for modelling global 21-cm signal from dark ages to reionization

Resumo Técnico de "ECHO21: Uma ferramenta para modelar o sinal global de 21 cm das eras escuras à reionização"

Problema e Motivação
As eras escuras e o amanhecer cósmico são períodos pouco compreendidos da linha do tempo cósmica devido à falta de sondas observacionais diretas. A linha de transição hiperfina de 21 cm do hidrogênio neutro oferece uma potencial sonda para desvendar a física das primeiras estrelas, do meio intergaláctico (IGM) e dos parâmetros cosmológicos. Embora códigos numéricos e semi-numéricos como 21cmFAST e BEORN forneçam previsões precisas e fisicamente ricas, eles são computacionalmente lentos, dificultando a inferência astrofísica ou cosmológica eficiente a partir de grandes conjuntos de dados. Por outro lado, códigos analíticos existentes como ARES e Zeus21 oferecem velocidade, mas carecem de certos componentes físicos, especificamente os efeitos de aquecimento de fótons Lyman-$\alpha$ (Ly$\alpha$). Os autores identificam que a omissão do aquecimento por Ly$\alpha$ pode levar a vieses fortes na modelagem do sinal. Para abordar a necessidade de uma ferramenta rápida, suficientemente realista e flexível, capaz de lidar com grandes espaços de parâmetros para inferência, os autores introduzem o ECHO21 (Exploring the Cosmos with Hydrogen Observation).

Metodologia
O ECHO21 é um pacote Python que utiliza um arcabouço analítico para modelar o sinal global de 21 cm de redshift $z \approx 1500$ até o presente dia. O código opera em um modelo de duas zonas do IGM, distinguindo entre o IGM em massa e as regiões H II totalmente ionizadas.

• Arcabouço Físico: O código resolve equações diferenciais acopladas para a temperatura cinética do gás ($T_k$), a fração eletrônica do IGM em massa ($x_e$) e o fator de preenchimento de volume das regiões H II ($Q$).

  • Histórico Térmico: A evolução de $T_k$ considera o resfriamento adiabático, o aquecimento de Compton, o aquecimento por raios-X e, crucialmente, o aquecimento por Ly$\alpha$. Os autores observam que, enquanto códigos analíticos anteriores ignoravam o aquecimento por Ly$\alpha$, o ECHO21 o inclui, calculando a taxa de aquecimento com base na intensidade específica de Ly$\alpha$ ($J_{Ly}$) derivada da densidade de taxa de formação estelar (SFRD).
  • Histórico de Ionização: O código rastreia a recombinação do IGM em massa e a ionização impulsionada por raios-X e fótons UV. Ele modela o fator de preenchimento de volume das regiões H II comparando a taxa de emissão de fótons ionizantes contra as recombinações radiativas.
  • Densidade de Taxa de Formação Estelar (SFRD): Uma característica inovadora do ECHO21 é a inclusão de três modelos distintos de SFRD:
    1. Fisicamente motivado: Baseado no colapso de bárions em halos de matéria escura (usando funções de massa de halo Press-Schechter, Sheth-Tormen ou Tinker) com uma eficiência de formação estelar constante.
    2. Semi-empírico: Uma variante do modelo físico incorporando uma escala de tempo de formação estelar característica e um corte abrupto em uma temperatura virial mínima.
    3. Empiricamente motivado: Um modelo restrito pelos dados da função de luminosidade UV (forma Madau-Dickinson para baixo $z$, queda exponencial para alto $z$).

• Parâmetros: O modelo é controlado por seis parâmetros astrofísicos livres: normalização de Ly$\alpha$ ($f_{Ly}$), índice espectral da série de Lyman ($s$), normalização de raios-X ($f_X$), índice espectral de raios-X ($w$), fração de escape de fótons ionizantes ($f_{esc}$) e temperatura virial mínima ($min(T_{vir})$). Adicionalmente, o código permite a variação dos parâmetros cosmológicos padrão ($H_0, \Omega_m, \Omega_b, \sigma_8, n_s, T_{\gamma0}, Y_p$).

• Desempenho: A abordagem analítica permite que o ECHO21 gere uma única instância do sinal de 21 cm global em aproximadamente $O(1)$ segundo em uma única CPU. O código foi projetado para paralelização MPI, permitindo a geração de milhares de modelos em minutos.

Principais Contribuições

1. Inclusão do Aquecimento por Ly$\alpha$: O ECHO21 é o primeiro código analítico a incluir explicitamente o aquecimento por Ly$\alpha$ na evolução térmica do IGM. Os autores demonstram que negligenciar este termo pode introduzir vieses significativos, particularmente em cenários de baixa emissividade de raios-X e forte emissão de Ly$\alpha$.
2. Diversidade de Modelos de SFRD: O código oferece unicamente três abordagens diferentes para modelar a SFRD (física, semi-empírica e empírica), permitindo que os usuários avaliem como a escolha de uma prescrição de formação estelar impacta o sinal de 21 cm resultante.
3. Flexibilidade Cosmológica: Ao contrário de muitos códigos especializados, o ECHO21 trata parâmetros cosmológicos com a mesma facilidade que parâmetros astrofísicos, facilitando o estudo das incertezas cosmológicas (ex: a tensão de Hubble) no sinal de 21 cm.
4. Velocidade e Escalabilidade: O código é otimizado para inferência rápida, capaz de gerar 15.625 modelos em um espaço de parâmetros astrofísicos de seis dimensões em ~17 minutos usando 50 CPUs.

Resultados

• Modelo Fiducial: Usando uma SFRD fisicamente motivada e a função de massa de halo Press-Schechter, o modelo fiducial prevê um histórico de reionização completo até $z \approx 5.3$ com um redshift de ionização de 50% de $z_{re} = 7.42$. A espessura óptica de espalhamento de elétrons é calculada como $\tau_e \approx 0.0592$, consistente com as restrições do Planck 2018.
• Características do Sinal: O sinal global de 21 cm fiducial exibe um primeiro pico de absorção de $-41.1$ mK em $z=87.5$ (eras escuras) e um segundo pico de absorção mais profundo de $-142.4$ mK em $z=19.2$ (amanhecer cósmico), seguido por um pico de emissão de $22.2$ mK em $z=11.2$.
• Impacto do Aquecimento por Ly$\alpha$: Comparações com modelos que excluem o aquecimento por Ly$\alpha$ mostram que a omissão leva a um sinal de absorção de 21 cm mais forte (mais profundo). Os autores enfatizam que o aquecimento por Ly$\alpha$ é um componente crucial para uma modelagem precisa, pois pode contrapor efeitos de resfriamento e alterar a força do sinal de forma não monotônica dependendo dos parâmetros de raios-X.
• Impacto dos Modelos de SFRD: A escolha do modelo de SFRD altera significativamente o sinal. Modelos empiricamente motivados, que preveem um início de formação estelar mais tardio em comparação com modelos fisicamente motivados, resultam em sinais que são atrasados em $\sim 8.8$ em redshift e mais rasos em $\sim 39.4$ mK. Da mesma forma, a escolha da Função de Massa de Halo (HMF) afeta o tempo e a profundidade das características do sinal.
• Sensibilidade de Parâmetros: O sinal é mais sensível a $f_{Ly}$, $f_X$, $w$ e $min(T_{vir})$. Variações nos parâmetros cosmológicos dentro das incertezas observacionais atuais resultam em variações relativamente pequenas na amplitude do sinal (ex: uma faixa de $\sim 18$ mK na profundidade do pico de absorção).

Significância e Alegações
Os autores posicionam o ECHO21 como uma ferramenta de código aberto flexível e extensível, projetada para realizar "inferências astrofísicas rápidas, porém suficientemente realistas" a partir de futuros conjuntos de dados de 21 cm (ex: de REACH, SARAS ou futuros experimentos espaciais). Eles alegam que o código é particularmente adequado para:

• Realizar inferência astrofísica e cosmológica robusta via exploração de grande escala de parâmetros.
• Testar pipelines de análise de dados (especificamente para a colaboração REACH) com sinais fisicamente motivados em vez de apenas ajustes empíricos.
• Construir sinergias entre observações de 21 cm e outras sondas cosmológicas, como previsões de BBN, anisotropias do CMB e a tensão de Hubble.

O artigo reconhece modestamente que, embora o ECHO21 seja menos preciso que modelos semi-numéricos devido às suas aproximações analíticas (ex: negligenciar múltiplos espalhamentos de Ly$\alpha$), sua velocidade e modularidade o tornam uma ferramenta ideal para inferências iniciais e para incorporar nova física em versões futuras. O código está publicamente disponível via PyPI e GitHub.