Formation of intermediate-mass black holes in young massive clusters detected with JWST: analytic mass estimates

Resumo Técnico: Formação de Buracos Negros de Massa Intermediária em Aglomerados Estelares Massivos Jovens Detectados pelo JWST

Declaração do Problema
O Telescópio Espacial James Webb (JWST) revelou uma população de sistemas estelares densos e massivos em alto desvio para o vermelho ($z \sim 4-10$), incluindo o arco "Cosmic Gems", o "Firefly Sparkle" e galáxias Little Red Dot (LRD). Esses sistemas exibem massas estelares variando de $10^5$ a $10^7 M_\odot$ e raios de meia-massa de $\sim 0,4$ a $15$ pc. Embora esses ambientes sejam teoricamente ideais para a formação de buracos negros de massa intermediária (BNMI) via colisões estelares descontroladas, simulações diretas de $N$-corpos são computacionalmente inviáveis para sistemas tão massivos e grandes populações. Além disso, a eficiência da formação de BNMI nesses ambientes específicos de alto desvio para o vermelho, particularmente no que diz respeito à interação entre taxas de colisão e perda de massa por ventos estelares, requer estimativa quantitativa para avaliar sua viabilidade como sementes para os buracos negros supermassivos (BNSM) observados no Universo primordial.

Metodologia
Para abordar a inviabilidade computacional de simulações diretas de $N$-corpos para esses aglomerados massivos, os autores empregam um quadro híbrido analítico e semianalítico:

1. Evolução Dinâmica (Fokker-Planck): Os autores utilizam o solver PhaseFlow (parte da biblioteca Agama) para resolver a equação de Fokker-Planck unidimensional média de órbita. O modelo assume aglomerados compostos por estrelas de massa igual distribuídas conforme um perfil de densidade de Plummer. A evolução é rastreada ao longo de uma escala de tempo característica de 100 Myr. Embora se reconheça que aglomerados reais evoluem mais rápido devido à segregação de massa (que não é explicitamente modelada), essa abordagem fornece uma estimativa conservadora do colapso do núcleo e da evolução da densidade central.
2. Modelo de Colisão Descontrolada: O crescimento de uma Estrela Muito Massiva (VMS) central é modelado usando um quadro analítico baseado em Portegies Zwart & McMillan (2002), estendido por Fujii et al. (2024) e Pacucci et al. (2025) para incluir perda de massa por ventos. A taxa de acreção de massa ($\dot{M}_{acc}$) é derivada da taxa de colisão no núcleo do aglomerado, dependendo da massa do núcleo ($M_{core}$), do tempo de relaxamento do núcleo ($t_{relax, core}$) e do logaritmo de Coulomb.
3. Equilíbrio de Perda de Massa: O crescimento líquido da VMS é determinado equilibrando o ganho de massa via colisões contra a perda de massa via ventos estelares. A taxa de perda de massa por ventos ($\dot{M}_{wind}$) é calculada usando a prescrição de Vink (2018), que depende fortemente da metalicidade ($Z$).
4. Validação: O quadro analítico é validado contra simulações diretas de $N$-corpos (Dragon II, Vergara et al. 2025b, e aglomerados FROST) onde tais comparações são possíveis. Os autores encontram que seu modelo reproduz a ordem de grandeza das massas de BNMI dentro de um fator de 2–5, consistente com as incertezas esperadas de modelos simplificados.

Resultados Principais
A aplicação desse quadro a uma amostra de aglomerados observados pelo JWST produz as seguintes descobertas:

• Faixa de Massa de BNMI: As massas finais estimadas de BNMI (assumindo colapso direto da VMS) variam tipicamente de $\sim 10^2 M_\odot$ a $\sim 4 \times 10^3 M_\odot$.
• Resultados Específicos por Aglomerado:
  • Cosmic Gems (Adamo et al. 2024): Esses aglomerados são extremamente compactos ($R_h \sim 1$ pc) com altas massas estelares ($\sim 10^6 M_\odot$) e baixa metalicidade ($Z \approx 0,005 Z_\odot$). O modelo prevê a formação das sementes mais massivas da amostra, variando de $1,7 \times 10^3$ a $2,7 \times 10^3 M_\odot$.
  • Firefly Sparkle (Mowla 2024): Esses aglomerados são mais extensos ($R_h \sim 4-6$ pc) com menores densidades. As massas de BNMI previstas são menores, variando de $\sim 1,8 \times 10^2$ a $8,1 \times 10^2 M_\odot$.
  • Sunburst/Campos de Fronteira (Vanzella et al.): Apesar dos grandes raios, as altas massas totais ($>10^7 M_\odot$) permitem núcleos massivos, gerando BNMI de até $\sim 2,1 \times 10^3 M_\odot$.
• Eficiência de Formação: A eficiência de formação (definida como $M_{BH}/M_{cl}$) é encontrada na ordem de alguns por cento. Isso alinha-se com ajustes derivados de simulações diretas de $N$-corpos por Vergara et al. (2025a).
• Dependência da Metalicidade: Baixa metalicidade ($Z \lesssim 0,02 Z_\odot$) é identificada como um fator crítico. Metalicidades mais altas aumentam significativamente a perda de massa por ventos, suprimindo a massa final da semente. Por exemplo, aglomerados modelados com $Z = 0,1 Z_\odot$ (Messa et al. 2025) mostram eficiências de formação mais baixas em comparação com contrapartes pobres em metais de massa similar.

Significado e Alegações
O artigo afirma que aglomerados estelares massivos jovens, densos e pobres em metais detectados pelo JWST são fábricas viáveis para "sementes pesadas" de buracos negros de massa intermediária. As massas estimadas ($10^2 - 4 \times 10^3 M_\odot$) são suficientes para servir como sementes que podem crescer até os BNSM observados em alto desvio para o vermelho ($z \sim 6-10$) dentro do tempo cósmico disponível, aliviando as restrições temporais associadas a sementes mais leves (por exemplo, remanescentes da População III).

Especificamente, os autores sugerem que a extrema compactação dos aglomerados Cosmic Gems facilita a formação de sementes com massas de $1600-2700 M_\odot$. Essas sementes pesadas exigiriam significativamente menos tempos de dobra (e-folding times) para atingir $10^7-10^8 M_\odot$ por $z \sim 6$ em comparação com sementes mais leves. Além disso, os autores notam que essas sementes pesadas poderiam desempenhar um papel relevante na emissão de ondas gravitacionais de fusões de BNMI, potencialmente detectáveis por futuros observatórios como o LISA e o Einstein Telescope.

O estudo enfatiza que, embora o modelo seja simplificado e dependa de suposições (como estrelas de massa igual e negligência de binários primordiais), os resultados são consistentes com simulações numéricas mais complexas e fornecem uma estimativa robusta de ordem de grandeza para a formação de BNMI no Universo primordial. Os autores reconhecem incertezas quanto à perda de massa durante colisões e ao potencial de funções de massa inicial pesadas no topo, que poderiam aumentar ainda mais a formação de BNMI, mas não foram explicitamente modeladas.