La PizarraMatemática.

Máster y doctorado

Sistemas integrables

Conservar energía ayuda a seguir una trayectoria, pero en un sistema grande no suele bastar. La integrabilidad de Liouville exige suficientes cantidades conservadas, independientes y compatibles. Veremos cada palabra con osciladores antes de llegar a toros y coordenadas acción–ángulo.

1. Un grado de libertad usa una posición y un momento

En coordenadas canónicas (q₁,…,qₙ,p₁,…,pₙ), el espacio de fases tiene dimensión 2n. Un Hamiltoniano H determina \dot q_i=\partial_{p_i}H,\quad\dot p_i=-\partial_{q_i}H. Trabajamos con funciones suaves y estructura simpléctica.

Una trayectoria es una curva en ese espacio, no solo el movimiento de q. Para H autónomo la energía se conserva. Un nivel de energía restringe la curva, pero puede contener muchas trayectorias y no determina toda una dinámica de varios grados.

Un oscilador de frecuencia dos

H=(p²+4q²)/2.

  1. Ecuaciones

    q′=p y p′=−4q.
  2. Candidato

    Con q(0)=1, p(0)=0: q=cos(2t), p=−2sin(2t).
  3. Comprobar

    q′=p, p′=−4q y H=2 para todo t.
Resultado

Posición y momento recorren juntos una elipse.

Lo comprobamos

El signo negativo de p′ produce oscilación, no crecimiento exponencial.

Contar dimensiones

Tres posiciones y tres momentos.

  1. Fase

    Hay seis coordenadas.

  2. Grados

    Se habla de tres grados de libertad.

  3. Integrabilidad

    Buscar tres integrales independientes es el conteo de Liouville, no seis.

Resultado

Grados de libertad y dimensión del espacio de fases son números distintos.

Lo comprobamos

El número de fórmulas escritas no garantiza independencia.

2. El corchete de Poisson comprueba conservación

Definimos \{F,G\}=\sum_i(F_{q_i}G_{p_i}-F_{p_i}G_{q_i}).

La regla de la cadena junto a Hamilton da dF/dt={F,H} cuando F no depende explícitamente de t.

Si {F,H}=0, F es constante a lo largo de ese flujo.
Si F(t,q,p) depende explícitamente de t, dF/dt=\partial_tF+\{F,H\}.
La antisimetría da {H,H}=0 para H autónomo.

No se omite la derivada temporal en una cantidad que sí contiene t.

El corchete básico

F=q, G=p, en un grado.

  1. Primer producto

    F_{qG}_p=1·1=1.
  2. Segundo

    F_{pG}_q=0·0=0.
  3. Restar

    {q,p}=1; al intercambiar, {p,q}=−1.
Resultado

El orden de F y G importa.

Lo comprobamos

No es un producto ordinario conmutativo.

Una cantidad conservada con tiempo explícito

H=p²/2, F=q−tp.

  1. Derivada explícita

    ∂ₜF=−p.
  2. Corchete

    {F,H}=1·p−(−t)·0=p.
  3. Total

    dF/dt=−p+p=0.
Resultado

F recupera la posición inicial de una partícula libre.

Lo comprobamos

Mirar solo su corchete no habría reconocido la conservación.

3. Liouville exige número, independencia e involución

Para n grados se requieren n integrales F₁,…,Fₙ independientes en la región regular y en involución: \{F_i,F_j\}=0. Habitualmente una es H o H depende de ellas de manera pertinente. Independencia significa rango n de sus diferenciales.

Dos expresiones distintas pueden describir la misma restricción.

Si F₂=F₁², dF₂=2F₁dF₁ y no añade rango.

La involución permite flujos compatibles; conservar cantidades sin comprobar sus corchetes no basta.

Dos osciladores separados

H₁=(p₁²+q₁²)/2 y H₂=(p₂²+q₂²)/2; H=H₁+H₂.

  1. Compatibilidad

    Cada H usa variables diferentes, luego {H₁,H₂}=0.
  2. Conservación

    {H₁,H}={H₂,H}=0.
  3. Independencia

    Si ambos pares (qᵢ,pᵢ) son no nulos, sus diferenciales independientes tienen soportes distintos.

Resultado

Hay dos integrales adecuadas en esa región regular.

Lo comprobamos

En un oscilador detenido en el origen baja el rango.

Una lista engañosa

Se proponen H₁, H₁² y 7H₁ para tres restricciones.

  1. Diferenciar

    d(H₁²)=2H₁dH₁ y d(7H₁)=7dH₁.
  2. Rango

    Todos están en la misma dirección de covectores.

  3. Contar

    El rango es como máximo uno, no tres.

Resultado

Reescribir una cantidad no crea nuevas integrales independientes.

Lo comprobamos

Puede bajar incluso a cero en un punto crítico de H₁.

4. Liouville–Arnold describe niveles regulares, compactos y conexos

Si un nivel conjunto no vacío de n integrales en involución es regular, compacto y conexo, el teorema de Liouville–Arnold lo identifica con un toro Tⁿ. En una vecindad adecuada existen coordenadas acción–ángulo.

Regular exige que los diferenciales sean independientes en ese nivel. Compacto excluye escapar indefinidamente dentro del nivel. Conexo selecciona una pieza. No se concluye que todo nivel, ni todo el espacio de fases, sea un único toro con coordenadas globales.

Un toro como producto de dos círculos

H₁=1/2 y H₂=2 para dos osciladores unitarios.

  1. Primero

    q₁²+p₁²=1: círculo de radio uno.
  2. Segundo

    q₂²+p₂²=4: círculo de radio dos.
  3. Producto

    Ambos positivos: nivel compacto, conexo y regular de dimensión dos, un T².

Resultado

El toro de nivel contiene todas las parejas de fases posibles.

Lo comprobamos

No es necesariamente la imagen recorrida por una sola trayectoria.

Integrable sin nivel compacto

Una partícula libre en q∈R, H=p²/2, nivel H=1/2.

  1. Resolver

    p=1 o p=−1.
  2. Cada pieza

    q puede ser cualquier real: dos rectas no compactas.

  3. Teorema

    No se deduce un círculo porque falta compacidad.

Resultado

Ser integrable no obliga a que todas las órbitas sean cerradas.

Lo comprobamos

La hipótesis global no se sustituye por conservar H.

5. Las acciones fijan el nivel y los ángulos avanzan linealmente

En acción–ángulo, H=H(I), \dot I=0,\quad\dot\theta=\omega(I)=\partial_IH. Cada ángulo se toma módulo 2π. Las acciones son constantes y la dinámica consiste en girar con frecuencias que dependen del nivel.

La acción mide el área de una órbita de fase
H=(p²+4q²)/2=2(0,2)(0,−2)(1,0)(−1,0)qpÁrea=π·1·2=2πI=área/(2π)=1
Plano de fase con posición q horizontal y momento p vertical, usando la misma escala en ambos ejes. Para H igual p al cuadrado más cuatro q al cuadrado dividido por dos, el nivel de energía dos es una elipse con semiejes uno en q y dos en p. Se indican los cuatro extremos. La flecha en el extremo q uno p cero apunta hacia momento decreciente, según las ecuaciones q prima igual p y p prima igual menos cuatro q. El área encerrada es dos pi y la acción integral de p dq sobre dos pi vale uno. El nivel no es un círculo con esta misma escala; el origen es crítico y queda fuera de esta órbita.
Para un oscilador de frecuencia ω>0, H=(p²+ω²q²)/2.

Una acción es I=(2\pi)^{-1}\oint p\,dq=E/\omega en un nivel E>0, con el ciclo orientado por el flujo. La integral mide el área encerrada dividida por 2π.

Medir una elipse

ω=2, E=2.

  1. Semiejes

    q máximo=√(2E)/ω=1; p máximo=√(2E)=2.
  2. Área

    π·1·2=2π.
  3. Acción

    I=2π/(2π)=1, también E/ω=1.
Resultado

La acción integra una vuelta completa, no solo una semicircunferencia.

Lo comprobamos

La orientación elegida hace positiva esta integral.

Reconstruir el movimiento

Para I>0, q=√(2I/ω)sinθ y p=√(2Iω)cosθ.

  1. Hamiltoniano

    Al sustituir, H=ωI.
  2. Dinámica

    I′=0 y θ′=ω.
  3. Verificar

    q′=√(2I/ω)cosθ·ω=p y p′=−ω²q.
Resultado

La transformación explica la forma sinusoidal y su frecuencia.

Lo comprobamos

En I=0 el ángulo deja de identificar una posición única.

6. Integrabilidad permite órbitas periódicas y cuasiperiódicas

Con dos ángulos, θ(t)=θ₀+ωt módulo 2π.

Si ambas frecuencias no son cero y su razón es racional, hay un período común. Si su razón es irracional, no hay período positivo y la órbita es densa en el toro bidimensional.

Una razón racional cierra y una irracional no tiene período
Ángulos módulo 2πω=(1,2): cierraθ₁θ₂ω=(1,√2): no cierraθ₁θ₂Bordes opuestos identificados.
Dos cuadrados representan el mismo tipo de toro de dos ángulos, con lados opuestos identificados. Arriba las frecuencias son uno y dos: tras tiempo dos pi la trayectoria vuelve a la fase inicial. Abajo son uno y raíz de dos: se muestra solo un tramo finito de cuatro vueltas del primer ángulo, sin afirmar que ese tramo haya llenado el toro. Las líneas que llegan a un borde reaparecen por el opuesto; los cortes del dibujo no son saltos de la trayectoria en el toro. La razón irracional impide un período y hace la órbita completa densa. Las escalas angulares son iguales en ambos cuadrados.

En más dimensiones importan todas las relaciones enteras k·ω=0. Pueden confinar la órbita a un subtoro. Densidad no significa recorrer cada punto en un tiempo finito ni hace que la trayectoria se vuelva aleatoria.

Encontrar un período común

ω=(1,2).

  1. Primer ángulo

    En T=2π avanza 2π.
  2. Segundo

    Avanza 4π, dos vueltas.

  3. Cerrar

    Ambos vuelven a la fase inicial y la órbita es periódica.

Resultado

El toro completo tiene más puntos que esa órbita cerrada.

Lo comprobamos

Los dos ángulos se interpretan módulo 2π.

Probar que no cierra

ω=(1,√2).

  1. Si cerrara

    T=2πm y √2T=2πn, con enteros y m≠0.
  2. Dividir

    √2=n/m sería racional.
  3. Contradicción

    No existe tal T; la órbita es cuasiperiódica y densa en T².

Resultado

El sistema puede ser perfectamente integrable sin repetir su posición de fase.

Lo comprobamos

En tres frecuencias hay que comprobar relaciones enteras, no solo una pareja.

7. Un nivel crítico puede cambiar de dimensión o contener una separatriz

Para H=(p²+q²)/2, dH=(q,p).

En niveles positivos no se anula; en E=0 solo queda el origen y el diferencial es cero. No se aplican allí las coordenadas angulares de los círculos positivos como si nada cambiara.

En sistemas más ricos, los niveles críticos pueden contener equilibrios y separatrices. El teorema local regular no describe automáticamente esas piezas. La topología global también puede impedir una única elección de acciones y ángulos en todo el sistema.

Un círculo que colapsa

Dos osciladores: H₁=0 y H₂=1/2.

  1. Primero

    q₁=p₁=0: no queda círculo.
  2. Segundo

    q₂²+p₂²=1: queda un círculo.
  3. Rango

    dH₁=0; el nivel no cumple rango dos.
Resultado

El producto es un punto por un círculo, no el toro regular T².

Lo comprobamos

Que aún haya movimiento no restaura la hipótesis perdida.

El péndulo en su punto alto

H=p²/2+1−cos q, con q módulo 2π.

  1. Equilibrio

    q=π, p=0.
  2. Energía

    H=2 y dH=(sin q,p)=(0,0).
  3. Consecuencia

    Ese nivel incluye una pieza crítica; no es un nivel regular al que aplicar directamente Liouville–Arnold.

Resultado

El comportamiento próximo a una separatriz exige otro análisis.

Lo comprobamos

No basta comprobar solo que el valor de energía sea positivo.

8. Lax y las perturbaciones añaden herramientas, no garantías automáticas

Si una matriz cumple \dot L=[B,L], entonces tr(Lᵏ) se conserva por la propiedad cíclica de la traza. Son candidatas a integrales. Para deducir integrabilidad de Liouville aún hay que comprobar independencia, involución y su relación con la estructura de Poisson.

KAM estudia persistencia de ciertos toros bajo pequeñas perturbaciones con condiciones de regularidad, no degeneración y frecuencias. No afirma que todo toro sobreviva ni que toda perturbación destruya integrabilidad. En EDP, “integrable” puede referirse a estructuras de dimensión infinita con otro marco.

Una evolución por conjugación

L(t)=R(t)diag(1,3)R(t)ᵀ, donde R es una rotación y B=R′Rᵀ.

  1. Derivar

    Como (Rᵀ)′=−RᵀB, queda L′=BL−LB.
  2. Traza

    tr L=1+3=4.
  3. Otras

    det L=3 y tr(L²)=1+9=10 para todo t.
Resultado

La ecuación conserva autovalores aunque cambien entradas.

Lo comprobamos

Esos números por sí solos no verifican todas las hipótesis de Liouville.

Una perturbación que sigue siendo integrable

H₀(I)=(I₁²+I₂²)/2; añade εI₁I₂.

  1. Ángulos

    El nuevo H sigue sin depender de θ.

  2. Acciones

    I₁ e I₂ continúan constantes y en involución.

  3. Frecuencias

    Pasan a (I₁+εI₂,I₂+εI₁).

Resultado

Perturbar un Hamiltoniano no implica automáticamente caos.

Lo comprobamos

Para una perturbación que dependa de ángulos el análisis puede ser muy distinto.

Ahora lo resolvemos de principio a fin.

EJEMPLO RESUELTO 1

Un nivel de energía

Para H=(p²+q²)/2 y H=2, describe el nivel.

  1. Hamiltoniano

    H=(p²+q²)/2.
  2. Nivel

    H=2 equivale a p²+q²=4.
  3. Geometría

    Es un círculo centrado en el origen.

  4. Radio

    r²=4, luego r=2.
  5. Regularidad

    dH=(q,p) no se anula en ese círculo.
  6. Comprobar

    El origen tiene energía cero y no pertenece al nivel.

Resultado

El nivel regular es un círculo de radio 2.

Lo comprobamos

El origen H=0 es crítico y se excluye de ese nivel regular.

EJEMPLO RESUELTO 2

Dos energías compatibles

Con H₁=(p₁²+q₁²)/2,H₂=(p₂²+q₂²)/2, calcula {H₁,H₂}.

  1. Variables

    H₁ usa solo q₁,p₁; H₂ solo q₂,p₂.

  2. Definir

    El corchete suma un producto cruzado por cada índice.

  3. Índice uno

    Las dos derivadas de H₂ respecto a q₁,p₁ son cero.

  4. Índice dos

    Las dos derivadas de H₁ respecto a q₂,p₂ son cero.

  5. Sumar

    Todos los términos son cero: {H₁,H₂}=0.
  6. Comprobar

    La involución está probada; independencia requiere además revisar los diferenciales en el nivel.

Resultado

El corchete es 0.

Lo comprobamos

La separación de variables garantiza involución.

Los símbolos que hemos usado

n grados de libertad
Espacio de fases de dimensión2n en el marco hamiltoniano.
Primera integral F
Función constante a lo largo de las trayectorias.
{F,G}
Corchete de Poisson que mide compatibilidad dinámica.
Independencia funcional
Los diferenciales de las integrales son linealmente independientes.
Involución
Todos los corchetes de un conjunto de integrales son cero.
Nivel regular compacto
Nivel conjunto con independencia de diferenciales y compacidad.
I,θ
Acciones constantes y ángulos periódicos de una descripción local.

Ahora te toca a ti.

Tienes cinco ejercicios. Inténtalos; las pistas te ayudan si te atascas.

Practicar este tema →

Para ampliar