In the present paper, which is a development of an earlier study by the author \cite{Sosnitskii08}, we consider the stability of triangular libration points in the spatial circular restricted three-body problem and improve the result of author's work \cite{Sosnitskii08}. Unlike \cite{Sosnitskii08}, where the instability of libration points was established on the base of reduced approximate equations, we succeeded in this paper to use a new approach that made it possible to prove the instability of triangular libration points on the base of a closed complete system of equations. The relationship between the Lyapunov stability and Birkhoff stability (formal stability) is also discussed
ON the stability of triangular Lagrangian points in the spatial restricted three-body problem 论文ID : 2510.13388标题 : ON the stability of triangular Lagrangian points in the spatial restricted three-body problem作者 : S. P. Sosnitskii (Institute of Mathematics of National Academy of Sciences of Ukraine)分类 : astro-ph.SR (Solar and Stellar Astrophysics)发表时间 : 2025年10月15日 (arXiv预印本)论文链接 : https://arxiv.org/abs/2510.13388 本文是作者早期研究的发展,考虑了空间圆形限制性三体问题中三角拉格朗日点的稳定性,并改进了作者之前工作的结果。与之前基于简化近似方程建立拉格朗日点不稳定性的工作不同,本文成功使用了一种新方法,基于封闭完整的方程组证明了三角拉格朗日点的不稳定性。文中还讨论了李雅普诺夫稳定性与伯克霍夫稳定性(形式稳定性)之间的关系。
历史悠久的经典问题 :三角拉格朗日点的稳定性问题可追溯到19世纪,至今仍未完全解决理论局限性 :KAM理论虽然能解决平面圆形限制性问题,但应用于空间圆形问题时无法得到期望结果方法学挑战 :现有方法主要依赖于简化的近似方程,缺乏基于完整方程组的严格证明作者旨在通过挖掘雅可比积分的深层特性,找到新的研究途径,建立扰动运动的特殊形式方程,从而在现有稳定性理论框架内有效解决稳定性问题。
方法学创新 :提出了基于封闭完整方程组的新方法,克服了以往依赖简化近似方程的局限性理论证明 :严格证明了空间圆形限制性三体问题中三角拉格朗日点L₄和L₅的李雅普诺夫不稳定性积分利用 :创新性地利用雅可比积分的内在结构,构建了更精确的扰动运动方程稳定性理论 :阐明了李雅普诺夫稳定性与伯克霍夫稳定性之间的关系研究空间圆形限制性三体问题中,质量为m₁和m₂的两个主体在圆轨道上运动时,第三个质量可忽略的粒子在三角拉格朗日点L₄和L₅附近的稳定性。
在旋转坐标系中,运动方程为:
x'' - 2y' = x - (1-μ)(x-μ)/ρ₁₃³ - μ(x+1-μ)/ρ₂₃³
y'' + 2x' = y - (1-μ)y/ρ₁₃³ - μy/ρ₂₃³
z'' = -(1-μ)z/ρ₁₃³ - μz/ρ₂₃³
其中μ = m₂/(m₁+m₂),ρ₁₃和ρ₂₃分别表示小粒子到两个主体的距离。
作者引入了创新的距离方程系统(10),包含:
距离平方的二阶导数方程 能量变量E₁₃, E₂₃的演化方程 y坐标的一阶关系式 定义扰动变量:
x₁ = ρ₁₃² - 1, x₂ = ρ₂₃² - 1
y₁ = E₁₃ + 1, y₂ = E₂₃ + 1
η = y - y₀
作者将雅可比积分表示为三种等价形式(26)-(28),其中最紧凑的形式为:
(1-μ)y₁ + μy₂ + (1-2μ)η' + 2y₀ξ' + 2(1-2μ)ξ - 4y₀η - 2(ξη'-ηξ'+ξ²+η²) = 2h*
引理1 :将系统(19)转化为包含集体变量v = (1-μ)y₁ + μy₂的形式引理2 :建立关键等式ζ'' = 4(h-ζ²) + O(||q*||²) + O(||q||³)引入新变量ζ̃ = ζ² - h*,将系统转化为更适合稳定性分析的形式(42)。
本文采用纯理论分析方法,不涉及数值实验,主要基于:
特征方程分析 :研究线性化系统的特征根李雅普诺夫稳定性理论 :分析非零特征李雅普诺夫指数保守系统性质 :利用原系统的保守性质分析在约束条件27μ(1-μ) ≤ 1下进行,涵盖了实际天体力学问题的参数范围。
系统(42)的线性化特征方程可分解为:
[3/4·λ³ + y₀(1-2μ)λ² - (5-8μ+8μ²)λ + 8(1-μ+μ²) + 4y₀(1-2μ)] × [λ⁴ + λ² + 27/4·μ(1-μ)] = 0
定理 :在空间圆形限制性三体问题框架下,拉格朗日点L₄和L₅在李雅普诺夫意义下是不稳定的。
证明要点:
三次方程(44)的所有三个根都具有非零实部 系统(42)具有非零特征李雅普诺夫指数 由于原系统的保守性,平衡态不稳定 稳定性与共振的分离 :不稳定性与共振问题无关,主要反映KAM理论本身的特征双重稳定性现象 :空间圆形限制性三体问题中,形式稳定性(伯克霍夫稳定性)与李雅普诺夫不稳定性可以共存吸引轨道的存在 :存在当τ→∞和τ→-∞时都被拉格朗日点吸引的轨道19世纪基础 :Gascheau (1843)和Routh (1875)的早期工作20世纪理论发展 :Kolmogorov-Arnold-Moser理论的建立现代研究 :Deprit & Deprit-Bartholomé (1967), Leontovich (1962), Markeev (1978)等的贡献相比已有工作,本文的优势在于:
不依赖于简化近似方程 对雅可比积分常数h*无限制条件 提供了基于完整方程组的严格证明 不稳定性的确立 :严格证明了三角拉格朗日点L₄和L₅的李雅普诺夫不稳定性方法的有效性 :通过巧妙选择相关变量和利用雅可比积分的内在资源,成功简化了稳定性分析理论的完整性 :基于封闭完整的方程组,避免了近似方法的局限性参数限制 :分析假设27μ(1-μ) ≤ 1,虽然涵盖实际情况但仍有限制理论性质 :纯理论分析,缺乏数值验证特定系统 :结果特定于圆形限制性三体问题扩展到更一般的参数范围 研究椭圆限制性三体问题中的类似现象 探索数值方法验证理论结果 理论严谨性 :提供了基于完整方程组的严格数学证明方法创新性 :创新性地利用雅可比积分的深层结构问题重要性 :解决了天体力学中的经典难题写作清晰度 :数学推导逻辑清晰,论证严密应用范围 :结果仅适用于特定的参数范围和问题设置数值验证缺失 :纯理论分析缺乏数值实验支持物理解释 :对不稳定性的物理机制解释有待深入理论贡献 :为天体力学稳定性理论提供了新的分析工具方法学价值 :所提出的方法可能适用于其他类似的动力系统问题实用意义 :对空间任务设计和天体力学应用具有指导意义天体力学研究 :拉格朗日点附近的轨道设计和稳定性分析动力系统理论 :具有多个自由度的哈密顿系统的稳定性研究空间工程 :利用拉格朗日点的航天器任务规划论文引用了该领域的重要文献,包括:
Arnold, V. I. (1963): KAM理论的基础工作 Kolmogorov, A. N. (1954): 条件周期运动保存的开创性研究 Szebehely, V. (1967): 限制性三体问题的经典教材 作者自己的前期工作 Sosnitskii (2008, 2009) 这篇论文在天体力学的经典问题上取得了重要理论突破,通过创新的数学方法严格证明了空间圆形限制性三体问题中三角拉格朗日点的不稳定性,为该领域的理论发展做出了重要贡献。