ON the Stability of Triangular Lagrangian Points in the Spatial Restricted Three-Body Problem

Basic Information

• Paper ID: 2510.13388
• Title: ON the stability of triangular Lagrangian points in the spatial restricted three-body problem
• Author: S. P. Sosnitskii (Institute of Mathematics of National Academy of Sciences of Ukraine)
• Classification: astro-ph.SR (Solar and Stellar Astrophysics)
• Publication Date: October 15, 2025 (arXiv preprint)
• Paper Link: https://arxiv.org/abs/2510.13388

Abstract

This paper represents a development of the author's earlier research, investigating the stability of triangular Lagrangian points in the spatial circular restricted three-body problem and improving upon the results of previous work. Unlike earlier work that established the instability of Lagrangian points based on simplified approximate equations, this paper successfully employs a novel methodology to prove the instability of triangular Lagrangian points based on a closed, complete system of equations. The paper also discusses the relationship between Lyapunov stability and Birkhoff stability (formal stability).

Research Background and Motivation

Importance of the Problem

1. A Classical Problem with Historical Significance: The stability problem of triangular Lagrangian points dates back to the 19th century and remains incompletely resolved
2. Theoretical Limitations: Although KAM theory can address the planar circular restricted problem, it fails to yield expected results when applied to the spatial circular problem
3. Methodological Challenges: Existing methods primarily rely on simplified approximate equations, lacking rigorous proofs based on complete equation systems

Research Motivation

The author aims to discover new research pathways by exploiting the deeper properties of the Jacobi integral, establishing special-form equations for perturbed motion, and thereby effectively resolving stability problems within the existing framework of stability theory.

Core Contributions

1. Methodological Innovation: Proposes a novel approach based on closed, complete equation systems, overcoming the limitations of previous methods that relied on simplified approximate equations
2. Theoretical Proof: Rigorously proves the Lyapunov instability of the triangular Lagrangian points L₄ and L₅ in the spatial circular restricted three-body problem
3. Integral Utilization: Innovatively exploits the intrinsic structure of the Jacobi integral to construct more precise equations of perturbed motion
4. Stability Theory: Clarifies the relationship between Lyapunov stability and Birkhoff stability

Detailed Methodology

Problem Definition

Investigation of the stability near the triangular Lagrangian points L₄ and L₅ of a third particle with negligible mass when two primary bodies with masses m₁ and m₂ move in circular orbits in the spatial circular restricted three-body problem.

Model Architecture

1. Basic Equations of Motion

In the rotating coordinate system, the equations of motion are:

x'' - 2y' = x - (1-μ)(x-μ)/ρ₁₃³ - μ(x+1-μ)/ρ₂₃³
y'' + 2x' = y - (1-μ)y/ρ₁₃³ - μy/ρ₂₃³  
z'' = -(1-μ)z/ρ₁₃³ - μz/ρ₂₃³

where μ = m₂/(m₁+m₂), and ρ₁₃ and ρ₂₃ denote the distances from the small particle to the two primary bodies, respectively.

2. Distance Equation System

The author introduces an innovative distance equation system (10), comprising:

• Second-order derivative equations for squared distances
• Evolution equations for energy variables E₁₃, E₂₃
• First-order relation for the y-coordinate

3. Perturbation Variable Transformation

Perturbation variables are defined as:

x₁ = ρ₁₃² - 1, x₂ = ρ₂₃² - 1
y₁ = E₁₃ + 1, y₂ = E₂₃ + 1
η = y - y₀

Technical Innovations

1. Innovative Utilization of the Jacobi Integral

The author expresses the Jacobi integral in three equivalent forms (26)-(28), with the most compact form being:

(1-μ)y₁ + μy₂ + (1-2μ)η' + 2y₀ξ' + 2(1-2μ)ξ - 4y₀η - 2(ξη'-ηξ'+ξ²+η²) = 2h*

2. Establishment of Key Lemmas

• Lemma 1: Transforms system (19) into a form containing the collective variable v = (1-μ)y₁ + μy₂
• Lemma 2: Establishes the key equality ζ'' = 4(h-ζ²) + O(||q*||²) + O(||q||³)

3. Variable Substitution Strategy

Introduces a new variable ζ̃ = ζ² - h*, transforming the system into a form more suitable for stability analysis (42).

Experimental Setup

Theoretical Analysis Framework

The paper employs a purely theoretical analytical approach without numerical experiments, primarily based on:

1. Characteristic Equation Analysis: Investigation of characteristic roots of the linearized system
2. Lyapunov Stability Theory: Analysis of nonzero characteristic Lyapunov exponents
3. Conservative System Properties: Utilization of the conservative nature of the original system

Parameter Constraints

The analysis is conducted under the constraint 27μ(1-μ) ≤ 1, which covers the parameter ranges of practical celestial mechanics problems.

Experimental Results

Main Theoretical Results

Structure of the Characteristic Equation

The characteristic equation of the linearized system (42) can be factored as:

[3/4·λ³ + y₀(1-2μ)λ² - (5-8μ+8μ²)λ + 8(1-μ+μ²) + 4y₀(1-2μ)] × [λ⁴ + λ² + 27/4·μ(1-μ)] = 0

Proof of Instability

Theorem: Within the framework of the spatial circular restricted three-body problem, the Lagrangian points L₄ and L₅ are unstable in the Lyapunov sense.

Key points of the proof:

1. All three roots of the cubic equation (44) possess nonzero real parts
2. System (42) exhibits nonzero characteristic Lyapunov exponents
3. Due to the conservative nature of the original system, the equilibrium is unstable

Theoretical Findings

1. Separation of Stability and Resonance: Instability is unrelated to resonance phenomena and primarily reflects the characteristics of KAM theory itself
2. Dual Stability Phenomenon: In the spatial circular restricted three-body problem, formal stability (Birkhoff stability) and Lyapunov instability can coexist
3. Existence of Attracting Trajectories: There exist trajectories that are attracted to the Lagrangian points as τ→∞ and τ→-∞

Related Work

Historical Development

1. 19th Century Foundations: Early work by Gascheau (1843) and Routh (1875)
2. 20th Century Theoretical Development: Establishment of Kolmogorov-Arnold-Moser theory
3. Modern Research: Contributions by Deprit & Deprit-Bartholomé (1967), Leontovich (1962), Markeev (1978), and others

Uniqueness of This Work

Compared to existing work, this paper's advantages include:

1. Independence from simplified approximate equations
2. No restrictions on the Jacobi integral constant h*
3. Provision of rigorous proofs based on complete equation systems

Conclusions and Discussion

Main Conclusions

1. Establishment of Instability: Rigorously proves the Lyapunov instability of triangular Lagrangian points L₄ and L₅
2. Effectiveness of the Method: Successfully simplifies stability analysis through clever selection of relevant variables and exploitation of the intrinsic resources of the Jacobi integral
3. Completeness of Theory: Based on closed, complete equation systems, avoiding the limitations of approximate methods

Limitations

1. Parameter Restrictions: The analysis assumes 27μ(1-μ) ≤ 1, which, while covering practical cases, remains restrictive
2. Theoretical Nature: Pure theoretical analysis lacking numerical verification
3. System-Specific: Results are specific to the circular restricted three-body problem

Future Directions

1. Extension to more general parameter ranges
2. Investigation of analogous phenomena in the elliptic restricted three-body problem
3. Exploration of numerical methods to verify theoretical results

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Provides rigorous mathematical proofs based on complete equation systems
2. Methodological Innovation: Innovatively exploits the deep structure of the Jacobi integral
3. Problem Importance: Addresses a classical problem in celestial mechanics
4. Clarity of Exposition: Mathematical derivations are logically clear and arguments are rigorous

Weaknesses

1. Limited Scope of Application: Results apply only to specific parameter ranges and problem settings
2. Absence of Numerical Verification: Pure theoretical analysis lacks support from numerical experiments
3. Physical Interpretation: Physical mechanisms underlying instability require deeper explanation

Impact

1. Theoretical Contribution: Provides new analytical tools for stability theory in celestial mechanics
2. Methodological Value: The proposed method may be applicable to other similar dynamical systems problems
3. Practical Significance: Offers guidance for space mission design and celestial mechanics applications

Applicable Scenarios

1. Celestial Mechanics Research: Orbital design and stability analysis near Lagrangian points
2. Dynamical Systems Theory: Stability studies of Hamiltonian systems with multiple degrees of freedom
3. Space Engineering: Spacecraft mission planning utilizing Lagrangian points

References

The paper cites important literature in the field, including:

• Arnold, V. I. (1963): Foundational work on KAM theory
• Kolmogorov, A. N. (1954): Pioneering research on preservation of conditionally periodic motion
• Szebehely, V. (1967): Classical textbook on the restricted three-body problem
• Author's own earlier work: Sosnitskii (2008, 2009)

This paper achieves an important theoretical breakthrough on a classical problem in celestial mechanics, rigorously proving the instability of triangular Lagrangian points in the spatial circular restricted three-body problem through innovative mathematical methods, making a significant contribution to theoretical development in this field.