Impact of subhalo dynamical friction heating on the formation of the first structures in the universe

1. O Problema

A história térmica do gás em halos de matéria escura regula fundamentalmente a formação estelar. Embora processos de feedback (como supernovas e AGNs) e resfriamento sejam bem estudados, o papel do aquecimento gravitacional via fricção dinâmica (DF) de subhalos orbitantes em diferentes épocas cósmicas e massas de halo permanece menos explorado, especialmente no universo primordial.

• Contexto: Em halos massivos ($z \approx 0$), o aquecimento gravitacional é conhecido por suprimir o resfriamento do gás e "apagar" (quench) galáxias. No entanto, a evolução em redshift e os limites de massa para esse efeito não foram totalmente quantificados.
• Questão Central: O aquecimento por fricção dinâmica de subhalos é significativo o suficiente no universo primordial ($z \sim 15$) para competir com o resfriamento molecular de hidrogênio ($H_2$) em minihalos? Se sim, isso poderia suprimir a formação das primeiras estrelas (População III) e favorecer a formação de sementes de Buracos Negros de Colapso Direto (DCBH).

2. Metodologia

Os autores desenvolveram um modelo analítico para calcular as taxas de aquecimento por fricção dinâmica, validado e calibrado com dados de simulações cosmológicas de alta resolução.

• Dados de Simulação: Utilizaram o catálogo de halos e subhalos da simulação TNG50 (parte do projeto IllustrisTNG) para derivar a Função de Massa de Subhalos (SHMF) em uma ampla faixa de redshifts ($z=15$ a $z=0$). Para minihalos de baixa massa não resolvidos no TNG50, utilizaram simulações de alta resolução de Smith et al. (2015).
• Modelo Analítico:
  • Basearam-se no modelo de perturbação linear dependente do tempo de Ostriker (1999) para calcular a força de arrasto e a taxa de aquecimento.
  • A taxa de aquecimento total ($\dot{E}_{DF}$) é obtida integrando a SHMF sobre as massas dos subhalos, considerando a densidade do gás, a velocidade do subhalo e o fator de eficiência da fricção dinâmica.
  • Correções Não-Lineares: Investigaram correções para perturbações não-lineares (baseadas em Kim & Kim 2009) e a dependência do número de Mach, utilizando dados do TNG50 para estimar distribuições de velocidade e parâmetros de perturbação.
• Análise Estatística: Realizaram amostragem de Monte Carlo (500 realizações) para quantificar a dispersão (scatter) na taxa de aquecimento devido à variância intrínseca no número de subhalos (flutuações de Poisson).
• Comparação com Resfriamento: Compararam as taxas de aquecimento DF com as taxas de resfriamento (molecular e atômico) calculadas usando a biblioteca química Grackle, considerando diferentes metalicidades e frações de $H_2$.

3. Principais Contribuições e Resultados

A. Evolução da Função de Massa de Subhalos (SHMF)

• Derivaram um ajuste analítico para a SHMF evoluída em função do redshift.
• Encontraram que, em altos redshifts ($z > 6$), a SHMF se assemelha mais à SHMF "não evoluída" (antes do arrasto de maré), com uma inclinação mais rasa ($\alpha$ menor) e amplitude diferente em comparação com $z=0$.
• A dispersão no número de subhalos é essencialmente Poissoniana, levando a uma variação de até duas ordens de magnitude na taxa de aquecimento DF para uma massa de halo hospedeiro fixa.

B. Aquecimento em Halos Massivos de Baixo Redshift ($z \le 6$)

• Confirmaram que a fricção dinâmica é um mecanismo importante para o apagamento de galáxias em halos massivos ($M_{halo} \gtrsim 10^{13} M_\odot$) em redshifts baixos.
• A taxa de aquecimento DF média é menor do que o aquecimento gravitacional total estimado em estudos anteriores (como DB07), provavelmente porque estudos anteriores ignoraram o arrasto de maré (que reduz a massa dos subhalos) e trataram o aquecimento como um limite superior.
• O aquecimento DF é altamente variável; em halos com poucos subhalos, ele pode ser insuficiente para suprimir a formação estelar, mesmo em halos massivos.

C. Impacto no Universo Primordial ($z \sim 15$) e Formação de Pop III/DCBH

• Competição Térmica: Em minihalos ($10^6 - 10^8 M_\odot$) em $z=15$, as taxas de aquecimento DF podem ser comparáveis às taxas de resfriamento de $H_2$ quando a fração de $H_2$ é baixa ($\sim 10^{-5}$).
• Supressão de Pop III: O aquecimento DF aumenta a fração crítica de $H_2$ necessária para o colapso e formação de estrelas da População III em até uma ordem de magnitude. Isso torna a formação de estrelas Pop III mais difícil do que o previsto por modelos que ignoram esse aquecimento.
• Favorecimento de DCBH: Em combinação com o fundo de radiação Lyman-Werner (LW), o aquecimento DF favorece a formação de sementes de DCBH em halos de resfriamento atômico. Mesmo com uma fração média de $H_2$ de $\sim 10^{-5}$ durante a fase de progenitor do minihalo, o aquecimento DF pode suprimir o resfriamento molecular suficiente para permitir o colapso isotérmico via resfriamento atômico, levando a um DCBH.
• Resolução Necessária: Para capturar corretamente esse efeito em simulações, é necessário resolver subhalos com uma razão de massa subhalo/halo hospedeiro de $\psi \gtrsim 0.05$.

4. Significado e Conclusões

O trabalho estabelece que a fricção dinâmica de subhalos é um mecanismo de aquecimento gravitacional contínuo e ubíquo no modelo $\Lambda$CDM, que não depende de eventos explosivos ou de feedback estelar/AGN.

• Relevância Cosmológica: O mecanismo é crucial para entender a transição entre a formação de estrelas Pop III e a formação de DCBHs. Ele oferece uma explicação alternativa ou complementar para a supressão de estrelas Pop III em ambientes onde a radiação LW sozinha não seria suficiente.
• Implicações para Simulações: Os resultados destacam a necessidade crítica de alta resolução em simulações cosmológicas hidrodinâmicas. Simulações que não resolvem subhalos acima de $\psi \approx 0.05$ subestimam significativamente o aquecimento do gás e, consequentemente, podem prever incorretamente a formação de estrelas primordiais ou a abundância de buracos negros supermassivos.
• Modelagem Futura: O modelo analítico desenvolvido serve como uma ferramenta eficiente para incorporar o aquecimento DF em modelos semi-analíticos de formação de galáxias e sementes de buracos negros, onde simulações de N-corpos de alta resolução são computacionalmente proibitivas.

Em suma, o aquecimento por fricção dinâmica é um fator chave, frequentemente negligenciado, que regula o balanço térmico do gás no universo primordial e influencia diretamente o destino das primeiras estruturas cósmicas.