Cascada de duplicacion de periodo
La ruta universal al caos y la constante de Feigenbaum
La cascada universal
En el diagrama de bifurcacion del mapa logistico observamos una secuencia notable de transiciones. Cuando \(r\) aumenta desde 3, el punto fijo se desdobla en un ciclo de periodo 2. Poco despues, el periodo 2 se desdobla en periodo 4. Luego 4 se convierte en 8, 8 en 16, 16 en 32, y asi sucesivamente: \(1 \to 2 \to 4 \to 8 \to 16 \to 32 \to \ldots \to \infty \to \text{caos}\).
Esta secuencia se llama cascada de duplicacion de periodo (period-doubling cascade), y es una de las rutas mas comunes por las que un sistema pasa del orden al caos. Lo extraordinario es que esta cascada tiene una estructura matematica precisa y universal.
Las bifurcaciones se aceleran: cada nueva duplicacion ocurre en un intervalo de parametro mas pequeno que la anterior. La cascada completa ocupa un intervalo finito de \(r\) y culmina en un valor critico \(r_\infty\) donde el periodo se vuelve infinito. Mas alla de \(r_\infty\), el sistema es caotico (con las interrupciones de las ventanas de periodicidad).
Feigenbaum y la universalidad
En 1975, Mitchell Feigenbaum, un joven fisico del Laboratorio Nacional de Los Alamos, hizo un descubrimiento que le valio fama inmortal en la matematica. Mientras estudiaba la cascada de duplicacion de periodo en el mapa logistico usando una calculadora de bolsillo HP-65, noto que la razon entre intervalos sucesivos de bifurcacion convergia a un numero fijo.
Sea \(r_n\) el valor del parametro donde ocurre la bifurcacion que crea un ciclo de periodo \(2^n\). Feigenbaum definio la razon:
Este numero, \(\delta \approx 4.669\), es la constante de Feigenbaum. Significa que cada intervalo de bifurcacion es aproximadamente 4.669 veces mas pequeno que el anterior. La cascada se comprime geometricamente, convergiendo al punto critico \(r_\infty\).
Lo verdaderamente asombroso, lo que convirtio un resultado numerico curioso en un descubrimiento de primer orden, fue lo que Feigenbaum hizo a continuacion: verifico que el mismo numero \(\delta = 4.669...\) aparecia en el mapa seno (\(x_{n+1} = A \sin(\pi x_n)\)), en el mapa cuadratico general, en el mapa exponencial, y en cualquier otro mapa unimodal (un maximo, continuo, diferenciable). El numero era universal.
La constante de Feigenbaum
La primera constante de Feigenbaum mide la contraccion de los intervalos de parametro. Pero existe una segunda constante que mide la contraccion de las ramas del diagrama de bifurcacion en la direccion vertical. Si \(d_n\) es la "anchura" de la bifurcacion de orden \(n\) (la distancia entre los valores extremos del ciclo), entonces:
Ambas constantes, \(\delta\) y \(\alpha\), son numeros trascendentes que aparecen en todos los sistemas con ruta de duplicacion de periodo al caos. No importa si el sistema es un mapa logistico, una reaccion quimica, un circuito electronico, o un modelo de crecimiento poblacional. Si la ruta al caos es por duplicacion de periodo, las constantes de Feigenbaum estan ahi.
Universalidad
La universalidad de las constantes de Feigenbaum fue demostrada rigurosamente por Oscar Lanford en 1982, usando la teoria del operador de renormalizacion. La idea es que la duplicacion de periodo corresponde a un punto fijo de un operador funcional (el operador de composicion y reescalado), y las constantes \(\delta\) y \(\alpha\) son los autovalores de ese operador en el punto fijo.
Esta es la misma idea de universalidad que aparece en la fisica de transiciones de fase: los detalles microscopicos del sistema son irrelevantes; solo importa la clase de universalidad. Para la duplicacion de periodo, la clase esta determinada por el tipo de maximo del mapa (cuadratico, cubico, etc.). Todos los mapas con maximo cuadratico comparten las mismas \(\delta\) y \(\alpha\).
La universalidad de Feigenbaum conecta la teoria del caos con la fisica estadistica de una manera profunda e inesperada. No es una coincidencia superficial; es la manifestacion de que ambos fenomenos (transiciones de fase y ruta al caos) son gobernados por la misma estructura matematica: el grupo de renormalizacion.
Puntos de bifurcacion del mapa logistico
| n | Periodo 2ⁿ | rₙ | rₙ - rₙ₋₁ | Ratio (→ δ) |
|---|---|---|---|---|
| 1 | 2 | 3.000000 | — | — |
| 2 | 4 | 3.449490 | 0.449490 | — |
| 3 | 8 | 3.544090 | 0.094600 | 4.7514 |
| 4 | 16 | 3.564407 | 0.020317 | 4.6562 |
| 5 | 32 | 3.568759 | 0.004352 | 4.6684 |
| 6 | 64 | 3.569692 | 0.000933 | 4.6686 |
| 7 | 128 | 3.569891 | 0.000200 | 4.6692 |
| ∞ | ∞ | 3.569946... | 0 | 4.6692... |
La convergencia del ratio hacia \(\delta = 4.6692...\) es visible incluso con solo las primeras bifurcaciones. Con cada bifurcacion adicional, ganamos aproximadamente un digito de precision. El valor acumulado \(r_\infty = 3.569946...\) marca el "borde del caos": mas alla de este valor, el periodo es infinito y el sistema es caotico (con las interrupciones periodicas).
"La constante de Feigenbaum es tan universal como π: aparece en sistemas completamente diferentes. Es una constante de la naturaleza, no de un sistema particular. Su universalidad conecta la teoria del caos con la fisica de transiciones de fase a traves del grupo de renormalizacion."
Verificacion experimental
La universalidad de Feigenbaum ha sido verificada experimentalmente en sistemas fisicos reales. En 1981, Albert Libchaber y Jean Maurer midieron la cascada de duplicacion de periodo en una celda de conveccion de Rayleigh-Benard (el mismo fenomeno que inspiro a Lorenz). Encontraron que la razon entre los intervalos de temperatura de las bifurcaciones sucesivas convergia a \(\delta \approx 4.4 \pm 0.1\), consistente con el valor teorico.
Otros experimentos han verificado la constante en diodos Josephson, reacciones quimicas oscilantes (Belousov-Zhabotinsky), oscilaciones acusticas, y circuitos electronicos no lineales. En todos los casos, la misma constante \(\delta \approx 4.669\) emerge, confirmando la prediccion teorica de universalidad.
Ejercicios
-
Calcular bifurcaciones
Usando el mapa logistico \(x_{n+1} = rx_n(1-x_n)\), encuentra numericamente los primeros 4 puntos de bifurcacion \(r_1, r_2, r_3, r_4\). Pista: para encontrar \(r_1\), busca donde el punto fijo \(x^* = 1 - 1/r\) pierde estabilidad (cuando \(|f'(x^*)| = 1\)).
-
Verificar la convergencia
Con los valores de la tabla, calcula las razones \(\frac{r_3 - r_2}{r_4 - r_3}\), \(\frac{r_4 - r_3}{r_5 - r_4}\), etc. ¿A que velocidad convergen a \(\delta\)? ¿Cuantas bifurcaciones necesitas para tener 3 digitos significativos?
-
¿Por que es universal?
Reflexiona: ¿por que la misma constante aparece en el mapa logistico, el mapa seno, y un circuito electronico? ¿Que tienen en comun? ¿Que propiedad del mapa es relevante y cuales son irrelevantes?
LAB: Cascada de Feigenbaum
Diagrama de bifurcacion enfocado en la region de duplicacion de periodo, con anotaciones que marcan cada punto de bifurcacion \(r_n\) y las razones entre intervalos sucesivos convergiendo a \(\delta = 4.669...\). Usa el boton "Animar zoom" para ver como la cascada se repite a escalas cada vez mas finas.