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

Resumen Técnico: Formación de Agujeros Negros de Masa Intermedia en Cúmulos Estelares Jóvenes y Masivos Detectados por el JWST

Enunciado del Problema
El Telescopio Espacial James Webb (JWST) ha revelado una población de sistemas estelares densos y masivos a alto corrimiento al rojo ($z \sim 4-10$), incluyendo el arco "Cosmic Gems", el "Firefly Sparkle" y las galaxias Little Red Dot (LRD). Estos sistemas exhiben masas estelares que oscilan entre $10^5$ y $10^7 M_\odot$ y radios de media masa desde $\sim 0.4$ hasta $15$ pc. Si bien estos entornos son teóricamente ideales para la formación de agujeros negros de masa intermedia (IMBH) mediante colisiones estelares descontroladas, las simulaciones directas de $N$-cuerpos son computacionalmente inviables para sistemas tan masivos y grandes poblaciones. Además, la eficiencia de la formación de IMBH en estos entornos específicos de alto corrimiento al rojo, particularmente en lo que respecta a la interacción entre las tasas de colisión y la pérdida de masa por vientos estelares, requiere una estimación cuantitativa para evaluar su viabilidad como semillas para los agujeros negros supermasivos (SMBH) observados en el Universo temprano.

Metodología
Para abordar la inviabilidad computacional de las simulaciones directas de $N$-cuerpos para estos cúmulos masivos, los autores emplean un marco híbrido analítico y semi-analítico:

1. Evolución Dinámica (Fokker-Planck): Los autores utilizan el solucionador PhaseFlow (parte de la biblioteca Agama) para resolver la ecuación de Fokker-Planck promediada en órbitas unidimensional. El modelo asume cúmulos compuestos por estrellas de masa igual distribuidas según un perfil de densidad de Plummer. La evolución se rastrea durante una escala de tiempo característica de 100 Myr. Si bien se reconoce que los cúmulos reales evolucionan más rápido debido a la segregación de masas (que no se modela explícitamente), este enfoque proporciona una estimación conservadora del colapso del núcleo y la evolución de la densidad central.
2. Modelo de Colisión Descontrolada: El crecimiento de una Estrella Muy Masiva (VMS) central se modela utilizando un marco analítico basado en Portegies Zwart & McMillan (2002), extendido por Fujii et al. (2024) y Pacucci et al. (2025) para incluir la pérdida de masa por vientos. La tasa de acreción de masa ($\dot{M}_{acc}$) se deriva de la tasa de colisión en el núcleo del cúmulo, dependiendo de la masa del núcleo ($M_{core}$), el tiempo de relajación del núcleo ($t_{relax, core}$) y el logaritmo de Coulomb.
3. Balance de Pérdida de Masa: El crecimiento neto de la VMS se determina equilibrando la ganancia de masa mediante colisiones contra la pérdida de masa mediante vientos estelares. La tasa de pérdida de masa por vientos ($\dot{M}_{wind}$) se calcula utilizando la prescripción de Vink (2018), que depende fuertemente de la metalicidad ($Z$).
4. Validación: El marco analítico se valida frente a simulaciones directas de $N$-cuerpos (Dragon II, Vergara et al. 2025b, y cúmulos FROST) donde tales comparaciones son posibles. Los autores encuentran que su modelo reproduce el orden de magnitud de las masas de IMBH dentro de un factor de 2–5, consistente con las incertidumbres esperadas de los modelos simplificados.

Resultados Clave
La aplicación de este marco a una muestra de cúmulos observados por el JWST arroja los siguientes hallazgos:

• Rango de Masa de IMBH: Las masas finales estimadas de IMBH (asumiendo el colapso directo de la VMS) típicamente oscilan entre $\sim 10^2 M_\odot$ y $\sim 4 \times 10^3 M_\odot$.
• Resultados Específicos por Cúmulo:
  • Cosmic Gems (Adamo et al. 2024): Estos cúmulos son extremadamente compactos ($R_h \sim 1$ pc) con altas masas estelares ($\sim 10^6 M_\odot$) y baja metalicidad ($Z \approx 0.005 Z_\odot$). El modelo predice la formación de las semillas más masivas de la muestra, con rangos de $1.7 \times 10^3$ a $2.7 \times 10^3 M_\odot$.
  • Firefly Sparkle (Mowla 2024): Estos cúmulos son más extendidos ($R_h \sim 4-6$ pc) con densidades más bajas. Las masas de IMBH predichas son menores, oscilando entre $\sim 1.8 \times 10^2$ y $8.1 \times 10^2 M_\odot$.
  • Sunburst/Campos Fronterizos (Vanzella et al.): A pesar de sus grandes radios, las altas masas totales ($>10^7 M_\odot$) permiten núcleos masivos, produciendo IMBH de hasta $\sim 2.1 \times 10^3 M_\odot$.
• Eficiencia de Formación: La eficiencia de formación (definida como $M_{BH}/M_{cl}$) se encuentra en el orden de unos pocos por ciento. Esto se alinea con los ajustes derivados de simulaciones directas de $N$-cuerpos por Vergara et al. (2025a).
• Dependencia de la Metalicidad: Se identifica la baja metalicidad ($Z \lesssim 0.02 Z_\odot$) como un factor crítico. Las metalicidades más altas aumentan significativamente la pérdida de masa por vientos, suprimiendo la masa final de la semilla. Por ejemplo, los cúmulos modelados con $Z = 0.1 Z_\odot$ (Messa et al. 2025) muestran eficiencias de formación más bajas en comparación con contrapartes pobres en metales de masa similar.

Significado y Afirmaciones
El artículo afirma que los cúmulos estelares jóvenes, densos y pobres en metales detectados por el JWST son fábricas viables para "semillas pesadas" de agujeros negros de masa intermedia. Las masas estimadas ($10^2 - 4 \times 10^3 M_\odot$) son suficientes para servir como semillas que pueden crecer hasta convertirse en los SMBH observados a alto corrimiento al rojo ($z \sim 6-10$) dentro del tiempo cósmico disponible, aliviando las restricciones de tiempo asociadas con semillas más ligeras (por ejemplo, remanentes de la Población III).

Específicamente, los autores sugieren que la extrema compacidad de los cúmulos Cosmic Gems facilita la formación de semillas con masas de $1600-2700 M_\odot$. Estas semillas pesadas requerirían significativamente menos tiempos de e-folding para alcanzar $10^7-10^8 M_\odot$ para $z \sim 6$ en comparación con semillas más ligeras. Además, los autores señalan que estas semillas pesadas podrían desempeñar un papel relevante en la emisión de ondas gravitacionales provenientes de fusiones de IMBH, potencialmente detectables por futuros observatorios como LISA y el Telescopio Einstein.

El estudio enfatiza que, aunque el modelo es simplificado y se basa en suposiciones (como estrellas de masa igual y la omisión de binarias primordiales), los resultados son consistentes con simulaciones numéricas más complejas y proporcionan una estimación robusta del orden de magnitud para la formación de IMBH en el Universo temprano. Los autores reconocen incertidumbres respecto a la pérdida de masa durante las colisiones y el potencial de funciones de masa inicial con exceso de masas altas, lo cual podría mejorar aún más la formación de IMBH pero no fue modelado explícitamente.