Contenido: Tema 10
10. Teoría canónica de perturbaciones
10.1 Transformaciones canónicas y teoría de perturbaciones
10.2 Teoría de perturbaciones para sistemas unidimensionales
10.3 Teoría de perturbaciones para sistemas de n-dimensiones
10.4 Invariantes adiabáticos
1
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 1/35
/35
Contenido: Tema 10
10. Teoría canónica de perturbaciones
10.1 Transformaciones canónicas y teoría de perturbaciones
10.2 Teoría de perturbaciones para sistemas unidimensionales
10.3 Teoría de perturbaciones para sistemas de n-dimensiones
10.4 Invariantes adiabáticos
2
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 2/35
/35
Transformaciones canónicas y teoría de perturbaciones
Fundamentos
Supongamos que tenemos el Hamiltoniano,
H(q, p, t, ) = H0 (q, p, t) + H1 (q, p, t),
en donde:
• H0 : es el Hamiltoniano del sistema no-perturbado.
• H1 : es la perturbación al Hamiltoniano, siendo la diferencia H −H0
muy pequeña, controlada por el parámetro .
Entonces, se puede considerar que existe una transformación canónica
dependiente de tal que:
(q, p, t) → (Q, P, t) , ∀ H(q, p, t) → K(Q, P, t)
en donde K(Q, P, t) sea mas sencillo de resolver que H(q, p, t).
En este enfoque se tendrá:
• = 0: la transformación canónica es la identidad.
• 6= 0: una transformación canónica que nos lleva a K(Q, P, t).
3
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 3/35
/35
Transformaciones canónicas y teoría de perturbaciones
Transformación canónica via S(q, P, t)
Consideremos una transformación canónica tipo S(q, P, t; ) en donde
es fijo y muy pequeño, entonces:
∂S(q, P, t) 1
K(Q, P, t) = H(q, p, t) + ,
∂t
en donde q = q(Q, P, t), p = p(Q, P, t).
Ahora, expandiendo q, p, K, y S en potencias de :
K = K0 + K1 + 2 K2 + . . .
pi = Pi + p1,i + 2 p2,i + . . .
qi = Qi + q1,i + 2 q2,i + . . .
S = qi Pi + S1 + 2 S2 + . . .
en donde S(q, P, t) = qi Pi representa la función generatriz relativa a
la transformación identidad.
1
en donde se ha escrito S(q, P, t; ) = S(q, P, t) por simplicidad. 4
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 4/35
/35
Transformaciones canónicas y teoría de perturbaciones
Transformación canónica via S(q, P, t)
Pero de la transformación canónica se sabe:
∂S ∂S
pi = & Qi = ,
∂qi ∂Pi
entonces, sustituyendo S por su expansión en pi ,
∂S ∂
pi = = qi Pi + S1 + 2 S2 + . . . ,
∂qi ∂qi
∂S1 ∂S2
= Pi + + 2 + ...,
∂qi ∂qi
pero,
pi = Pi + p1,i + 2 p2,i + . . . ,
comparando término a término las últimas dos ecuaciones se tiene:
∂Sk
pk,i = .
∂qi
5
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 5/35
/35
Transformaciones canónicas y teoría de perturbaciones
Transformación canónica via S(q, P, t)
Realizando el mismo análisis con Qi ,
∂S ∂
Qi = = qi Pi + S1 + 2 S2 + . . . ,
∂Pi ∂Pi
∂S1 ∂S2
= qi + + 2 + ...,
∂Pi ∂Pi
pero,
qi = Qi + q1,i + 2 q2,i + . . . ,
sustituyendo la ecuación anterior en Qi :
∂S1 ∂S2
2
Qi = Qi + q1,i + + q2,i + + ...,
∂Pi ∂Pi
para que la ecuación anterior se mantenga, se debe cumplir:
∂Sk
qk,i = − .
∂Pi
6
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 6/35
/35
Transformaciones canónicas y teoría de perturbaciones
Teoría perturbativa
Expandiendo ahora al Hamiltoniano H(q, p, t) a primer orden,
H(q, p, t) = H0 (q, p, t) + H1 (q, p, t),
para ello se aplica Taylor en H0 respecto (Q, P ),2
∂H0 ∂H0
H0 (p, q, t) = H0 (Q, P, t) + (qi − Qi ) + (pi − Pi ) ,
∂qi =0
∂pi =0
donde, qi = Qi + q1,i , pi = Pi + p1,i ,
∂H0 ∂H0
⇒ H0 (p, q, t) = H0 (Q, P, t) + q1,i + p1,i ,
∂Qi ∂Pi
∂Sk ∂Sk
pero: qk,i = − & pk,i = ,
∂Pi ∂qi
∂S1 ∂H0 ∂S1 ∂H0
⇒ H0 (p, q, t) = H0 (Q, P, t) + − + ,
∂Pi ∂Qi ∂Qi ∂Pi
∴ H0 (p, q, t) = H0 (Q, P, t) + [S1 , H0 ].
2
f (x, y) = f (x0 , y0 ) + (x − x0 )∂f /∂x|x0 ,y0 + (y − y0 )∂f /∂y|x0 ,y0 + O(2) 7
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física
/35
7/35
Transformaciones canónicas y teoría de perturbaciones
Teoría perturbativa
Para H1 (q, p, t) se realiza un desarrollo similar,
H1 (q, p, t) = H1 (Q, P, t)+O(), ⇒ H1 (q, p, t) = H1 (Q, P, t)+O(2 ).
Por tanto, sustituyendo en H(q, p, t),
H(q, p, t) = H0 (q, p, t) + H1 (q, p, t),
∴ H(q, p, t) = H0 (Q, P, t) + (H1 (Q, P, t) + [S1 , H0 ]) .
Ahora, recordando la transformación canónica para el Hamiltoniano,
∂S
K(Q, P, t) = H(q, p, t) + ,
∂t
entonces se sustituye la expresión para H(q, p, t),
∂S
K(Q, P, t) = H0 (Q, P, t) + (H1 (Q, P, t) + [S1 , H0 ]) + ,
∂t
∂S1
= H0 (Q, P, t) + H1 (Q, P, t) + [S1 , H0 ] + ,
∂t
donde: S = qi Pi + S1 → ∂S/∂t = ∂S1 /∂t.
8
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 8/35
/35
Transformaciones canónicas y teoría de perturbaciones
Teoía perturbativa
Recordando finalmente la expresión de K(Q, P, t) en potencias de ,
K(Q, P, t) = K0 + K1 + 2 K2 + . . . ,
y comparando con el resultado anterior,
∂S1
K(Q, P, t) = H0 (Q, P, t)+ H1 (Q, P, t) + [S1 , H0 ] + +O(2 ),
∂t
se llega al siguiente resultado:
K0 (Q, P, t) = H0 (Q, P, t),
∂S1
K1 (Q, P, t) = H1 + [S1 , H0 ] + ,
∂t
La ecuación anterior muestra como se obtiene K1 en función de S1 ,
sin embargo en problemas de dinámica ambas cantidades necesitan ser
encontradas de manera simultánea.
9
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 9/35
/35
Contenido: Tema 10
10. Teoría canónica de perturbaciones
10.1 Transformaciones canónicas y teoría de perturbaciones
10.2 Teoría de perturbaciones para sistemas unidimensionales
10.3 Teoría de perturbaciones para sistemas de n-dimensiones
10.4 Invariantes adiabáticos
10
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 10/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Consideremos que el Hamiltoniano es una cantidad conservada ⇒
H(q, p, t) = H(q, p), en donde el H perturbado viene dado como:
H(q, p) = H0 (q, p) + H1 (q, p).
Asumiendo además que se tratan sistemas de naturaleza periódica,
entonces un set adecuado de variables canónicas para describir a H
son las variables angulares y de acción:
(q, p) → (w, J) ∀ H(q, p) ≡ K(J) = E(J).
Para hallar tal transformación canónica, primero se procede a obtener:
(q, p) → (w0 , J0 ) ∀ K0 (w0 , J0 ) = H0 [q(w0 , J0 ), p(w0 , J0 )] = K0 (J0 ),
definiendo:
K1 (w0 , J0 ) = H1 [q(w0 , J0 ), p(w0 , J0 )] ,
∀ K(w0 , J0 ) = K0 (J0 ) + K1 (w0 , J0 ),
siendo (w0 , J0 ) las variables angulares y de acción de H0 . 11
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 11/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Del resultado anterior:
K(w0 , J0 ) = K0 (J0 ) + K1 (w0 , J0 ),
es imporante enfatizar que K(w0 , J0 ), K0 (J0 ), y K1 (w0 , J0 ) son can-
tidades conocidas:
• H(p, q) es conocido del propio problema a tratar.
• La transformación canónica (q, p) → (w0 , J0 ) se obtiene al aplicar
el procedimiento de VAA al sistema de H0 .
Ahora, recordando que K es integrable y con movimiento periódico, y
por tanto que posee un set de variables (w, J) propias ⇒ existe una
transformación canónica tal que:
(w0 , J0 ) → (w, J) ∀ K [w0 (w, J), J0 (w, J)] ≡ K(J) = E(J),
encontrando tal transformación canónica, y relacionándola con la transf.
(q, p) → (w0 , J0 ) el problema puede ser resuelto, es decir, es posible
obtener E(J). 12
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física
/35
12/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Para resolver (w0 , J0 ) → (w, J), se aplica la siguiente transf.:3
S(w0 , J) = w0 J + S1 (w0 , J) + 2 S2 (w0 , J) + . . . ,
en donde w0 J es la transformación identidad.
De la transformación canónica se conocen las siguientes relaciones:
∂S ∂S
J0 = , w= ,
∂w0 ∂J
en las cuales se sustituye S expresada como una expansión:
∂S ∂S1 ∂S2
J0 = =J + + 2 + ...,
∂w0 ∂w0 ∂w0
∂S ∂S1 ∂S2
w= = w0 + + 2 + ...,
∂J ∂J ∂J
3
del tipo 2: F2 = S(q, P ). 13
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 13/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
De igual manera se expande E(J) en potencias de ,
E(J) = E0 (J) + E1 (J) + 2 E2 (J) + . . . ,
= K(w0 , J0 ) ≡ K0 (J0 ) + K1 (w0 , J0 ).
Para describir cada término de E(J), también se expande en potencias
de (J0 − J) a K(w0 , J0 ):
K(w0 , J0 ) = K0 (J0 ) + K1 (w0 , J0 ),
∂K0
2
1 2 ∂ K0
= K0 (J) + (J0 − J) + (J0 − J) + ...
∂J0 J 2 ∂J02 J
∂K1
. . . + K1 (w0 , J) + (J0 − J) + ...
∂J0 J
en donde:
2
∂S1 ∂S2 ∂S1
J0 − J = + 2 + . . . & (J0 − J)2 = 2 + ...
∂w0 ∂w0 ∂w0
14
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 14/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Sustituyendo las expresiones de (J0 − J),
∂S1 ∂S2 ∂K0
K(w0 , J0 ) = K0 (J) + + 2 + ... + ...
∂w0 ∂w0 ∂J
" #
1 2 ∂S1 2 ∂ 2 K0
... + + ... + ...
2 ∂w0 ∂J 2
∂S1 ∂S2 ∂K1
. . . + K1 (w0 , J) + + 2 + ... + ...
∂w0 ∂w0 ∂J
agrupando en potencias de ,
∂S1 ∂K0
K(w0 , J0 ) = K0 (J) + + K1 (w0 , J) + . . .
∂w0 ∂J
" #
2 ∂S2 ∂K0 1 ∂S1 2 ∂ 2 K0 ∂S1 ∂K1
... + + + + ...
∂w0 ∂J 2 ∂w0 ∂J 2 ∂w0 ∂J
15
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 15/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Comparando el resultado anterior con la expansión de E(J) = K(w0 , J0 ),
E(J) = E0 (J) + E1 (J) + 2 E2 (J) + . . .
∂S1 ∂K0
K(w0 , J0 ) = K0 (J) + + K1 (w0 , J) + . . .
∂w0 ∂J
" #
2 ∂S2 ∂K0 1 ∂S1 2 ∂ 2 K0 ∂S1 ∂K1
... + + + + ...
∂w0 ∂J 2 ∂w0 ∂J 2 ∂w0 ∂J
e igualando término a término hasta segundo grado de , se obtiene:
E0 (J) = K0 (J),
∂S1 ∂K0
E1 (J) = K1 (w0 , J) + = K1 (w0 , J) + [S1 , K0 ] ,
∂w0 ∂J
∂S2 ∂K0 1 ∂S1 2 ∂ 2 K0 ∂S1 ∂K1
E2 (J) = + 2
+ ,
∂w0 ∂J 2 ∂w0 ∂J ∂w0 ∂J
donde es necesario conocer simultaneamente Ek (J) y ∂Sk /∂w0 , debi-
endo determinar cuatro cantidades, de solo dos ecuaciones. 16
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física
/35
16/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Para obtener las expresiones de Ek (J) se evalúan los promedios en un
periodo del movimiento, respecto a w0 :
Z 2π
dw0
hf (w0 )i =
f (w0 ),
2π 0
en donde el promedio es realizado a un J fijo, lo cual conlleva que cada
una de las derivadas ∂Sk /∂w0 sean periódicas en w0 .
Debido a esta periodicidad en las derivadas de Sk , es que se puede
concluir que Sk es igual periódica, por lo que es posible expresarlas
(función y derivadas) en términos de series de Fourier,
∞
X
Sk (w0 , J) = Sk (J; m)eimw0 ,
m=−∞
con lo cual, los promedios de ∂Sk /∂w0 se anularán:
Z 2π
∂Sk 1 ∂Sk
= dw0 = 0.
∂w0 2π 0 ∂w0
17
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 17/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Calculando ahora el promedio de las expresiones para Ek (J), comen-
zando con E1 (J):
∂S1 ∂K0
E1 (J) = K1 (w0 , J) + ,
∂w0 ∂J
∂S1 ∂K0 4
⇒ E1 (J) = hK1 (w0 , J)i + ,
∂w0 ∂J
∴ E1 (J) = hK1 (w0 , J)i , 5
despejando de E1 (J) a ∂S1 /∂w0 ,
∂S1 hK1 i − K1 ∂K0
= ∀ ν0 (J) = ,
∂w0 ν0 (J) ∂J
La solución de la ec. anterior es S1 (w0 , J) + σ(J), donde σ(J) es una
función arbitraria cuyo papel es añadir una fase a w, ∴ σ(J) = 0.
4
hf (J)i = f (J).
5
h∂Sk /∂w0 i = 0. 18
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 18/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformaciót canónica a variables de ángulo y acción
Aplicando ahora el promedio para E2 (J),
∂S2 ∂K0 1 ∂S1 2 ∂ 2 K0 ∂S1 ∂K1
E2 (J) = + 2
+ ,
∂w0 ∂J 2 ∂w0 ∂J ∂w0 ∂J
* +
∂S2 1 ∂ν0 ∂S1 2 ∂S1 ∂K1
⇒ E2 (J) = ν0 + + ,
∂w0 2 ∂J ∂w0 ∂w0 ∂J
pero se sabe que:
∂Sk ∂S1 hK1 i − K1
=0 & = ,
∂w0 ∂w0 ν0
* 2 +
∂S1 1 D E
por tanto: = hK1 i2 − 2K1 hK1 i + K12 ,
∂w0 ν02
1 h D Ei
= 2 hK1 i2 − 2 hK1 i2 + K12 ,
ν0
1 hD E i
= 2 K12 − hK1 i2 .
ν0 19
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 19/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
También:
∂S1 ∂K1 1 ∂K1
= (hK1 i − K1 ) ,
∂w0 ∂J ν0 ∂J
1 ∂K1 ∂K1
= hK1 i − K1 ,
ν0 ∂J ∂J
sustituyendo ambos resultados en E2 (J),
1 ∂ν0 hD 2 E i 1 ∂K ∂K1
2 1
E2 (J) = 2 K1 − hK1 i + hK1 i − K1 .
2ν0 ∂J ν0 ∂J ∂J
Entonces, de la expresión original de E2 (J):
2 2
∂S2 ∂K0 1 ∂S1 ∂ K0 ∂S1 ∂K1
E2 (J) = + + ,
∂w0 ∂J 2 ∂w0 ∂J 2 ∂w0 ∂J
y con el resultado anterior, es posible despejar ∂S2 /∂w0 .
20
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 20/35
/35
Teoría de perturbaciones para sistemas unidimensionales
Transformación canónica a variables de ángulo y acción
Realizando el despeje e insertando la expresión para E2 (J) se obtiene:
∂S2 1 ∂K1 ∂K1 ∂K1 ∂K1
= 2 hK1 i − K1 − hK1 i + K1 + . . .
∂w0 ν0 ∂J ∂J ∂J ∂J
1 ∂ν0 hD 2 E i
... + K1 − 2 hK1 i2 + 2K1 hK1 i − K12 .
2ν0 ∂J
Se observa, por tanto, que el requerimiento de que E = E(J) nos lleva
a ecuaciones diferenciales parciales para S1 y S2 .
Con las expresiones para Ek (J) obtenidas anteriormente se obtiene:
E(J) = E0 (J) + E1 (J) + 2 E2 (J) + . . .
2 1 1 ∂ν0 hD 2 E
i
E(J) = K0 (J) + hK1 i + K1 − hK1 i2 + . . .
ν 2ν0 ∂J
0
∂K1 ∂K1
... + hK1 i − K1 ,
∂J ∂J
en donde las frecuencias del sistema perturbado son ν = ∂E/∂J. 21
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 21/35
/35
Contenido: Tema 10
10. Teoría canónica de perturbaciones
10.1 Transformaciones canónicas y teoría de perturbaciones
10.2 Teoría de perturbaciones para sistemas unidimensionales
10.3 Teoría de perturbaciones para sistemas de n-dimensiones
10.4 Invariantes adiabáticos
22
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 22/35
/35
Teoría de perturbaciones para sistemas de n-dimensiones
Generalización en n-dimensiones
Generalizando las consideraciones anteriores para n grados de libertad,
en donde w0 ≡ {w01 , w02 , . . . , w0n } y J0 ≡ {J01 , J02 , . . . , J0n }, se tiene:
∂S ∂S1 ∂S2
J0α =
α = J α + α + 2 α + . . .
∂w0 ∂w0 ∂w0
∂S ∂S1 ∂S2
wα = α
= w0α + α + 2 α + . . .
∂J ∂J ∂J
Con lo cual, aplicando el mismo formalismo desarrollado anteriormente,
se obtiene para la expansión de E(J),
E0 (J) = K0 (J),
X∂S1 ∂K0
E1 (J) = K1 (w0 , J) + αν0α (J)
∀ ν0α (J) = ,
α ∂w0 ∂J α
X ∂S2 ∂S1 ∂K1 1 X ∂ν0α ∂S1 ∂S1
E2 (J) = ν0α α + α + ,
α ∂w0 ∂w0 ∂J α 2 αβ ∂J β ∂w0α ∂w0β
en donde J representa al set {J 1 , J 2 , . . . , J n }. 23
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 23/35
/35
Teoría de perturbaciones para sistemas de n-dimensiones
Generalización en n-dimensiones
Para calcular los promedios en la expansión de E(J) usamos,
Z 2π Z 2π Z 2π
1
hf (w0 , J)i = dw01 dw02 . . . dw0n f (w0 , J).
(2π)n 0 0 0
Como en el caso unidimensional, E1 (J) = hK1 i, con lo cual se tiene:
X ∂S1
ν0α = hK1 i − K1 .
α ∂w0α
Tal como se realizó anteriormente, el lado derecho de la ecuación an-
terior se puede expandir en series de Fourier:
X X
hK1 i − K1 = Θm (J)eim·w0 ∀ m · w0 = mα w0α ,
m α
en donde los coeficientes Θm vienen dados por,
Z 2π Z 2π Z 2π
1
Θm = dw01 dw02 . . . dw0n [hK1 i − K1 ] e−im·w0 .
(2π)n 0 0 0
24
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 24/35
/35
Teoría de perturbaciones para sistemas de n-dimensiones
Generalización en n-dimensiones
De igual manera expandiendo S1 ,
X
S1 = S1 (J; m)eim·w0 ,
m
por tanto,
X ∂S1 X
ν0α = i (ν0 · m)S1 (J; m)eim·w0 ,
α ∂w0α m
comparando ambas expansiones se obtiene:
Θm (J)
S1 (J; m) = .
i(ν0 · m)
Esta expresión provee una solución al problema a primer orden (S1 )
siempre y cuando el denominador sea diferente de cero.
Cuando ν0 ·m es cero se tiene el caso de resonancias, conocido también
como sistemas degenerados, en los cuales la teoría perturbativa falla
al describir el sistema. 25
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física
/35
25/35
Contenido: Tema 10
10. Teoría canónica de perturbaciones
10.1 Transformaciones canónicas y teoría de perturbaciones
10.2 Teoría de perturbaciones para sistemas unidimensionales
10.3 Teoría de perturbaciones para sistemas de n-dimensiones
10.4 Invariantes adiabáticos
26
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 26/35
/35
Invariantes adiabáticos
Condición adiabática
Cuando H depende explícitamente del tiempo, su solución es compli-
cada, debido a la dificultad de encontrar ctes. de movimiento.
Sin embargo, cuando la dependencia temporal del sistema es lo sufi-
cientemente lenta, entonces:
• Existirán variables dinámicas que varían poco, por tanto, es posible
considerarlas practicamente ctes.
• Tales cantidades resultan ser las variables de acción de sistemas
completamente integrables.
Consideremos un sist. unidimensional, completamente integrable, que
realiza un movimiento finito, caracterizado por un parámetro λ:
H = H(q, p, λ),
y en donde λ varía con el tiempo, pero lentamente, en un periodo T
del movimiento, tal que:
T dλ/dt = T λ̇ λ ⇒ λ̇ ≈ cte.
obteniendo así lo que se conoce como condición adiabática. 27
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 27/35
/35
Invariantes adiabáticos
Condición adiabática
Para analizar como afecta tal consideración al sistema, se calcula la
razón de cambio de E,
dH ∂H ∂H ∂H ∂H dE
= q̇ + ṗ + λ̇ = λ̇ = ,
dt ∂q ∂p ∂λ ∂λ dt
en donde los primeros dos términos se cancelan debido a las ecuaciones
canónicas de Hamilton.
Por tanto, para analizar la variación de E se toma el promedio en un
periodo de movimiento,
dE ∂H
= λ̇ ,
dt ∂λ
siendo λ̇ = cte. en el caso de la condición adiabática, y donde el
promedio viene dado como:
Z T
∂H 1 ∂H
= dt.
∂λ T 0 ∂λ
28
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 28/35
/35
Invariantes adiabáticos
Condición adiabática
Para calcular el promedio anterior, es posible utilizar una de las ecua-
ciones canónicas de Hamilton,
−1
∂H 1 ∂H
q̇ = ⇒ dt = dq = dq,
∂p ∂H/∂p ∂p
donde el periodo T se puede calcular en términos de q como:
Z T −1
∂H
I
T = dt = dq.
0 ∂p
Sustituyendo los resultados anteriores en la expresión de ∂H/∂λ,
Z T
dE ∂H λ̇ ∂H
= λ̇ = dt,
dt ∂λ T 0 ∂λ
−1
λ̇ ∂H ∂H
I
= H −1 dq.
∂H
dq ∂λ ∂p
∂p
29
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 29/35
/35
Invariantes adiabáticos
Condición adiabática
Debido a que las integrales de la ec. anterior se realizan en trayectorias
donde λ = cte. ⇒ E = cte., pudiendo definir p = p(q, E, λ), en donde
E y λ son independientes por tanto calculando:
d dE
H(q, p, λ) = ,
dλ dλ
∂H ∂p ∂H
⇒ + = 0 ∀ E 6= E(λ),
∂p ∂λ ∂λ
∂H ∂H −1 ∂p
∴ =− .
∂λ ∂p ∂λ
Sustituyendo el resultado anterior en la expresión de hdE/dti,
−1
dE λ̇ ∂H ∂H
I
= H −1 dq,
dt ∂H
dq ∂λ ∂p
∂p
λ̇ ∂p
I
= −H dq.
∂p
dq ∂λ
∂H
30
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 30/35
/35
Invariantes adiabáticos
Condición adiabática
Reordenando el resultado anterior,
dE ∂p ∂p
I I
dq + λ̇ dq = 0,
dt ∂H ∂λ
dE ∂p ∂p
I
∴ + λ̇ dq = 0.6
dt ∂E ∂λ
Ahora, por otro lado, analizando la acción del sistema,
I
J= pdq,
y calculando la diferencial total de J respecto al tiempo, considerando
p = p(q, E, λ),
dJ dp ∂p ∂E ∂p ∂p ∂p ∂E ∂p
I I I
= dq = + λ̇ + q̇ dq = + λ̇ dq.
dt dt ∂E ∂t ∂λ ∂q ∂E ∂t ∂λ
6
donde
H hdE/dti es
H independiente Hde q y haciendo H → E dentro la integral.
7
(∂p/∂q)q̇dq = (∂H/∂q)dq = dH = 0. 31
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 31/35
/35
Invariantes adiabáticos
Condición adiabática
Calculando el promedio temporal del resultado anterior:
dJ ∂p dE ∂p
I
= + λ̇ dq,
dt ∂E dt ∂λ
sin embargo, se había encontrado que para sistemas en donde aplica la
condición adiabática se cumple que:
dE ∂p ∂p
I
+ λ̇ dq = 0,
dt ∂E ∂λ
dJ
⇒ = 0,
dt
por tanto, se tiene que la acción es un invariante adiabático, es decir,
J permanece constante cuando el parámetro λ varía lentamente tal
que λ̇T λ ⇒ λ̇ ≈ cte.
32
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 32/35
/35
Invariantes adiabáticos
Ecuaciones de movimiento
λ = cte. ⇒ S = S(q, J) λ = λ(t) ⇒ S = S(q, J, λ(t))
Transformación canónica Transformación canónica
∂S ∂S ∂S ∂S
p= , w= , p= , w= ,
∂q ∂J ∂q ∂J
Ecuaciones de movimiento de Ecuaciones de movimiento de
Hamilton Hamilton
∂S ∂S
K=H+ ⇒ K = H(J), K=H+ ,
∂t ∂t
∂K
J˙ = − = 0 ⇒ J = cte., en donde:
∂w
∂K
Z
ẇ = = ν ⇒ w = νt + β. S= pdq.
∂J
33
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 33/35
/35
Invariantes adiabáticos
Ecuaciones de movimiento
Se observa que es necesario hallar ∂S/∂t para obtener las ecuaciones
de movimiento en el caso λ = λ(t):
∂K ∂K
J˙ = − , ẇ = .
∂w ∂J
Para ello se usa:
Z
S= pdq, ∀ p = p(q, E, λ(t)),
∂S ∂p ∂p ∂
Z Z Z
⇒ = dq = λ̇ dq = λ̇ pdq ,
∂t ∂t ∂λ ∂λ
∂S
= λ̇ = λ̇Λ(w, J),
∂λ
en donde se ha definido:
∂S
= Λ(w, J).
∂λ
34
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 34/35
/35
Invariantes adiabáticos
Ecuaciones de movimiento
Por tanto, de la expresión para el Hamiltoniano transformado,
∂S
K=H+ ,
∂t
= H(J) + λ̇Λ(w, J),
se hallan las ecuaciones de movimiento cuando λ = λ(t):
∂K ∂Λ
J˙ = − =− λ̇,
∂w ∂w
∂K ∂Λ
ẇ = = ν(J) + λ̇,
∂J ∂J
en donde ν(J) = ∂H/∂J es la frecuencia de oscilación calculada
considerando λ = cte.
35
Omar De la Peña-Seaman | IFUAP Mecánica Clásica − M.C. Física 35/35
/35
