Impact of initial mass function on the chemical evolution of high-redshift galaxies

Título: Impacto da Função de Massa Inicial na Evolução Química de Galáxias de Alto Redshift

1. Problema e Contexto

Observações recentes do Telescópio Espacial James Webb (JWST) revelaram uma relação invariante entre massa estelar, metalicidade e taxa de formação estelar (relação M-Z-SFR) em galáxias até $z \sim 8$, seguida de uma ruptura em redshifts mais altos. Além disso, há um "excesso" de galáxias luminosas no UV e assinaturas químicas peculiares (como enriquecimento de N e C) em galáxias de alto redshift ($z \gtrsim 5$).
O problema central é entender a física subjacente que regula essas correlações. Especificamente, o artigo investiga se a Função de Massa Inicial (IMF) das estrelas, particularmente sua extensão para massas muito elevadas e a presença de estrelas supermassivas ($\gtrsim 100 M_\odot$), é o fator determinante para explicar a evolução química e a história de formação estelar nessas épocas cósmicas.

2. Metodologia

Os autores utilizaram o código semi-analítico de evolução galáctica a-sloth, aplicado a árvores de fusão derivadas de uma simulação cosmológica de alta resolução (Ishiyama et al. 2016) no universo $\Lambda$CDM.

Principais Atualizações e Configurações do Modelo:

• Prescrição de Formação Estelar: Atualização do raio característico da região de formação estelar ($R_*$). Em vez de usar o raio de escala do halo ($R_s$), o modelo agora utiliza $R_* = 2.5 R_{50}$ (onde $R_{50}$ é o raio de meia-massa estelar), baseado em simulações "zoom-in" e observações, resultando em uma dispersão na sequência principal de formação estelar mais consistente com dados observacionais.
• Populações Estelares e IMF:
  • População III (Pop III): Fixada com parâmetros de melhor ajuste anteriores ($m_{min}=13.6 M_\odot$, $m_{max}=197 M_\odot$, $\alpha=1.77$).
  • População II (Pop II): A IMF segue uma forma tipo Kroupa, mas o limite superior de massa ($m_{max}$) é variado entre 100 e 600 $M_\odot$. O foco é testar como a presença de estrelas muito massivas afeta o enriquecimento químico.
• Feedback Estelar e Evolução:
  • Uso de novos tracks de evolução estelar (PARSEC v2.0) cobrindo uma ampla faixa de metalicidade ($Z \sim 10^{-11} - 0.03$) e massa ($2 - 600 M_\odot$).
  • Inclusão de rendimentos metálicos de ventos estelares, supernovas de colapso do núcleo (CCSN), supernovas de instabilidade de par (PISNe), supernovas pulsacionais de instabilidade de par (PPISNe) e supernovas do Tipo Ia.
  • Modelagem de ventos galácticos (outflows) com eficiência dependente da massa do halo ($M_h$) e da energia das supernovas.
• Espaço de Parâmetros: Foram executados 165 modelos variando três parâmetros chave:
  • $m_{max}$ (limite superior da IMF Pop II): $100 - 600 M_\odot$.
  • $\alpha_{out}$ (inclinação da eficiência de outflow): $0, 0.5, 1$.
  • $M_{out0}$ (escala de massa de normalização do outflow): $1 - 5 \times 10^9 M_\odot$.

3. Restrições Observacionais

Os modelos foram confrontados com três conjuntos de dados observacionais para calcular uma verossimilhança global ($L_{all}$):

1. História de Formação Estelar Cósmica (CSFH): Densidade de taxa de formação estelar (SFRD) e densidade de massa estelar cósmica (CSMD) em $z \sim 4-16$.
2. Evolução Química (MZSFR): A relação invariante Massa-Metalicidade-Taxa de Formação Estelar (FMR e MZR) em $z \sim 4-8$ e sua ruptura em $z > 8$.
3. Conexão Galáxia-Halo: Relação Massa Estelar-Massa do Halo (SHMR).

4. Resultados Principais

• Necessidade de Estrelas Muito Massivas ($m_{max} \gtrsim 200 M_\odot$):
  • Modelos com $m_{max} \le 100 M_\odot$ falham em reproduzir simultaneamente a CSFH e a relação MZSFR observada. Eles produzem rendimentos metálicos insuficientes.
  • O melhor ajuste é alcançado com $m_{max} \approx 250 M_\odot$. Estrelas com massas superiores a $200 M_\odot$ produzem Supernovas de Instabilidade de Par (PISNe), que liberam quantidades massivas de metais, elevando o rendimento metálico médio por unidade de massa estelar para $\sim 0.025$.
• Parâmetros de Outflow:
  • O melhor ajuste ocorre com $\alpha_{out} = 0$ (eficiência de outflow independente da massa do halo) e $M_{out0} = 3 \times 10^9 M_\odot$.
  • Isso sugere que galáxias de alto redshift (especialmente as anãs) são menos vulneráveis a feedbacks de outflow do que suas contrapartes locais, permitindo maior retenção de metais e formação estelar eficiente.
• Reprodução de Observações:
  • O modelo de melhor ajuste ($m_{max}=250 M_\odot, \alpha_{out}=0$) reproduz a CSFH e a MZSFR dentro de $1\sigma$ para $z \sim 4-13$.
  • O modelo explica a ruptura da relação MZSFR em $z > 8$ (galáxias mais pobres em metais) devido ao influxo de gás primordial não enriquecido.
  • O modelo prevê uma população significativa de galáxias fracas e de baixa massa ($M_* \sim 10^3 M_\odot$) que contribuem com $>90\%$ da SFRD total em $z \gtrsim 12$, mas que são invisíveis nas observações atuais do JWST.

5. Contribuições e Significância

• Papel Crucial das PISNe: O estudo fornece evidência teórica robusta de que estrelas muito massivas ($\gtrsim 200 M_\odot$) e suas subsequentes PISNe são essenciais para o enriquecimento químico do Universo primordial, resolvendo a discrepância entre modelos tradicionais e os dados do JWST.
• Implicações para Transientes e Ondas Gravitacionais:
  • Taxa de PISNe: O modelo prevê uma taxa de PISNe de População II muito alta, picando em $z \sim 5$ ($\sim 10 \text{ eventos/ano/deg}^2$), muito superior à prevista apenas por Pop III.
  • Fusões de Buracos Negros Binários: A alta taxa de formação estelar em galáxias pobres em metais (predominantes em $z \gtrsim 8$) implica uma taxa de fusão de buracos negros binários significativamente maior do que estimativas anteriores baseadas apenas em galáxias luminosas. Isso tem implicações diretas para detectores de ondas gravitacionais de 3ª geração.
• Reionização: A população de galáxias fracas e pobres em metais, dominantes na SFRD em altos redshifts, é identificada como provável motor principal da reionização cósmica.

6. Conclusão

O trabalho conclui que a suposição de uma IMF universal e "padrão" (limitada a $\sim 100 M_\odot$) é insuficiente para explicar a evolução química de galáxias em $z \gtrsim 4$. A existência de uma IMF "top-heavy" estendida até $\sim 250 M_\odot$, impulsionando PISNe frequentes, é necessária para reconciliar a história de formação estelar cósmica com as relações de escalas químicas observadas pelo JWST. Isso redefine nossa compreensão sobre a formação estelar primordial, a reionização e a origem de eventos transientes e de ondas gravitacionais no Universo jovem.