The THESAN project: Lyman-alpha emitters as probes of ionized bubble sizes

1. Problema e Contexto

Durante a Época da Reionização (EoR), ocorrida entre $z \approx 5-20$, o hidrogênio neutro do meio intergaláctico (IGM) foi ionizado pelas primeiras gerações de galáxias. Este processo não foi uniforme, mas sim "em manchas", formando bolhas ionizadas que cresceram e se fundiram.

• O Desafio: Medir diretamente o tamanho dessas bolhas ionizadas é extremamente difícil. A linha de emissão Lyman-$\alpha$ (Ly$\alpha$) é uma das sondas observacionais mais promissoras, pois sua detecção depende criticamente do tamanho da bolha ionizada ao redor da galáxia hospedeira. Em bolhas pequenas, os fótons Ly$\alpha$ sofrem espalhamento ressonante e são absorvidos pelo hidrogênio neutro circundante; em bolhas grandes, a transmissão é alta.
• A Lacuna: Interpretar observações de Ly$\alpha$ para inferir tamanhos de bolhas é complexo devido a degenerescências físicas (poeira, ventos galácticos, transferência radiativa no meio interestelar e circumgaláctico) e à simplificação excessiva de modelos analíticos esféricos ideais. É necessário um modelo teórico robusto e calibrado para conectar as propriedades observáveis das galáxias à topologia da reionização.

2. Metodologia

O estudo utiliza a suite de simulações cosmológicas de hidrodinâmica radiativa thesan, que combina a física de formação de galáxias do modelo IllustrisTNG com transporte de radiação ionizante em tempo real (usando o código arepo-rt).

• Simulações: A simulação principal (thesan-1) possui uma caixa de $95.5$ cMpc e alta resolução, permitindo resolver halos de resfriamento atômico e capturar a evolução de bolhas ionizadas de forma autoconsistente.
• Modelagem de Ly$\alpha$:
  • Emissão Intrínseca: Calculada a partir de recombinações, excitação colisional e emissão estelar, considerando a fração de fuga de fótons ionizantes ($f_{esc}$).
  • Transmissão no IGM: Utiliza rastreamento de raios (768 direções) a partir do raio virial das galáxias para calcular a transmissão dependente da frequência através do IGM, incluindo efeitos de amortecimento de asa (damping wing).
  • Calibração Empírica: Para corrigir processos em pequena escala não resolvidos (poeira ISM/CGM e ventos), os autores aplicaram um modelo empírico calibrado. Este modelo ajusta o perfil espectral (desvio de velocidade, largura) e a atenuação por poeira, calibrado contra as funções de luminosidade (LF) observadas de Ly$\alpha$ em $z=5.5$ e $z=6.6$ (dados de Ouchi et al. e Konno et al.).
• Análise de Bolhas: O tamanho das bolhas é definido pelo caminho livre médio (MFP) dos raios ionizantes a partir da galáxia. Os autores correlacionam as propriedades das galáxias (magnitude UV, luminosidade Ly$\alpha$, Equivalente de Largura - EW) com o tamanho da bolha ao longo da linha de visada.

3. Principais Contribuições

1. Catálogo Calibrado de Emissores Ly$\alpha$ (LAEs): Desenvolvimento de um modelo atualizado que gera catálogos realistas de LAEs a partir das simulações thesan, incorporando efeitos de poeira e ventos de forma consistente com observações.
2. Mapeamento Estatístico: Estabelecimento de uma ligação estatística robusta entre a fração de LAes observáveis e a distribuição de tamanhos de bolhas ionizadas, validando o uso de LAEs como sondas da topologia da reionização.
3. Identificação de Vieses: Demonstração de como modelos analíticos simplificados podem falhar ao não considerar a estrutura 3D inhomogênea do IGM e os efeitos de transferência radiativa internos às galáxias.
4. Disponibilidade de Dados: Os catálogos calibrados serão disponibilizados publicamente para auxiliar na interpretação de levantamentos do JWST e de banda estreita.

4. Resultados Chave

• Correlação com Tamanho da Bolha: Antes do ponto médio da reionização ($z \gtrsim 7$), galáxias dentro de bolhas maiores ($\gtrsim 10$ cMpc) exibem luminosidades Ly$\alpha$ observadas e Equivalentes de Largura (EW) significativamente maiores. Isso confirma que a fração de LAEs é um traçador prático do tamanho das bolhas.
• Seleção de LAEs: Em amostras selecionadas por $L_{Ly\alpha} > 10^{41.5}$ erg s$^{-1}$, as propriedades do Ly$\alpha$ (especialmente a luminosidade e o EW) correlacionam-se mais fortemente com o tamanho da bolha do que a magnitude UV, especialmente em $z \gtrsim 7$.
• Evolução Temporal:
  • A correlação entre o tamanho da bolha e a visibilidade do Ly$\alpha$ enfraquece à medida que a reionização avança e as bolhas se fundem (percolação), tornando o IGM cada vez mais ionizado.
  • Galáxias mais brilhantes em UV tendem a ter menores frações de fuga e EWs devido à maior atenuação por poeira e ventos, mas residem preferencialmente em bolhas maiores.
• Tensões com Observações: O modelo reproduz bem a queda na fração de LAEs em altos redshifts ($z \gtrsim 8$), mas tende a superestimar a fração de LAEs em redshifts mais baixos ($z < 7$) em comparação com dados recentes do JWST. Isso sugere que a história de reionização nas simulações thesan pode terminar um pouco cedo demais ou que há incertezas na modelagem do continuum e da poeira.
• Densidade de LAEs: A densidade de LAEs em bolhas grandes supera a de bolhas pequenas a partir de $z \approx 7.5$, e a razão entre a densidade em bolhas grandes e a densidade total segue de perto a história de ionização global.

5. Significado e Conclusões

Este trabalho fornece um framework baseado em simulações de alta fidelidade para interpretar observações de Ly$\alpha$ no contexto da reionização.

• Validação de Sondas: Confirma que a fração de LAEs e suas propriedades espectrais são sondas viáveis para inferir o tamanho das bolhas ionizadas, especialmente na primeira metade da reionização.
• Ferramenta para o JWST: O framework e os catálogos gerados são essenciais para interpretar os dados atuais e futuros do Telescópio Espacial James Webb (JWST) e levantamentos de banda estreita, permitindo extrair informações sobre a topologia da reionização a partir de estatísticas de galáxias individuais.
• Correção de Vieses: Ao quantificar como a poeira e a transferência radiativa afetam a visibilidade, o estudo ajuda a corrigir vieses em modelos mais idealizados, oferecendo uma estimativa mais precisa da evolução da reionização.

Em suma, o artigo conecta a física de formação de galáxias e a transferência radiativa com a cosmologia de larga escala, demonstrando que a observação de emissores Ly$\alpha$ pode revelar a estrutura "em manchas" do Universo primordial.