Bifurcaciones en Lorenz
De la calma al caos: el parametro \(\rho\) decide todo
¿Que es una bifurcacion?
Una bifurcacion es un cambio cualitativo en el comportamiento de un sistema dinamico que ocurre cuando un parametro cruza un valor critico. No es un cambio gradual en la amplitud o frecuencia del movimiento, sino una transformacion estructural: el tipo de atractor cambia, nuevos equilibrios nacen o mueren, orbitas periodicas aparecen o desaparecen. La palabra proviene del latin bifurcus ("dos horquillas"), porque en el caso mas simple, un camino unico se divide en dos.
Los principales tipos de bifurcacion son:
- Saddle-node (nodo-silla): Dos equilibrios (uno estable, uno inestable) colisionan y se aniquilan mutuamente.
- Pitchfork (horquilla): Un equilibrio pierde estabilidad y nacen dos nuevos equilibrios simetricos.
- Hopf: Un equilibrio pierde estabilidad y nace (o muere) un ciclo limite, una orbita periodica cerrada.
- Duplicacion de periodo: Una orbita periodica duplica su periodo, generando una cascada que conduce al caos.
El sistema de Lorenz exhibe al menos dos de estos tipos al variar \(\rho\), y el resultado es una de las secuencias de bifurcaciones mas ricas de la dinamica no lineal.
El viaje de \(\rho\)
Al variar \(\rho\) desde 0 hasta valores grandes (manteniendo \(\sigma = 10\) y \(\beta = 8/3\) fijos), el sistema de Lorenz atraviesa una secuencia fascinante de regimenes dinamicos. Cada transicion es una bifurcacion.
\(\rho < 1\): el reposo absoluto
Para \(\rho < 1\), el unico equilibrio es el origen y es globalmente estable. Toda trayectoria, independientemente de su condicion inicial, converge asintoticamente a \((0, 0, 0)\). No hay conveccion, no hay oscilaciones, no hay caos. El gradiente termico es insuficiente para vencer la disipacion viscosa. En el diagrama de bifurcacion, este regimen aparece como un unico punto en \(z = 0\).
\(\rho = 1\): bifurcacion de horquilla
En \(\rho = 1\) exactamente, el origen sufre una bifurcacion de horquilla supercritica (supercritical pitchfork). Uno de los eigenvalores del origen cruza cero: el origen se vuelve inestable y nacen simultaneamente dos nuevos equilibrios, \(C^+\) y \(C^-\).
Condicion de la bifurcacion de horquilla
$$\rho = 1: \quad \lambda = \frac{-(\sigma+1) + \sqrt{(\sigma+1)^2 + 4\sigma(\rho-1)}}{2} = 0$$Los nuevos equilibrios \(C^{\pm} = (\pm\sqrt{\beta(\rho-1)}, \pm\sqrt{\beta(\rho-1)}, \rho-1)\) nacen con tamano cero en \(\rho = 1\) y crecen continuamente con \(\rho\). Son estables: las trayectorias convergen a uno u otro en espiral, representando conveccion estacionaria en una de las dos direcciones.
\(1 < \rho < \rho_H \approx 24.74\): espirales estables
En este rango, \(C^+\) y \(C^-\) son puntos fijos estables con eigenvalores complejos: las trayectorias espiralan hacia ellos con oscilaciones decrecientes. A medida que \(\rho\) aumenta, la parte imaginaria de los eigenvalores crece (oscilaciones mas rapidas) y la parte real negativa se acerca a cero (oscilaciones que decaen mas lentamente). El sistema se aproxima peligrosamente al umbral de inestabilidad.
\(\rho = \rho_H \approx 24.74\): bifurcacion de Hopf
En el valor critico \(\rho_H\), la parte real de los eigenvalores complejos de \(C^{\pm}\) cruza cero. Esta es una bifurcacion de Hopf subcritica: los equilibrios pierden estabilidad, pero no nace un ciclo limite estable. En cambio, las trayectorias son expulsadas directamente hacia el atractor extrano que ya existe como estructura global.
Valor critico de Hopf
$$\rho_H = \sigma\,\frac{\sigma + \beta + 3}{\sigma - \beta - 1} = 10 \times \frac{10 + 8/3 + 3}{10 - 8/3 - 1} \approx 24.74$$Hay un detalle sutil: el atractor caotico ya existe para valores de \(\rho\) ligeramente menores que \(\rho_H\) (desde \(\rho \approx 24.06\)), pero coexiste con los puntos fijos estables. En ese rango, la cuenca de atraccion de \(C^{\pm}\) se encoge progresivamente, y perturbaciones grandes pueden enviar la trayectoria al atractor caotico. Este es el regimen de caos transitorio y metaestabilidad.
\(\rho > \rho_H\): el atractor extrano
Para \(\rho > \rho_H\), los tres equilibrios son inestables. No hay puntos fijos ni ciclos limites estables. Las trayectorias viven sobre el atractor extrano de Lorenz: una estructura fractal que atrae todas las trayectorias pero sobre la cual el movimiento es caotico. La trayectoria espiralea alrededor de \(C^+\) y \(C^-\), saltando entre ambos de forma aperiodica y sensible a las condiciones iniciales.
Diagrama de bifurcacion
El diagrama de bifurcacion es una herramienta visual poderosa: para cada valor de \(\rho\), integramos el sistema, descartamos un transitorio largo, y registramos los maximos locales de \(z(t)\). Si el movimiento converge a un punto fijo, obtenemos un unico punto. Si es periodico, obtenemos unos pocos puntos discretos. Si es caotico, obtenemos una nube densa de puntos.
El eje horizontal del diagrama es \(\rho\) (de 0 a 200) y el eje vertical muestra los valores de \(z\) en los maximos locales. El resultado revela toda la estructura de bifurcaciones del sistema de una sola vista: la transicion del punto fijo al caos, las ventanas periodicas, las cascadas de duplicacion de periodo.
Ventanas de periodicidad
Dentro del regimen caotico, existen islas de orden. Para ciertos rangos estrechos de \(\rho\), el movimiento caotico da paso a orbitas periodicas estables. Estas son las ventanas periodicas, analogas a las que se observan en el mapa logistico.
La ventana periodica mas notable ocurre alrededor de \(\rho \approx 100\), donde aparece una orbita periodica estable de periodo largo. Tambien hay ventanas prominentes cerca de \(\rho \approx 145\) y \(\rho \approx 166\). En el diagrama de bifurcacion, las ventanas aparecen como bandas estrechas donde los puntos se condensan en unas pocas lineas discretas, en contraste con las bandas anchas del caos donde los puntos llenan regiones densas.
La alternancia entre caos y periodicidad demuestra que el caos no es un estado permanente e irreversible: el sistema puede entrar y salir del caos al variar un parametro de forma continua. Esta estructura intrincada, con ventanas periodicas anidadas dentro del caos, es una de las huellas universales de los sistemas dinamicos no lineales.
"Las bifurcaciones son las puertas entre mundos dinamicos. Cada transicion revela una nueva faceta del sistema. El parametro rho es la llave que abre todas esas puertas, desde la calma del punto fijo hasta la tormenta del caos."
| Rango de \(\rho\) | Bifurcacion | Comportamiento | Tipo |
|---|---|---|---|
| \(0 < \rho < 1\) | — | Convergencia al origen | Estable |
| \(\rho = 1\) | Pitchfork supercritica | Nacen \(C^+\) y \(C^-\) | Bifurcacion |
| \(1 < \rho < 24.06\) | — | Espirales estables en \(C^{\pm}\) | Estable |
| \(24.06 < \rho < 24.74\) | — | Coexistencia: caos transitorio + \(C^{\pm}\) | Metaestable |
| \(\rho \approx 24.74\) | Hopf subcritica | \(C^{\pm}\) pierden estabilidad | Bifurcacion |
| \(24.74 < \rho \lesssim 313\) | — | Caos con ventanas periodicas | Caotico |
| \(\rho \approx 100, 145, 166\) | Saddle-node de orbitas | Ventanas periodicas | Periodico |
Ejercicios
-
Identifica las bifurcaciones en el diagrama.
En el LAB, pasa el cursor lentamente por el diagrama de bifurcacion. ¿Puedes ver la transicion en \(\rho = 1\) donde aparecen los puntos fijos? ¿Y la explosion de caos en \(\rho \approx 24.74\)?
-
Encuentra una ventana periodica.
Busca una region dentro del caos donde los puntos se condensan en pocas lineas discretas. Haz clic en esa region para ver el atractor correspondiente en el panel derecho. ¿Cuantas ramas tiene la orbita periodica?
-
Calcula \(\rho_H\) para otros parametros.
Usa la formula \(\rho_H = \sigma(\sigma + \beta + 3)/(\sigma - \beta - 1)\). ¿Que pasa si \(\sigma = 5\)? ¿Y si \(\sigma = 20\)? ¿El caos aparece antes o despues?
LAB: Diagrama de bifurcacion de Lorenz
Diagrama de bifurcacion: maximos locales de \(z(t)\) vs \(\rho\), barriendo de 0 a 200. El panel izquierdo muestra el diagrama completo. Haz clic en cualquier punto del diagrama para ver la proyeccion \(x\text{-}z\) del atractor correspondiente a ese valor de \(\rho\) en el panel derecho.