AGN-driven BBH mergers: Black hole populations and hierarchical growth across the AGN parameter space

Resumen Técnico: Fusiones de BBH impulsadas por AGN: Poblaciones de agujeros negros y crecimiento jerárquico a través del espacio de parámetros de los AGN

Problema y Motivación
Los discos de núcleos galácticos activos (AGN) han emergido como un canal dinámico prometedor para la formación y fusión de sistemas binarios de agujeros negros (BBH). Si bien los agujeros negros (BH) de masa estelar incrustados en estos discos pueden migrar, emparejarse y fusionarse mediante torques de gas y emisión de ondas gravitacionales (GW), las propiedades predichas de esta población siguen siendo altamente inciertas. Estudios previos han destacado una fuerte dependencia de ingredientes físicos como la estructura del disco, los torques de migración, la eficiencia del endurecimiento por gas y la vida útil de los AGN, adoptando a menudo diferentes suposiciones que dificultan la comparación directa. Las restricciones observacionales recientes sugieren que el canal de los AGN no puede explicar la totalidad de la población de fusiones BBH detectada por los observatorios de ondas gravitacionales, pero sus contribuciones específicas y firmas distintivas requieren una exploración sistemática del espacio de parámetros de los AGN para identificar características observacionales robustas y fuentes de incertidumbre.

Metodología
Los autores presentan un marco semi-analítico actualizado implementado dentro del entorno de software fastcluster para modelar la formación, evolución y crecimiento jerárquico de los BBH en los discos de AGN. El enfoque rastrea el ciclo de vida de los agujeros negros de origen estelar (primera generación, 1g) a través de la captura, migración radial, formación de binarias, endurecimiento impulsado por gas, perturbaciones dinámicas y la eventual fusión.

Componentes metodológicos clave incluyen:

• Modelado del Disco: El disco de acreción se modela utilizando el modelo analítico de estado estacionario de Sirko & Goodman (2003) (SG). Los autores resuelven de forma autoconsistente las ecuaciones de disco unidimensional para variaciones en las masas de los agujeros negros supermasivos ($M_\bullet$), razones de Eddington ($\lambda_\bullet$) y parámetros de viscosidad del disco ($\alpha$). Se prueban dos valores de viscosidad: $\alpha = 0.01$ (baja viscosidad) y $\alpha = 0.1$ (alta viscosidad).
• Condiciones Iniciales: Las simulaciones inicializan las poblaciones de BH basándose en las relaciones de escala de los cúmulos estelares nucleares (NSC). Los BH de 1g se extraen de un catálogo de síntesis de población (sevn) asumiendo metalicidad solar. Los espines se asignan basándose en una distribución de Maxwell, con orientaciones alineadas al momento angular del disco.
• Física Evolutiva:
  • Migración: Los BH migran a través de los regímenes de Tipo I y Tipo II, con prescripciones de Grishin et al. (2024) y Kanagawa et al. (2018).
  • Emparejamiento: La formación de binarias sigue un criterio asistido por gas, considerando la fracción de órbitas prógradas ($f_{progr}$).
  • Endurecimiento y Encuentros: Tras la formación, las binarias evolucionan bajo torques de gas combinados (interacción con el disco circumbinario) y encuentros de tres cuerpos con BH individuales. Los autores integran las ecuaciones de evolución orbital combinando el endurecimiento impulsado por gas (usando Ishibashi & Gröbner 2024) y la emisión de GW (Peters 1964). Los encuentros se modelan utilizando una rejilla precomputada de experimentos de dispersión post-newtonianos (tsunami).
  • Crecimiento Jerárquico: Los remanentes de fusión se calculan utilizando prescripciones de espín alineado y velocidades de retroceso (kick). Los remanentes que permanecen ligados al disco participan en generaciones de fusiones subsecuentes (generación $N$), sujetos a límites de vida útil del AGN y criterios de eyección.
• Síntesis de Población: Para construir una población intrínseca de BBH representativa del Universo local ($z \lesssim 0.3$), los autores ponderan los resultados de la simulación de acuerdo con las distribuciones observadas de AGN de $M_\bullet$ y $\lambda_\bullet$ (Greene & Ho 2007).

Resultados Clave
El estudio explora sistemáticamente la dependencia de las poblaciones de fusión de BBH con respecto a $M_\bullet$, $\lambda_\bullet$ y $\alpha$, arrojando los siguientes hallazgos:

1. Crecimiento Jerárquico y Distribución de Masas:

   • Los discos de AGN pueden producir fusiones repetidas, construyendo una cola de alta masa que se extiende más allá de la brecha de la supernova de inestabilidad de par (PISN) ($50 \lesssim m_1/M_\odot \lesssim 130$) hacia el rango de los agujeros negros de masa intermedia (IMBH).
   • Dependencia de la Viscosidad: El crecimiento jerárquico es significativamente más eficiente en discos de baja viscosidad ($\alpha = 0.01$), donde las densidades superficiales de gas más altas acortan los tiempos de escala dinámica, permitiendo un mayor número de generaciones de fusión ($N_{gen} > 100$). Los discos de alta viscosidad ($\alpha = 0.1$) suprimen la formación de remanentes masivos, saturándose típicamente en generaciones más bajas ($N_{gen} \approx 20$).
   • Razones de Masa: La población evoluciona hacia razones de masa cada vez más desiguales en masas primarias altas. Mientras que las masas primarias crecen mediante fusiones sucesivas, las masas secundarias se saturan en valores más bajos (típicamente por debajo de la brecha PISN), ya que las primarias de alta generación se fusionan con secundarias de baja generación.

2. Eficiencia de Fusión y Espacio de Parámetros:

   • La eficiencia de fusión ($\eta$) generalmente aumenta con la razón de Eddington ($\lambda_\bullet$).
   • La dependencia de $M_\bullet$ es no trivial; la eficiencia y la masa máxima del remanente tienden a alcanzar su máximo en masas de SMBH intermedias a grandes, pero ningún comportamiento característico único abarca todo el espacio de parámetros. Algunas configuraciones de alta masa no muestran fusiones, lo que indica que una gran $M_\bullet$ por sí sola no garantiza un ensamblaje eficiente.

3. Firmas de Espín y Excentricidad:

   • Espín Efectivo ($\chi_{eff}$): La distribución de $\chi_{eff}$ depende fuertmente de la fracción prógrada inicial ($f_{progr}$). Las configuraciones prógradas producen $\chi_{eff}$ positivos, mientras que las retrógradas producen valores negativos. Las fracciones intermedias producen una distribución bimodal. Existe una correlación entre la masa primaria y $|\chi_{eff}|$, ya que las fusiones repetidas refuerzan las orientaciones de espín coherentes.
   • Excentricidad: La mayoría de las fusiones entran en la banda de los detectores terrestres (a $f_{peak} = 10$ Hz) con excentricidades bajas ($e \sim 10^{-3}–10^{-2}$). Sin embargo, una subpoblación retiene una excentricidad alta ($e \gtrsim 0.8$), habiéndose fusionado antes de que ocurriera una circularización significativa.

4. Comparación con Eventos de GW:

   • El canal de los AGN puede reproducir sistemas consistentes con los eventos masivos GW190521 y GW231123, aunque estos ocupan las colas de las distribuciones predichas.
   • Los sistemas que coinciden con los parámetros de masa de GW190521 ocurren a una tasa de $\sim 2 \times 10^{-3}$ en modelos de baja viscosidad, mientras que los sistemas similares a GW231123 son más raros ($\sim 3 \times 10^{-4}$).
   • El modelo tiende a subpredecir la masa secundaria ($m_2$) para estos eventos específicos en comparación con los parámetros inferidos, lo que sugiere que, si bien el canal de los AGN puede producir la masa primaria y las características de espín, igualar la distribución completa de masas de los componentes es un desafío.

5. Robustez del Modelo:

   • Las propiedades generales de la población (cola de alta masa, correlaciones de espín) son robustas ante variaciones en los formalismos de torque de migración (por ejemplo, correcciones térmicas) y las prescripciones de endurecimiento por gas (IG20 vs. IG24 vs. C24).
   • Dependencia Crítica: La formación de remanentes masivos y la cola de alta masa se suprimen completamente cuando se desactiva el endurecimiento por gas, confirmando que la evolución de las binarias impulsada por el gas es un ingrediente clave para el crecimiento jerárquico en los discos de AGN.

Significado y Reivindicaciones
El artículo afirma proporcionar un marco exhaustivo y actualizado que mapea sistemáticamente el espacio de parámetros de los AGN para evaluar la robustez de las predicciones de fusión de BBH. Los autores sostienen que el canal de los AGN no produce una población característica única, sino una amplia diversidad de resultados. El estudio identifica firmas distintivas a nivel de población —específicamente la correlación entre la masa primaria y $|\chi_{eff}|$, la tendencia hacia razones de masa desiguales en masas altas, y la retención de memoria respecto a las configuraciones orbitales iniciales en la distribución de espín.

El trabajo enfatiza que, si bien el canal de los AGN puede explicar eventos de BBH masivos como GW190521 y GW231123, estos eventos probablemente representan ocurrencias raras dentro del canal. Los autores concluyen que la evolución de las binarias impulsada por el gas es esencial para un crecimiento jerárquico eficiente, y que las suposiciones físicas específicas sobre la estructura del disco y el endurecimiento influyen significámente en la eficiencia de la formación de remanentes masivos, aunque la existencia cualitativa de una cola de alta masa sigue siendo una característica robusta del canal.