Atractores: la forma del caos

La paradoja que hace interesante el tema es que las dos propiedades conviven: el volumen se contrae —todo cabe en un conjunto de medida nula— y al mismo tiempo las distancias dentro de ese conjunto crecen. La única manera de que ambas cosas pasen a la vez es que el atractor se pliegue sobre sí mismo infinitas veces. De ahí sale su dimensión no entera, y de ahí sale que predecir a largo plazo sea imposible aunque las ecuaciones sean exactas y tengan tres líneas.

Este artículo arma el aparato de medición completo y después lo aplica: cada capítulo trae una figura que integra el sistema que cita, con Runge-Kutta de cuarto orden y paso fijo, y el catálogo entero está validado en Node contra las cifras publicadas. En el medio está la vitrina: 23 sistemas en tres dimensiones, girables desde cualquier ángulo, con sus parámetros en vivo.

Dado un flujo φ_t en el espacio de estados, un atractor es un conjunto compacto A que cumple tres cosas: es invariante (una órbita que empieza en A se queda en A), atrae a un entorno suyo —el conjunto de puntos cuyas órbitas terminan en A es su cuenca de atracción— y es mínimo: ningún subconjunto propio suyo cumple lo anterior. Esa minimalidad es lo que impide llamar «atractor» a cualquier bolsa que contenga uno.

Lo que distingue a los tipos es la geometría del conjunto, no la complejidad del dibujo. Un punto fijo tiene dimensión cero; un ciclo límite, uno; un toro, dos. Un atractor extraño tiene dimensión no entera, y esa es su definición estricta. La confusión habitual es tratarlo como sinónimo de caótico: no lo es. Existen atractores extraños no caóticos —fractales con λ₁ < 0, descritos por Grebogi y colaboradores en 1984— y el orden inverso también se da. En los 23 sistemas de este catálogo las dos propiedades coinciden, pero conviene saber que son propiedades distintas.

En el sistema de Lorenz la divergencia no depende del punto: ∇·f = −σ − 1 − β = −13,667. El volumen se multiplica por e^{−13,667 t}, así que a los dos segundos ya se redujo a una billonésima, y en el límite es exactamente cero. El atractor tiene medida nula en el espacio de estados.

Y sin embargo la nube no se junta en un punto: se estira en una dirección mientras se aplasta en las otras dos. Ese es el reparto que después leeremos como espectro de Lyapunov, y la suma de los tres exponentes tiene que dar exactamente la divergencia. En nuestra validación esa igualdad se cumple hasta la cuarta cifra en todos los sistemas polinomiales del catálogo: es el control más barato de que la medición está bien hecha.

Un sistema conservativo tiene ∇·f = 0 y por el teorema de Liouville conserva el volumen: no puede tener atractores. Su caos existe igual, pero vive en un mar de órbitas donde flotan islas regulares. El catálogo incluye uno a propósito —Nosé-Hoover— porque el contraste aclara qué parte del fenómeno se debe a la disipación y qué parte no.

Esa dependencia logarítmica es la mala noticia del asunto: mejorar la medición inicial en un factor mil sólo compra ln(1000)/λ₁ de tiempo extra —siete unidades y media en el caso de Lorenz—. No es un problema de instrumentos; es la estructura del sistema.

El espectro completo se mide con el algoritmo de Benettin: se integra el sistema junto a una base de vectores tangentes, se re-ortonormaliza con Gram-Schmidt cada cierto número de pasos y se acumula el logaritmo de cuánto se estiró cada dirección. Un detalle que cuesta caro: la parte tangente tiene que integrarse con el mismo orden que el estado. Con el estado en RK4 y los vectores tangentes en Euler, nuestro Lorenz daba λ₁ = 1,079 —19 % de más— y λ₂ no bajaba a cero. Integrando el sistema aumentado completo con RK4, λ₁ = 0,9072 contra 0,9056 publicado y λ₂ = −0,0005.

Ese λ₂ ≈ 0 no es un detalle cosmético: en un flujo continuo acotado que no termina en un punto fijo, la dirección a lo largo de la trayectoria no se estira ni se encoge, así que uno de los exponentes tiene que ser nulo. Es el testigo más honesto de que el método funciona, porque no se le impone: aparece solo.

La manera directa de medirla es contar: se cubre el conjunto con cajas de lado ε y se mira cómo crece el número necesario cuando ε se achica. Si N(ε) ~ ε^−D, ese D es la dimensión de caja. Para el atractor de Hénon la literatura da D₀ = 1,261 ± 0,003; nuestro conteo sobre medio millón de puntos da 1,28. La diferencia no es un desacuerdo: el conteo de cajas converge muy despacio y a resolución finita sobreestima sistemáticamente, algo conocido desde los primeros trabajos de Grassberger.

La vía indirecta es más barata y se usa en todo el catálogo: la conjetura de Kaplan-Yorke reconstruye la dimensión desde los exponentes, sumándolos hasta que la suma deja de ser positiva e interpolando la fracción que falta.

Seguir una órbita en tres dimensiones es incómodo. Poincaré propuso mirar sólo los cruces por un plano: cada vez que la trayectoria lo atraviesa en un sentido se anota el punto. El flujo continuo se convierte así en un mapa discreto de una dimensión menos, que conserva la estabilidad de las órbitas periódicas y la estructura del atractor.

En el sistema de Rössler el resultado es notable: los cruces no forman una nube sino una curva, y graficar cada cruce contra el anterior da una función unimodal, con una sola joroba. Es decir, un objeto tridimensional se comporta, en lo esencial, como el mapa logístico de una sola variable que veremos en el capítulo siguiente. Esa reducción es la que permite aplicar toda la teoría de mapas unidimensionales a sistemas de ecuaciones diferenciales.

Lo que hace de esto un resultado y no una curiosidad es la palabra cualquiera. δ no depende de la función: sale igual para el logístico, para el seno, para un circuito electrónico o para una llave de agua goteando. Es universalidad en el sentido fuerte, y fue lo que convirtió el caos en una disciplina con predicciones contrastables en vez de una colección de anécdotas numéricas.

Las otras dos rutas clásicas son la cuasiperiódica de Ruelle-Takens —un toro que se rompe tras tres bifurcaciones de Hopf— y la intermitencia de Pomeau y Manneville, donde el sistema alterna tramos casi periódicos con estallidos irregulares cada vez más frecuentes.

Difeomorfo significa que las propiedades que nos importan se conservan: los exponentes de Lyapunov, la dimensión y la estructura de órbitas periódicas son las mismas. La forma se ve deformada; los números, no. Sobre esa base funcionan la estimación de dimensión en series experimentales, la predicción a corto plazo por vecinos cercanos y los métodos de causalidad de Sugihara en ecología.

La advertencia va junta con el teorema: aplicado sin cuidado produce dimensiones fractales para cualquier cosa, incluido el ruido coloreado. La prueba mínima es la de datos subrogados de Theiler (1992): se generan series con el mismo espectro de potencia pero con las fases aleatorizadas, se les aplica el mismo análisis y se verifica que el resultado cambie. Buena parte de la literatura de los años noventa sobre «caos» en series fisiológicas y financieras no sobrevive a ese control.

Sincronizar. Dos sistemas caóticos idénticos pueden ponerse de acuerdo si uno recibe del otro una sola variable. Pecora y Carroll lo mostraron en 1990: el receptor descarta su propia x y usa la del emisor, y el error entre ambos cae exponencialmente aunque ninguno de los dos sea predecible por separado. Es contraintuitivo y es la base de las propuestas de comunicación por caos.

Controlar. Un atractor extraño contiene infinitas órbitas periódicas inestables. El método OGY —Ott, Grebogi y Yorke, 1990— aprovecha eso: con perturbaciones diminutas y bien cronometradas en un parámetro, el sistema se estabiliza sobre la órbita periódica que uno elija. Justamente la sensibilidad que impide predecir es la que permite controlar barato.

Predecir mientras se pueda. El horizonte es finito, pero no es cero. En meteorología ese límite es del orden de dos semanas para la escala sinóptica, y el trabajo actual mide cuánto falta para alcanzarlo. En aprendizaje automático, los reservoir computers de Pathak y colaboradores (2018) reproducen la dinámica de sistemas espacio-temporales caóticos durante varios tiempos de Lyapunov, y —más interesante que la predicción— reproducen después el atractor correcto en el largo plazo, aunque la trayectoria puntual ya se haya perdido.

Reconocer el mecanismo. El modelo de Rikitake produce inversiones de polaridad magnética con la estadística irregular que muestra el registro geológico; el circuito de Chua pone un atractor extraño en un protoboard con un osciloscopio al lado; los láseres de cavidad en anillo siguen el mapa de Ikeda. Cuando un sistema real exhibe una cascada de Feigenbaum medida, la universalidad deja de ser una metáfora y pasa a ser una predicción cumplida.

Lo que no está resuelto conviene decirlo con la misma claridad. El cifrado por caos de baja dimensión fue roto por ataques de reconstrucción ya en los años noventa y no es criptografía en el sentido moderno. La «dimensión fractal» como biomarcador clínico sigue sin un resultado replicado con controles adecuados. Y afirmar que una serie económica es caótica exige, como mínimo, pasar la prueba de datos subrogados, cosa que casi nunca se hace.

El archivo validar.mjs corre en Node las mismas funciones que corren en pantalla y somete cada sistema a cuatro controles fijados de antemano: que la órbita esté acotada y no colapse; que sea disipativo (Σλ < 0, salvo el conservativo); que sea caótico (λ₁ > 0 con margen); y que λ₂ ≈ 0, el exponente a lo largo del flujo, que el método tiene que recuperar solo. Donde hay una cifra publicada se compara contra ella con una tolerancia del 5 % en λ₁ y del 3 % en la dimensión.

Dos resultados de la corrida vale la pena declarar. El primero: en todos los sistemas polinomiales, Σλ coincide con la divergencia del campo hasta la cuarta cifra —−13,6667 en Lorenz, −19,0000 en Lü, −41,0000 en Qi—, que es una identidad exacta y por lo tanto un control duro. El segundo: el atractor de Halvorsen con el valor que suele citarse, a = 1,89, no nos dio caos sino un ciclo límite —λ₂ = −0,36, sin exponente nulo— con dos convenciones de signo y varias semillas. El catálogo lo publica con a = 1,4, donde el caos es inequívoco, y deja dicho lo anterior en vez de esconderlo.