Expansion of the planetary Hamiltonian for small eccentricities and inclinations: a vector formalism

Resumo Técnico: Expansão do Hamiltoniano Planetário para Pequenas Excentricidades e Inclinações: Uma Formalização Vetorial

Enunciado do Problema
A evolução de longo prazo das órbitas planetárias é governada pela função perturbadora, que encapsula as interações gravitacionais entre os planetas. Desenvolvimentos clássicos na dinâmica planetária tipicamente expandem esta função como uma série de potências em elementos orbitais (excentricidades e inclinações). Embora expansões em razões de semieixo maior sejam frequentemente tratadas usando formalismos vetoriais (por exemplo, polinômios de Legendre), expansões em excentricidades e inclinações tradicionalmente dependem de quantidades escalares. Embora trabalhos anteriores de Abdullah (posteriormente estendido por Boué e Fabrycky) tenham reconstruído com sucesso a parte secular da função perturbadora utilizando variáveis do tipo spin derivadas do momento angular e vetores de excentricidade, uma derivação geral para o Hamiltoniano completo — incluindo harmônicos de Fourier não seculares (desigualdades) — utilizando uma abordagem puramente baseada em vetores desde o início, era inexistente.

Metodologia
O autor apresenta uma derivação da expansão do Hamiltoniano planetário que emprega uma formalização vetorial durante todo o processo, evitando a dependência inicial de elementos orbitais escalares. A metodologia procede da seguinte forma:

1. Representação Vetorial de Posição e Elementos: O vetor posição $\rho$ de um planeta é expresso como uma função da longitude excêntrica $F$ e dos elementos orbitais. Ao utilizar a conservação do vetor de Laplace–Runge–Lenz e a fórmula de rotação de Rodrigues, o vetor posição é decomposto em um componente $\rho_0$ (representando o limite circular e coplanar) e uma perturbação $\delta$.
2. Decomposição do Plano Invariante: A derivação é conduzida dentro de um referencial inercial onde o eixo vertical $\hat{k}$ alinha-se com o momento angular total do sistema (o plano invariante). O vetor posição $\rho$ é mostrado como dependente de três vetores primários situados neste plano invariante:
   • $\rho_0$: O vetor de posição adimensional na posição circular.
   • $z$: Um vetor que codifica a excentricidade e a longitude do periélio.
   • $\zeta$: Um vetor que codifica a inclinação e a longitude do nó ascendente.
3. Expansão da Função Perturbadora: A parte direta da função perturbadora é expandida usando polinômios de Legendre. Ao tratar a perturbação $\epsilon$ (derivada de $\delta$) como pequena em relação ao vetor de separação $\nu$, a expansão é sistematicamente reescrita inteiramente em termos de produtos escalares envolvendo o conjunto de vetores $\mathcal{V} = \{\rho_0, \rho'_0, z, z', \zeta, \zeta'\}$.
4. Média Secular e Não Secular:
   • Parte Secular: A função perturbadora secular é obtida através da média sobre as longitudes médias. O autor deriva uma proposição geral (Proposição 1) para a média de produtos de produtos escalares envolvendo $\rho_0$. Isso permite que os termos seculares sejam expressos unicamente como produtos escalares dos vetores de variação lenta ($z, \zeta$).
   • Desigualdades Não Seculares: Para harmônicos de Fourier arbitrários (desigualdades), o método utiliza expansões de Jacobi–Anger para expressar os fatores exponenciais das longitudes médias em termos da longitude excêntrica. Isso permite que qualquer desigualdade específica seja computada diretamente como uma função de produtos escalares envolvendo os vetores $\rho_0, z, \zeta$ e os vetores de longitude média $\ell, \ell'$.
5. Transformação para Vetores Físicos: Finalmente, os vetores $z$ e $\zeta$ são invertidos para expressar as expansões em termos do vetor de momento angular adimensional $\mathbf{h}$ e do vetor de excentricidade $\mathbf{e}$. Este passo recupera a forma compacta do Hamiltoniano secular encontrada na literatura anterior e estende-a para os termos não seculares.

Principais Contribuições e Resultados

• Formalismo Vetorial Unificado: O artigo demonstra que tanto a parte secular quanto qualquer harmônico de Fourier da função perturbadora podem ser expressos inteiramente através de produtos escalares de vetores situados no plano invariante do sistema.
• Algoritmo Geral para Desigualdades: Um algoritmo sistemático é fornecido para computar qualquer desigualdade específica (harmônico de Fourier) sem a necessidade de determinar toda a expansão de Fourier da função perturbadora. Isso é alcançado expressando as integrais em termos de coeficientes de Laplace e produtos escalares de vetores.
• Expansões Explícitas: O autor fornece expansões explícitas para a função perturbadora secular até o grau 4 e para ressonâncias de movimento médio específicas (1:1, 2:1, 3:2) até o grau 2, todas expressas em termos de produtos escalares de $\mathbf{h}$ e $\mathbf{e}$.
• Recuperação de Resultados Conhecidos: A expansão secular derivada coincide com a forma compacta estabelecida anteriormente por Abdullah e Boué & Fabrycky, validando a abordagem vetorial.
• Independência de Referencial: As expressões resultantes são independentes da escolha específica do referencial, desde que o eixo vertical esteja alinhado com o momento angular total. O artigo observa que, embora o vetor $\hat{k}$ desapareça das expressões seculares finais, ele permanece implícito na definição dos vetores de longitude média $\ell$ e $\ell'$.

Significância e Alegações
O artigo alega que a principal significância desta formalização baseada em vetores reside em sua invariância rotacional. Como as expansões são construídas inteiramente a partir de produtos escalares de vetores no plano invariante, elas podem ser avaliadas em qualquer sistema de referência sem a necessidade de rederivar os coeficientes. Isso é particularmente vantajoso para modelar sistemas planetários que estão inclinados em relação a um plano de referência ou que interagem com companheiros inclinados, onde o vetor de momento angular total pode não ser constante ou não estar alinhado com um eixo inercial fixo.

O autor enfatiza que esta abordagem facilita a construção de modelos dinâmicos simplificados, como os de ressonância de movimento médio, ao reduzir significativamente o número de termos que precisam ser armazenados e processados. O trabalho baseia-se em sistemas de álgebra computacional (especificamente o TRIP) para lidar com a manipulação simbólica de vetores, sugerindo que o método é computacionalmente robusto para expansões de alta ordem.