Desvios das leis de Kepler são causados por exemplo pela influência gravitacional de Júpiter na órbita dos outros planetas do sistema solar, assim como por correções na energia potencial efetiva do sistema advindas da relatividade geral.
Considere uma partícula de massa m sob a ação da força central \vec{F}=F(r) \hat{r}, de modo que a energia potencial do sistema seja dada por
V(r)=-\frac{k}{r}+\frac{\alpha}{2 r^{2}},
representando um desvio do modelo kepleriano.
a) Determine a força agindo na partícula, em função de r.
b) Calcule o momento angular (com relação à origem) de uma órbita circular de raio a, para a>\alpha / k.
c) Calcule a frequência angular de rotação \Omega(a) da órbita circular de raio a.
d) Calcule a frequência \omega_{0} de pequenas oscilações radiais em torno da órbita circular de raio a.
e) Se \Delta \varphi é o ângulo [em coordenadas polares (r, \varphi) ] varrido pela partícula durante uma oscilação radial completa, o ângulo de precessão \Delta \psi devido à passagem da partícula por dois pericentros consecutivos é \Delta \psi=\Delta \varphi-2 \pi. Calcule \Delta \psi para pequenas oscilações ao redor de uma órbita circular de raio a, em primeira ordem em \frac{\alpha}{k a}. Explicite o passo-a-passo de sua dedução.
Passo 1
O enunciado dá a energia potencial do sistema e pede várias quantidades relacionadas a órbitas circulares e pequenas perturbações. Vou listar claramente os dados e as perguntas (peguei todas estas informações do enunciado):
Considera-se uma órbita circular de raio a, com a condição a>\dfrac{\alpha}{k} (isso garante certas quantidades positivas, que usaremos).
O que o enunciado pede:
(a) determinar a força \vec{F}(r) em função de r;
(b) calcular o momento angular (relativo à origem) de uma órbita circular de raio a;
(c) calcular a frequência angular de rotação \Omega(a) dessa órbita circular;
(d) calcular a frequência \omega_{0} de pequenas oscilações radiais em torno da órbita circular de raio a;
(e) calcular o ângulo de precessão \Delta\psi=\Delta\varphi-2\pi para pequenas oscilações ao redor da órbita circular de raio a, em primeira ordem em \dfrac{\alpha}{k a}.
Passo 2
(a) Determinar a força a partir do potencial.
Lembrete: para uma força central derivada de um potencial escalar V(r), a componente radial da força é F_r=-\dfrac{dV}{dr} e o vetor força é \vec F = F_r \,\hat r. Vou calcular a derivada explicitamente.
\vec{F}(r)=\left(-\frac{k}{r^{2}}+\frac{\alpha}{r^{3}}\right)\hat r.
Passo 3
(b) Momento angular da órbita circular de raio a.
Método usado (básico e prático para provas): usar o potencial efetivo e a condição de equilíbrio (mínimo) para órbita circular.
Definição: potencial efetivo para movimento radial com momento angular L (módulo do momento angular relativo à origem) é
V_{\text{eff}}(r)=V(r)+\frac{L^{2}}{2 m r^{2}}.
Uma órbita circular de raio a corresponde a um ponto de equilíbrio radial, logo
(e) Ângulo de precessão \Delta\psi para pequenas oscilações, em primeira ordem em \dfrac{\alpha}{k a}.
Definição e relação útil: seja T_{r} o período da oscilação radial (tempo entre dois pericentros consecutivos). Durante esse tempo a partícula gira um ângulo
\Delta\varphi=\Omega\,T_{r}.
Para uma oscilação harmônica radial pequena, T_{r}=2\pi/\omega_{0}. Assim
Voltando a \varepsilon=\dfrac{\alpha}{k a}, obtemos, em primeira ordem em \dfrac{\alpha}{k a},
\Delta\psi \approx -\pi\,\frac{\alpha}{k a}.
Isto significa que, para \alpha>0, o pericentro sofre uma precessão retrograda (ângulo negativo) de magnitude \pi\alpha/(k a) por oscilação radial — a direção (sinal) depende do sinal de \alpha.