Exploring the High-Redshift 21-cm Signal via Self-Consistent Simulations using Artificial Neural Network Emulation

Resumo Técnico: Explorando o Sinal de 21 cm de Alto Desvio para o Vermelho via Simulações Autoconsistentes usando Emulação por Redes Neurais Artificiais

Declaração do Problema
Compreender o Amanhecer Cósmico (AC), a era das primeiras estrelas e galáxias ($z \sim 15-60$), é crucial para rastrear a evolução cósmica. Embora sondas observacionais como o JWST e arranjos de rádio (HERA, LOFAR, MWA) estejam avançando, a interpretação de dados da linha de 21 cm exige modelos teóricos robustos. As abordagens existentes enfrentam uma dicotomia: modelos analíticos simplificados são computacionalmente eficientes, mas dependem de pressupostos restritivos, enquanto simulações hidrodinâmicas complexas oferecem alta fidelidade, mas são limitadas a volumes pequenos e espaços de parâmetros estreitos devido ao custo computacional. Além disso, muitos modelos seminumericos carecem de tratamentos autoconsistentes da formação estelar (FE) em pequena escala calibrados à hidrodinâmica, histórias estocásticas de fusão de halos de matéria escura (ME) e feedback radiativo em grande escala. Este artigo aborda a necessidade de um framework seminumerico, em grande escala e autoconsistente, que preencha essa lacuna para fornecer previsões precisas do sinal global de 21 cm e do espectro de potência.

Metodologia
Os autores apresentam um framework de simulação seminumerico aprimorado, construindo sobre o trabalho de Feathers et al. (2025, F25). A metodologia central envolve:

1. Volume de Simulação: Um volume em grande escala de $192^3$ Mpc$^3$ (composto por $64^3$ células de $(3 \text{ Mpc})^3$) é gerado usando 21cmvFAST. Este volume inclui flutuações espaciais na densidade de matéria ($\delta$) e velocidades de fluxo bárion-ME ($v_{bc}$), que são críticas para previsões precisas de FE.
2. Evolução de Halos de Matéria Escura: O framework utiliza o formalismo Press-Schechter Estendido (EPS) para construir árvores estocásticas de fusão de halos de ME. Uma biblioteca de $\sim 10^4$ árvores de fusão é gerada para 40 faixas de massa de halo ($10^{5.6} - 10^{9.5} M_\odot$). Essas árvores são atribuídas às células de simulação com base na sobre-densidade local, servindo como andaime para a FE.
3. Emulação por Rede Neural (RN): Para superar as proibições computacionais de executar um Modelo Semi-Analítico de Subgrade (SAM) complexo para cada célula, os autores empregam emuladores de Redes Neurais Artificiais.
   • Dados de Treinamento: As RNs são treinadas em um SAM calibrado às simulações hidrodinâmicas mais avançadas: as simulações Renaissance (O'Shea et al. 2015) e AEOS (Brauer et al. 2025).
   • Prescrição de Formação Estelar (FE): O modelo fiducial incorpora prescrições de FE de Hazlett et al. (2025a,b). Isso inclui:
     • FE PopIII: Modelada com uma função de distribuição cumulativa para massas estelares, resultando em uma média de quatro estrelas de $35.6 M_\odot$ por halo ao cruzar o limiar de massa crítica ($M_{crit}$).
     • FE PopII: Modelada como um processo de duas fases: uma fase inicial "explosiva" seguida por uma fase "estável" assim que a temperatura virial do halo atinge $T_{vir, steady} \approx 1.76 \times 10^4$ K. Isso inclui um período de quiescência ($t_{quies}$) após as explosões.
     • Massa Crítica ($M_{crit}$): Calculada considerando o feedback de Lyman-Werner (LW) e o auto-escudo de H$_2$ (Kulkarni et al. 2021).
   • Eficiência: A emulação por RN reduz o custo computacional de gerar previsões de Densidade de Taxa de Formação Estelar (SFRD) por um fator de $\sim 200$ em comparação com o SAM completo, permitindo que o volume completo seja simulado em $\sim 21$ horas-CPU.
4. Cálculo do Sinal de 21 cm: O SFRD simulado é pós-processado para calcular a temperatura de brilho de 21 cm ($\delta T_b$) e o espectro de potência ($\Delta^2(k)$). Isso envolve:
   • Calcular o fundo de Lyman-$\alpha$ ($J_\alpha$) para determinar o coeficiente de acoplamento Wouthuysen-Field (W-F) ($x_\alpha$).
   • Resolver para a temperatura de spin ($T_S$) e a temperatura cinética ($T_K$) do meio intergaláctico (MIG), contabilizando o aquecimento por Compton e o aquecimento por raios-X (assumido para rastrear o SFR).
   • Resolver iterativamente para $T_S$ para determinar $\delta T_b$ e seu espectro de potência.

Principais Contribuições

• Calibração Novel: Esta é a primeira simulação seminumerica em grande escala a unir autoconsistentemente árvores estocásticas de fusão de halos de ME com prescrições de FE explicitamente calibradas às simulações hidrodinâmicas modernas (Renaissance e AEOS).
• Modelagem PopII de Duas Fases: O framework introduz uma transição de "explosiva" para "estável" para a FE PopII, calibrada aos resultados hidrodinâmicos, o que difere significativamente dos modelos anteriores de eficiência constante.
• Histórias Estocásticas de Fusão: A inclusão de histórias estocásticas de fusão de halos de ME permite um tratamento realista da transição da FE PopIII para PopII, incluindo um período de atraso ($t_{delay}$) e os efeitos de $v_{bc}$ sobre $M_{crit}$.
• Eficiência Computacional: A arquitetura otimizada de RN e o conjunto de treinamento reduzido permitem múltiplas realizações com diferentes pressupostos astrofísicos (por exemplo, diferentes massas PopIII, $t_{delay}$, $M_{crit}$) serem executadas eficientemente.

Resultados
Os autores apresentam resultados de um modelo fiducial e várias realizações alternativas (variando massa PopIII, $M_{crit}$ e $t_{delay}$):

• Sinal Global ($\delta T_b$):
  • O modelo fiducial prevê um mínimo de vale de absorção de $\delta T_b \approx -146$ mK em $z \approx 20.25$.
  • Impacto PopIII vs. PopII: A FE PopIII domina o sinal de 21 cm em $z \gtrsim 34$ e em grandes escalas. A FE PopII domina em $z \lesssim 20$ e em escalas menores ($k \gtrsim 0.2$ Mpc$^{-1}$ em $z \sim 20-34$).
  • Período de Atraso ($t_{delay}$): O atraso entre a FE PopIII e PopII impacta significativamente o sinal. Um atraso mais longo (por exemplo, 100 Myr) resulta em um vale de absorção mais tardio, mais profundo e mais estreito em comparação com o atraso fiducial de 10 Myr.
  • Comparação: O modelo fiducial produz vales de absorção mais precoces e profundos do que muitos modelos anteriores (por exemplo, Muñoz et al. 2022) devido à combinação de árvores de fusão de ME, auto-escudo de H$_2$ e eficiências de FE calibradas.
• Espectro de Potência ($\Delta^2(k)$):
  • Os espectros de potência são detectáveis pelo HERA a partir de $z \approx 25$ para $k \sim 0.2$ Mpc$^{-1}$, expandindo-se para $k \lesssim 0.9$ Mpc$^{-1}$ até $z=15$ sob suposições moderadas de foreground (1080 horas de observação).
  • Modelos com histórias estocásticas de fusão mostram potência em grande escala aumentada em altos desvios para o vermelho em comparação com modelos sem elas (como o modelo de Muñoz et al. 2022), impulsionado por flutuações de $v_{bc}$ afetando $M_{crit}$.
  • Precisão da RN: Os emuladores reproduzem os resultados do SAM com uma precisão média de $\sim 8.3\%$ para PopII e $\sim 5.0\%$ para PopIII em todos os desvios para o vermelho e modelos.

Significado e Alegações
O artigo afirma que a combinação de calibração de FE hidrodinâmica, um modelo de massa crítica de halo que considera o auto-escudo de H$_2$ e histórias estocásticas de fusão de halos de ME resulta em formação estelar mais precoce e SFRDs mais altos durante o Amanhecer Cósmico em comparação com trabalhos seminumericos anteriores.

Os autores afirmam que:

1. Observabilidade: O sinal de 21 cm fiducial é potencialmente detectável pelo HERA em $z \lesssim 25$.
2. Restrições à Física: A falta de detecção em $z \gtrsim 20$ sugeriria um período de atraso ($t_{delay}$) entre a FE PopIII e PopII de $\gtrsim 30$ Myr.
3. Necessidade Metodológica: Contabilizar adequadamente a transição entre a FE PopIII e PopII requer a inclusão de histórias de fusão de halos de ME, pois modelos analíticos simplificados (como integrais de HMF) falham em capturar a estocasticidade e o timing dessa transição.
4. Diferenciação: Embora distinguir entre variações sutis na eficiência de FE possa ser difícil, a sensibilidade do espectro de potência do HERA é suficiente para distinguir entre modelos com diferentes valores de $t_{delay}$ e para diferenciar o modelo dos autores (com árvores de fusão e calibração hidrodinâmica) de modelos que carecem desses recursos (como a realização de Muñoz et al. 2022).

O trabalho posiciona-se como uma ferramenta novel e complementar aos códigos seminumericos existentes (por exemplo, 21cmFAST, Zeus21), oferecendo uma abordagem mais fundamentada fisicamente para prever o sinal de 21 cm para campanhas observacionais atuais e futuras.