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

Technische Samenvatting: Expansie van de Planetaire Hamiltoniaan voor Kleine Excentriciteiten en Inclinaties: Een Vectoreenformalisme

Probleemstelling
De langetermijn evolutie van planetaire banen wordt beheerst door de verstorende functie, die de gravitationele interacties tussen planeten inkapselt. Klassieke ontwikkelingen in de planetaire dynamica brengen deze functie doorgaans tot een machtsreeks uit in orbitale elementen (excentriciteiten en inclinaties). Hoewel expansies in de verhoudingen van de halve lange as vaak worden afgehandeld met vectoreenformalismen (bijv. Legendre-polynomen), vertrouwen expansies in excentriciteiten en inclinaties traditioneel op scalaire grootheden. Hoewel eerder werk van Abdullah (later uitgebreid door Boué en Fabrycky) er succesvol in is geslaerd de seculaire del van de verstorende functie te herformuleren met behulp van spin-achtige variabelen afgeleid van het impulsmoment en excentriciteitsvectoren, ontbrak een algemene afleiding voor de volledige Hamiltoniaan—inclusclusief niet-seculaire Fourier-harmonieken (ongelijkheden)—met behulp van een puur vectorgebaseerde benadering vanaf het begin.

Methodologie
De auteur presenteert een afleiding van de planetaire Hamiltoniaanse expansie die een vectoreenformalisme gebruikt gedurende het gehele proces, waarbij het initiële vertrouwen op scalaire orbitale elementen wordt vermeden. De methodologie verloopt als volgt:

1. Vectorepresentatie van Positie en Elementen: De positievector $\rho$ van een planeet wordt uitgedrukt als een functie van de excentrische lengte $F$ en orbitale elementen. Door gebruik te maken van de behoudswet van de Laplace–Runge–Lenz-vector en de Rodrigues-rotatieformule, wordt de positievector ontleed in een component $\rho_0$ (die de circulaire, coplanaire limiet vertegenwoordigt) en een perturbatie $\delta$.
2. Invariante Vlak Decompositie: De afleiding vindt plaats binnen een inertieel referentiekader waarin de verticale as $\hat{k}$ samenvalt met het totale impulsmoment van het systeem (het invariante vlak). De positievector $\rho$ blijkt afhankelijk te zijn van drie primaire vectoren die in dit invariante vlak liggen:
   • $\rho_0$: De dimensieloze positievector in de circulaire limiet.
   • $z$: Een vector die de excentriciteit en de lengte van het perihelium codeert.
   • $\zeta$: Een vector die de inclinatie en de lengte van de opgaande knoop codeert.
3. Expansie van de Verstorende Functie: Het directe deel van de verstorende functie wordt geëxpandeerd met behulp van Legendre-polynomen. Door de perturbatie $\epsilon$ (afgeleid van $\delta$) als klein te behandelen ten opzichte van de scheidingsvector $\nu$, wordt de expansie systematisch herschreven in termen van scalaire producten betrokken bij de verzameling vectoren $\mathcal{V} = \{\rho_0, \rho'_0, z, z', \zeta, \zeta'\}$.
4. Seculaire en Niet-Seculaire Middeling:
   • Seculaire Del: De seculaire verstorende functie wordt verkregen door te middelen over de gemiddelde lengtes. De auteur leidt een algemeen propositie af (Propositie 1) voor het gemiddelde van producten van dot-producten betrokken bij $\rho_0$. Dit maakt het mogelijk om de seculaire termen uitsluitend uit te drukken in termen van de scalaire producten van de traag variërende vectoren ($z, \zeta$).
   • Niet-Seculaire Ongelijkheden: Voor willekeurige Fourier-harmonieken (ongelijkheden) maakt de methode gebruik van Jacobi–Anger-expansies om de exponentiële factoren van de gemiddelde lengtes uit te drukken in termen van de excentrische lengte. Dit maakt het mogelijk om elke specifieke ongelijkheid direct te berekenen als een functie van de scalaire producten van de vectoren $\rho_0, z, \zeta$ en de gemiddelde lengtevectoren $\ell, \ell'$.
5. Transformatie naar Fysieke Vectoren: Ten slotte worden de vectoren $z$ en $\zeta$ geïnverteerd om de expansies uit te drukken in termen van de dimensieloze impulsmomentvector $\mathbf{h}$ en de excentriciteitsvector $\mathbf{e}$. Deze stap herstelt de compacte seculaire Hamiltoniaan-vorm gevonden in eerdere literatuur en breidt deze uit naar niet-seculaire termen.

Kernbijdragen en Resultaten

• Verenigd Vectoreenformalisme: Het artikel toont aan dat zowel de seculaire del als elke Fourier-harmoniek van de verstorende functie volledig kan worden uitgedrukt via de scalaire producten van vectoren die in het invariante vlak van het systeem liggen.
• Algemeen Algoritme voor Ongelijkheden: Er wordt een systematisch algoritme geleverd om elke specifieke ongelijkheid (Fourier-harmoniek) te berekenen zonder de gehele Fourier-expansie van de verstorende functie te hoeven bepalen. Dit wordt bereikt door de integralen uit te drukken in termen van Laplace-coëfficiënten en scalaire producten van vectoren.
• Expliciete Expansies: De auteur biedt expliciete expansies voor de seculaire verstorende functie tot graad 4 en voor specifieke gemiddelde bewegingsresonanties (1:1, 2:1, 3:2) tot graad 2, allen uitgedrukt in termen van scalaire producten van $\mathbf{h}$ en $\mathbf{e}$.
• Herstel van Bekende Resultaten: De afgeleide seculaire expansie komt overeen met de compacte vorm die eerder door Abdullah en Boué & Fabrycky is vastgesteld, wat de vectorbenadering valideert.
• Referentiekaderonafhankelijkheid: De resulterende uitdrukkingen zijn onafhankelijk van de specifieke keuze van het referentiekader, mits de verticale as samenvalt met het totale impulsmoment. Het artikel merkt op dat hoewel de vector $\hat{k}$ uit de uiteindelijke seculaire uitdrukkingen verdwijnt, deze impliciet aanwezig blijft in de definitie van de gemiddelde lengtevectoren $\ell$ en $\ell'$.

Betekenis en Claims
De auteur stelt dat de primaire betekenis van dit vectorgebaseerde formalisme ligt in de rotationele invariantie. Omdat de expansies volledig zijn geconstrueerd uit scalaire producten van vectoren in het invariante vlak, kunnen ze in elk referentiekader worden geëvalueerd zonder de coëfficiënten opnieuw te moeten afleiden. Dit is bijzonder voordelig voor het modelleren van planetaire systemen die gekanteld zijn ten opzichte van een referentievlak of interageren met gekantelde begeleiders, waarbij de totale impulsmomentvector mogelijk niet constant is of niet is uitgelijnd met een vaste inertiële as.

De auteur benadrukt dat deze benadering de constructie van vereenvoudigde dynamische modellen faciliteert, zoals die voor gemiddelde bewegingsresonanties, door het aantal termen dat moet worden opgeslagen en verwerkt aanzienlijk te verminderen. Het werk steunt op computer-algebra-systemen (specifiek TRIP) om de symbolische manipulatie van vectoren af te handelen, wat suggereert dat de methode computationeel robuust is voor expansies van hoge orde.