Cien mil estrellas en el denso núcleo de un cúmulo globular se mueven como bolas de acero en una gigantesca máquina de pinball. Cada «empujón» —un encuentro gravitacional entre un sistema binario y una estrella solitaria— desencadena una cascada de eventos impredecibles. La pareja puede separarse, el intruso puede expulsar a uno de los compañeros y ocupar su lugar, y a veces los tres quedan atrapados en una enmarañada danza jerárquica. Estos encuentros duran días, pero las consecuencias determinan el destino de millones de astros y engendran objetos exóticos: desde agujeros negros hasta fuentes de ondas gravitacionales.
El clásico problema de los tres cuerpos es caótico: una pequeña desviación en la velocidad inicial y el resultado cambia radicalmente. Henri Poincaré fue el primero en comprenderlo a finales del siglo XIX, demostrando que no existe una solución analítica general. Sin embargo, un grupo de astrofísicos ha demostrado que el caos no es total. Creando casi un millón de simulaciones numéricas con el código REBOUND, entrenaron un modelo de aprendizaje automático XGBoost para adivinar el resultado. Como entrada se usaron cantidades físicas adimensionales: la rigidez del sistema binario, el factor de enfoque y el número de Safronov. La rigidez, cociente entre la energía de enlace del par y la energía cinética de la estrella incidente, resultó ser el director principal: si la pareja es «dura», es casi seguro que sobrevive al impacto o intercambia compañeros.
Los errores no están distribuidos uniformemente. Se agrupan en una estrecha frontera entre el intercambio y los sistemas triples jerárquicos. Aquí, la redistribución de energía y momento es tan sensible que la captura temporal de las estrellas genera una resonancia caótica. El resultado depende de fluctuaciones microscópicas, como una bola que se detiene en la cresta entre dos hoyos. El análisis de la entropía de las cuencas confirmó que esta frontera tiene una estructura fractal. Es justo allí donde el modelo vacila y arroja baja confianza. No es un defecto, sino un reflejo de la impredecibilidad fundamental descubierta por Poincaré.
La ganancia práctica es colosal. Antes, simular un solo encuentro de este tipo requería minutos de tiempo de cómputo. Ahora, el sustituto lo resuelve en milisegundos. Una aceleración de cientos de miles de veces permite incorporar el algoritmo en simulaciones globales de cúmulos, donde antes cada evento se reemplazaba por una aproximación simplificada. Ahora, la dinámica de resonancias y la formación de sistemas triples se reproducirán con precisión física. El siguiente paso es añadir fuerzas de marea, tamaños finitos de las estrellas y correcciones relativistas para predecir fusiones de objetos compactos con un detalle sin precedentes. Esto impactará directamente en los estudios poblacionales de fuentes de ondas gravitacionales.
El pinball estelar gigantesco encuentra sus reglas. El caos, que parecía absoluto, se somete a un orden oculto que un algoritmo entrenado es capaz de leer. Allí donde las matemáticas capitulan, la física deja pistas, y el modelo aprende a interpretarlas. Por lo tanto, el estruendo de las bolas en colisión en la oscuridad de los cúmulos ya no será ruido aleatorio.
En algunas configuraciones, tres estrellas forman un tándem jerárquico temporal que vaga por el cúmulo durante millones de años antes de disolverse.
En la famosa novela de Liu Cixin «El problema de los tres cuerpos», la civilización de Trisolaris sufre las órbitas impredecibles de tres soles. Los astrónomos se enfrentan a un problema similar a escala de cúmulos, y ahora parte de esa impredecibilidad ha sido domeñada por un algoritmo.