Sistemas dinámicos y caos · Artículo 04

KAM: el orden que resiste

El caos hamiltoniano no invade a sangre y fuego: los toros de acción-ángulo se defienden, y los más irracionales caen los últimos. El teorema KAM, los huecos de Kirkwood en el cinturón de asteroides — y el veredicto final sobre el sistema solar. Cierre del módulo.

Queda el caos SIN fricción — el de los sistemas hamiltonianos, donde no hay atractores (Liouville lo prohibió en III.2) y la pregunta es otra: los toros de acción-ángulo de III.3, ese orden perfecto de relojes… ¿sobreviven cuando el sistema no es exactamente separable? La respuesta — el teorema KAM (Kolmogórov 1954, Arnold, Moser) — es el resultado más profundo de la mecánica del siglo XX: la mayoría de los toros sobrevive, deformada pero intacta. El caos no conquista: se infiltra por las grietas — y las grietas son las resonancias.

Prerrequisitos: los toros de III.3 (el módulo entero cobra aquí su segunda factura) y la separatriz de III.2.

Quién cae y quién resiste: la aritmética decide

En un toro con dos frecuencias, todo depende del cociente ω₁/ω₂. Racional (órbita que se cierra: resonancia): el toro es frágil — la perturbación lo golpea en fase, como el columpio de II.3, y lo rompe en un collar de islas (teorema de Poincaré-Birkhoff) rodeadas de una fina capa caótica nacida de su separatriz (la grieta anunciada en la hoja de III.3). Irracional: los empujones nunca se sincronizan y el toro aguanta — y aguanta más cuanto peor aproximable por racionales es su cociente. El campeón de la resistencia es el número más irracional que existe, la razón áurea (√5−1)/2: su toro es el último en morir en medio mundo de modelos. Que la frontera entre orden y caos la dicte la teoría de números — fracciones continuas decidiendo destinos dinámicos — es de esas conexiones que justifican una carrera entera.

Ejemplo resuelto 1 · Los huecos de Kirkwood

Problema. El histograma de semiejes del cinturón de asteroides tiene huecos tajantes: casi nadie orbita donde el período es 1/3, 2/5 o 1/2 del de Júpiter. Calcula dónde cae el hueco 3:1 y explícalo con el vocabulario del artículo.

Solución. Kepler III con T = T_J/3:

a=aJ(13)2/3=5,20×0,481=2,50 UAa = a_J \left(\tfrac{1}{3}\right)^{2/3} = 5{,}20 \times 0{,}481 = 2{,}50\ \mathrm{UA}

— y ahí está el hueco, medido desde 1866. El mecanismo es este artículo: en la resonancia 3:1, los tirones de Júpiter llegan sincronizados (mismo punto de la órbita cada tres vueltas), el toro del asteroide se rompe, y su capa caótica bombea la excentricidad a saltos hasta cruzar la órbita de Marte o la Tierra — donde un encuentro lo expulsa (o lo entrega: buena parte de los meteoritos que caen aquí salió por la puerta 3:1, cosa que confirman sus órbitas medidas).

Resultado. Kirkwood vio los huecos en 1866; Wisdom calculó en 1983 la difusión caótica que los vacía — resolviendo de paso por qué los huecos son tan nítidos. El cinturón de asteroides es un experimento KAM de 4500 millones de años: donde había toros racionales hay vacío; donde irracionales, asteroides. La aritmética, esculpida en el cielo a escala de unidades astronómicas.

El veredicto del sistema solar — y el módulo en limpio. Ahora sí se puede enunciar I.5 con todas las palabras: el sistema solar es un sistema KAM — casi integrable, con sus planetas en toros supervivientes y sus capas caóticas reales (Lyapunov de 5 millones de años, el 1 % de Mercurio, la oblicuidad de Marte de III.4, los huecos de Kirkwood). Ni reloj perfecto ni ruleta: un archipiélago de orden en un mar fino de caos, cartografiado por el teorema más difícil del curso. Y el módulo entero en tres líneas: la disipación fabrica atractores (01), la no linealidad duplica períodos con ritmo universal (02), tres dimensiones bastan para el hojaldre extraño (03), y sin fricción el orden resiste por irracionalidad (04). El curso de mecánica está completo salvo un estado de la materia: el que fluye. III.7.

Ejercicios

Ejercicio 1

La sección de Poincaré — mirar el sistema solo cada vez que cruza un plano — convierte los toros en curvas y el caos en polvo de puntos. ¿Por qué es LA herramienta para «ver» KAM, y qué esperas ver al aumentar la perturbación?

Solución

Baja una dimensión: el toro (2D) corta el plano en una curva cerrada; una órbita caótica lo salpica todo. La película KAM estándar al subir la perturbación: curvas concéntricas (toros) → collares de islas donde había racionales (Poincaré-Birkhoff) → mares caóticos creciendo desde las separatrices → las últimas curvas doradas resistiendo → caos conectado cuando caen (el criterio de Chirikov: resonancias que se solapan). Cada retrato de sección en los libros de Strogatz o en los papers de dinámica planetaria es esta película congelada — y ya sabes leer todos sus fotogramas.

Ejercicio 2

Los anillos de Saturno tienen su Kirkwood: la división de Cassini está en resonancia 2:1 con la luna Mimas. Compruébalo con Kepler III (T_Mimas = 22,6 h; la división está a ~117 500 km; Saturno: GM tal que un anillo a 117 500 km orbita en ~11,3 h).

Solución

11,3 h = 22,6/2 exacto: las partículas de la división recibirían el tirón de Mimas sincronizado cada dos vueltas — toro roto, excentricidad bombeada, choques con las vecinas, y la franja se vacía: 4700 km de hueco visibles con telescopio de aficionado. Los anillos son el laboratorio KAM más fotogénico del sistema solar (con estructura fina esculpida por decenas de resonancias más — cada surco, una fracción). Kirkwood, Cassini y el ejercicio 1: el mismo teorema en tres decorados.

Ejercicio 3

¿Por qué el teorema KAM fue también la salvación de la teoría de perturbaciones de III.3 — la que mantiene el GPS y las efemérides — tras dos siglos de «pequeños divisores» que amenazaban con invalidarla?

Solución

Las series de perturbaciones dividen por combinaciones n₁ω₁ − n₂ω₂ que cerca de las resonancias se hacen minúsculas: términos gigantes, series sospechosas — Poincaré llegó a temer que TODA la mecánica celeste perturbativa fuera formal pero falsa. KAM demostró dónde converge de verdad: sobre los toros suficientemente irracionales (condición diofántica), que son la mayoría — la teoría de perturbaciones es legítima exactamente donde el orden sobrevive, y peligrosa exactamente en las grietas resonantes. Dos siglos de cálculos de Le Verrier a los GPS quedaron cimentados — con la advertencia cartografiada. Un teorema que a la vez salva y delimita una disciplina entera: así se cierra la mecánica clásica del siglo XX.

Resumen en frío · Módulo III.6 · Sistemas dinámicos y caos

Todo lo que el módulo deja utilizable, con dónde se dedujo cada cosa — y, cuando toca, qué deuda de un módulo anterior salda cada fila. Pensado para leerse dentro de seis meses sin releer nada: si para resolver un problema típico del módulo hay que volver al texto a buscar una constante, esta tabla ha fallado.

QuéFórmula o valorDónde
Qué cambia con disipaciónLiouville deja de valer: el volumen de fases se contrae y nacen los atractores, conjuntos donde las trayectorias terminanart. 01 · rompe deliberadamente III.2 art. 04
Catálogo de atractoresPunto fijo · ciclo límite · atractor extrañoart. 01 · art. 03
Ciclo límiteÓrbita periódica atractora: el sistema bombea energía cuando la amplitud es pequeña y la disipa cuando es grande. Ningún sistema lineal oscila soloart. 01
Oscilador de van der Pol (1926)ẍ − μ(1 − x²)ẋ + x = 0: fricción negativa mientras |x| < 1, positiva fuera. Modelo de relojes, láseres, cuerda de violín, neuronas y corazónart. 01
Amplitud rígida, fase blandaMoverse a lo largo del ciclo no cuesta nada, así que un acoplo débil solo puede empujar fases: modelo de Kuramoto (1975) — luciérnagas, marcapasos, aplausos, generadores de una red eléctricaart. 01, ej. 1 · generaliza los relojes de Huygens de II.3
Teorema de Poincaré-BendixsonEn el plano, toda trayectoria acotada acaba en (a) un punto fijo o (b) un ciclo límite. Nada más cabeart. 01, ej. res. 1
El caos exige tres variablesPéndulo simple (θ, p): jamás. Péndulo forzado (θ, p, fase del forzamiento): puede. Lorenz necesitó exactamente 3; el péndulo doble (4D) sobraba de sitioart. 01, ej. res. 1 · cobra la regla anunciada en III.2 art. 02
El reloj como ciclo límiteEl escape no impone ritmo: empuja cuando el péndulo pasa. La frecuencia la pone el péndulo (√(g/L)), no la cuerda — por eso da igual el peso de la pesaart. 01, ej. 2
FibrilaciónOtro atractor del mismo sistema —ondas espirales autosostenidas—, no ausencia de orden. Por eso el desfibrilador saca al corazón de todas las cuencas de golpe en vez de empujar suaveart. 01, ej. 3
Mapa logísticoxn+1 = r·xn(1 − xn)art. 02
Punto fijo y su estabilidadx* = 1 − 1/r; la perturbación se multiplica por f′(x*) = 2 − r en cada vuelta, luego es estable mientras |2 − r| < 1, es decir r < 3art. 02, ej. 1 · es III.5 (autovalores como test) en versión de bolsillo
Cascada de duplicacionesr₁ = 3 · r₂ = 3,4495 (= 1+√6) · r₃ = 3,5441 · r₄ = 3,5644 … acumulando en r∞ = 3,5699, donde empieza el caosart. 02
Constante de Feigenbaumδ = 4,669201…, el cociente límite de encogimiento entre ventanas de parámetro consecutivasart. 02, ej. res. 1
UniversalidadCualquier sistema que llegue al caos duplicando períodos lo hace con el mismo δ: cerca de la acumulación solo importa que el máximo del mapa sea cuadrático. Convección de helio (Libchaber, 1979): δ ≈ 4,4 ± 0,1; grifo que gotea (Shaw, 1984); láseres y circuitosart. 02, ej. res. 1
r = 4 tiene fórmula cerradaCon x = sen²(πθ) el mapa se vuelve θ → 2θ (mód 1): el error se dobla cada iteración — exponente de Lyapunov ln 2, exacto. Estirar y plegar, la receta del panaderoart. 02, ej. 2 · la gota filamentándose de III.2
Ventanas dentro del caosLa de período 3 en r ≈ 3,83 (= 1+√8), con réplicas en miniatura del árbol entero: autosimilaridad fractal. «Período tres implica caos» (Li y Yorke, 1975) y el orden de Sharkovskiiart. 02, callout
Alternans cardíacoLa primera duplicación del ciclo límite del corazón; su presencia en el electrocardiograma es predictor reconocido de riesgo de arritmia. La universalidad no distingue helio líquido de miocardioart. 02, ej. 3
Sistema de Lorenz (1963)ẋ = σ(y − x), ẏ = x(ρ − z) − y, ż = xy − βz, con σ = 10, ρ = 28, β = 8/3art. 03
Contracción del volumendiv = −σ − 1 − β = −13,67: el volumen de fases cae unos 6 órdenes de magnitud por unidad de tiempoart. 03, ej. 1
Espectro de Lyapunovλ₁ ≈ +0,9 · λ₂ = 0 (a lo largo del flujo) · λ₃ ≈ −14,6, y su suma es exactamente la divergencia. Estirar en una dirección mientras se aplasta en otra: el hojaldreart. 03, ej. 1
El ritmo de la impredecibilidadEl error se duplica cada ln2/0,9 ≈ 0,77 unidades de tiempo, y una vuelta del atractor dura ~0,75: casi exactamente una duplicación por vuelta, un bit de información perdido por paseoart. 03
Presupuesto de la mariposaDe un error de 10⁻⁴ a orden 1: ln(10⁴)/0,9 ≈ 10 unidades de tiempo, unas 13 vueltas — media página de la impresora de 1961art. 03, ej. 2 · pone número a la anécdota de I.5
Dimensión fractal del atractor≈ 2,06: más que una sábana, menos que un sólido. Su sección transversal es un conjunto de Cantorart. 03
Qué se puede hacer con el caosUsarlo (cifrado con circuitos de Lorenz o Chua; sincronización de caóticos — Pecora y Carroll, 1990) y detectarlo: reconstruido el espacio de fases (Takens), el determinismo cae sobre un objeto de dimensión baja mientras el ruido lo llena todoart. 03, ej. res. 1
Honestidad sobre la mariposaEl atractor es acotado: caos no es «cualquier cosa puede pasar». Y Lorenz no demuestra que la atmósfera sea caótica, sino un principio de posibilidad — modelo que ilustra un mecanismo, no que describe un sistemaart. 03, callout y ej. 3
Teorema KAM (Kolmogórov 1954, Arnold, Moser)Sin fricción no hay atractores, y la pregunta es otra: los toros de acción-ángulo, ¿sobreviven a una perturbación? La mayoría sí, deformados pero intactos. El caos no conquista: se infiltra por las resonanciasart. 04 · segunda factura de III.3
Quién cae y quién resisteCociente ω₁/ω₂ racional (resonancia): la perturbación golpea en fase y rompe el toro en un collar de islas (Poincaré-Birkhoff) con una capa caótica nacida de su separatriz. Irracional: los empujones nunca se sincronizan y el toro aguantaart. 04
El último toro en morirEl de cociente igual a la razón áurea, (√5−1)/2 = 0,618…, el número peor aproximable por racionales que existe. La frontera entre orden y caos la dicta la teoría de númerosart. 04
Sección de PoincaréMirar el sistema solo cuando cruza un plano: un toro corta en curva cerrada, una órbita caótica lo salpica de polvo. La película al subir la perturbación: curvas → collares de islas → mares caóticos creciendo desde las separatrices → últimas curvas doradas → caos conectado (criterio de Chirikov: resonancias que se solapan)art. 04, ej. 1
Huecos de KirkwoodResonancia 3:1 con Júpiter por Kepler III: a = a_J·(1/3)^(2/3) = 5,20 × 0,481 = 2,50 UA. Vistos en 1866; Wisdom calculó en 1983 la difusión caótica que bombea la excentricidad y los vacíaart. 04, ej. res. 1
División de CassiniA ~117 500 km el período es 11,3 h, la mitad de las 22,6 h de Mimas: resonancia 2:1 y 4700 km de huecoart. 04, ej. 2
Veredicto sobre el sistema solarEs un sistema KAM: casi integrable, con planetas en toros supervivientes y capas caóticas reales — Lyapunov de 5 millones de años, el 1 % de probabilidad de Mercurio, la oblicuidad de Marte, los huecos de Kirkwood. Ni reloj ni ruletaart. 04, callout · cierra I.5 y III.4
KAM salva la teoría de perturbacionesLos «pequeños divisores» n₁ω₁ − n₂ω₂ hacían sospechosas las series de la mecánica celeste. KAM demuestra que convergen sobre los toros suficientemente irracionales (condición diofántica): legítima donde el orden sobrevive, peligrosa en las grietas resonantesart. 04, ej. 3 · cimenta la teoría canónica de III.3 art. 02