Bifurcaciones y organización paramétrica

Cambios cualitativos en una bifurcación

Un parámetro puede modificar la posición de una órbita sin alterar la organización del sistema. Una bifurcación ocurre cuando esa organización sí cambia: aparecen o desaparecen equilibrios, una rama pierde estabilidad, aparece un ciclo o se transforma la estructura de las órbitas recurrentes. El valor crítico no se identifica solo porque una gráfica “se vea distinta”. Hay que localizar la solución implicada, calcular su estabilidad y comprobar que el cambio persiste al mejorar el muestreo.

Consideremos una familia autónoma

$$\dot{\mathbf{x}}=\mathbf{f}(\mathbf{x},\mu).$$

Una rama de equilibrios satisface $\mathbf{f}(\mathbf{x}^*(\mu),\mu)=0$. Si el jacobiano $D_{\mathbf{x}}\mathbf{f}$ es invertible, el teorema de la función implícita permite continuar localmente esa rama: pequeños cambios de $\mu$ producen pequeños cambios de $\mathbf{x}^*$. Una bifurcación local exige que esa invertibilidad falle; en una dimensión, que $f_x(x^*,\mu)=0$. Esta condición señala un candidato, no completa el análisis. También se revisan los términos no lineales que determinan el mecanismo.

Las formas normales conservan los términos mínimos necesarios después de cambios suaves de coordenadas y parámetros. No pretenden reproducir toda la física lejos del punto crítico. Strogatz las utiliza para separar la geometría universal del detalle particular del modelo Strogatz, 2015 Steven H. Strogatz (2015) Nonlinear Dynamics and Chaos: With Applications to Physics, Biology, Chemistry, and Engineering 2.ª ed. · Westview Press Ubicación consultada: cap. 3, §§3.1–3.6 .

Silla-nodo, transcrítica y horquilla

La forma normal silla-nodo es

$$\dot x=\mu-x^2.$$

Los equilibrios cumplen $x_\pm^*=\pm\sqrt{\mu}$ cuando $\mu>0$. Para $\mu<0$ no hay equilibrios reales y en $\mu=0$ ambos coinciden en el origen. Como $f_x=-2x$, la rama positiva es estable y la negativa es inestable. La distancia

$$x_+^*-x_-^*=2\sqrt{\mu}$$

obedece una ley de raíz cuadrada. Esta ley es una firma local: si una medición se ajusta a ella en una ventana pequeña, apoya el mecanismo silla-nodo, pero no demuestra que el modelo completo sea exactamente la forma normal.

En la transcrítica,

$$\dot x=\mu x-x^2=x(\mu-x),$$

las ramas $x^*=0$ y $x^*=\mu$ existen a ambos lados del valor crítico y intercambian estabilidad. Los multiplicadores lineales son $\mu$ en el origen y $-\mu$ en la segunda rama. En la horquilla supercrítica,

$$\dot x=\mu x-x^3=x(\mu-x^2),$$

el origen es estable para $\mu<0$ e inestable para $\mu>0$; entonces aparecen dos ramas estables $x_\pm^*=\pm\sqrt{\mu}$. La simetría $x\mapsto-x$ obliga a que nazcan juntas. Si el sistema físico rompe esa simetría, un término constante o cuadrático deforma la horquilla y puede separar los puntos de giro.

El potencial permite verificar la estabilidad sin memorizar el diagrama. Para

$$V_\mu(x)=-\frac{\mu}{2}x^2+\frac14x^4,$$

se tiene $\dot x=-V_\mu'(x)$ y

$$\frac{dV_\mu}{dt}=V_\mu'(x)\dot x=-[V_\mu'(x)]^2\leq0.$$

Los mínimos de $V_\mu$ son atractores del sistema gradiente. Esta derivación explica por qué el origen deja de ser preferido y aparecen dos estados estables; la animación hace visible el mismo cambio.

Pérdida de una rama estable en catástrofes mecánicas

En dinámica no lineal, catástrofe no significa que la ecuación explote. Significa que la respuesta estable que el sistema seguía termina en un pliegue y ya no existe para el siguiente valor del control. El salto observado puede ser grande aunque el parámetro cambie muy poco. Thompson y Stewart distinguen este cambio discontinuo de las bifurcaciones en las que una rama atractora continúa sin salto, y lo sitúan en el espacio conjunto de estado y controles. Thompson, 2002 J. M. T. Thompson y H. B. Stewart (2002) Nonlinear Dynamics and Chaos 2.ª ed. · John Wiley & Sons Ubicación consultada: cap. 13, §§13.1–13.3, pp. 249–284

El potencial de cúspide concentra el mecanismo en una ecuación:

$$V(x;a,b)=\frac{x^4}{4}+\frac{a x^2}{2}+bx, \qquad \dot x=-\frac{\partial V}{\partial x}=-(x^3+ax+b).$$

La coordenada $x$ puede representar la deflexión de una pieza elástica; $a$ modifica su rigidez efectiva y $b$ actúa como carga o imperfección. Los equilibrios satisfacen

$$x^3+ax+b=0,$$

y un pliegue exige además $V''=3x^2+a=0$. Eliminando $x$ se obtiene la curva de valores críticos

$$a=-3x^2,\qquad b=2x^3,\qquad 4a^3+27b^2=0.$$

Cuando $4a^3+27b^2<0$ hay tres equilibrios: dos mínimos del potencial separados por un máximo. Fuera de esa región queda uno. La cúspide no es una silueta decorativa; es la frontera entre esas dos organizaciones.

Fijemos $a=-3$. Los pliegues aparecen en $x=\pm1$ y $b=\pm2$. Para $b=0$, los equilibrios son $x=0,\pm\sqrt3$: los laterales son mínimos y el central es máximo. Si $b$ aumenta lentamente desde cero mientras el sistema ocupa el mínimo positivo, ese mínimo se aproxima a $x=1$ y desaparece en $b=2$ junto con el máximo. La trayectoria ya no puede permanecer allí y cae hacia el mínimo restante, cercano a $x=-2$. Al invertir el barrido, el salto de regreso ocurre en el otro pliegue; de ahí surge una histéresis mecánica.

La figura de la horquilla corresponde a la rebanada simétrica $b=0$ con $a=-\mu$. Añadir $b\neq0$ inclina el potencial: los dos pozos dejan de ser equivalentes y la horquilla perfecta se despliega en pliegues. Esta lectura conecta la animación del potencial con una pieza real: una imperfección geométrica pequeña no se limita a mover el umbral, también puede seleccionar qué rama sobrevive.

Pérdida de estabilidad de ciclos

Para un mapa $x_{n+1}=F_\mu(x_n)$, un punto fijo cambia de estabilidad cuando su multiplicador cruza el círculo unidad. En una bifurcación de duplicación de periodo, el multiplicador atraviesa $-1$: la perturbación alterna de signo y deja de contraerse. Entonces aparece un ciclo de periodo dos. Si ese ciclo vuelve a cruzar $-1$ como punto fijo de $F_\mu^2$, aparece un ciclo de periodo cuatro.

El diagrama logístico muestra esta secuencia. Para cada $r$ se descarta un transitorio y se dibujan varios estados posteriores. Una rama, dos ramas y cuatro ramas corresponden a respuestas que deben verificarse calculando su periodo y su multiplicador. Las bandas densas contienen ventanas periódicas y estructuras más finas que la resolución de la figura puede ocultar. Alligood, Sauer y Yorke presentan el diagrama como síntesis de muchas órbitas y conectan sus ramas con bifurcaciones genéricas Alligood, 1996 Kathleen T. Alligood, Tim D. Sauer y James A. Yorke (1996) Chaos: An Introduction to Dynamical Systems Springer Ubicación consultada: cap. 11, §§11.1–11.4, y cap. 12 .

La parte inferior de una figura combinada suele mostrar un exponente de Lyapunov a tiempo finito. Un valor negativo es compatible con atracción periódica; uno positivo, con expansión promedio. Cerca de una bifurcación la convergencia es lenta y el exponente puede quedar próximo a cero. Un color o una curva no convierte automáticamente esa estimación en un resultado asintótico.

Amplitud y estabilidad en la forma normal de Hopf

Consideremos el oscilador radial reducido

$$\dot r=\mu r-r^3,\qquad r\geq0.$$

Además de $r=0$, para $\mu>0$ existe $r^*=\sqrt{\mu}$. La derivada $g_r=\mu-3r^2$ vale $\mu$ en el origen y $-2\mu$ en la rama no nula. Así, el origen atrae para $\mu<0$ y la nueva amplitud atrae para $\mu>0$. Si este modelo radial proviene de una bifurcación de Hopf supercrítica, $r^*$ representa la amplitud de un ciclo límite.

La ley $r^*\sim\sqrt{\mu}$ ofrece una comprobación cuantitativa. Supongamos que un barrido numérico produce amplitudes máximas $A(\mu)$. Para probar la ley se grafica $A^2$ contra $\mu$ cerca del umbral; una relación aproximadamente lineal es más informativa que observar que “la oscilación crece”. También se verifica que el periodo permanezca finito y que el resultado no dependa del transitorio usado para medir la amplitud.

Una Hopf no puede justificarse solo porque aparece una oscilación. En el equilibrio debe cruzar el eje imaginario un par complejo conjugado, mientras los demás autovalores mantienen parte real distinta de cero. La dirección supercrítica o subcrítica depende del término no lineal efectivo. Esta distinción importa: en el caso subcrítico puede existir una región de coexistencia y saltos, con consecuencias visibles al barrer el parámetro en sentidos opuestos.

Construcción de barridos paramétricos

Antes de ejecutar, se fija qué objeto ocupará el eje vertical:

• estados retenidos de un mapa;
• máximos locales de una componente;
• cruces orientados de una sección;
• amplitud de un ciclo;
• valor de un observable promedio.

Estas elecciones no son equivalentes. Muestrear una serie a tiempos equiespaciados puede llenar una curva periódica con muchos puntos; guardar máximos produce una rama por pico distinto. Para cada parámetro también se decide si la condición inicial se reinicia o si se usa el estado final anterior como continuación. La continuación detecta ramas que un reinicio podría no alcanzar, pero puede conservar memoria del sentido del barrido cuando coexisten atractores.

Un protocolo básico para el mapa logístico usa $r\in[2.8,4]$ con al menos $1200$ valores, $1500$ iteraciones por valor, descarte de $1000$ y $200$ estados retenidos. Luego se repite una ventana pequeña con malla más fina. Para un flujo se declaran integrador, paso, tolerancias, duración, transitorio y algoritmo de detección de máximos o cruces. Se prueban al menos dos pasos y se aumenta el tiempo cerca de regiones con convergencia lenta.

El diagrama resultante es evidencia finita. Una rama que cambia al duplicar el transitorio era probablemente transitoria; una rama que aparece solo al usar continuación puede indicar coexistencia, pero se confirma con condiciones iniciales independientes. La discretización del parámetro también importa: una ventana estrecha puede desaparecer entre dos valores de la malla.

Interpretación de diagramas de bifurcación

Una explicación completa enlaza cuatro piezas:

1. Objeto: qué equilibrio, ciclo u observable se sigue.
2. Mecanismo: qué autovalor o multiplicador cruza el umbral.
3. Escala: cómo crece la rama cerca del valor crítico.
4. Protocolo: cómo se obtuvo la figura y qué controles se repitieron.

Devaney desarrolla la bifurcación como una combinación de experimento y argumento, especialmente en mapas cuadráticos Devaney, 2020 Robert L. Devaney (2020) A First Course in Chaotic Dynamical Systems: Theory and Experiment 2.ª ed. · CRC Press Ubicación consultada: caps. 6 y 8 . La figura ayuda a localizar valores interesantes; el cálculo de estabilidad explica por qué cambia la órbita; el refinamiento numérico delimita cuánto podemos afirmar sobre el modelo implementado.

Ejercicios sobre formas normales y diagramas paramétricos

1. Dibuja las ramas de silla-nodo, transcrítica y horquilla usando sus ecuaciones, y anota la estabilidad calculada en cada tramo.
2. Deriva la ley de raíz cuadrada para la amplitud de la horquilla y explica en qué intervalo esperarías verla en un sistema físico.
3. Para la forma $\dot x=\mu x+x^3$, clasifica las ramas y compárala con el caso supercrítico.
4. Localiza numéricamente el primer desdoblamiento del mapa logístico usando el multiplicador, no solo el diagrama.
5. Diseña una comparación entre reinicio y continuación en un sistema con coexistencia. Indica cómo registrarás el sentido del barrido.
6. Elige una región densa de un diagrama y redacta tres pruebas que separen periodo largo, transitorio y respuesta aperiódica.

La caracterización de una bifurcación incluye la solución, el cambio de estabilidad, la amplitud de las ramas y el protocolo numérico que respalda esa lectura.