The SEEDZ Simulations: Methodology and First Results on Massive Black Hole Seeding and Early Galaxy Growth

1. El Problema Científico

El campo de la astrofísica de alto desplazamiento al rojo (high-z) está experimentando una transformación debido a las observaciones del Telescopio Espacial James Webb (JWST). Se han detectado poblaciones de agujeros negros masivos (MBHs) con masas de hasta $10^6 - 10^8 M_\odot$ en los primeros mil millones de años del universo ($z > 6$). Estas detecciones presentan dos desafíos principales para los modelos teóricos actuales:

• Crecimiento Rápido: Explicar cómo los agujeros negros crecieron tan rápidamente desde sus semillas iniciales hasta alcanzar masas tan grandes en tan poco tiempo.
• Densidad y Entorno: La alta densidad numérica de estos MBHs y su presencia en galaxias extremadamente compactas (como los sistemas "Little Red Dots" o LRD), donde la masa del agujero negro puede ser comparable a la masa estelar total de la galaxia anfitriona.

Los modelos existentes a menudo luchan por reconciliar las vías de formación de "semillas ligeras" (restos de estrellas de Población III) y "semillas pesadas" (colapso directo de estrellas supermasivas o cúmulos estelares) con las observaciones actuales, especialmente considerando los efectos de retroalimentación (feedback) y la metalicidad.

2. Metodología: Las Simulaciones SEEDZ

El equipo ha desarrollado la suite de simulaciones SEEDZ, un conjunto de simulaciones hidrodinámicas cosmológicas de gran escala diseñadas para explorar la formación y crecimiento de los primeros agujeros negros.

• Código Numérico: Se utiliza AREPO2, un código de malla móvil de última generación que resuelve las ecuaciones de Euler en una malla Voronoi no estructurada.
• Resolución y Configuración:
  • Se inicia con una caja "padre" de 40 $h^{-1}$ Mpc (solo materia oscura) hasta $z=10$.
  • Se seleccionan tres regiones de interés para simulaciones "zoom-in" con hidrodinámica completa: Rarepeak (centrada en el halo más masivo, $M \sim 1.8 \times 10^{11} M_\odot$) y dos regiones Normales (con halos de $\sim 10^{10} M_\odot$).
  • Resolución espacial mínima: $\sim 6$ pc físicos en $z=15$.
  • Resolución de masa de partículas de materia oscura: $\sim 6.9 \times 10^4 M_\odot h^{-1}$.
• Modelado de Estrellas y Agujeros Negros (SmartStar):
  • Se implementa un esquema de partículas "sink" (sumidero) llamado SmartStar.
  • Semillas Ligeras: Se forman a partir del final de la vida de estrellas de Población III (PopIII). Si la tasa de acreción es baja y la metalicidad es alta, se forman cúmulos estelares de Población II (PopII).
  • Semillas Pesadas: Se forman estocásticamente cuando se cumple una condición crítica de tasa de flujo de masa inminente ($> 1 M_\odot \text{yr}^{-1}$) en una celda gravitacionalmente inestable. Este método es agnóstico respecto al mecanismo subyacente (colapso de estrellas supermasivas o cúmulos densos) y no depende de umbrales de metalicidad o fracción de $H_2$.
  • Rango de Masas Iniciales: Las semillas pesadas se generan entre $5 \times 10^3$ y $10^5 M_\odot$.
• Física Subgrid:
  • Química y Enfriamiento: Solucionador de química no equilibrada que rastrea H, D, He y $H_2$, incluyendo el fondo de radiación Lyman-Werner (LW).
  • Explosiones de Supernova (SN): Modelado de retroalimentación mediante inyección de momento y enriquecimiento de metales. Las estrellas PopIII de diferentes masas evolucionan a agujeros negros ligeros, estallan o colapsan directamente.
  • Acreción y Feedback: Se utiliza el modelo de acreción de Bondi-Hoyle-Lyttleton ajustado por vorticidad. Se incluye retroalimentación térmica tanto para tasas sub-Eddington como super-Eddington (usando parametrizaciones de Madau et al. 2014).
  • Fricción Dinámica: Se implementa un modelo de fricción dinámica para corregir la falta de resolución en la interacción de agujeros negros con el fondo de materia oscura, permitiendo que los agujeros negros migren hacia el centro de los halos.

3. Contribuciones Clave

1. Marco Unificado: SEEDZ integra en un solo marco coherente la formación de semillas ligeras (PopIII) y pesadas (vía inestabilidad gravitacional y alto flujo de masa), sin depender de prescripciones fijas de metalicidad para las semillas pesadas.
2. Modelado de Crecimiento Super-Eddington: La simulación permite tasas de acreción super-Eddington sostenidas durante periodos de 1 a 15 millones de años, crucial para explicar el crecimiento rápido observado.
3. Comparación con Observaciones JWST: El estudio está diseñado específicamente para conectar los mecanismos de formación temprana con las observaciones de alto desplazamiento al rojo ($z \sim 10-15$) y futuras comparaciones con $z=10$.

4. Resultados Principales (hasta $z=15$)

• Crecimiento de Agujeros Negros vs. Galaxias:
  • Los agujeros negros crecen inicialmente más rápido que sus galaxias anfitrionas. Esto resulta en agujeros negros sobremasivos en relación con las relaciones de escalado locales ($z=0$).
  • En el extremo superior de masas, las relaciones $M_{BH} - M_*$ simuladas coinciden con las tendencias observadas en alto desplazamiento al rojo ($M_{BH}/M_* \sim 10^{-2}$), sugiriendo que las relaciones fundamentales observadas hoy ($z=0$) aún no se han establecido en estas galaxias jóvenes.
• Dinámica de Crecimiento:
  • El crecimiento es dominado por la acreción ($>99\%$), con una contribución mínima de fusiones ($<1\%$) debido a limitaciones de resolución en la dinámica de binarias de agujeros negros.
  • Las semillas pesadas crecen desde masas iniciales de $\sim 10^4 M_\odot$ hasta $\sim 10^6 M_\odot$ (un aumento de más de dos órdenes de magnitud) en un periodo corto.
  • La acreción super-Eddington ocurre durante los primeros 1-15 Myr tras la formación, antes de que la retroalimentación térmica frene el crecimiento.
• Entorno de Formación:
  • Aunque los criterios de formación de semillas pesadas son independientes de la metalicidad, la mayoría de los agujeros negros masivos se forman en entornos de baja metalicidad ($Z \ll 10^{-4} Z_\odot$).
  • Algunos de los agujeros negros más masivos permanecen en entornos de baja metalicidad incluso a $z=15$, lo que tiene paralelismos con observaciones recientes de SMBHs en galaxias casi prístinas.
• Densidad Numérica:
  • Se produce una densidad numérica de agujeros negros con $M > 10^5 M_\odot$ del orden de $10^{-2} \text{cMpc}^{-3} h^3$ a $z=15$. Solo una pequeña fracción ($\sim 0.2\%$) necesita continuar creciendo para ser consistente con las observaciones de LRDs y quasares.

5. Significado e Implicaciones

• Resolución de la Tensión Observacional: Los resultados sugieren que la aparente discrepancia entre las masas de los agujeros negros y sus galaxias anfitrionas en alto desplazamiento al rojo es un fenómeno transitorio. Las relaciones de escalado actuales ($M_{BH} - M_*$) son el resultado de una evolución posterior en galaxias maduras, no una condición inicial.
• Viabilidad de Semillas Pesadas: La capacidad de las semillas pesadas para crecer rápidamente mediante acreción super-Eddington en entornos de baja metalicidad ofrece una vía viable para explicar la existencia de MBHs masivos en la época de reionización, sin necesidad de mecanismos de crecimiento extremadamente eficientes a largo plazo.
• Futuro del Proyecto: La suite SEEDZ continuará evolucionando hasta $z=10$, momento en el cual se realizará una comparación exhaustiva con los datos del JWST. Se planean futuras iteraciones con mayor resolución para capturar adecuadamente el crecimiento de semillas ligeras y la dinámica de fusiones de agujeros negros.

En resumen, las simulaciones SEEDZ proporcionan un marco físico autoconsistente que demuestra cómo la combinación de semillas pesadas y acreción super-Eddington puede explicar la población de agujeros negros masivos observada en el universo temprano, desafiando la idea de que las relaciones galaxia-agujero negro son universales desde el principio.