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

Resumen Técnico: Expansión del Hamiltoniano Planetario para Pequeñas Excentricidades e Inclinaciones: Un Formalismo Vectorial

Planteamiento del Problema
La evolución a largo plazo de las órbitas planetarias está gobernada por la función perturbadora, la cual encapsula las interacciones gravitatorias entre los planetas. Los desarrollos clásicos en la dinámica planetaria suelen expandir esta función como una serie de potencias en los elementos orbitales (excentricidades e inclinaciones). Si bien las expansiones en razones de semieje mayor suelen manejarse utilizando formalismos vectoriales (por ejemplo, polinomios de Legendre), las expansiones en excentricidades e inclinaciones dependen tradicionalmente de cantidades escalares. Aunque el trabajo previo de Abdullah (posteriormente extendido por Boué y Fabrycky) logró reformular la parte secular de la función perturbadora utilizando variables de tipo espín derivadas de los vectores de momento angular y excentricidad, ha faltado una derivación general para el Hamiltoniano completo —incluyendo los armónicos de Fourier no seculares (inequalities)— utilizando un enfoque puramente basado en vectores desde el inicio.

Metodología
El autor presenta una derivación de la expansión del Hamiltoniano planetario que emplea un formalismo vectorial en todo momento, evitando la dependencia inicial de los elementos orbitales escalares. La metodología procede de la siguiente manera:

1. Representación Vectorial de la Posición y los Elementos: El vector de posición $\rho$ de un planeta se expresa en función de la longitud excéntrica $F$ y los elementos orbitales. Al utilizar la conservación del vector de Laplace–Runge–Lenz y la fórmula de rotación de Rodrigues, el vector de posición se descompone en un componente $\rho_0$ (que representa el límite circular y coplanar) y una perturbación $\delta$.
2. Descomposición del Plano Invariante: La derivación se lleva a cabo dentro de un marco de referencia inercial donde el eje vertical $\hat{k}$ se alinea con el momento angular total del sistema (el plano invariante). Se demuestra que el vector de posición $\rho$ depende de tres vectores primarios que yacen en este plano invariante:
   • $\rho_0$: El vector de posición adimensional en el límite circular.
   • $z$: Un vector que codifica la excentricidad y la longitud del perihelio.
   • $\zeta$: Un vector que codifica la inclinación y la longitud del nodo ascendente.
3. Expansión de la Función Perturbadora: La parte directa de la función perturbadora se expande utilizando polinomios de Legendre. Al tratar la perturbación $\epsilon$ (derivada de $\delta$) como pequeña en relación con el vector de separación $\nu$, la expansión se reescribe sistemáticamente utilizando únicamente productos escalares que involucran al conjunto de vectores $\mathcal{V} = \{\rho_0, \rho'_0, z, z', \zeta, \zeta'\}$.
4. Promedio Secular y No Secular:
   • Parte Secular: La función perturbadora secular se obtiene promediando sobre las longitudes medias. El autor deriva una proposición general (Proposición 1) para el promedio de productos de productos escalares que involucran a $\rho_0$. Esto permite que los términos seculares se expresen únicamente mediante productos escalares de los vectores de variación lenta ($z, \zeta$).
   • Inequalities No Seculares: Para armónicos de Fourier arbitrarios (inequalities), el método utiliza expansiones de Jacobi–Anger para expresar los factores exponenciales de las longitudes medias en términos de la longitud excéntrica. Esto permite que cualquier inequality específica se compute directamente como una función de los productos escalares de los vectores $\rho_0, z, \zeta$ y los vectores de longitud media $\ell, \ell'$.
5. Transformación a Vectores Físicos: Finalmente, los vectores $z$ y $\zeta$ se invierten para expresar las expansiones en términos del vector de momento angular adimensional $\mathbf{h}$ y el vector de excentricidad $\mathbf{e}$. Este paso recupera la forma compacta del Hamiltoniano secular encontrada en la literatura previa y la extiende a los términos no seculares.

Contribuciones Clave y Resultados

• Formalismo Vectorial Unificado: El artículo demuestra que tanto la parte secular como cualquier armónico de Fourier de la función perturbadora pueden expresarse enteramente a través de productos escalares de vectores que yacen en el plano invariante del sistema.
• Algoritmo General para Inequalities: Se proporciona un algoritmo sistemático para computar cualquier inequality específica (armónico de Fourier) sin necesidad de determinar toda la expansión de Fourier de la función perturbadora. Esto se logra expresando las integrales en términos de coeficientes de Laplace y productos escalares de vectores.
• Expansiones Explícitas: El autor proporciona expansiones explícitas para la función perturbadora secular hasta grado 4 y para resonancias de movimiento medio específicas (1:1, 2:1, 3:2) hasta grado 2, todas expresadas en términos de productos escalares de $\mathbf{h}$ y $\mathbf{e}$.
• Recuperación de Resultados Conocidos: La expansión secular derivada coincide con la forma compacta establecida previamente por Abdullah y Boué & Fabrycky, validando el enfoque vectorial.
• Independencia de Marco de Referencia: Las expresiones resultantes son independientes de la elección específica del marco de referencia, siempre que el eje vertical esté alineado con el momento angular total. El artículo señala que, si bien el vector $\hat{k}$ desaparece de las expresiones seculares finales, permanece implícito en la definición de los vectores de longitud media $\ell$ y $\ell'$.

Significancia y Reivindicaciones
El autor sostiene que la principal significancia de este formalismo basado en vectores radica en su invariancia de rotación. Debido a que las expansiones se construyen enteramente a partir de productos escalares de vectores en el plano invariante, pueden evaluarse en cualquier sistema de referencia sin necesidad de re-derivar los coeficientes. Esto es particularmente ventajoso para modelar sistemas planetarios que están inclinados respecto a un plano de referencia o que interactúan con compañeros inclinados, donde el vector de momento angular total puede no ser constante o no estar alineado con un eje inercial fijo.

El autor enfatiza que este enfoque facilita la construcción de modelos dinámicos simplificados, tales como los de resonancia de movimiento medio, al reducir significamente el número de términos que deben almacenarse y procesarse. El trabajo depende de sistemas de álgebra computacional (específicamente TRIP) para manejar la manipulación simbólica de vectores, lo que sugiere que el método es computacionalmente robusto para expansiones de alto orden.