No confunda la falta de solución con el caos absoluto
En la novela de ciencia ficción El problema de los tres cuerpos de Liu Cixin y en el imaginario colectivo, el problema de los tres cuerpos suele considerarse el arquetipo de lo intrínsecamente impredecible. Sin embargo, se trata de un malentendido conceptual que se remonta a más de un siglo. En la década de 1880, Henri Poincaré demostró que no existe una solución analítica general en términos de integrales primeras algebraicas para este sistema, un hallazgo que sentó las bases de la teoría del caos moderna. No obstante, que no exista una fórmula cerrada universal no implica en absoluto que dentro de un sistema caótico no existan soluciones periódicas exactas.
Un proyecto web llamado Three Body Orbits ha reunido 3.915 soluciones periódicas exactas procedentes de 13 artículos de investigación en un atlas interactivo. Tres cuerpos de masa idéntica, bajo condiciones iniciales rigurosamente ajustadas, orbitan entre sí completando ciclos completos y regresando con precisión milimétrica a sus posiciones y velocidades de partida. Que las matemáticas teóricas carezcan de una fórmula general cerrada es algo radicalmente distinto de la impredecibilidad en el cálculo de ingeniería. Cuando se despliega la potencia de cálculo suficiente, es posible cartografiar islas de orden en medio del océano del caos.
Figura: Vista general de la página principal de Three Body Orbits, donde cada órbita periódica ocupa un lugar en el mapa celeste y las trayectorias similares se agrupan en islas. Fuente: Captura de pantalla de Three Body Orbits
Desafiando el límite de 32 dígitos de precisión: por qué se descartó la simulación en tiempo real en el navegador
Detrás de las 3.915 órbitas recopiladas en el atlas existe una exigente labor de reconstrucción ingenieril. El autor no se limitó a transcribir las condiciones iniciales publicadas en la literatura académica. Por motivos de maquetación y espacio, los artículos científicos suelen redondear los valores a solo 10 o 15 dígitos significativos. En un sistema caótico como el problema de los tres cuerpos, un error de truncamiento minúsculo del orden de 10⁻¹⁵ se amplifica exponencialmente en apenas unas pocas órbitas, provocando que la trayectoria diverja con rapidez y el sistema se desintegre.
Para lograr que cada órbita cierre de manera continua de principio a fin, el autor volvió a hacer converger las aplicaciones de sección de Poincaré (Poincaré section maps) de las ecuaciones diferenciales mediante librerías de aritmética de precisión arbitraria y el método de Newton-Raphson. Calibrando iterativamente las posiciones y velocidades iniciales hasta 32 dígitos decimales de precisión, consiguieron reducir la mediana del error de cierre a un asombroso 1,4 × 10⁻²⁴. En esta escala, el error de posición residual tras completar un periodo entero es mucho menor que el diámetro de un solo átomo. Un total de 3.889 órbitas alcanzan un cierre superior a 10⁻⁹.
Estas trayectorias, calculadas con una precisión extrema, plantearon un dilema técnico en el desarrollo frontend. Si se hubiera optado por una simulación en tiempo real en el navegador —permitiendo que JavaScript integrara numéricamente las ecuaciones mediante aritmética nativa de punto flotante de 64 bits y doble precisión—, las órbitas teóricamente cerradas habrían derivado rápidamente por la acumulación de errores de redondeo.
Por ello, el sitio web adoptó la arquitectura de “splines precalculados + reproducción en el cliente”. El servidor comprime las coordenadas reales de alta precisión en más de 9,7 millones de nodos de splines polinómicos de orden superior y los transmite al navegador. El cliente solo tiene que interpolar y renderizar con suavidad dichos puntos en función de las marcas de tiempo. Al cambiar espacio de almacenamiento precalculado por fidelidad numérica en tiempo de ejecución y aislar la computación pesada en el backend, el navegador puede dedicar todo su rendimiento a una renderización visual fluida a 60 fotogramas por segundo.
Figura: Vista de la órbita Figure-8 (figura en ocho), donde tres cuerpos se entrelazan en una trayectoria de larga exposición. Fuente: Captura de pantalla de Three Body Orbits
Tres generaciones de investigadores logran 3.915 soluciones exactas
En estas 3.915 órbitas se condensa una auténtica odisea académica en dinámica no lineal desarrollada a lo largo de décadas. No es el fruto de un único laboratorio, sino un mapa de clasificación articulado a través de generaciones por equipos de investigación internacionales que trabajaron en supercomputadores.
El equipo liderado por Shijun Liao en la Universidad Jiao Tong de Shanghái (grupo Li–Liao) aportó la inmensa mayoría: 1.964 órbitas. Empleando supercomputadores y su estrategia de simulación numérica limpia (Clean Numerical Simulation, CNS), exploraron sistemáticamente un vasto universo de estructuras topológicas desconocidas hasta la fecha. El equipo de la Universidad de Sofía en Bulgaria (Sofia) remitió 709 órbitas complejas, enriqueciendo notablemente las familias de armónicos superiores. Las célebres coreografías (Choreographies) de las matemáticas representan 583 órbitas, en las cuales los tres cuerpos parecen hacer malabares sobre una misma cuerda floja, compartiendo una única curva geométrica idéntica.
Asimismo, se catalogan 311 soluciones especiales con características de caída libre (Free-fall), y 230 soluciones de gran belleza geométrica descubiertas por el grupo de investigación de Belgrado en Serbia (Šuvakov y Dmitrašinović). La lista incluye incluso 52 órbitas clásicas deducidas por Roger Broucke en 1975. Junto con 66 órbitas recién descubiertas mediante búsqueda colaborativa distribuida en la propia web, configuran un compendio exhaustivo de la dinámica de tres cuerpos. Cada conjunto de coordenadas representa cientos de miles de horas de procesamiento consumidas en clústeres de supercomputación.
Un toque con el deslizador para poner a prueba la estabilidad de Floquet
En el atlas, 1.593 órbitas están explícitamente clasificadas como “linealmente estables”, lo que representa el 40,7 % del total. En los sistemas dinámicos, el parámetro de referencia para determinar si una órbita periódica resiste perturbaciones externas es el multiplicador de Floquet. Este valor mide cómo evoluciona una pequeña perturbación tras completarse un periodo entero.
Si el valor absoluto (módulo) de un multiplicador de Floquet es mayor que 1, el error se acumula exponencialmente y el sistema colapsa de inmediato. En cambio, cuando el multiplicador máximo de Floquet es estrictamente igual a 1,000, el sistema se sitúa en la frontera de estabilidad. Este 1,000 significa que las perturbaciones solo provocan oscilaciones periódicas acotadas, conservando intacta la topología global. Un multiplicador de 1,000 es el anclaje matemático que certifica la persistencia de estas órbitas en el mar del caos.
Para que la diferencia física entre órbitas estables e inestables resulte intuitiva, la plataforma incorporó un deslizador “Kick” en la parte derecha de la interfaz. Al desplazarlo en el navegador para aplicar un ligero vector de velocidad a cualquiera de los cuerpos, el integrador numérico en tiempo real del frontend toma el control de inmediato, apartándose de los splines precalculados.
En órbitas linealmente estables como Broucke R12, una pequeña perturbación solo provoca que los astros oscilen brevemente alrededor de su trayectoria original; gracias a su propia resiliencia dinámica, el sistema recupera el equilibrio y prosigue su danza a tres bandas. En cambio, en una órbita inestable, un leve toque del deslizador hace que el sistema triple se desintegre en unos pocos ciclos, expulsando los cuerpos al espacio profundo de forma caótica. Este diseño interactivo minimalista traduce el abstracto análisis de autovalores de ecuaciones diferenciales en una experiencia física palpable y directa. En los foros de debate de Hacker News, físicos e ingenieros mantuvieron intensas discusiones sobre los límites y umbrales tolerables de estas perturbaciones.
Convertir el barrido de parámetros en un proyecto de computación distribuida
Recientemente, la plataforma incorporó además una función de computación colaborativa en el navegador: el modo “Hunt”. Los usuarios pueden activar esta modalidad para donar ciclos de CPU inactivos y rastrear órbitas aún no registradas dentro de rangos paramétricos determinados. Cuando el dispositivo de un usuario converge en una nueva solución periódica linealmente estable, el descubridor tiene el privilegio de dar nombre a ese sistema planetario algorítmico. Una tarea tediosa de barrido numérico de condiciones iniciales se ha transformado así en un juego de computación distribuida para geeks.
El público tiende a concebir la ausencia de solución analítica en el problema de los tres cuerpos como una condena de caos perpetuo. Sin embargo, este atlas pone de manifiesto que, cuando los algoritmos permiten verificar soluciones cerradas con una precisión de 24 decimales, surgen del caos cerca de cuatro mil islas de riguroso orden matemático. El problema de los tres cuerpos ha dejado atrás el callejón sin salida del fatalismo filosófico para convertirse en un fértil laboratorio de ingeniería impulsado por datos concretos y una colosal potencia computacional.
Enlaces de referencia:
- Sitio web Three Body Orbits
- Physical Review Letters 110, 114301
- Science China Physics, Mechanics & Astronomy 60, 129511
- Sofia University Research Group Publications
- Celestial Mechanics 12, 439-462
- Debate en Hacker News (item?id=49670852)
