El método RAMCOAL se ha ampliado para modelar directamente la dinámica por debajo de la escala de la cuadrícula de los sistemas triples de agujeros negros supermasivos en simulaciones hidrodinámicas de galaxias. Los agujeros negros pasan por una fase de fricción dinámica, forman un par ligado que se compacta mediante la dispersión estelar, los momentos gaseosos y la vinculación con el disco circunuclear; cuando el triple jerárquico se vuelve caótico, la interacción se contrasta con una biblioteca de resultados de experimentos numéricos de N-cuerpos, lo que permite tener en cuenta fusiones, intercambios y eyecciones. Las pruebas en galaxias aisladas con diferentes geometrías han mostrado que la geometría de acercamiento altera qué par acaba fusionándose y al cabo de cuánto tiempo. Por primera vez se demuestra la evolución dinámica completa de un triple de agujeros negros hasta la fusión dentro de una simulación hidrodinámica viva. Esto proporciona un método integral para predecir coalescencias de agujeros negros causadas por interacciones triples, vinculado de forma autoconsistente con la evolución de la galaxia progenitora. Como resultado, cada evento de fusión conserva la historia de su entorno, lo que abre el camino a catálogos de fusiones que conectan directamente las firmas de ondas gravitacionales con las galaxias donde se forman.
Los agujeros negros masivos (ANM) en los centros de las galaxias inevitablemente forman pares y tripletes cuando sus hogares estelares colisionan. Sus fusiones son una fuente clave de ondas gravitacionales de baja frecuencia, tras las cuales cazan las redes de púlsares (PTA), heredando las ideas del descubrimiento de Jocelyn Bell Burnell. Sin embargo, el camino hacia la colisión final está lleno de obstáculos: después de acercarse a escalas de parsec, un par de agujeros negros puede «atascarse» durante miles de millones de años — este es el famoso problema del último parsec. La situación es aún más compleja cuando en el núcleo de una galaxia se encuentran tres ANM. Su interacción caótica puede tanto acelerar como cambiar por completo el desenlace. El nuevo trabajo incorpora estos escenarios triples directamente en las simulaciones hidrodinámicas, utilizando los pioneros cálculos de fricción dinámica de Subrahmanyan Chandrasekhar y la métrica de Karl Schwarzschild.
Los científicos modificaron el código RAMCOAL, que funciona en el entorno de malla adaptativa RAMSES. En él, los agujeros negros se representan como «partículas sumidero», y su evolución por debajo de la resolución de la malla (menos de ~0,4 kpc) se describe mediante recetas semianalíticas. Para los sistemas binarios se establecieron tres etapas: fricción dinámica sobre el fondo de estrellas y gas, ligadura gravitacional y, finalmente, endurecimiento de la órbita por dispersión estelar, interacción con el disco de gas y emisión de ondas gravitacionales, que se propagan a la velocidad de la luz. Para los sistemas triples se introdujo un clasificador de siete estados que reconoce el momento de transición de un sistema triple jerárquico a un régimen caótico según el criterio de estabilidad orbital. El desenlace de la danza caótica se extrae de una biblioteca construida a partir de cálculos newtonianos directos de miles de configuraciones triples. Todos los aspectos — acreción de gas, crecimiento de espines, retroceso por fusión — están vinculados al estado instantáneo de la galaxia anfitriona.
Las pruebas realizadas en una galaxia aislada revelaron una fuerte dependencia de la geometría orbital. Cuando el tercer agujero negro se acercaba en el plano del disco, el par original se fusionó solo después de 3,51 mil millones de años. Pero cuando el tercer objeto siguió una órbita inclinada, interceptó al agujero primario y la fusión final ocurrió en solo 0,29 mil millones de años, es decir, 12 veces más rápido. En el tercer escenario, el más compacto, tres ANM de igual masa formaron un sistema caótico. Se activó el mecanismo triple y, 0,9 mil millones de años después de la activación del caos, el agujero externo se fusionó con el primario, mientras que el segundo componente del par original permaneció en una órbita amplia. El modelo también mostró que, en un disco de gas común, el agujero secundario a menudo recibe más materia que el primario, lo que modifica la distribución de espines antes de la generación de ondas gravitacionales.
La nueva herramienta permite, por primera vez en simulaciones cosmológicas, diferenciar los destinos de los agujeros negros: fusión rápida, estancamiento por siglos o expulsión de la galaxia tras una interacción triple. Esto influye directamente en las predicciones del fondo estocástico de ondas gravitacionales para PTA y en las fuentes individuales para LISA (lanzamiento previsto en 2035). Mientras tanto, las observaciones del James Webb ya encuentran un exceso de núcleos activos dobles a altos desplazamientos al rojo, y la interpretación correcta de esos datos requiere modelos triples como este. Dado que cada agujero lleva consigo la historia de acreción y dinámica, surge la posibilidad de vincular los parámetros de la señal — masas, espines, desplazamientos — con la composición de la galaxia anfitriona, incluido su halo de materia oscura.
La siguiente etapa será el lanzamiento del proyecto cosmológico a gran escala COSCOAL, donde millones de agujeros negros evolucionarán según las nuevas reglas en tiempo real. Esto permitirá construir catálogos sintéticos de fusiones con retrasos físicamente fundamentados, de importancia crítica para la preparación de las misiones LISA y TianQin. Un campo de pruebas especial serán los «pequeños puntos rojos» — galaxias compactas de la época del universo temprano, descubiertas por el James Webb, donde las interacciones triples pueden ser especialmente intensas.
El método influirá en la astronomía multicanal, al combinar el análisis de ondas gravitacionales con los relevos electromagnéticos de núcleos activos de galaxias. También es importante para comprender el crecimiento de los agujeros negros supermasivos desde el universo temprano hasta nuestros días.
Los planes inmediatos incluyen la calibración del tiempo de interacción caótica y la ampliación de la biblioteca de desenlaces para el caso en que el objeto externo sea más masivo que el par interno.
El trabajo está directamente relacionado con el problema del último parsec: por qué los pares de agujeros negros logran fusionarse dentro del tiempo de Hubble. La inyección de un tercer cuerpo y la consideración del entorno gaseoso ofrecen una solución natural, al tiempo que plantean nuevas preguntas sobre la jerarquía de fusiones y el papel de la materia oscura en la fricción dinámica.
🎯 Los sistemas triples de agujeros negros pueden funcionar como aceleradores naturales de fusiones: en algunas configuraciones, la excentricidad aumenta tanto que las ondas gravitacionales disipan la energía en miles de años en lugar de miles de millones. Es como un billar cósmico con carambolas impredecibles.
🎬 El tema de los sistemas triples y las interacciones caóticas resuena con la trama de la famosa novela de Liu Cixin «El problema de los tres cuerpos», donde la inestabilidad orbital de tres estrellas determina el destino de la civilización trisolariana. En la realidad, nuestros modelos predicen desenlaces igualmente dramáticos: desde fusiones ultrarrápidas hasta la expulsión de un agujero negro al espacio intergaláctico.