Lyapunov Orbit

Author: CislunarSpace
Website: https://cislunarspace.cn

Definition

A Lyapunov orbit is a family of periodic orbits lying in the plane near a libration point, named after Russian mathematician Aleksandr Lyapunov. Lyapunov orbits are the in-plane counterparts of Halo orbits -- when the $z$ -direction amplitude of a Halo orbit approaches zero, the three-dimensional Halo orbit degenerates into a planar Lyapunov orbit. Lyapunov orbits serve as the foundation for studying libration point dynamics, providing the theoretical starting point for understanding more complex three-dimensional orbits. Lyapunov orbits belong to the Libration Point Orbit (LPO) family, together with Halo orbits, vertical orbits, and axial orbits.

Key Elements

Dynamic Characteristics of Lyapunov Orbits

Key characteristics of Lyapunov orbits in the CR3BP framework include:

Planar motion: Lyapunov orbits lie strictly in the $xOy$ plane with no $z$ -direction motion component
Periodicity: The orbits are precisely closed periodic orbits, forming closed curves in the synodic reference frame
Symmetry: Lyapunov orbits are symmetric about the $x$ -axis; when crossing the $x$ -axis, the $y$ -direction velocity is zero
Orbit shape: Near the libration point, the shape is approximately elliptical; as amplitude increases, the shape gradually distorts, with the side far from the libration point becoming pointed or twisted

Lyapunov orbit families are parameterized by the initial displacement $x_0$ on the $x$ -axis (relative to the libration point). When $x_0$ is small, the orbit approximates linearized simple harmonic oscillation; as $x_0$ increases, nonlinear effects become significant and the orbit shape deviates from elliptical.

Classification of Lyapunov Orbits

Based on libration point locations in the Earth-Moon system, Lyapunov orbits can be classified into the following families:

Orbit Family	Libration Point	Characteristics
Ly1 (Lyapunov L1)	L1	Located between Earth and Moon along the line connecting them
Ly2 (Lyapunov L2)	L2	Located on the far side of the Moon, away from Earth
Ly3 (Lyapunov L3)	L3	Located on the far side of Earth, away from the Moon

Orbital Parameter Characteristics

For the Earth-Moon system, the main parameter ranges of Lyapunov orbit families (based on the dynamic catalog statistics by Guzzetti et al.) are as follows:

Orbit Family	Jacobi Constant Range	Period Range (days)	Stability Index
Ly1	~3.188 – 2.124	Moderate	Relatively high
Ly2	~3.172 – 2.881	Moderate	Relatively high
Ly3	~2.5	Relatively long	Moderate

The Jacobi constant ranges of L1 and L2 Lyapunov orbits serve as important reference benchmarks for orbit comparison and selection in cislunar space. When the Jacobi constants of other orbit families (such as DRO, DPO, axial orbits, etc.) overlap with those of Lyapunov orbits, it suggests the existence of low-cost transfer corridors.

Linearized Analysis of Lyapunov Orbits

Near the libration point, the linearized CR3BP equations of motion have the following eigenvalue structure in the plane:

\lambda_{1,2} = \pm \sigma, \quad \lambda_{3,4} = \pm i\omega

where $\sigma$ is a real eigenvalue (corresponding to stable/unstable manifolds) and $\omega$ is an imaginary eigenvalue (corresponding to periodic oscillation). Lyapunov orbits correspond to motion that excites only the imaginary eigenvalue mode:

\mathbf{x}(t) = A \cos(\omega t + \phi) \mathbf{e}_{\text{center}} + \text{nonlinear corrections}

where $\mathbf{e}_{\text{center}}$ is the direction vector of the center manifold.

Relationship Between Lyapunov and Halo Orbits

There is a profound connection between Lyapunov and Halo orbits:

Degeneration relationship: Halo orbits degenerate into Lyapunov orbits as the $z$ -direction amplitude $A_z \to 0$
Bifurcation structure: In the parameter space of orbit families, Lyapunov orbit families generate Halo orbit families through pitchfork bifurcation
Frequency relationship: Lyapunov orbits involve only the in-plane oscillation frequency $\omega_{xy}$ , while Halo orbits require $\omega_{xy} = \omega_z$
Stability differences: Both are unstable, but the unstable mode structure of Lyapunov orbits is simpler (in-plane only)

Numerical Computation of Lyapunov Orbits

Precise computation of Lyapunov orbits typically employs the following methods:

Linearized initial guess: Using linearized analysis to obtain an approximate analytical solution
Differential correction: Using a shooting method to correct initial conditions so the orbit precisely closes
Parameter continuation: Starting from small-amplitude orbits, gradually increasing amplitude, using each orbit as the initial guess for the next

Application Value

Lyapunov orbits have value in both theoretical research and practical missions:

Foundation for dynamics research: Lyapunov orbits are the foundation for understanding the phase space structure near libration points, and a prerequisite for learning about more complex orbits like Halo and Lissajous
Invariant manifold analysis: The stable and unstable manifolds of Lyapunov orbits form the skeleton of low-energy transfer channels near libration points
Low-energy transfer design: Using the invariant manifolds of Lyapunov orbits, low-energy transfer trajectories connecting different libration point regions can be designed
Poincaré section analysis: Lyapunov orbits are commonly used as reference orbits in Poincaré sections for analyzing the global structure of phase space
Mission design reference: In the dynamic catalog framework proposed by Guzzetti et al., the Jacobi constant ranges of L1/L2 Lyapunov orbits serve as key benchmarks for evaluating the accessibility of other orbit families
Education and introduction: As the simplest periodic orbit family at libration points, Lyapunov orbits are an ideal starting point for orbital mechanics education

References

Richardson D L. Analytic construction of periodic orbits about the collinear points[J]. Celestial Mechanics, 1980, 22(3): 241-253.
Szebehely V. Theory of Orbits: The Restricted Problem of Three Bodies[M]. Academic Press, 1967.
Gomez G, Masdemont J, Simo C. Quasihalo orbits associated with libration points[J]. Journal of the Astronautical Sciences, 1998.
Guzzetti D, Bosanac N, Howell K C. A framework for efficient trajectory comparisons in the Earth-Moon design space[C]. AAS/AIAA Space Flight Mechanics Meeting, 2014.