Galaxy Proximate Damped Lyman-Alpha Systems and HI Reionization Topology in TECHNICOLOR DAWN

Resumo Técnico: Sistemas Proximos Damped Lyman-Alpha em Galáxias e Topologia da Reionização do HI em TECHNICOLOR DAWN

Declaração do Problema
Observações recentes do Telescópio Espacial James Webb (JWST) identificaram sistemas proximos damped Lyman-$\alpha$ (PDLAs) no primeiro plano de galáxias de alto desvio para o vermelho ($z > 5$). Estes sistemas são interpretados como meios circumgalácticos (CGM) neutros. Embora os PDLAs tenham sido historicamente utilizados para estudar ambientes de quasares, a sua presença em torno de galáxias de alto desvio para o vermelho oferece uma nova sonda potencial para rastrear o progresso da reionização cósmica. No entanto, a utilidade deste método é incerta. Permanece incerto se a história de ionização do CGM espelha a do meio intergaláctico (IGM), dada a topologia complexa da reionização. Além disso, é desconhecido se a densidade de coluna observada de hidrogénio neutro ($N_{\text{HI}}$) nos PDLAs correlaciona-se primariamente com a fração neutra global (um proxy para o progresso da reionização) ou se é impulsionada por outros fatores, como massa do halo, fração de massa gasosa ou fração neutra local.

Metodologia
Os autores utilizam a simulação hidrodinâmica cosmológica TECHNICOLOR DAWN para investigar estas questões. A simulação emprega um modelo $\Lambda$CDM com parâmetros cosmológicos do Planck, utilizando hidrodinâmica de partículas suavizadas (SPH) e transferência radiativa em tempo real via o método dos momentos. A caixa de simulação tem comprimentos de lado de 16,5 cMpc/$h$ e inclui $2 \times 704^3$ partículas.

A análise prossegue em duas etapas:

1. Análise de Topologia: Os autores calculam a sobre-densidade ($\Delta$) e a fração de hidrogénio neutro ($x_{\text{HI}}$) para partículas de gás individuais através de instantâneos de desvio para o vermelho ($z = 10, 8, 6, 5.5$). Analisam a razão entre as frações neutras médias em massa e médias em volume para caracterizar a topologia da reionização.
2. Análise Halo por Halo: As galáxias são identificadas usando o algoritmo SKID, e os seus halos de matéria escura hospedeiros são definidos via sobre-densidade esférica. O estudo foca em halos resolvidos com massas $M > 10^{8.5} M_\odot$. Os autores calculam a densidade de coluna de hidrogénio neutro ($N_{\text{HI}}$) no primeiro plano destes halos, ponderando o cálculo pela taxa de formação estelar local (SFR) para imitar linhas de visada observacionais através de regiões de formação estelar. Em seguida, correlacionam $N_{\text{HI}}$ com a massa do halo, massa estelar, fração de massa gasosa ($f_g$) e fração neutra local ponderada em massa ($\langle f_{\text{HI}} \rangle$).

Finalmente, os autores ajustam uma relação funcional entre $\log_{10}(N_{\text{HI}})$, massa do halo/estelar ($M$) e a fração neutra média em volume ($\langle x_{\text{HI}} \rangle$) para determinar a viabilidade de usar PDLAs para inferir o progresso da reionização.

Principais Resultados

• Topologia da Reionização: A simulação revela uma topologia de reionização "de dentro para fora para o meio" (IOM). Inicialmente, a reionização procede de dentro para fora, ionizando regiões densas imediatamente circundantes às fontes. Isto é seguido por uma fase dominante de fora para dentro, onde o fluxo ionizante escapa para regiões de baixa densidade, ionizando-as rapidamente devido a baixas taxas de recombinação. Finalmente, as frentes de ionização propagam-se de volta para regiões de densidade intermédia (filamentos). Consequentemente, o CGM reioniza após o IGM e permanece parcialmente neutro mesmo em $z = 5.5$.
• Impulsionadores da Densidade de Coluna: O estudo encontra que a densidade de coluna de hidrogénio neutro ($N_{\text{HI}}$) não é impulsionada primariamente pela fração de massa gasosa ou pela fração neutra local. Em vez disso, $N_{\text{HI}}$ mostra uma correlação clara e dominante com a massa do halo (e massa estelar). Halos massivos exibem $N_{\text{HI}}$ mais elevado mesmo quando as suas frações neutras locais são baixas, simplesmente porque contêm mais massa gasosa total.
• Correlação com o Progresso da Reionização: Apesar da dominância da massa, existe uma relação estatisticamente significativa entre $N_{\text{HI}}$ e a fração neutra global. Os autores estabelecem a seguinte forma funcional:
  $$\log_{10}(N_{\text{HI}}) = a + b \log_{10}(M) + c \langle x_{\text{HI}} \rangle$$
  onde $b$ (dependência da massa) é o termo dominante, mas $c$ (dependência da fração neutra) é não nulo.
• Viabilidade Observacional: Os autores realizam uma análise de propagação de erros para estimar o número de galáxias necessário para restringir $\langle x_{\text{HI}} \rangle$ a uma precisão fracionária específica. Encontram que o método é mais sensível em altos desvios para o vermelho, onde $\langle x_{\text{HI}} \rangle$ é elevado. À medida que a reionização se aproxima da conclusão e $\langle x_{\text{HI}} \rangle$ diminui, o número necessário de observações escala como $\langle x_{\text{HI}} \rangle^{-2}$, tornando o método menos sensível em desvios para o vermelho mais baixos.

Significado e Alegações
O artigo alega que os PDLAs de galáxias representam um método observacional viável e novo para rastrear o progresso da reionização cósmica, desde que existam estimativas precisas da massa do halo ou estelar para desvendar a dependência da massa do estado de ionização.

Os autores enfatizam que as suas descobertas suportam a interpretação de que altas densidades de coluna nos PDLAs não são meramente indicadores de altas frações gasosas ou ambientes localmente neutros, mas estão fortemente ligadas à massa da galáxia hospedeira. Ao contabilizar esta dependência da massa, os PDLAs podem servir como um traçador para a fração neutra global. O estudo sugere que este método é particularmente bem adequado para restringir a reionização precoce (alto $z$) usando dados atuais e projetados do JWST, embora note que variações de linha de visada para linha de visada e tamanhos limitados da caixa de simulação possam introduzir dispersão que requer investigação adicional. Os resultados alinham-se com simulações anteriores (por exemplo, Finlator et al. 2009) relativamente à topologia IOM e oferecem consistência com tendências observacionais relatadas na literatura recente (por exemplo, Heintz et al. 2024; Umeda et al. 2024).