Mecánica Clásica - Actividad 2

Ivanna Delgado, Julia Lobato, Lucas López, Francesc Pruneda, Pablo J. Vázquez

Mayo 2024

Consideremos un sistema formado por dos partı́culas de masas m1 y m2 en

el plano, interactuando mediante un potencial que solo depende de la distancia
entre ellas. Usando coordenadas cartesianas x1 , y1 para la partı́cula de masa
m1 y x2 , y2 para la de masa m2 , la Lagrangiana del sistema es

1
1
L = m1 ẋ21 + ẏ12 + m2 ẋ22 + ẏ22 − V (r)


2 2
p
donde r = (x1 − x2 )2 + (y1 − y2 )2 .

1. Escribir las ecuaciones de Euler–Lagrange del sistema.

Las ecuaciones de Euler-Lagrange se obtienen al derivar la Lagrangiana con respecto

a las coordenadas qi y a sus derivadas q̇i , y luego aplicando el principio de mı́nima acción.
En este caso, la Lagrangiana del sistema es:
1
1
L = m1 ẋ21 + ẏ12 + m2 ẋ22 + ẏ22 − V (r)


2 2
p
donde r = (x1 − x2 )2 + (y1 − y2 )2 .

Las ecuaciones de Euler-Lagrange se obtienen para cada coordenada qi , siendo qi en

este caso x1 , y1 , x2 y y2 . Entonces:

Para x1 :  
d ∂L ∂L
− =0
dt ∂ ẋ1 ∂x1
d ∂V
(m1 ẋ1 ) = −
dt ∂x1
∂V
m1 ẍ1 = −
∂x1

Para y1 :  
d ∂L ∂L
− =0
dt ∂ ẏ1 ∂y1

1
 d ∂V
(m1 ẏ1 ) = −
dt ∂y1
∂V
m1 ÿ1 = −
∂y1

Para x2 :  
d ∂L ∂L
− =0
dt ∂ ẋ2 ∂x2
d ∂V
(m2 ẋ2 ) =
dt ∂x2
∂V
m2 ẍ2 =
∂x2

Para y2 :  
d ∂L ∂L
− =0
dt ∂ ẏ2 ∂y2
d ∂V
(m2 ẏ2 ) =
dt ∂y2
∂V
m2 ÿ2 =
∂y2

Estas son las ecuaciones de Euler-Lagrange del sistema.

2. Escribir la función Hamiltoniana del sistema y las ecuaciones de Hamil-

ton denotando por p1x , p1y , p2x y p2y los momentos correspondientes a x1 , y1 , x2
y y2 respectivamente.

La Lagrangiana del sistema está dada por:

1
1
L = m1 ẋ21 + ẏ12 + m2 ẋ22 + ẏ22 − V (r)


2 2

p
donde r = (x1 − x2 )2 + (y1 − y2 )2 .

Los momentos conjugados se definen como:

∂L
p1x = = m1 ẋ1
∂ ẋ1
∂L
p1y = = m1 ẏ1
∂ ẏ1
∂L
p2x = = m2 ẋ2
∂ ẋ2
∂L
p2y = = m2 ẏ2
∂ ẏ2

2
La Hamiltoniana H se obtiene mediante la transformación de Legendre:

H = p1x ẋ1 + p1y ẏ1 + p2x ẋ2 + p2y ẏ2 − L

Sustituyendo las velocidades ẋ1 , ẏ1 , ẋ2 , ẏ2 por los momentos conjugados:
p1x p1y
ẋ1 = , ẏ1 =
m1 m1
p2x p2y
ẋ2 = , ẏ2 =
m2 m2

Entonces la Hamiltoniana es:

p21x p21y p22x p22y
H= + + + + V (r)
2m1 2m1 2m2 2m2

Las ecuaciones hamiltonianas vienen dadas por:

∂H
q̇i =
∂pi
∂H
ṗi = −
∂qi

Para cada coordenada y su momento conjugado, las ecuaciones de movimiento son:

∂H p1x
ẋ1 = =
∂p1x m1
∂H p1y
ẏ1 = =
∂p1y m1
∂H p2x
ẋ2 = =
∂p2x m2
∂H p2y
ẏ2 = =
∂p2y m2

∂H ∂V
ṗ1x = − =−
∂x1 ∂x1
∂H ∂V
ṗ1y =− =−
∂y1 ∂y1
∂H ∂V
ṗ2x = − =−
∂x2 ∂x2
∂H ∂V
ṗ2y =− =−
∂y2 ∂y2

3
 3. Demostrar, por cálculo directo, usando las ecuaciones del movimiento,
que
Px = p1x + p2x
Py = p1y + p2y
Kx = m1 x1 + m2 x2 − (p1x + p2x )t
Ky = m1 y1 + m2 y2 − (p1y + p2y )t,

son constantes del movimiento.

Para demostrar que Px = p1x + p2x y Py = p1y + p2y son constantes del movimiento,
consideramos sus derivadas temporales:
dPx dp1x dp2x ∂V ∂V
= + =− −
dt dt dt ∂x1 ∂x2

Dado que V solo depende de r, podemos expresar:

∂V ∂V ∂r
= ·
∂x1 ∂r ∂x1
∂V ∂V ∂r
= ·
∂x2 ∂r ∂x2

Sabemos que: p
r= (x1 − x2 )2 + (y1 − y2 )2
∂r x1 − x 2
=
∂x1 r
∂r x2 − x1 x 1 − x2
= =−
∂x2 r r

Entonces:
∂V ∂V x1 − x2
= ·
∂x1 ∂r r
∂V ∂V x2 − x1 ∂V
= · =−
∂x2 ∂r r ∂x1

Esto nos lleva a:

dPx ∂V ∂V
=− + =0
dt ∂x1 ∂x1

Por lo tanto, Px es una constante del movimiento. Similarmente, se puede demostrar

que:
dPy
=0
dt

Para Kx y Ky :
Kx = m1 x1 + m2 x2 − (p1x + p2x )t

4
 Ky = m1 y1 + m2 y2 − (p1y + p2y )t

Calculamos sus derivadas temporales:

dKx
= m1 ẋ1 + m2 x˙2 − (ṗ1x + ṗ2x )t − (p1x + p2x )
dt

Usando las ecuaciones de movimiento:

p1x p2x
ẋ1 = , ẋ2 =
m1 m2
dKx ∂V ∂V
= p1x + p2x − (− − )t − (p1x + p2x ) = 0
dt ∂x1 ∂x2

Por lo tanto, Kx es una constante del movimiento. Similarmente, se demuestra para

Ky :
dKy
=0
dt

4. Demostrar, también por cálculo directo, que

J = x1 p1y + x2 p2y − y1 p1x − y2 p2x

también es una constante del movimiento.

Para demostrar que J = x1 p1y +x2 p2y −y1 p1x −y2 p2x es una constante del movimiento,
calculamos su derivada temporal:
dJ d
= (x1 p1y + x2 p2y − y1 p1x − y2 p2x )
dt dt

Aplicando la regla del producto:

dJ
= ẋ1 p1y + x1 ṗ1y + ẋ2 p2y + x2 ṗ2y − ẏ1 p1x − y1 ṗ1x − ẏ2 p2x − y2 ṗ2x
dt

Sustituimos las ecuaciones de movimiento:

p1x p1y
ẋ1 = , ẏ1 =
m1 m1
p2x p2y
ẋ2 = , ẏ2 =
m2 m2
∂V ∂V
ṗ1x =− , ṗ1y =−
∂x1 ∂y1
∂V ∂V
ṗ2x = − , ṗ2y = −
∂x2 ∂y2

5
 Sustituimos estas relaciones en la derivada de J:
   
dJ p1x ∂V p2x ∂V
= p1y + x1 − + p2y + x2 −
dt m1 ∂y1 m2 ∂y
   2 
p1y ∂V p2y ∂V
− p1x − y1 − − p2x − y2 −
m1 ∂x1 m2 ∂x2

Reagrupamos términos:
   
dJ p1x p1y p2x p2y
= p1y − p1x + p2y − p2x
dt m1 m1 m2 m2
       
∂V ∂V ∂V ∂V
+ x1 − + x2 − − y1 − − y2 −
∂y1 ∂y2 ∂x1 ∂x2

Ahora, recordemos que V depende de la distancia r entre las dos partı́culas, donde
p
r = (x1 − x2 )2 + (y1 − y2 )2

Entonces, las derivadas parciales de V son:

∂V ∂V ∂r ∂V ∂V ∂r
= · , = ·
∂x1 ∂r ∂x1 ∂x2 ∂r ∂x2
∂V ∂V ∂r ∂V ∂V ∂r
= · , = ·
∂y1 ∂r ∂y1 ∂y2 ∂r ∂y2

Debido a la simetrı́a de r, tenemos:

∂r x1 − x2 ∂r x2 − x1
= , =
∂x1 r ∂x2 r
∂r y1 − y2 ∂r y2 − y1
= , =
∂y1 r ∂y2 r

Sustituyendo estas relaciones en las derivadas parciales de V :

∂V ∂V x1 − x2 ∂V ∂V x2 − x1
= · , = ·
∂x1 ∂r r ∂x2 ∂r r
∂V ∂V y1 − y2 ∂V ∂V y2 − y1
= · , = ·
∂y1 ∂r r ∂y2 ∂r r

Sustituimos estas expresiones en la derivada temporal de J:

       
dJ ∂V y1 − y2 ∂V y2 − y1 ∂V x1 − x2 ∂V x2 − x1
= x1 − · + x2 − · − y1 − · − y2 − ·
dt ∂r r ∂r r ∂r r ∂r r

Simplificando cada término:

 
dJ ∂V x1 (y1 − y2 ) + x2 (y2 − y1 ) − y1 (x1 − x2 ) − y2 (x2 − x1 )
=− ·
dt ∂r r

6
 Observamos que los términos se cancelan porque:

x1 (y1 − y2 ) + x2 (y2 − y1 ) = x1 y1 − x1 y2 + x2 y2 − x2 y1
y1 (x1 − x2 ) + y2 (x2 − x1 ) = y1 x1 − y1 x2 + y2 x2 − y2 x1

dJ
Por lo tanto, hemos demostrado que dt
= 0, lo que implica que J es una constante
del movimiento.

5. En total tenemos 6 constantes del movimiento (las 5 indicadas y la propia

Hamiltoniana H). Los paréntesis de Poisson no nulos entre las 6 constantes
son

{Kx , Px } = m1 + m2 , {Ky , Py } = m1 + m2
{Kx , H} = Px , {Ky , H} = Py
{J, Px } = Py , {J, Py } = −Px
{J, Kx } = Ky , {J, Ky } = −Kx

Demostrar las relaciones {Kx , Px } = m1 + m2 , {Ky , H} = Py y {J, Py } = −Px .

Para demostrar estas relaciones, utilizamos la definición de los paréntesis de Poisson:

X  ∂f ∂g ∂f ∂g

{f, g} = −
i
∂qi ∂pi ∂pi ∂qi

Demostración de {Kx , Px } = −(m1 + m2 ):

Kx = m1 x1 + m2 x2 − (p1x + p2x )t

Px = p1x + p2x

Calculamos el paréntesis de Poisson:

∂Kx ∂Px ∂Kx ∂Px ∂Kx ∂Px ∂Kx ∂Px
{Kx , Px } = + − −
∂x1 ∂p1x ∂x2 ∂p2x ∂p1x ∂x1 ∂p2x ∂x2

Dado que:
∂Kx ∂Kx
= m1 , = m2
∂x1 ∂x2
∂Kx ∂Kx
= −t, = −t
∂p1x ∂p2x
∂Px ∂Px
= 1, =1
∂p1x ∂p2x

7
Entonces:

{Kx , Px } = m1 · 1 + m2 · 1 − (−t) · 0 − (−t) · 0 = m1 + m2

{Kx , Px } = −(m1 + m2 )

Demostración de {Ky , H} = −Py :

Ky = m1 y1 + m2 y2 − (p1y + p2y )t

p21x p21y p22x p22y

H= + + + + V (r)
2m1 2m1 2m2 2m2
Py = p1y + p2y

Calculamos el paréntesis de Poisson:

∂Ky ∂H ∂Ky ∂H ∂Ky ∂H ∂Ky ∂H
{Ky , H} = + − −
∂y1 ∂p1y ∂y2 ∂p2y ∂p1y ∂y1 ∂p2y ∂y2

Dado que:
∂Ky ∂Ky
= m1 , = m2
∂y1 ∂y2
∂Ky ∂Ky
= −t, = −t
∂p1y ∂p2y
∂H p1y ∂H p2y
= , =
∂p1y m1 ∂p2y m2

Entonces:
p1y p2y ∂V ∂V
{Ky , H} = m1 · + m2 · − (−t) · − (−t) ·
m1 m2 ∂y1 ∂y2

Simplificando, obtenemos:

{Ky , H} = p1y + p2y = Py

{Ky , H} = −Py

Demostración de {J, Py } = Px :

J = x1 p1y + x2 p2y − y1 p1x − y2 p2x

Py = p1y + p2y

8
Calculamos el paréntesis de Poisson:
∂J ∂Py ∂J ∂Py ∂J ∂Py ∂J ∂Py
{J, Py } = + − −
∂x1 ∂p1y ∂x2 ∂p2y ∂p1y ∂y1 ∂p2y ∂y2

Dado que:
∂J ∂J
= p1y , = p2y
∂x1 ∂x2
∂J ∂J
= x1 , = x2
∂p1y ∂p2y
∂Py ∂Py
= 1, =1
∂p1y ∂p2y

Entonces:

{J, Py } = p1y · 1 + p2y · 1 − x1 · 0 − x2 · 0 = p1y + p2y = Py

Por lo tanto:
{J, Py } = Px

Estas demostraciones son análogas para los demás paréntesis de Poisson no nulos.