Hamilton Jacobi
Hamilton Jacobi
Hamilton Jacobi
ECUACIÓN DE HAMILTON-JACOBI
que no es otra cosa mas que el momento angular. Es razonable pensar que el potencial
es un potencial central,
V = V (q12 + q 22 ) ,
enteramente consistente con la cuarta ecuación en derivadas parciales.
--------------------------------
32
∂D
= [H , D ]
∂t
lo que, utilizando la forma de D, lleva a la ecuación
∂H ∂f ∂H ∂g ∂H ∂f ∂H ∂g
g= + t− − t
∂q ∂p ∂q ∂p ∂p ∂q ∂p ∂q
Igualando los coeficientes de t, se llega a las ecuaciones
∂H ∂g ∂H ∂g
− =0
∂q ∂p ∂p ∂q
∂H ∂f ∂H ∂f
− =g
∂q ∂p ∂p ∂q
La primer ecuación muestra que g = H es solución. Utilizando este resultado en la
segunda, ésta queda
1 ∂f ∂f p 2 1
−p = − 2
q ∂p
3
∂q 2 2 q
con lo que
∂f q ∂f p
=− ; =−
∂p 2 ∂q 2
es decir:
pq
f =−
2
--------------------------------------
33
Si P = p + q, la condición anterior queda como:
∂Q ∂Q
− = 1.
∂q ∂p
La solución de la ecuación anterior es suma de la correspondiente a la ecuación
homogénea y de una solución particular de la ecuación completa. La integral de la
ecuación homogénea es de la forma:
Q hom = f ( p + q )
siendo f una función arbitraria. Por otra parte, Q = q es una solución particular de la
ecuación inhomogénea, luego la solución general es:
Q = q + f ( p + q)
Para el caso general en que P es una función cualquiera de p y q, la solución a la
ecuación homogénea, [P,Q] = 0, es de la forma Q hom f ( P ) , como puede verificarse por
simple sustitución.
-----------------------------------
4 2
34
1 2 1
p + mω 2 q 2 = E
2m 2
con un área encerrada dada por 2πE/ω. Por lo tanto el invariante adiabático es I = E/ω.
La condición de invariancia adiabática viene dada por:
E1 E2
=
ω1 ω2
Si escribimos la expresión anterior en términos de amplitudes de oscilación, la amplitud
final puede expresarse como:
1 4
k
A2 = A1 1
k2
-----------------------------------
35
∂K P &
Q& = = −D
∂P m
∂
P& = − V (Q + D)
∂Q
Combinándolas:
&& = − ∂ &&
mQ V (Q + D ) − m D
∂Q
-------------------------------------------
(
F1 = − 2m 2 g (Q − q ) )3 2
3m 2 g .
p2
Deseamos pasar de H = + mgq a K = mgQ . De la identificación de ambas
2m
expresiones, obtenemos:
∂F1
p = (2m) 1 2 (mgQ − mgq )
1 2
= ,
∂q
en donde F1 es la función generatriz problema. Integrando la ecuación anterior,
36
siendo ∆ ≡ Qq Pp − Q p Pq . Es pues necesario y suficiente que el determinante ∆ sea la
unidad.
--------------------------------
b) Pi = q i ; Qi = p i
c) Pi = − a ij q j ; Qi = bij p j
d) Pi = − Aij p j ; Qi = Bij q j
Q = q α cos β p
P = q α sin β p
una transformación canónica? ¿Cuál es la forma de la función generatriz F3 ?
37
1
α= ; β =2
2
La transformación es:
Q = q cos 2 p
P = q sin 2 p
∂F3 ∂F3
además, q=− ; p=− y F3 = F3 ( p, Q)
∂p ∂q
Rescribiendo las variables P y q en función de p y Q:
Q2
q= ; P = Qtg 2 p
cos 2 p
Q2 ∂F ∂F
de modo que = − 3 y Qtg 2 p = − 3
cos 2 p ∂p ∂Q
Las dos ultimas ecuaciones se integran inmediatamente a
Q2
F3 = − tg 2 p + ϕ ( p)
2
donde ϕ ( p ) =constante de modo que,salvo una constante irrelevante se verifica
finalmente que
Q2
F3 (Q, p ) = − tg2 p .
2
-------------------------------------------
1
L= mx& 2 (1 + 4a 2 x 2 ) − mgax 2
2
Solo un grado de libertad, x . Para obtener el Hamiltoniano ,
∂L
p= = mx& (1 + 4a 2 x 2 ) , y:
∂x
1 p2
H= + mgax 2 .
2m 1 + 4a 2 x 2
38
∂F2
La ecuación de Hamilton-Jacobi se obtiene de substituir p i por en la expresión
∂q i
∂F2
H ( pi , qi , t ) + = 0 , de modo que:
∂t
2
1 ∂F2 ∂F
+ mgax 2 + 2 = 0 (Ecuación de Hamilton-Jacobi).
2m(1 + 4a x ) ∂x
2 2
∂t
Puesto que H no depende explícitamente de t, F2 = W − αt donde α = const , y
2
1 ∂W
+ mgax = α ,
2
2m(1 + 4a x ) ∂x
2 2
39
donde G (Q) es en principio una función arbitraria.
Ahora bien, para este función generadora resulta que al Hamiltoniano transformado K
es:
K=H (5)
K=
1
2
1
[
+ P 2 Q 4 = x 2 + p 2 = H = 2 + − PQ 2 + G (Q) ]
2
Q Q
↑
en virtud de (4)
de este modo,
[
P 2 Q 4 = − PQ 2 + G (Q) ]2
⇒ G =0.
Finalmente,
p = − PQ 2
1
x=
Q
es la transformación.
-------------------------------------
12. Sean las variables q y p y sean Q y P las nuevas variables que resultan de
p df (q )
una transformación canónica. Sabiendo que Q = , donde f ′(q) =
f ′(q) dq
,y f es una función dada que sólo depende de q ,
a. Obtener la función generadora de la transformación canónica.
b. Determinar completamente la nueva variable P en función de q y p .
c. ¿Qué puede decirse del Hamiltoniano que gobierna la evolución de q y p ?
d. Si no se conociese la función generadora, para el caso particular
Q = pq ; P = − ln(q) , ¿cómo podría averiguarse directamente si la
transformación es canónica ó no?
40
∂F1
p = Qf ′(q ) = ,
∂q
e integrando,
F1 = Qf (q ) + ϕ (Q)
siendo ϕ una función arbitraria.
b).- De (1,b) ⇒ P = − f (q) − ϕ ′(Q) y en función de p y q:
p
Q=
f ′(q)
p
P = − f (q ) − ϕ ′( ) (2)
f ′(q )
c).- La condición para que la transformación sea canónica es independiente del
Hamiltoniano, Nada puede pues decirse de H .
d).-
Q = pq
(3)
P = − ln q
Está claro que este es un caso particular de (2) con ϕ ′ = 0; f = ln q . Por tanto la
transformación es canónica. Si no supiésemos la forma de F1 , podríamos hacer una
comprobación directa de que se cumple la condición general
∂Q ∂Q
∂q ∂p 0 1
MJM T = J , con M= , y J = .
∂P ∂P −1 0
∂q ∂p
De hecho
p q − q p 1
p − 0 1
M= 1 ; MJM = 1 q =
T
− q 0 0 − q q 0 − 1 0
Se verifica pues que MJM T = J , de modo que la transformación (3) es canónica, por
verificación directa. Nótese, por cierto, que para problemas con solo grado de libertad
como éste, la condición
MJM T = J ⇔ Det M = 1 , que es aún más trivial de probar.
-------------------------------------------
41
La ecuación, una vez hecha la separación de variables, queda como
dW ( z )
2 2 2
1 dW1( x ) dW ( y )
2 3
+ + + mgz = E
2m dx dy dz
2 2
1 dW1( x ) 1 dW2 ( x )
=α , =α
2m dx 1 2m dy 2
2
1 dW3 ( z )
+ mgz = α
2m dz 3
Y las soluciones
W3 ( z ) = ±
8
2
(α3 − mgz )32
9mg
∂S m ∂S m
β = = ±x − t, β = = ±x −t
1 ∂α 2α 2 ∂α 2α
1 1 2 2
β =
∂S
=± 3
(
2 α − mgz
−t
)
3 ∂α 2
3 mg
2α 2α
x=± 1 β +t,
m 1
( ) y=±
m
(
2 β +t
2
)
α g
z= 3 − β +t 2
mg 2 3
( )
-------------------------------------------
42
H =
1 3 2
[( )] [( )]
∑ p + exp − φ1 − φ3 + exp − φ2 − φ1 + exp − φ3 − φ2 − 3
2i = 1 i
[( )]
Aparte del hamiltoniano, existe una integral del movimiento obvia ¿Cuál es?
Genere una transformación canónica que `ponga en evidencia esta nueva
integral en el hamiltoniano transformado. (Nota: en la segunda parte del
problema utilice la función generatriz F2 , definiendo como nuevo momento esa
integral de movimiento adicional)
p1
φ
1
φ2
p3
p2 φ 3
H′=
2
[
1 2
P1 + P22 + (P3 − P1 − P2 )2 + ]
exp(- Φ 1 ) + exp[- (Φ 2 − Φ )] + exp(Φ 2 ) − 3.
43
Se tiene que
∂Q q ∂Q
= ; = tan p,
∂p cos 2 p ∂q
Luego
P = lnsen p + f (q tan p )
44
εq
F2 = f (q, P ) − cos ω t
ω
45
p 2 − 2γ t 1
H = e + mω 2 q 2e 2γ t .
2m 2
Nos piden ahora una función generatriz de la forma F2 (q, P, t ) = f (t ) qP para que
transforme H en un nuevo hamiltoniano K que sea constante del movimiento. Nos
bastará para ello asegurar que K no depende explícitamente del tiempo. Para ello,
aplicando las ecuaciones de transformación
∂f
p = f (t )P ; Q = f (t )q ; K = H + qP
∂t
tenemos
P 2 2 − 2γ t 1 ∂f
K= f e + mω 2Q 2 f − 2e 2γ t + QP f −1
2m 2 ∂t
Si hacemos la identificación
f (t ) = eγ t
obtenemos el resultado exigido.
-------------------------------------------
∂F2 ∂F
a) Buscaremos una función F2 (q, P, t ) tal que Qi = y p i = 2 , de :
∂Pi ∂q i
Q1 = q12 y Q2 = q1 + q 2 → F2 = q12 P1 + (q1 + q 2 ) P2 + g (q1 , q 2 , t ) .
Por tanto:
∂g ∂g
p1 = 2q1 P1 + P2 + ; p 2 = P2 +
∂q1 ∂q 2
con lo que
∂g 1 ∂g 1 ∂g ∂g
P2 = p 2 − ; P1 = p1 − P2 − = p1 − p 2 + −
∂q 2 2q1 ∂q1 2q1 ∂q 2 ∂q1
b) Si suponemos que g es independiente del tiempo, K = H y si sustituimos las
expresiones anteriores de P1 y P2 en K tenemos que
2
p − p2 1 ∂g ∂g ∂g
H = K = 1 + − + p 2 −
2q1 2q1 ∂q 2 ∂q1 ∂q 2
que será igual a la hamiltoniana original si
∂g ∂g ∂g
= −(q1 + q 2 ) = −q12 − q 22 − 2q1 q 2
2
= y
∂q1 ∂q 2 ∂q 2
46
y esta última expresión es fácilmente integrable a
1 ∂g ∂g
g = −q12 q 2 − q 23 − q1 q 22 + f (q1 ) y de la condición = = −q12 − q 22 − 2q1 q 2
3 ∂q1 ∂q 2
∂f
tenemos que = −q12 y por tanto
∂q1
1 1
F2 = q12 P1 + (q1 + q 2 )P2 − q12 q 2 − q 23 − q1 q 22 − q13 .
3 3
c) En el nuevo Hamiltoniano no aparecen ni Q1 ni Q2 por lo que P1 y P2 son
constantes y
∂H
Q& 1 = = 2P1 → Q1 = 2 P1t + α ,
∂P1
∂H
Q& 2 = = 1 → Q2 = t + β .
∂P2
-------------------------------------------
p1 − aq1
F1 = q1 q 2 y F2 =
q2
sean constantes del movimiento.
b) ¿Existen más constantes del movimiento algebraicas independientes? En
caso de respuesta afirmativa, ¿ podría dar alguna sugerencia de constante?
c) ¿Puede construirse alguna constante más a partir de la identidad de
Jacobi?. (Nota: identidad de Jacobi para las funciones H, v, w es
[H , [v, w]] + [v, [w, H ]] + [v, [H , w]] = 0 )
∂W ∂S
d) Resolver la ecuación de Hamilton-Jacobi: H q, + = 0 , con
∂q ∂t
S (q, P, t ) = W (q, P ) − α1t ¿Qué constantes del movimiento obtenemos?
a) Para que Fi sea constante del movimiento tiene que cumplir que [H , Fi ] = 0 , de la
definición de corchete de Poisson tenemos:
[H , F1 ] = ∂H ∂F1 − ∂H ∂F1 + ∂H ∂F1 − ∂H ∂F1 = − g 2 (q1 )q 2 + g 3 (q 2 )q1 = 0
∂q1 ∂p1 ∂p1 ∂q1 ∂q 2 ∂p 2 ∂p 2 ∂q 2
que tiene como solución g 2 (q1 ) = q1 y g 3 (q 2 ) = q 2 (o con los signos cambiados).
Introduciendo esto en el corchete de la segunda constante:
47
1 aq − p 1 ∂g1
[H , F2 ] = − ∂g1 + p1 ∂g 2 + g2
a
+ g 3 1 2 1 = − + p1 + aq1 + aq1 − p1 = 0
∂q1 ∂q1 q 2 q2 q 2 q 2 ∂q1
que se cumple si g1 (q1 ) = aq 1 .
2
48
-------------------------------------------
49
c) Como la transformación es independiente del tiempo la nueva hamiltoniana se
obtiene de la substitución de las nuevas variables, q = − P y p = Q − aP 2 en la antigua
hamiltoniana:
H ' = H = P 4 (a 2 − 4a + 4) + P 2 Q(4 − 2a) + P 2 + Q 2
en la que tomando a = 2 nos lleva a
H ' = P 2 + Q 2 , formalmente la del oscilador armónico.
~ ~ ~ ~
Para esta hamiltoniana la transformación P = 2 P cos(Q ) Q = 2 P sen(Q ) , nos lleva
~
a H ' = 2 P . Sencillamente podemos igualar las expresiones de P y Q en ambas
coordenadas para obtener:
~ ~ ~ ~
P = 2 P cos(Q ) = −q y Q = 2 P sen(Q ) = 2q 2 + p
y de aquí despejar para cualquier conjunto:
~ −q
~ ~ Q = tn −1 2
q = − 2 P cos Q, 2q + p
o
~ ~ ~ 2 ~
p = 2 P senQ − 4 P cos Q ~ q 2 + (2q 2 + p) 2
P=
2
La composición de 2 transformaciones canónicas siempre es canónica.
-------------------------------------------
a) No siempre será posible encontrar funciones de los 4 tipos para una misma
transformación canónica. Partiendo de las ecuaciones de las nuevas variables, Q(q, p, t )
y P(q, p, t ) , hay que comprobar qué parejas formadas por una variable antigua y una
nueva son expresables en función de las otras dos. Como ejemplo simple de este caso
tenemos la transformación: F1 (q, Q, t ) = qQ . Esta genera el cambio
∂F ∂F
p = 1 = Q , P = 1 = − q , para el que no es posible expresar p como función de q y
∂q ∂Q
P, ni q como función de p y Q. Por tanto existirán sólo las funciones F1 (q, Q, t ) y
F4 ( p, P, t ) .
Sin embargo sí hay casos en los que es posible, un ejemplo simple sería:
p = 2Q 2 q
, generada por la función F1 = q 2 Q 2 . Para este caso tenemos todas las
P = −2Qq 2
combinaciones:
50
1
p 2
p P2 3
p= q= q =
2q 3 2Q 2 2p
, , y podremos calcular todas las funciones
p 2p2 1
Q=− 2 P=− p2 3
2q 4Q 3 Q =
2P
generatrices:
2
P 2
p2 p5 3
F2 = − 2 , F3 = − , F4 = 2 .
4q 4Q 2 2P
Estas funciones generatrices, obtenidas unas de otras a través de F1 = F2 − QP ,
F1 = F3 − qp y F1 = F4 + qp − QP y del oportuno cambio de coordenadas (ver sección
9-1 del Goldstein) siempre generan la misma hamiltoniana.
1 2
24.- Sean un sistema hamiltoniano H ( p, q) = ( p + q 2 ) , la variable dinámica
2
g ( p, q ) = q p y las condiciones iniciales q (t = 0) = 0, p (t = 0) = p 0 . Calcule la
evolución temporal de la variable, g (t , p 0 ) siguiendo los siguientes procedimientos
y comente brevemente las diferencias:
a) Integración directa de las ecuaciones de Hamilton.
b) A partir del formalismo de Poisson obtenga la ecuación diferencial de
segundo orden de g (t , p 0 ) .
c) Encuentre una transformación canónica cualquiera tal que la nueva variable
sea Q = g ( p, q ) . Calcule el nuevo hamiltoniano y obtenga g (t , p 0 ) a partir de
las nuevas ecuaciones de Hamilton.
NOTA: las siguientes fórmulas pueden serle de utilidad:
∂F ∂F ∂F ∂F
F1 = F1 (q, Q, t ), p = 1 , P = − 1 F2 = F2 (q, P, t ), p = 2 , Q = 2
∂q ∂Q ∂q ∂P
∂F ∂F ∂F ∂F
F3 = F3 ( p, Q, t ), q = − 3 , P = − 3 F4 = F4 ( p, P, t ), q = − 4 , Q = 4
∂p ∂Q ∂p ∂P
51
∂H
q& = =p
∂p q (t ) = A sin t + B cos t
a) , q&& = p& = −q y por tanto , que para las
∂H p (t ) = A cos t − B sin t
p& = − = −q
∂q
q (t = 0) = B = 0
condiciones iniciales dadas: , tenemos q(t ) = q 0 sin t .
p (t = 0) = A = p0
p 02
g (t ) = p (t )q (t ) = p02 sin t cos t = sin 2t
2
dg ∂g ∂H ∂g ∂H
b) g& = = [g , H ] = − = p2 + q2 y
dt ∂q ∂p ∂p ∂q
dg& ∂g& ∂H ∂g& ∂H
g&& = = [g& , H ] = − = 2qp − 2 p (−q ) = 4 g
dt ∂q ∂p ∂p ∂q
g (t ) = A sin(2t ) + B cos(2t ) , con las condiciones iniciales g (t = 0) = 0 , g& (t ) = p02 , nos
p02
lleva a g (t ) = sin 2t .
2
c) Buscamos una transformación canónica desde el sistema de coordenadas p, q a otro
Q
P, Q donde, por ejemplo, Q = g = pq . Si despejamos q = , podemos considerar la
p
función generatriz F3 ( p, Q) tal que
∂F Q
q=− 3 =
∂p p
, → F3 = Q ln p + f (Q) → P = ln p + f ' (Q)
∂F3
P=−
∂Q
siguiendo el enunciado del problema podemos tomar el caso más simple f (Q) = 0 :
p = e P , q = Qe − P
y el nuevo hamiltoniano será
1
H = (Q 2e −2 P + e 2 P )
2
Las nuevas ecuaciones del movimiento serán:
∂H
Q& = = −Q 2 e −2 P + e 2 P
∂ P
y para la ecuación de g (t ) = Q(t ) será
& ∂H −2 P
P=− = Qe
∂Q
Q&& = −2Q 2 e −2 P P& + 2Qe −2 P Q& + 2e 2 P P& , dónde por substitución de las anteriores
ecuaciones obtenemos Q && = 2(−Q 2 e −2 P + e 2 P )Qe −2 P + 2e −2 P Q(−Q 2 e −2 P + e 2 P ) = 4Q , es
decir la misma ecuación que en el apartado anterior.
--------------------------------------------
52
b.- Utilícese esto para encontrar una constante de movimiento para una partícula
r r
r a⋅r r
en dos dimensiones bajo el potencial V (r ) = 3 (siendo a un vector constante
r
dado).
df
a) Para demostrar que = 0 utilizaremos el formalismo de Poisson:
dt
s
∂f ∂H ∂f ∂H ∂f ∂H ∂f ∂H ∂f ∂H ∂f ∂f ∂H ∂f
f& = [ f , H ] = ∑ − = − = − =0
i =1 ∂q i ∂p i ∂pi ∂qi ∂q1 ∂p1 ∂p1 ∂q1 ∂q1 ∂f ∂p1 ∂p1 ∂f ∂q1
b) El hamiltoniano de una particula sometida al potencial del problema la podemos
escribir en coordenadas polares como:
r r
pr2 pθ2 a⋅r r r
H = + 2 + 3 , si elegimos el eje x coincidiendo con el vector a podemos
2m 2r m r
p 2
p2 a ⋅ r cosθ p r2 pθ2 + 2ma ⋅ cosθ
escribir H = r + 2θ + = + de modo que podemos
2m 2 r m r3 2m 2r 2 m
p 2 f (θ , pθ )
escribir que H (r , p r , f (θ , pθ )) = r + 2
, con f (θ , pθ ) = pθ2 + 2ma ⋅ cosθ ,
2m 2r m
constante del movimiento.
53