🔢 mathematics

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

This paper presents a derivation of the planetary Hamiltonian for small eccentricities and inclinations using a vector formalism, demonstrating that secular and Fourier terms can be expressed through scalar products of angular momentum and eccentricity vectors within the system's invariant plane.

Original authors: Federico Mogavero

Published 2026-07-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Federico Mogavero

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the solar system not as a collection of lonely rocks, but as a grand, chaotic dance floor where planets are constantly bumping into each other's gravitational vibes. For a long time, scientists trying to predict how these planets move over millions of years have had to do some very messy math. They usually break the problem down into two ways: either looking at how far apart the planets are (like measuring the size of the dance floor), or looking at how wobbly their paths are (how eccentric they are) and how tilted they are relative to each other.

The paper you're asking about, written by Federico Mogavero, tackles that second, wobbly way. Specifically, it looks at what happens when planets have small eccentricities (nearly circular orbits) and small mutual inclinations (orbits that are nearly flat, like a stack of plates).

The Old Way vs. The New Vector Way

Traditionally, when scientists tried to expand the "disturbing function" (a fancy term for the gravitational tug-of-war between planets), they used a lot of scalar quantities. Think of scalars as just numbers—like saying "the planet is 5 units away" or "the tilt is 10 degrees." It works, but it's like trying to describe a 3D dance move using only a list of numbers. You lose the sense of direction and how the pieces fit together in space.

Mogavero's main finding is a fresh, cleaner way to do this math using vectors. Instead of just numbers, he treats the orbits as arrows pointing in specific directions. He shows that the entire gravitational interaction can be described using scalar products (a specific way of multiplying vectors that tells you how much they point in the same direction) of arrows lying on a special flat surface called the invariant plane.

Think of the invariant plane as the "dance floor" that the whole system naturally wants to stay on. Mogavero proves that you don't need to worry about the messy 3D angles anymore; you can just look at how these arrows (representing the planet's speed and wobble) interact with each other on that flat floor.

What This Paper Rules Out (and What It Doesn't)

It's important to know what this paper doesn't do. Mogavero isn't saying the old scalar math is "wrong" or that it can't be done. He isn't arguing against the existence of the planets or the laws of gravity.

However, he is explicitly moving away from the standard method of using complex scalar orbital elements (like specific angles and distances) for this particular type of expansion. He argues that the vector approach is superior because it is frame-independent. This means the math looks the same no matter how you rotate your head or your telescope. If you use the old scalar way, you often have to pick a specific "up" direction, and if you change that direction, the math gets messy. Mogavero's vector method stays clean and simple regardless of how you look at it.

He also clarifies that this specific expansion is only for small eccentricities and inclinations. If a planet is flying around in a crazy, highly tilted, or very stretched-out orbit, this specific "small wobble" math doesn't apply. He isn't trying to solve the problem for wild, chaotic orbits; he's perfecting the math for the calm, nearly flat ones.

The "Secret Sauce": Angular Momentum and Eccentricity Vectors

The paper goes a step further. It shows how to translate these vector arrows into two very specific, physical concepts: the angular momentum vector (which points in the direction the planet is spinning) and the eccentricity vector (which points toward the closest part of the orbit).

Mogavero demonstrates that you can write the entire gravitational interaction using just the dot products of these two vectors. This is a big deal because it connects the math directly to the physical shape and spin of the orbit. It's like realizing that instead of calculating the dance moves based on the music's tempo and the dancer's shoe size, you can just look at the dancer's momentum and their specific wobble.

How Sure Are We?

The paper is very confident in its results, but it's important to understand how that confidence was built. Mogavero didn't just guess or simulate a few examples; he provided a mathematical derivation. He built a logical proof from the ground up, starting with the basic laws of motion and using vector algebra to show that the expansion must work this way.

He also used a computer algebra system called TRIP to check his work. He didn't just run a simulation to see if it looked right; he used the computer to manipulate the symbolic vectors (the arrows with letters on them) to ensure the math held up to high degrees of complexity (up to degree 4 and 2 in different sections). So, when he says the secular part (the long-term average effect) can be expressed this way, it's a proven mathematical fact within the limits of the small eccentricity assumption.

The Takeaway

In short, Mogavero has taken a notoriously messy problem in planetary dynamics—the math of how planets tug on each other when they are nearly circular and flat—and rewritten it in a language of arrows. This new language is simpler, more elegant, and doesn't care about which way you are facing. It proves that the long-term dance of the planets can be understood entirely by looking at how their spin and wobble vectors interact on the system's invariant plane. It's a tool that makes the math of the solar system's slow, steady evolution much easier to handle, especially for those who want to build simplified models of how these cosmic dancers move over eons.

Technical Summary: Expansion of the Planetary Hamiltonian for Small Eccentricities and Inclinations: A Vector Formalism

Problem Statement
The long-term evolution of planetary orbits is governed by the disturbing function, which encapsulates gravitational interactions between planets. Classical developments in planetary dynamics typically expand this function as a power series in orbital elements (eccentricities and inclinations). While expansions in semi-major axis ratios are often handled using vector formalisms (e.g., Legendre polynomials), expansions in eccentricities and inclinations traditionally rely on scalar quantities. Although previous work by Abdullah (later extended by Boué and Fabrycky) successfully restated the secular part of the disturbing function using spin-like variables derived from angular momentum and eccentricity vectors, a general derivation for the full Hamiltonian—including non-secular Fourier harmonics (inequalities)—using a purely vector-based approach from the outset has been lacking.

Methodology
The author presents a derivation of the planetary Hamiltonian expansion that employs a vector formalism throughout, avoiding the initial reliance on scalar orbital elements. The methodology proceeds as follows:

Vector Representation of Position and Elements: The position vector $\rho$ of a planet is expressed as a function of the eccentric longitude $F$ and orbital elements. By utilizing the conservation of the Laplace–Runge–Lenz vector and Rodrigues' rotation formula, the position vector is decomposed into a component $\rho_0$ (representing the circular, coplanar limit) and a perturbation $\delta$ .
Invariant Plane Decomposition: The derivation is conducted within an inertial reference frame where the vertical axis $\hat{k}$ $\hat{k}$ aligns with the system's total angular momentum (the invariant plane). The position vector $\rho$ $ρ$ is shown to depend on three primary vectors lying in this invariant plane:
- $\rho_0$ : The dimensionless position vector in the circular limit.
- $z$ : A vector encoding the eccentricity and longitude of perihelion.
- $\zeta$ : A vector encoding the inclination and longitude of the ascending node.
Expansion of the Disturbing Function: The direct part of the disturbing function is expanded using Legendre polynomials. By treating the perturbation $\epsilon$ (derived from $\delta$ ) as small relative to the separation vector $\nu$ , the expansion is systematically rewritten entirely in terms of scalar products involving the set of vectors $\mathcal{V} = \{\rho_0, \rho'_0, z, z', \zeta, \zeta'\}$ .
Secular and Non-Secular Averaging:
- Secular Part: The secular disturbing function is obtained by averaging over the mean longitudes. The author derives a general proposition (Proposition 1) for the average of products of dot products involving $\rho_0$ . This allows the secular terms to be expressed solely as scalar products of the slow-varying vectors ( $z, \zeta$ ).
- Non-Secular Inequalities: For arbitrary Fourier harmonics (inequalities), the method utilizes Jacobi–Anger expansions to express the exponential factors of the mean longitudes in terms of the eccentric longitude. This allows any specific inequality to be computed directly as a function of scalar products involving the vectors $\rho_0, z, \zeta$ and the mean longitude vectors $\ell, \ell'$ .
Transformation to Physical Vectors: Finally, the vectors $z$ and $\zeta$ are inverted to express the expansions in terms of the dimensionless angular momentum vector $\mathbf{h}$ and the eccentricity vector $\mathbf{e}$ . This step recovers the compact secular Hamiltonian form found in previous literature and extends it to non-secular terms.

Key Contributions and Results

Unified Vector Formalism: The paper demonstrates that both the secular part and any Fourier harmonic of the disturbing function can be expressed entirely through scalar products of vectors lying in the system's invariant plane.
General Algorithm for Inequalities: A systematic algorithm is provided to compute any specific inequality (Fourier harmonic) without needing to determine the entire Fourier expansion of the disturbing function. This is achieved by expressing the integrals in terms of Laplace coefficients and scalar products of vectors.
Explicit Expansions: The author provides explicit expansions for the secular disturbing function up to degree 4 and for specific mean-motion resonances (1:1, 2:1, 3:2) up to degree 2, all expressed in terms of scalar products of $\mathbf{h}$ and $\mathbf{e}$ .
Recovery of Known Results: The derived secular expansion coincides with the compact form previously established by Abdullah and Boué & Fabrycky, validating the vector approach.
Frame Independence: The resulting expressions are independent of the specific choice of the reference frame, provided the vertical axis is aligned with the total angular momentum. The paper notes that while the vector $\hat{k}$ disappears from the final secular expressions, it remains implicit in the definition of the mean longitude vectors $\ell$ and $\ell'$ .

Significance and Claims
The paper claims that the primary significance of this vector-based formalism lies in its rotational invariance. Because the expansions are constructed entirely from scalar products of vectors in the invariant plane, they can be evaluated in any reference system without re-deriving the coefficients. This is particularly advantageous for modeling planetary systems that are tilted relative to a reference plane or interact with inclined companions, where the total angular momentum vector may not be constant or aligned with a fixed inertial axis.

The author emphasizes that this approach facilitates the construction of simplified dynamical models, such as those for mean-motion resonances, by significantly reducing the number of terms that need to be stored and processed. The work relies on computer algebra systems (specifically TRIP) to handle the symbolic manipulation of vectors, suggesting that the method is computationally robust for high-order expansions.

The Old Way vs. The New Vector Way

What This Paper Rules Out (and What It Doesn't)

The "Secret Sauce": Angular Momentum and Eccentricity Vectors

How Sure Are We?

The Takeaway

More like this