Projective Transformations for Regularized Central-Force Dynamics: Hamiltonian Formulation

This paper presents a Hamiltonian framework for regularizing and linearizing central-force dynamics by introducing a new canonical extension of projective decomposition, which yields closed-form solutions for inverse-square and inverse-cubic forces and is validated numerically for the two-body problem with J2 perturbations.

The Gist

Imagine you are trying to predict the path of a satellite orbiting a planet. In the real world, gravity pulls the satellite in a curve, and if you try to write down the math for this, the equations get messy, non-linear, and very difficult to solve, especially if the satellite gets very close to the planet (where the math can "break" or become infinite).

This paper introduces a new mathematical "magic trick" to make these difficult orbital problems easy to solve. Here is how the authors do it, using simple analogies:

1. The Problem: The Tangled Knot

Think of the standard way of describing a satellite's orbit as a tangled knot of string. The string represents the satellite's position and speed. As the satellite moves, the string twists and turns in complex ways because gravity changes strength depending on how close the satellite is. Solving the motion means untangling this knot, which is hard work.

2. The Solution: A New Perspective (Projective Transformation)

The authors propose changing the "lens" through which we look at the satellite. Instead of looking at the satellite directly in 3D space, they project its position onto a new, slightly larger set of coordinates (4 dimensions instead of 3).

• The Analogy: Imagine you are trying to draw a perfect circle on a piece of paper, but your hand is shaking, making the lines wobbly and hard to control. The authors suggest stepping back and looking at the drawing from a different angle, or perhaps using a special projector that turns that wobbly circle into a perfect, straight line on a wall.
• The "Projective" Part: They use a specific type of math called "projective transformation." Think of this like a camera lens that can stretch and shrink space. By stretching the space in a very specific way, the curved, twisting path of the satellite turns into a simple, straight, or perfectly oscillating line (like a pendulum swinging back and forth).

3. The "Hamiltonian" Twist: Keeping the Rules

In physics, there are strict rules about how energy and momentum behave (called the "Hamiltonian" framework). Many previous methods that simplified the math broke these rules, making the results physically inaccurate.

• The Analogy: Imagine you are rearranging a deck of cards to make a game easier to play. Some people just throw the cards on the floor (breaking the rules). The authors, however, rearrange the cards inside the deck so that the game is easier, but the rules of the deck remain perfectly intact. They created a "canonical transformation," which is a fancy way of saying they rearranged the math without breaking the fundamental laws of physics.

4. The "Knobs" and the Best Setting

The authors didn't just find one way to do this; they found a whole family of ways, controlled by "knobs" (mathematical parameters).

• They tested different settings and found one specific combination (where the knobs are set to -1) that works best.
• Why it's special: This specific setting connects the math directly to the satellite's "local view" (what the satellite sees as up, down, and forward). It separates the satellite's spinning motion (rotation) from its moving-in-and-out motion (radial distance).
  • Rotation: The spinning part becomes a simple, constant rotation (like a clock hand).
  • Distance: The moving-in-and-out part becomes a simple spring-like motion (like a weight on a spring).

5. What This Solves

By using this new method, the authors show that:

• Linearization: The messy, curved equations turn into simple, straight-line equations (linear equations). This is like turning a complex maze into a straight hallway.
• Closed-Form Solutions: Because the equations are now simple, they can write down the exact answer for where the satellite will be at any time without needing a computer to guess step-by-step. It's like having a direct formula instead of a long list of instructions.
• More Than Just Gravity: This trick works not just for standard gravity (Kepler dynamics), but also for slightly more complex gravity models (Manev dynamics) that include tiny relativistic effects.
• Perturbations: They even tested it with a real-world complication: the Earth isn't a perfect sphere; it's slightly squashed (oblate). They showed their method can handle this "squash" (called the $J_2$ perturbation) while keeping the math clean.

Summary

The paper presents a new mathematical tool that takes the difficult, curved problem of satellite orbits and "flattens" it into a simple, straight-line problem. It does this by changing the coordinate system (the map we use) and the time parameter (the clock we use) in a way that respects all the laws of physics. The result is a set of simple equations that can be solved instantly and exactly, offering a clearer and more intuitive way to understand and calculate orbital motion than previous methods.

Technical Summary

Technical Summary: Projective Transformations for Regularized Central-Force Dynamics

Problem Statement
The paper addresses the challenge of regularizing and linearizing the equations of motion for central-force dynamics within a Hamiltonian framework. While the classical two-body problem (Kepler dynamics) is well-understood, its non-linear differential equations present singularities at the origin ($r=0$) and complicate analytical and numerical treatments, particularly when perturbations are introduced. Existing methods, such as the Kustaanheimo-Stiefel (KS) transformation, successfully linearize Kepler dynamics using quaternion-based redundant coordinates and evolution parameter transformations. However, the authors note that projective transformations, which offer a different geometric perspective, have been less established in the Hamiltonian context. Furthermore, standard projective point transformations (e.g., Burdet-Ferrándiz) often require specific choices of parameters that may not yield the most intuitive connection to orbital reference frames or may fail to linearize broader classes of central-force potentials, such as the Manev potential.

Methodology
The authors develop a systematic method to extend point transformations from minimal to redundant coordinates into canonical transformations compatible with Hamiltonian mechanics.

1. Canonical Extension of Redundant Coordinates:
   Starting from Hamilton's principle, the authors derive a general procedure to lift a point transformation $\mathbf{x} = \mathbf{x}(\bar{\mathbf{q}}, t)$ (where $\bar{\mathbf{q}}$ are redundant coordinates subject to constraints) to a canonical transformation in phase space. This involves introducing Lagrange multipliers to handle the constraints, ensuring the preservation of the symplectic structure. The resulting transformation maps the original $2n$-dimensional phase space to a $2(n+m)$-dimensional redundant phase space while maintaining Hamilton's equations.

2. Family of Projective Transformations:
   The authors define a generalized family of projective point transformations for the central-force problem:
   $$\mathbf{r} = u^n q^m \mathbf{q}$$
   subject to the constraint $q = \|\mathbf{q}\| = 1$. Here, $\mathbf{r}$ is the inertial position, $\mathbf{q}$ is a unit vector, and $u$ is a scalar related to the radial distance. The parameters $n$ and $m$ act as "knobs" to generate different canonical extensions.

   • The classic Burdet-Ferrándiz (BF) transformation corresponds to $m=0, n=-1$.
   • The authors identify a preferred transformation where $m = n = -1$, leading to $\mathbf{r} = \frac{1}{u} \hat{\mathbf{q}}$.

3. Preferred Coordinate Set and Physical Interpretation:
   The choice of $m=n=-1$ is motivated by its intuitive physical interpretation. In this set:

   • The coordinate $\mathbf{q}$ corresponds to the radial unit vector $\hat{\mathbf{r}}$.
   • The conjugate momentum $\mathbf{p}$ is orthogonal to $\mathbf{q}$ and relates directly to the angular momentum.
   • The triplet $(\mathbf{q}, \mathbf{p}, \boldsymbol{\ell})$ forms an orthonormal basis that aligns with the Local Vertical, Local Horizontal (LVLH) frame.
   • Crucially, the radial distance $r$ depends only on $u$ ($r = 1/u$), decoupling the central-force potential from the angular dynamics in the Hamiltonian.

4. Time Reparameterization and Linearization:
   To achieve full linearity, the evolution parameter is transformed from time $t$ to new parameters $s$ and $\tau$ (where $dt = r^2 ds = r^2/\ell d\tau$).

   • The parameter $s$ is related to the Burdet time.
   • The parameter $\tau$ corresponds to the true anomaly.
     Using these parameters, the authors demonstrate that the Hamiltonian equations of motion become fully linear for both Kepler and Manev potentials.

Key Contributions and Results

• Linearization of Manev Dynamics: Unlike the KS transformation, which is specific to inverse-square forces, the proposed projective transformation linearizes dynamics for the Manev potential ($V_0 = -k_1/r - k_2/2r^2$). This potential includes an inverse-cubic force term and is relevant for modeling relativistic corrections (perihelion precession) while maintaining integrability.
• Closed-Form Solutions: The paper derives explicit closed-form solutions for the phase space variables $(\mathbf{q}, \mathbf{p}, u, p_u)$ in terms of the evolution parameters. These solutions describe a 4-dimensional linear harmonic oscillator with constant driving forces.
• State Transition Matrices (STM): The authors construct the Kepler State Transition Matrix (STM) in the projective coordinates. They distinguish between the matrix exponential of the system matrix (often called the STM in linear systems) and the true Jacobian of the flow ($\partial \mathbf{x}_\tau / \partial \mathbf{x}_0$), providing the explicit form of the latter which accounts for the dependence on angular momentum.
• Decoupling of Dynamics: The transformation naturally decouples the rotational (angular) motion from the radial motion. The angular dynamics are linear and invariant for any central force, while the radial dynamics become linear specifically for Kepler and Manev potentials.
• Numerical Verification: The method is applied to the $J_2$-perturbed two-body problem (the "main satellite problem"). The authors formulate the perturbed equations in the new coordinates and verify the results numerically, demonstrating the method's capability to handle non-conservative and conservative perturbations.

Significance and Claims
The authors claim that their work provides a simpler and more intuitive alternative to the KS transformation for linearizing central-force dynamics.

• Intuition: The coordinates are directly linked to the LVLH basis and attitude dynamics, offering a clearer geometric interpretation than the quaternion-based KS coordinates.
• Generality: The method extends linearization to the Manev potential, a broader class of forces than those handled by standard KS regularization.
• Robustness: The formulation handles coordinate redundancy transparently via Hamilton's principle, and the linearization holds regardless of whether the extra integrals of motion (constraints) are explicitly enforced during propagation, provided initial conditions are set correctly.
• Utility: The resulting linear, Hamiltonian system is well-suited for symplectic integration algorithms and semi-analytic perturbation methods (e.g., Lie series), as the transformed system preserves the symplectic structure while eliminating non-linearities.

The paper concludes that this projective approach offers a powerful tool for celestial mechanics, particularly for problems requiring high-precision propagation or analytical treatment of perturbations in central-force fields.