异质三轴主星与次级天体修正牛顿力作用下的椭圆型限制性三体问题

《Symmetry》：Motion in the ER3BP Under a Heterogeneous Triaxial Primary and a Modified Newtonian Force of the Secondary

　　

本文研究了在椭圆型限制性三体问题（Elliptic Restricted Three-Body Problem，ER3BP）中，无穷小质量天体在异质三轴主星与产生修正牛顿引力场的球形次级天体引力作用下的动力学行为。两主星绕共同质心沿椭圆轨道运动，第三体质量足够小而不反作用于主星。研究人员建立了旋转—脉动坐标系下的无量纲运动方程，其中主星为含 N 层不同密度 ρi 与半长轴 a1、a2、a3 的三轴刚体，次级天体引力由连分式势参数 α 修正牛顿势。研究人员由两主星间势与力导出系统平均运动（mean motion）n，其退化为 a=1、e=0、三轴与修正项为零时得到经典圆型限制性三体问题（Circular Restricted Three-Body Problem，CR3BP）的 n=1。研究人员解析求出了共线平动点（collinear equilibrium points）与非共线平动点（non-collinear equilibrium points）：非共线点构成 Scalene 三角形，位置由质量参数 μ、三轴系数 A、B 与修正参数相关；共线点在区间划分下可得到多于经典 3 个的点，某些参数下共出现 4 个共线点。稳定性方面，非共线点满足特征方程 λ4+...=0，临界质量比 μc 受偏心率 e、半长轴 a、异质三轴项 A、B 与修正牛顿项 α 共同影响，异质三轴与 e、a 具 destabilizing 趋势，修正牛顿势具 stabilizing 趋势；共线点因特征方程积条件不满足而总是不稳定。数值上研究人员用 Mathematica 展示了势面（potential surfaces）、平动点位置、运动允许区（regions of motion，由 Jacobian 常数 C 判定）、周期轨道（periodic orbits）、牛顿迭代吸引盆（basins of attraction）与线性稳定性根。结果表明：修正牛顿效应使势面扩张，半长轴减小使势面收缩；平动点最多 8 个、最少 2 个；除特定三角形点外平动点多不稳定。该结果有助于理解天体力学与动力天文学中异质三轴天体和修正引力下的小天体运动。

研究背景方面，牛顿之后三体问题被证不可解，学者转向圆型限制性三体问题（CR3BP）与椭圆型限制性三体问题（ER3BP）。已有研究纳入扁率（oblateness）、辐射压（radiation pressure）、光引力（photo-gravitational）、Yarkovsky 效应等，但较少同时考虑主星为多层异质三轴刚体、次级天体采用修正牛顿势，并在椭圆轨道下统一处理平均运动、平动点与稳定性。因此研究人员在《Symmetry》发表该文，以弥补异质三轴主星与修正牛顿次级这一组合模型在 ER3BP 中的空白。

主要关键技术方法上，研究人员采用旋转—脉动参考系下的无量纲方程建模；由两主星间连分式修正势推导平均运动 n；用速度加速度同时为零的代数条件解析求解共线与非共线平动点；通过变分方程与特征方程做线性稳定性分析；数值上使用 Mathematica，基于 Jacobian 积分常数 C 作运动允许区，用相位空间初值积分作周期轨道，用牛顿—拉夫森（N-R）迭代作吸引盆。

研究结果如下。第 2 节 Equations of Motion：研究人员设主星为 N 层密度 ρ_i、半轴 a₁、a₂、a₃ 的三轴异质体，次级为参数 α 修正牛顿势球形体，在 x 轴两主星绕质心运动，给出旋转—脉动系下第三体方程与势函数 Ω，引入修正力参数 α。

第 3 节 Mean Motion：研究人员由主从势 V 与两体椭圆距离 r=a(1-e²)/(1+e cos f) 推导主星与次星运动方程，相加并略去 μ、α、三轴项高阶量，得平均运动 n=n(a,e,A,B,α)；当 a=1、e=0、A=B=α=0 时 n=1，回归 CR3BP。

第 4 节 Analytical Investigations：研究人员解析处理平动点与稳定性。4.1 Equilibrium Points 中由 u=v=w=0 得运动在平面 z=0；非共线点经式 (13)(14) 令 x=x₀+ξ、y=y₀+η，得三角形点 L₄、L₅ 成 scalene 三角形；共线点按 x 区间讨论，经典 3 个外还可出现 L_new1、L_new2，参数使某区间无共线点，总计可达 4 个共线点。4.2 Stability 中，非共线点偏导代入变分方程，z 向周期有界，平面内特征方程 λ⁴+pλ²+q=0，稳定需根纯虚，得 e、μ 条件与临界质量比 μ_c=μ_c^*+Δ(A,B,α,a,e)；异质三轴与 e、a destabilizing，修正牛顿势 stabilizing。共线点因 A_xxA_yy<0 使特征方程 (30) 根积为正，不可能全纯虚，故始终 unstable。

第 5 节 Numerical Explorations：5.1 Potential Surfaces 用势函数绘六类图，异质使主极分裂，修正牛顿使势面扩张，a 减小使势面收缩。5.2 Locations of Equilibrium Points 解 U_x=U_y=0，经典 5 点；均质 7 点；修正牛顿 8 点；e=0.3、0.6、0.9、0.96 分别得 7、8、6、2 点；a 变化亦改点数，说明参数强烈影响位置。5.3 Regions of Motion 由 Jacobian 积分 v²=2Ω-C，取 C 为平动点值，着色区禁止、白色区允许。5.4 Periodic Orbits 由相位空间初值积分得周期 6.32、10.6、9.95、20.1、48.815、16.27 单位轨道。5.5 Basins of Attracting Domain 用 N-R 迭代 (32) 作六例吸引盆，吸引区延至无穷。5.6 Linear Stability 数值解特征方程 (30)，表中至少一根为正实或正实部，故除图 3a 非共线点外均 unstable。

讨论与结论翻译方面，研究人员总结：以异质三轴刚体为主星、修正牛顿势球体为次星的 ER3BP 中，系统动力学显著受物理与轨道参数影响。解析平动点数与数值不一致，源于解析用一阶近似、数值计入极小参数。非共线点在 μ< />_c 稳定、μ>μ_c 不稳定；共线点全质量比范围不稳定。数值上修正牛顿使势面扩张、a 减小使势面收缩；平动点最多 8 最少 2，异于经典 5 点。运动区中彩色为禁区、白色为允区；存在多周期轨道；吸引盆延展至无穷；除部分三角形点外平动点不稳定。该成果为天体力学、动力天文学及小天体运动与稳定性研究提供综合理解。