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

Riepilogo Tecnico: Formazione di Buchi Neri di Massa Intermedia in Ammassi Stellari Giovani e Massicci Rilevati da JWST

Enunciato del Problema
Il James Webb Space Telescope (JWST) ha rivelato una popolazione di sistemi stellari densi e massicci ad alto redshift ($z \sim 4-10$), inclusi l'arco "Cosmic Gems", il "Firefly Sparkle" e le galassie Little Red Dot (LRD). Questi sistemi presentano masse stellari comprese tra $10^5$ e $10^7 M_\odot$ e raggi di mezza massa da $\sim 0.4$ a $15$ pc. Sebbene questi ambienti siano teoricamente ideali per la formazione di buchi neri di massa intermedia (IMBH) tramite collisioni stellari a catena, le simulazioni dirette $N$-body sono computazionalmente non fattibili per sistemi così massicci e per grandi popolazioni. Inoltre, l'efficienza della formazione di IMBH in questi specifici ambienti ad alto redshift, in particolare per quanto riguarda l'interazione tra tassi di collisione e perdita di massa tramite venti stellari, richiede una stima quantitativa per valutarne la fattibilità come semi per i buchi neri supermassicci (SMBH) osservati nell'Universo primordiale.

Metodologia
Per affrontare l'impraticabilità computazionale delle simulazioni dirette $N$-body per questi ammassi massicci, gli autori impiegano un quadro ibrido analitico e semi-analitico:

1. Evoluzione Dinamica (Fokker-Planck): Gli autori utilizzano il solver PhaseFlow (parte della libreria Agama) per risolvere l'equazione di Fokker-Planck mediata sulle orbite unidimensionale. Il modello assume ammassi composti da stelle di massa uguale distribuite secondo un profilo di densità di Plummer. L'evoluzione è tracciata su una scala temporale caratteristica di 100 Myr. Pur riconoscendo che gli ammassi reali evolvono più rapidamente a causa della segregazione di massa (che non è modellata esplicitamente), questo approccio fornisce una stima conservativa del collasso del nucleo e dell'evoluzione della densità centrale.
2. Modello di Collisione a Catena: La crescita di una Stella Molto Massiccia (VMS) centrale è modellata utilizzando un quadro analitico basato su Portegies Zwart & McMillan (2002), esteso da Fujii et al. (2024) e Pacucci et al. (2025) per includere la perdita di massa tramite venti. Il tasso di accrescimento di massa ($\dot{M}_{acc}$) è derivato dal tasso di collisione nel nucleo dell'ammasso, dipendente dalla massa del nucleo ($M_{core}$), dal tempo di rilassamento del nucleo ($t_{relax, core}$) e dal logaritmo di Coulomb.
3. Bilancio della Perdita di Massa: La crescita netta della VMS è determinata bilanciando l'acquisto di massa tramite collisioni contro la perdita di massa tramite venti stellari. Il tasso di perdita di massa tramite venti ($\dot{M}_{wind}$) è calcolato utilizzando la prescrizione di Vink (2018), che dipende fortemente dalla metallicità ($Z$).
4. Validazione: Il quadro analitico è validato contro simulazioni dirette $N$-body (Dragon II, Vergara et al. 2025b, e ammassi FROST) dove tali confronti sono possibili. Gli autori riscontrano che il loro modello riproduce l'ordine di grandezza delle masse degli IMBH entro un fattore di 2–5, coerente con le incertezze attese dei modelli semplificati.

Risultati Chiave
L'applicazione di questo quadro a un campione di ammassi osservati da JWST produce le seguenti scoperte:

• Intervallo di Massa degli IMBH: Le masse finali stimate degli IMBH (assumendo il collasso diretto della VMS) tipicamente variano da $\sim 10^2 M_\odot$ a $\sim 4 \times 10^3 M_\odot$.
• Risultati Specifici per Ammasso:
  • Cosmic Gems (Adamo et al. 2024): Questi ammassi sono estremamente compatti ($R_h \sim 1$ pc) con alte masse stellari ($\sim 10^6 M_\odot$) e bassa metallicità ($Z \approx 0.005 Z_\odot$). Il modello prevede la formazione dei semi più massicci nel campione, con un intervallo da $1.7 \times 10^3$ a $2.7 \times 10^3 M_\odot$.
  • Firefly Sparkle (Mowla 2024): Questi ammassi sono più estesi ($R_h \sim 4-6$ pc) con densità inferiori. Le masse IMBH previste sono più basse, con un intervallo da $\sim 1.8 \times 10^2$ a $8.1 \times 10^2 M_\odot$.
  • Sunburst/Frontier Fields (Vanzella et al.): Nonostante i grandi raggi, le alte masse totali ($>10^7 M_\odot$) permettono nuclei massicci, producendo IMBH fino a $\sim 2.1 \times 10^3 M_\odot$.
• Efficienza di Formazione: L'efficienza di formazione (definita come $M_{BH}/M_{cl}$) risulta essere dell'ordine di qualche percento. Questo si allinea con gli adattamenti derivati dalle simulazioni dirette $N$-body di Vergara et al. (2025a).
• Dipendenza dalla Metallicità: La bassa metallicità ($Z \lesssim 0.02 Z_\odot$) è identificata come un fattore critico. Metallicità più elevate aumentano significativamente la perdita di massa tramite venti, sopprimendo la massa finale del seme. Ad esempio, gli ammassi modellati con $Z = 0.1 Z_\odot$ (Messa et al. 2025) mostrano efficienze di formazione inferiori rispetto a controparti povere di metalli di massa simile.

Significato e Affermazioni
Il lavoro afferma che gli ammassi giovani, massicci e poveri di metalli rilevati da JWST sono fabbriche valide per "semi pesanti" di buchi neri di massa intermedia. Le masse stimate ($10^2 - 4 \times 10^3 M_\odot$) sono sufficienti a servire come semi in grado di crescere fino agli SMBH osservati ad alto redshift ($z \sim 6-10$) entro il tempo cosmico disponibile, alleviando i vincoli temporali associati a semi più leggeri (ad esempio, i resti di Popolazione III).

Nello specifico, gli autori suggeriscono che l'estrema compattezza degli ammassi Cosmic Gems facilita la formazione di semi con masse di $1600-2700 M_\odot$. Questi semi pesanti richiederebbero significativamente meno tempi di raddoppio (e-folding times) per raggiungere $10^7-10^8 M_\odot$ entro $z \sim 6$ rispetto a semi più leggeri. Inoltre, gli autori notano che questi semi pesanti potrebbero svolgere un ruolo rilevante nell'emissione di onde gravitazionali da fusioni di IMBH, potenzialmente rilevabili da futuri osservatori come LISA e l'Einstein Telescope.

Lo studio sottolinea che, sebbene il modello sia semplificato e si basi su assunzioni (come stelle di massa uguale e trascuramento di binarie primordiali), i risultati sono coerenti con simulazioni numeriche più complesse e forniscono una stima robusta dell'ordine di grandezza per la formazione di IMBH nell'Universo primordiale. Gli autori riconoscono le incertezze riguardanti la perdita di massa durante le collisioni e il potenziale di funzioni di massa iniziale pesate verso l'alto, che potrebbero ulteriormente potenziare la formazione di IMBH ma non sono state modellate esplicitamente.