Lista 13 - Equações de Hamilton-Jacobi

1
a) Construa a formulação equivalente a Equação de Hamilton-Jacobi, mas baseado em uma função geratriz do tipo 3.

b)Aplique a na solução do oscilador harmônico unidimensional de massa m e frequência angular $\omega$. Resolva completamente o problema.

Resolução:

a) A função do tipo 3, $F_{3} = F_{3}(p,Q,t)$ será reescrita como uma função do tipo $S=S(p,Q,t)$, e temos as relações que são conhecidas desta função geradora:

$q = -\dfrac{\partial S}{\partial p}\,\, e\,\, P = \dfrac{\partial S}{\partial Q}$

Para a transformação canônica  podemos escrever uma Kamiltoniana da seguinte forma:

$\displaystyle K = H(q,p,t) + \dfrac{\partial S(p,Q,t)}{\partial t}$ e a Kamiltoniana será nula por hipótese se encontrarmos a transformação canônica adequada. Assim:

$ H(q,p,t) + \dfrac{\partial S(p,Q,t)}{\partial t} = 0$                               (1)

Para separar as variáveis vamos supor uma função $S$ do tipo:

$S(p,Q,t) = W(p,Q)-V(Q)t$

E tomando a derivada parcial de S e igualando a H vemos que: $H = V = \alpha$, pois as duas funções são de variáveis independentes. E reescrevendo S temos:

$S = W(p,Q) - \alpha t$, esta equação fornece as soluções desejadas para cada problema.

b) A Hamiltoniana do oscilador unidimensional de massa $m$ e frequência angular $\omega$ é dado por:

$H = \dfrac{p^{2}}{2m} + \dfrac{1}{2}m\omega^{2}q^{2}$

Para aplicarmos a solução encontrada precisaremos reescrever $H= H(\dfrac{\partial S}{\partial p})$ o que fornece a seguinte equação:

$H = \dfrac{p^{2}}{2m}+ \dfrac{m\omega^{2}}{2}(\dfrac{\partial S}{\partial p})^{2}$

do item a) vimos que $S = W(p,Q) - V(Q)t é uma forma de separar as variáveis do problema, fazendo esta substituição e tomando as derivadas das funções podemos escrever:

$H = \dfrac{p^{2}}{2m} + \dfrac{m\omega^{2}}{2}(\dfrac{\partial W(p,Q)}{\partial p})^{2}$

Da expressão para S, e da equação (1) vemos que $H = \dfrac{\partial S}{\partial t} = \alpha$, onde $\alpha$ é uma constante. E podemos, então escrever:

$ \dfrac{\partial W(p,Q)}{\partial p} =  \dfrac{1}{m\omega}\sqrt{2m\alpha - p^{2}}$

Que pode ser integrado em $p$ fornecendo:

$W = \dfrac{\alpha}{\omega}(\frac{1}{2}\sin(2\theta) + \theta)$

Retornando para as variáveis originais temos:

$W = \dfrac{\alpha}{\omega}(\frac{p}{\sqrt{2m\alpha}}\sqrt{1-\frac{p^{2}}{2m\alpha}}+ arcsin(\frac{p}{\sqrt{2m\alpha}})$

E podemos escrever a função de Hamilton como: 

$S = \dfrac{\alpha}{\omega}(\frac{p}{\sqrt{2m\alpha}}\sqrt{1-\frac{p^{2}}{2m\alpha}}+ arcsin(\frac{p}{\sqrt{2m\alpha}})-\alpha t$

E como $q = -\dfrac{\partial S}{\partial p}$, podemos escrever $\alpha \omega = \beta$ o que fornece uma equação que permite determinar $p$ como:

$p(t) = \sqrt{2\alpha}sin(\omega t - \beta)$

e

$q(t) = -\sqrt{\frac{2\alpha}{m\omega^{2}}}\cos(\omega t - \beta)$

2


a) O que são variáveis de ação-ângulo? Qual sua relação com a teoria de Hamilton-Jacobi?

b) Faça uso destas variáveis e resolva o problema do oscilador harmônico simples.

Resolução

a) Como vimos a função principal de Hamilton-Jacobi, pode ser relacionada a integral de ação do sistema $S = \int Ldt$, e pelo princípio de Hamilton esta grandeza é estacionária para a evolução do sistema, ou seja pode ser considerada cíclica, o que possibilita usar essa grandeza como variável canônica, e neste caso esta grandeza é dita ser uma variável de ação. e o seu momento conjugado é chamado de variável ângulo.


b) Como vimos para o oscilador harmônico temos que:

$H = \dfrac{p^{2}}{2m} + \dfrac{1}{2}m\omega^{2}q^{2} = \alpha$

Que é constante do movimento, já que $H$ independe explicitamente do tempo. É possível reescrever esta equação como :

$(\frac{p}{\sqrt{2m\alpha}})^{2} + (\sqrt{\dfrac{m}{2\alpha}}\omega q)^{2} = 1$

Precisamos determinar a integral de ação, para um caso de um momento J tal que:

$J = \frac{1}{2\pi}\oint pdq$

Fazendo $a = 2m\alpha\, \, e \, \, b = \frac{m\omega^{2}}{2\alpha}$ podemos determinar de modo direto o valor da área desta integral que fornece:

$J = \frac{\alpha}{\omega}$, 

Ou seja:

$H = \omega J$. 

fornece a frequência deste movimento. e as equações horárias são as mesmas equações paramétricas da elipse: 

$p(t) = acos(X)\, \, e \, \, q(t) = bsin(Y)$ 

O que resulta na mesma solução encontrada na questão 1)

Lista 12 - Transformações canônicas e Parêntesis de Poisson.

1 - Sabe-se que o conjunto de transformações entre um conjunto de coordenadas é:

$Q = \ln (1 + \sqrt{q}\cos(p))$

$P = 2(1 + \sqrt{q}\cos(p))\sqrt{q}\sin(p)$

a) Mostre diretamente que Q e P são variáveis canônicas se p e q também o forem.

b) Mostre que a função geratriz $F_{3}$ que dá origem a esta transformação canônica é:

$F_{3} = -(e^{Q}-1)^2 \tan(p)$

Resolução:

As condições que permitem demonstrar diretamente que as variáveis são canônicas são: 

$\displaystyle (\dfrac{\partial Q_{i}}{\partial q_{j}})_{q,p} = (\dfrac{\partial p_{j}}{\partial P_{i}})_{Q,P}$

$\displaystyle (\dfrac{\partial P_{i}}{\partial q_{j}})_{q,p} = -(\dfrac{\partial p_{j}}{\partial Q_{i}})_{Q,P}$


$\displaystyle (\dfrac{\partial Q_{i}}{\partial p_{j}})_{q,p} = -(\dfrac{\partial q_{j}}{\partial P_{i}})_{Q,P}$

$\displaystyle (\dfrac{\partial P_{i}}{\partial p_{j}})_{q,p} = (\dfrac{\partial q_{j}}{\partial Q_{i}})_{Q,P}$

Para aplicar estas condições precisamos agora escrever q e p como função de Q e P.

As equações de partida são: 

$Q = \ln (1 + \sqrt{q}\cos(p))$                         (1)

$P = 2(1 + \sqrt{q}\cos(p))\sin(p)$                 (2)

Tomando a exponencial de ambos os lados da equação (1), e substituindo expressão equivalente a  na equação (2), eliminando também $\sqrt{q}$, pode-se escrever uma equação pra $p$ em função de Q e P dada por:

$$p = \arctan(\dfrac{Pe^{-Q}}{2(e^{Q}-1)})$$

Isolando $\sqrt{q}\cos(p)\, e\, \sqrt{q}\sin(p)$ nas expressões 1 e 2 e somando o quadrado destas expressões obtemos o ultima expressão que precisamos:

$$q = (e^{Q}-1)^{2} + \dfrac{P^{2}}{4e^{2Q}}$$

Agora, tendo em mãos as 4 equações para obtenção das relações diretas procedemos os cálculos das derivadas desejas: como só temos um grau de liberdade não precisamos nos preocuparmos com os índices:

$\displaystyle (\dfrac{\partial Q}{\partial q})_{q,p} = \dfrac{\cos(p)}{2(\sqrt{q}+q\cos(p))} = (\dfrac{\partial p_{j}}{\partial P_{i}})_{Q,P}$

$\displaystyle (\dfrac{\partial P}{\partial q})_{q,p} = \sin(p)(2\cos(p) + \dfrac{1}{\sqrt{q}}) = -(\dfrac{\partial p}{\partial Q})_{Q,P}$




$\displaystyle (\dfrac{\partial Q}{\partial p})_{q,p} = -\dfrac{\sqrt{q}\sin(p)}{1+\sqrt{q}\cos(p)} = -(\dfrac{\partial q}{\partial P})_{Q,P}$



$\displaystyle (\dfrac{\partial P}{\partial p})_{q,p} = -2q(\sin^{2}(p)-\cos^{2}(p)) -2\sqrt{q}\cos(p) = (\dfrac{\partial q}{\partial Q})_{Q,P}$

E observamos que as derivadas satisfazem as condições diretas. Outra maneira de desenvolver este exercício é partir dos colchetes de Poisson, e mostrar que são válidos em cada caso como segue:

$[Q,Q]=[P,P] = 0$

$[Q,P]=1$

$[P,Q]=-1$

Que são os critérios exigidos para que a transformação seja canônica.

b) Para determinar se a função F, dada por:  $F_{3} = -(e^{Q}-1)^2 \tan(p)$ é função geratriz desta transformação, e sabendo que F é do tipo 3, precisamos mostrar que as expressões:

$q = - \dfrac{\partial F_{3}}{\partial p}$

$P = \dfrac{\partial F_{3}}{\partial p}$

Fazendo a derivada temos:

$\dfrac{\partial F_{3}}{\partial p} = -(e^{Q}-1)\sec^{2}(p) \Rightarrow \sqrt{q}\cos(p) = e^{Q}-1$

$\dfrac{\partial F_{3}}{\partial p} = -\tan(p)2(e^{Q}-1)e^{Q} \Rightarrow P = 2q\sin(p)\cos(p) + 2\sqrt{q}\sin{p}$

Nesta ultima expressão é preciso reescrever $e^{Q}$ em função de q para podermos encontrar a solução desejada. Como partindo da função $F_{3}$ encontramos as grandezas esperadas, pode-se dizer que esta é uma função geratriz desta transformação.

2) Mostre que se a Hamiltoniana e uma quantidade $F$ são constantes de movimento, então, $\dfrac{\partial F}{\partial t}$ também deverá ser constante de movimento. Como ilustração desse resultado, considere o movimento uniforme de uma partícula livre de massa m. A Hamiltoniana, certamente se conserva e a quantidade $F =  x- \dfrac{pt}{m}$ é uma constante do movimento. Mostre então, que a constante de movimento  $\dfrac{\partial F}{\partial t}$ concorda com $[H, F]$.

Resolução: 

Podemos economizar na notação definindo que $\dfrac{\partial U}{\partial t} = \partial ^{t} U$. Usando esta notação e sabendo que H e F são constantes do movimento, podemos escrever que:

$\dfrac{d H}{d t} = [H,H] +\partial^{t} H = 0$                               (1)

$\dfrac{d F}{d t} = [F,H] +\partial^{t} F = 0$                               (2)

Então temos que: 

$\dfrac{d }{d t} \partial^{t} F = [\partial^{t} F,H] + \partial^{t}(\partial^{t} F) = 0$

Usando o resultado da equação (2) temos que:

$\dfrac{d }{d t} \partial^{t} F = [-[F,H],H] - \partial^{t}([F,H]) = 0$

Que usando a regra da cadeia e os resultados em 1 e 2, pode ser escrito como:

$\dfrac{d }{d t} \partial^{t} F = [[H,F],H] + [[F,H],H] +[[H,H],F]$ 

Que é a identidade de Jacobi, portanto:

$\dfrac{d }{d t} \partial^{t} F = 0$, ou seja,  $\dfrac{\partial F}{\partial t}$ concorda com $[H, F]$.

b) A hamiltoniana do sistema é dada por: $H = \dfrac{p^{2}}{2m}$, portanto $\dot{p} = 0$ e p é constante. Usando o fato de F ser constante do movimento.

$\partial^{t} F = -[F,H] = [H,F]$, então:

$\partial^{t} F = -\dfrac{p}{m}$

Ou seja,  $\dfrac{\partial F}{\partial t}$ é constante do movimento.

3) Considere a Lagrangeana $\displaystyle L = e^{\beta t} |\dfrac{m\dot{x}^{2}}{2}-V(x)|$

a) Escreva a equação do movimento para o sistema descrito por ela.

b) Obtenha o momento canonicamente conjugado à coordenada generalizada.

c) Escreva a Hamiltoniana do sistema e obtenha a sua variação temporal $\displaystyle \dfrac{d H}{d t}$.

Resolução:


a) Temos duas situações a considerar.

$\displaystyle \dfrac{m\dot{x}^{2}}{2}>V(x)$        (1)

e 

$\displaystyle \dfrac{m\dot{x}^{2}}{2}<V(x)$     (2)

Vamos primeiramente analisar o caso (1). 

$L = e^{\beta t}(\dfrac{m\dot{x}^{2}}{2}-V(x))$ 

Substituindo na equação de Euler-Lagrange tem-se que :

$\ddot{x} + \beta \dot{x} + m^{-1}\dfrac{\partial V(x)}{\partial x} = 0$

Que é uma das equações do movimento.

No caso (2): 

$L = -e^{\beta t}(\dfrac{m\dot{x}^{2}}{2}-V(x))$ 

E substituindo na equação de Euler-Lagrange tem-se que :

$\ddot{x} + \beta \dot{x} + m^{-1} \dfrac{\partial V(x)}{\partial x}$

b)  

No caso 1, temos que :

$p = \dfrac{\partial L}{\partial \dot{x}} = e^{\beta t}m\dot{x}$

No caso 2:

$p = \dfrac{\partial L}{\partial \dot{x}} = - e^{\beta t}m\dot{x}$

c) $H =p\dot{x} -L$ , vamos analisar cada um dos casos.

No caso 1 teremos:

$H = e^{\beta t}m\dot{x}^{2} - e^{\beta t}(\dfrac{m\dot{x}^{2}}{2}-V(x))$

Caso 1:

$H =   e^{- \beta t}\dfrac{p^{2}}{2m}+ e^{\beta t}V(x)$

Caso 2:

$H = - e^{- \beta t}\dfrac{p^{2}}{2m}- e^{\beta t}V(x)$

Lembrando que $\dfrac{d H}{dt} = - \dfrac{\partial L}{\partial t}$,  e escrevendo a lagrangeana em função do momento generalizado temos:

Caso 1:

$\dot{H} =   e^{- \beta t}\dfrac{p^{2}}{2m}\beta+ e^{\beta t}V(x)\beta$

Caso 2:

$\dot{H} = - e^{- \beta t}\dfrac{p^{2}}{2m}\beta- e^{\beta t}V(x)\beta$

Lista de Exercícios 11 - Colchetes de Poisson.

Resolução:

2)  

Resolução:

Lista 10 - Transformação Canônica.

1)
a) Ache a transformação canônica definida pela função geradora:  $$F(q,Q,t) =\dfrac{1}{2}m\omega q^{2}cot(Q)$$.

Escreva as equações do movimento nas variáveis Q e P para um oscilador harmônico de frequência angular $\omega$.

b) Ache a transformação canônica definida pela função geradora: $$F(q,Q,t) = \dfrac{1}{2}m\omega (q - \dfrac{f(t)}{m\omega^{2}})^{2}cot(Q)$$.

Escreva as equações do movimento nas variáveis Q e P para o mesmo oscilador do item anterior quando este é submetido a uma força externa $f(t)$.


Resolução:

a) $F = F(q,Q,t)$ é uma função geratriz do tipo 1. E tem-se que:

$p = \dfrac{\partial F}{\partial q}$, 

$P = -\dfrac{\partial F}{\partial Q}$, e

$K(Q,P,t) = H(q,p,t) + \dfrac{\partial F}{\partial t}$ 

logo:

$p = m\omega qcot(Q)$,  

$P = -(-\dfrac{1}{2}m\omega q^{2}cosec^{2}(Q) ) = \dfrac{m\omega q^{2}}{2sin^{2}(Q)}$.

O que permite escrever: 

$q = \sqrt(\dfrac{2P}{m\omega})sin(Q)$; 

e

 $p = \sqrt(2Pm\omega)cos(Q)$. 

Que são as duas equações de transformação canônica. Para um oscilador harmônico simples de frequência $\omega$ e massa $m$ a expressão do Hamiltoniano é dado por:

$H = \dfrac{1}{2m}(p^{2} + m^{2}\omega^{2} q^{2})$.

Onde $w = \sqrt(\dfrac{k}{m})$.  Realizando a transformação canônica, e lembrando que $\dfrac{\partial F}{\partial t} = 0$, ou seja,  $K=H$, então :

$$K = P\omega$$.

Usando as expressões para a Kamiltoniana temos: 

$\dot{Q} = \dfrac{\partial K}{\partial P}$ 

e

 $\dot{P} = -\dfrac{\partial K}{\partial Q}$. 

Então: 

$$\dot{Q} = \omega$$

$$\dot{P} = 0$$

Que são as equações do movimento.

b) Esta função geradora também é do tipo 1. Repetindo os mesmo procedimentos utilizados temos:

$p = \dfrac{\partial F}{\partial q} = m\omega (q-\dfrac{f(t)}{m\omega^{2}})cot(Q)$        (1 b)

$P = -\dfrac{\partial F}{\partial Q}  = \dfrac{m\omega(q-\dfrac{f(t)}{m\omega^{2}})^{2}}{2\sin^{2}     (Q)}$     (2 b)

$\dfrac{\partial F}{\partial t} = -\dfrac{\dot{f}(t)}{\omega}(q-\dfrac{f(t)}{m\omega^{2}})cot(Q)$.    (3b)

De (2 b) podemos escrever:

$q = \sqrt{\dfrac{2P}{m\omega}}sin(Q) + \dfrac{f(t)}{m\omega^{2}}$  (4 b)

e substituindo este resultado em (1 b) temos:

$p = \sqrt{2Pm\omega}cos(Q)$    (5 b)

4 e 5 b são as duas equações de transformação.

A Hamiltoniana do oscilador sob ação de uma força externa pode ser escrita como:

$H = \dfrac{p^{2}}{2m} + \dfrac{ m\omega^{2}}{2} (q-\dfrac{f(t)}{m\omega^{2}})^{2}$.

Onde o fator $\dfrac{f(t)}{m\omega^{2}}$  representa a distensão do oscilador de constante de mola $ k = m\omega^{2}$,  quando sujeito a força adicional  $ f(t) $.   A Kamiltoniana é obtida da expressão $K(q,Q,t) = H(q,p,t) + \dfrac{\partial F}{\partial t}$, portanto:

$K = \omega P -\dfrac{\dot{f}(t)}{\omega}(q-\dfrac{f(t)}{m\omega^{2}})cot(Q)$

ou, 

$K = \omega P -\dfrac{\dot{f}(t)}{\omega}\sqrt{\dfrac{2P}{m\omega}}cos(Q)$

portanto: 

$\dot{Q} = \dfrac{\partial K}{\partial P} = \omega - \dfrac{\dot{f}(t)}{2\omega}\sqrt{\dfrac{2}{m\omega P}}cos(Q)$

e

$\dot{P} = \dfrac{\partial K}{\partial Q} = \dfrac{\dot{f}(t)}{\omega}\sqrt{\dfrac{2P}{m\omega }}sin(Q)$

Que são as duas equações do movimento nas coordenadas P e Q.

2)
a) Resolva usando a formulação Hamiltoniana o problema de uma partícula de massa $m$ que sofre uma queda livre a partir do repouso


Consideremos uma partícula de massa $m$, em queda livre sob ação da força gravitacional de aceleração $\vec{a} = g$, energia cinética desta partícula é: $T = \dfrac{1}{2}m\dot{y}^{2}$ e a energia potencial é : $V = -mgy$, cuja lagrangeana é $L = \dfrac{1}{2}m\dot{y}^{2} + mgy$,  portanto a Hamiltoniana deste sistema é:
$H = \dfrac{1}{2}m\dot{y}^{2} - mgy$, que reescrito e considerando que $ p = \dfrac{\partial L}{\partial q}$, ou seja $ p = m\dot{y}$. Temos:

$H = \dfrac{p^{2}}{2m} - mgq$

E pode-se obter escrever:

$\dot{p} = -\dfrac{\partial H}{\partial q} = mg$  (2a)

e

$\dot{q} = \dfrac{\partial H}{\partial p} = mp$ (2b)

Substituindo 2a em 2b, temos: $\ddot{q} = g$, ou seja,

$q(t) = g\dfrac{t^{2}}{2} + at + b$, usando as condições de contorno: $q(0) = 0$, e $\dot{q}(0) = 0$, temos:

$q(t) = \dfrac{1}{2}gt^{2}$

b) É possível reescrever este problemas em termos de novas coordenadas e momentos de modo que única coordenada seja cíclica?

Resposta:

Sim, é possível. Observando que a energia mecânica se conserva, neste caso, pode-se simplesmente definir uma grandeza para ser a grandeza física idêntica a energia mecânica e então encontrar a variável cíclica associada, Como o momento linear total se conserva, vamos por exemplo supor que $P = \dfrac{p^{2}}{2m} - mgq$. Então a Kamiltoniana é $K = P$, e neste caso $Q$ é a variável cíclica associada. Mas, a forma de realizar esta transformação não é única, e outros casos são possíveis de ser encontrados.


c) Em caso afirmativo, determine F, a função geratriz de transformação canônica correspondente. Se necessário utilize a resposta do item a) para descobrir a transformação. usando a transformação $P = \dfrac{p^{2}}{2m} - mgq$, pode-se escrever : $q = \dfrac{1}{mg}(\dfrac{p^{2}}{2m}-P)$. Vamos escolher uma função do tipo $F_{4}$, assim: $q_{i} = -\dfrac{\partial F_{4}}{\partial p_{i}}$, utilizando esta expressão e a função conhecida para $q_{i}$, podemos então integrar em $p_{i}$ e obter $F_{4} = F(p,P,t)$.

$F_{4} = -\int \dfrac{1}{mg}(\dfrac{p^{2}}{2m}-P)dp$

$F_{4} = \dfrac{1}{mg}(Pp - \dfrac{p^{3}}{6m})$

Esta é a função desejada.

Lista 9 - Teorema de Noether

1) Demonstre o teorema de Noether.

Solução:


Vamos supor um sistema descrito pela Lagrangeana $L(q(t),\dot{q}(t),t)$ e que a ação associada é : $S = \int_{t_{1}}^{t_{2}}L(q(t),\dot{q}(t),t)dt$.

Supondo a seguinte transformação:

$t \mapsto t' = t + \epsilon X(q(t),t)$ e $q(t) \mapsto q(t') = q(t) + \epsilon \Psi_{i} (q(t),t)$. Onde $\epsilon$  é um parâmetro infinitesimal.


Agora vamos supor que a ação é invariante sob a transformação considerada. Agora, calculando a variação da ação temos:

$$\Delta S = \int_{t'_{1}}^{t'_{2}}L(q(t'),\frac{dq(t')}{dt'},t')dt' - \int_{t_{1}}^{t_{2}}L(q(t),\frac{dq(t)}{dt},t)dt$$


Fazendo uma mudança na integral da Lagrangeana em $dt'$, para a variável de integração $dt$, levando em conta as transformações temos que :

$\dfrac{dt'}{dt} = 1+\epsilon \dot{X}$, de onde $\frac{dt}{dt'} = (1+\epsilon \dot{X})^{-1}$, e expandindo este último termo usando a expansão binomial, considerando apenas os termos de primeira ordem tem-se:

$\dfrac{dt}{dt'} = 1-\epsilon \dot{X}$, portanto:

$\dfrac{dq'(t')}{dt'} = \dfrac{dq'(t')}{dt} \dfrac{dt}{dt'} $ tomando a derivada da transformação de $q'(t')$ e substituindo a expressão encontrada para $\dfrac{dt}{dt'}$ temos que: $\dfrac{dq'(t')}{dt'} = \dot{q}_{i}+\epsilon \xi_{i}$.

Sendo que $\xi = \dot{\Psi}_{i}-\dot{q}_{i}\dot{X}$, tendo descartados os termos de ordem maior que 1 em $\epsilon$.

Reescrevendo a variação na integral de ação em termos da variável $dt$ obtem-se:

$$\Delta S = \int_{t_{1}}^{t_{2}}L(q(t)+\epsilon \Psi, \dot{q}+\epsilon \xi,t + \epsilon X)(1+\epsilon \dot{X})dt - \int_{t_{1}}^{t_{2}}L(q(t),\frac{dq(t)}{dt},t)dt$$

Para calcular a variação na ação vamos primeiramente expandir a  primeira Lagrangeana numa série de Taylor de primeira ordem e então calcular a variação na ação.


$L(q(t)+\epsilon \Psi, \dot{q}+\epsilon \xi,t + \epsilon X) = L(q,\dot{q},t) + \sum_{i=1}^{n} (\dfrac{\partial L}{\partial q_{i}}\epsilon \Psi_{i} + \dfrac{\partial L}{\partial \dot{q}_{i}}\epsilon \xi_{i}) + \dfrac{\partial L}{\partial t}\epsilon X$

Então a variação na ação pode ser escrita como: $\Delta S = \epsilon \int_{t_{1}}^{t_{2}} [\sum_{i=1}^{n} (\dfrac{\partial L}{\partial q_{i}}\Psi_{i} + \dfrac{\partial L}{\partial \dot{q}_{i}}\xi_{i}) + \dfrac{\partial L}{\partial t}X + L\dot{X}]dt$

Como $\epsilon$ e os limites de integração são arbitrários, para que a ação seja nula, é necessário que o integrando seja nulo. Portanto : $\sum_{i=1}^{n} [\dfrac{\partial L}{\partial q_{i}}\Psi_{i} + \dfrac{\partial L}{\partial \dot{q}_{i}}(\dot{\Psi}_{i}-\dot{q}_{i}\dot{X})] + \dfrac{\partial L}{\partial t}X + L\dot{X} = 0$


Agora iremos utilizar a equação de Lagrange, $\dfrac{d}{dt}(\dfrac{\partial L}{\partial \dot{q}_{i}}) - \dfrac{\partial L}{\partial q_{i}} = 0$ para escrever a derivada parcial em $\dot{q}_{i}$. Então:

$\sum_{i=1}^{n} [\dfrac{d}{dt}(\dfrac{\partial L}{\partial \dot{q}_{i}})\Psi_{i} + \dfrac{\partial L}{\partial \dot{q}_{i}}\dot{\Psi}_{i}]-(\sum_{i=1}^{n}\dot{q}_{i}\dfrac{\partial L}{\partial \dot{q}_{i}}- L)\dot{X} + \dfrac{\partial L}{\partial t}X= 0$.

Como sabemos que : $\dfrac{dh}{dt}=-\dfrac{\partial L}{\partial t}$, escrevendo o primeiro termo como a derivada do produto, e lembrando que $ h = \sum_{i=1}^{n}\dot{q}_{i}\dfrac{\partial L}{\partial \dot{q}_{i}}- L$ então:

$$\sum_{i=1}^{n} [\dfrac{d}{dt}(\dfrac{\partial L}{\partial \dot{q}_{i}}\Psi_{i}) -hX - \dfrac{dh}{dt}X$$

O termos negativos podem ser reescritos como: $-\dfrac{d}{dt}[hX]$, e vemos que a constante do movimento pode ser escrita como:  $C = hX - \dfrac{d}{dt}(\dfrac{\partial L}{\partial \dot{q}_{i}})\Psi_{i}$, ou reescrita como:

$$\sum_{i=1}^{n}\dfrac{\partial L}{\partial \dot{q}_{i}}(\dot{q}_{i}X-\Psi_{i}) - LX = C$$

Onde C é uma constante do Movimento.


2) Apresente uma situação Física para cada um dos casos abaixo:

a) A função de Hamilton é conservada, mas não corresponde à energia mecânica.
b) A função de Hamilton corresponde a energia mecânica, mas não é conservada.

Resolução:
a) Se encontrarmos um termo que depende da velocidade na função de Hamilton, isto já é suficiente para podermos dizer que esta expressão não corresponde somente a energia mecânica do corpo. Um exemplo deste fato é o caso de uma partícula carregada num campo eletromagnético. Camo já vimos a Lagrangeana deste sistema é dado por: $L =\dfrac{m\dot{x}^{2}}{2}-q(\phi + \dfrac{\dot{\vec{x}}}{c}.\vec{A})$, onde q é a carga da partícula, c é a velocidade da luz,  $\vec{A}$ é o potencial vetor, e $\phi$ o potencial escalar. Estes potenciais não dependem do tempo, e a Lagrangeana não depende diretamente do tempo, logo $\dfrac{\partial L}{\partial t} = 0$ e imediatamente, $\dfrac{dH}{dt} = 0$, ou seja, H é conservado. Porém o potencial depende da velocidade, e portanto não corresponde somente a energia mecânica.

b) Neste caso basta que a Hamiltoniana dependa diretamente do tempo, e que o potêncial não seja função da velocidade, para que a Hamiltoniana corresponda a energia mecânica e que não seja conservada. Por exemplo um pêndulo de massa m tem seu comprimento l encurtado a uma taxa constante $\alpha$, ou seja, $l(t) = l_{0}-\alpha t$. O Lagrangeano do sistema é $L = \dfrac{1}{2}m\alpha^{2} + \dfrac{1}{2}ml(t)^{2}\dot{\theta}^{2} + mgl(t)cos\theta$. Logo como $\dfrac{\partial L}{\partial t}\neq 0 \mapsto \dfrac{dH}{dt} \neq 0$. A energia não é conservada, mas corresponde a energia mecânica.

3) Escreva a função de Hamilton para um sistema massa-mola cuja constante elástica vale k e a massa vale m. Escolha dois conjuntos distintos de condições iniciais e represente o movimento do sistema no espaço de fase em cada caso.

Resolução:

O sistema pode ser descrito pela lagrangeana $L = \dfrac{1}{2}m\dot{x}^{2} - \dfrac{1}{2}kx^{2}$.
Sistema com 1 grau de liberdade. O hamiltoniano do sistema é descrito por : $H = p\dot{q}-L$, onde as coordenadas generalizadas são $q = x$, e $p = \dfrac{\partial L}{\partial \dot{q}} = m\dot{x}$ que implica que $\dot{x} = \dfrac{p}{m}$, então:

$ H =  \dfrac{p^{2}}{2m} + \dfrac{1}{2}kx^{2}$.

Usando as equação de Hamilton temos que :

$\dot{p} = -\dfrac{\partial H}{\partial x} = -kx$ e $\dot{x} = \dfrac{\partial H}{\partial p} = \dfrac{p}{m}$

Se prestarmos atenção na expressão do Hamiltoniano, não precisamos resolver este sistema, já que a energia total E é conservada, ou seja: $p = \sqrt{2mE}$,  em $x = 0$, e $x = \sqrt{\dfrac{2E}{k}}$, para $p=0$. Então:

$$\dfrac{x^2}{(2E/k)} + \dfrac{p^2}{(2mE)}=1.$$

Para duas condições iniciais $(x_{0},p_{0})$ diferentes teremos duas energias diferentes. De modo geral o espaço de fase é dado por curvas do tipo:

Onde E e E', representam condições iniciais diferentes para o sistema dentro dos limites de oscilações periódicas.

Lista de Exercícios 8 - Hamiltoniana.


1) Suponha que a função de Hamilton para um sistema de N partículas não varie quando ocorre uma translação infinitesimal de todo o sistema. Discuta a Conservação de Momentum.

Resolução:

Dada a Hamiltoniana de um sistema temos que a equação de Hamilton que descreve a variação dessa Hamiltoniana com relação a uma variação translacional infinitesimal é: $\frac{\partial H}{\partial q_{j}} = -\dot{p_{j}}$, como não ocorreu variação em H, isto implica que o momentum conjugado a esta coordenada ciclica $q_{j}$ é conservado.

2) Obtenha a função de Hamilton para um oscilador anarmônico, cuja função de Lagrange é dada por:

$$L = \frac{\dot{x^{2}}}{2}-\frac{w^{2}x^{2}}{2}-\alpha x^{3} +\beta x\dot{x}^{2}$$

Resolução:

Dada a Lagrangeana, pode-se usando a transformação de Legendre escrever:

$H = p\dot{q}-L$, onde p é o momento conjugado a coordenada x, e é dado por: $p = \frac{\partial L}{\partial \dot{q}}$, sendo que para este sistema $q=x$. Assim temos que:

$p_{x} = \dot{x} + 2\beta x\dot{x}$, ou $\dot{x} = \frac{p_{x}}{1 + 2\beta x}$, e a Hamiltoniana fica : 

$$H= p_{x}(\frac{p_{x}}{1+2\beta x})-\frac{\dot{x}^{2}}{2}+ \frac{w^{2}x^{2}}{2}+ \alpha x^{3}-\beta x\dot{x}^{2}$$

que escrito como $H = H(x,p_{x})$ fica:

$$ H = \frac{1}{2} (\frac{p_{x}^{2}}{1+2\beta x}) + \frac{w^{2}x^{2}}{2}+ \alpha x^{3}$$ 

3) Ache a equação de movimento de uma partícula, cuja função de Hamilton é:

a) $H(x,p_{x}) = \frac{p_{x}^2}{2}+\frac{w^{2}x^{2}}{2}+ \lambda (\frac{p_{x}^{2}}{2}+\frac{w^{2}x^{2}}{2})^{2}$

b) $H(x,p_{x}) =A\sqrt{p}-xF$

Resolução: 

Usando as equações canônicas $ \dot{x} = \frac{\partial H}{\partial p_{x}}$ e $\dot{p_{x}}= -\frac{\partial H}{\partial x} $  temos:

a) $ \dot{x} = p_{x} + \lambda (p_{x}^{2}+w^{2}x^{2})p_{x}$, ou $$ \dot{x} = p_{x}(1 + \lambda (p_{x}^{2}+w^{2}x^{2}))$$ 

E,  $\dot{p}_{x} = -(w^{2}x + \lambda (p_{x}^{2}+w^{2}x^{2})w^{2}x)$, ou  $$\dot{p} = -w^{2}x(1+ \lambda(p_{x}^{2}+w^{2}x^{2}) )$$

 b) Se A e F forem constantes arbitrárias teremos  $\dot{x} = \frac{A}{2\sqrt{p_{x}}}$ e $\dot{p_{x}} = F$,

Ou seja, $p_{x} = F.t$, e $x = A\sqrt(\frac{t}{F})+C$, onde C é uma constante.