Planar Three-Body Problem: theoretical predictions and simulation results

La gravedad es la gran coreógrafa del cosmos, dictando cómo las estrellas, los planetas y los agujeros negros se mueven a través del vasto vacío del espacio. Cuando dos objetos masivos interactúan, su trayectoria es predecible y estable, como un planeta orbitando una estrella en un bucle perfecto y repetitivo. Pero si se introduce un tercer objeto, las reglas cambian por completo. El problema de los tres cuerpos es un desafío fundamental en la física donde tres cuerpos de peso similar se atraen entre sí, creando una danza caótica donde cambios diminutos al inicio conducen a resultados radicalmente diferentes después. Debido a que este caos hace imposible predecir la trayectoria exacta futura de cualquier grupo específico de tres, los científicos han recurrido a la estadística. En lugar de preguntar dónde terminará un trío específico, preguntan qué sucede en promedio a través de millones de escenarios. Este enfoque estadístico ayuda a explicar cómo se forman los planetas en discos giratorios de gas o cómo colisionan los agujeros negros en los centros de las galaxias, eventos que a menudo involucran arreglos casi planos, bidimensionales, en lugar del caos tridimensional completo del espacio profundo.

En un estudio reciente, los investigadores se propusieron probar una teoría estadística específica contra la realidad de estas interacciones planas y bidimensionales. Se centraron en un escenario donde tres cuerpos, partiendo del reposo, caen unos hacia otros bajo la influencia de la gravedad hasta que el sistema se rompe. En casi todos los casos, esta interacción caótica resulta en la formación de un par ligado y estrecho entre dos cuerpos, mientras el tercero es expulsado hacia el espacio. El equipo quería saber si una teoría moderna, que predice los resultados basándose en el flujo de estados posibles a través del tiempo, podría describir con precisión lo que sucede cuando el movimiento está confinado a un solo plano, como ocurre a menudo en los discos reales de materia. Para averiguarlo, realizaron un millón de simulaciones por computadora para ocho combinaciones diferentes de masas, rastreando cada detalle de cómo evolucionaban los sistemas y finalmente se fragmentaban.

Los investigadores comenzaron construyendo un laboratorio digital donde podían controlar las condiciones iniciales con precisión. Colocaron un par de cuerpos en una órbita circular y un tercer cuerpo, más ligero, a una distancia determinada, todos situados de forma plana sobre la misma superficie. Luego dejaron que la gravedad tomara el control, observando cómo el cuerpo distante caía hacia el par, desencadenando un revuelo caótico. La computadora rastreó el sistema hasta que se asentó en un estado final: un par binario y un viajero en fuga. Para asegurar que estaban estudiando la parte verdaderamente caótica del proceso, el equipo filtró las expulsiones raras y rápidas que ocurrían casi inmediatamente, centrándose solo en los sistemas que pasaron suficiente tiempo en un estado caótico y mezclado antes de separarse. Esto les permitió comparar sus resultados simulados directamente con las predicciones teóricas, que asumen que el sistema ha olvidado su punto de partida y ha alcanzado un estado de equilibrio estadístico.

Uno de los hallazgos más sorprendentes concernía la probabilidad de cuál sería el cuerpo que escaparía. La teoría predecía que, en un mundo plano y bidimensional, las probabilidades de escape seguirían un patrón específico basado en las masas involucradas. Sin embargo, las simulaciones revelaron algo inesperado. Los resultados de las simulaciones planas no coincidieron con las nuevas predicciones teóricas para dos dimensiones. En cambio, coincidieron con las predicciones para un mundo tridimensional con una exactitud casi perfecta, con una diferencia de solo alrededor del uno por ciento. Esto sugiere que, incluso cuando el movimiento está físicamente restringido a un plano, el comportamiento estadístico del sistema imita el caso tridimensional más complejo. Los investigadores señalaron que esto era una desviación significativa de lo que la teoría actual sugería para sistemas planares, insinuando que la teoría podría necesitar ajustes para dar cuenta de cómo se comparten la energía y el momento en estos entornos planos.

El estudio también examinó los detalles de la propia fuga, particularmente los momentos en que el tercer cuerpo apenas logra liberarse. En estas escapes "marginales", la teoría predijo una distribución específica para el momento angular del cuerpo que escapa. Las simulcciones confirmaron esta predicción con alta precisión, mostrando que los datos se alineaban estrechamente con la expectativa matemática. Del mismo modo, los investigadores observaron la forma de las órbitas dejadas por el par restante. Encontraron que la distribución de qué tan estiradas o circulares eran estas nuevas órbitas coincidía notablemente bien con el modelo teórico, a pesar de que el modelo dependía de una suposición que los investigadores no pudieron probar completamente para este caso específico. Este acuerdo sugiere que las ideas centrales de la teoría son robustas, incluso si algunos de los detalles más finos respecto a la "emisividad" del sistema siguen siendo un misterio.

Otra área clave de investigación fue cuánto duran estas interacciones caóticas. Los investigadores encontraron que la vida útil del sistema seguía un patrón predecible: para interacciones muy largas, el número de sistemas supervivientes disminuía según una ley de potencia específica, un resultado que coincidía con las expectativas teóricas derivadas de las leyes del movimiento orbital. Sin embargo, las etapas tempranas de la interacción mostraron una mezcla de comportamientos, con algunos sistemas rompiéndose rápidamente y otros permaneciendo en un estado de excursiones de sub-escape, donde los cuerpos se alejan pero son atraídos de nuevo antes de finalmente escapar. Esta mezcla de comportamientos creó una distribución compleja de vidas útiles que las simulaciones capturaron en detalle, proporcionando una imagen clara de cómo el caos y la estabilidad temporal coexisten en estos sistemas de tres cuerpos.

Finalmente, el equipo investigó si la dirección de la fuga importaba. En un sistema plano, el par restante puede orbitar en la misma dirección que el espín total del sistema (progrado) o en la dirección opuesta (retrógrado). La teoría predecía que estos dos resultados deberían ser igualmente probables, independientemente de las masas involucradas. Las simulaciones, sin embargo, contaron una historia diferente. Los investigadores observaron que, a medida que aumentaba la diferencia de masa entre los cuerpos, la probabilidad de una fuga prograda crecía significativamente por encima de una retrógrada. Este sesgo sistemático contradecía la predicción teórica de probabilidad igual, sugiriendo que la distribución de masa juega un papel más crítico en la determinación de la dirección de la órbita final de lo que el modelo actual contempla.

Los autores concluyen que, si bien su teoría estadística predijo con éxito muchos aspectos del problema de los tres cuerpos, como la distribución de las formas orbitales y el tiempo de las interacciones de larga duración, falló al predecir con exactitud las probabilidades de escape y la dirección de la órbita final en el caso planar. El hecho de que las simulaciones planas coincidieran tan estrechamente con las predicciones tridimensionales fue una gran sorpresa, indicando que las restricciones de un plano no simplifican el resultado estadístico tanto como se pensaba. Los investigadores sugieren que la ruptura de la aproximación "ciega a la emisividad", que asume que el sistema olvida su pasado rápidamente, podría ser la causa de estas discrepancias. Planean refinar su teoría para dar cuenta mejor de estos efectos, con el objetivo de crear una descripción estadística más completa que pueda aplicarse a entornos astrofísicos del mundo real, desde la formación de planetas en el pasado de nuestro propio sistema solar hasta las violentas colisiones de agujeros negros en los centros de galaxias distantes.