具有修正牛顿势与伸长主天体及次天体的微扰受限三体问题

《Symmetry》：Perturbed Restricted Three-Body Problem with Elongated Primary and Secondary Having Modified Newtonian Potential

　　

本研究考察了在伸长线段形主天体与服从修正牛顿势的次天体引力作用下，变质量测试粒子的动力学运动性质。研究人员利用Jeans定律与Meshcherskii时空变换，推导了系统运动方程、拟雅可比积分与微扰平均运动。此后给出静止点存在性与稳定性的解析研究，并通过数值手段展示势面投影、静止点位置、稳定性及吸引域。该研究为各类微扰对系统内测试粒子运动的影响提供了有价值的见解，增进了对天体动力学(celestial dynamics)的理解，并为天体力学与航天动力学(astrodynamics)的后续研究提供动机。

研究背景方面，受限三体问题(R3BP)是天体力学与空间动力学中的经典课题，两个有限质量主天体绕共同质心作圆或椭圆运动，第三体为测试粒子，其运动受主天体引力支配且不反作用于主天体。既有文献已引入多种微扰：主天体扁率、三轴性、光压、量子修正势、线段形天体、变质量（Jeans定律）等。然而同时考虑“伸长线段形主天体+服从修正牛顿势的次天体+按Jeans定律变质量的测试粒子”三者耦合微扰的R3BP尚未被系统研究。宇宙中某些小行星与小卫星呈伸长线段状，而牛顿逆平方律的修正（如引入参数α的修正势）亦源于低能量子结构等信息，因此构建此类非标准构型模型对理解真实不规则天体附近动力学具有现实意义。研究人员据此在《Symmetry》发文，旨在解析变质量测试粒子在上述耦合微扰下的运动性质。

主要关键技术方法包括：采用Jeans定律描述测试粒子质量演化m(t)=m₀e^-kt，结合Meshcherskii时空变换将变质量方程归化为常质量框架；推导含线段形主天体引力势（沿长度L积分）与次天体修正牛顿势（V～r^-1+α型）的合成势能函数Ω(x,y)；由Ω梯度为零解析求静止点，通过变分方程获得特征方程λ⁴+Aλ²+B=0判定稳定性（纯虚根稳定）；数值上用Mathematica绘制势面投影、曲线交点定静止点坐标、Newton–Raphson迭代构吸引域。样本并非生物样本，而是参数化数值实验：设质量比μ、线段长L、倾角θ、修正势参α、Jeans参k，分四情形（仅变质量、仅修正势、仅线段、全耦合）计算。

研究结果各小节浓缩如下。

第二节“运动方程”：研究人员建立旋转坐标系下变质量测试粒子方程，经Meshcherskii变换得等效常质量系统方程(4)与拟雅可比积分(6)C=2Ω(x,y)?(?²+?²)/e^2kt；由主从天体间修正势引力导出微扰平均运动n²=1+(3α/2)μ(1?μ)+…，无微扰时n=1。

第三节“解析研究—静止点”：共线静止点由y=0代入方程(12)得高次代数方程，确切个数留待数值；非共线点设x_L±δx, y_T±δy，微扰展开至一阶得偏移量δx,δy表达式，无微扰退化为经典三角点L₄/L₅。

第三节“解析研究—稳定性”：在静止点代入ξ=x?x₀, η=y?y₀得变分方程(14)，系数A_ij为Ω二阶导在平衡点取值，特征方程(16)λ⁴+Aλ²+B=0，由根全纯虚则稳定，否则不稳定。

第四节“数值研究—势面投影”：情形I（变质量）势面双极凹型；情形II（修正势）单极且等高线外扩；情形III（线段）双极凹型随θ变；情形IV（耦合）矩形单极外扩，显示微扰改变势拓扑。

第四节“数值研究—静止点位置”：情形I得5点（3共线+2非共线）同经典；情形II得6点（2共线+4非共线对称于x轴）；情形III中θ=0得5点，θ=π/2时1个共线点下移成非共线，仍5点；情形IV中θ=0得6点，θ=π/2得8点（含新对称点L₆/L₇），表明耦合微扰突破单微扰对称限制。

第四节“数值研究—静止点性质”：查表可知情形I非共线稳、共线不稳；情形II仅一对非共线(L₄型)稳；情形III同II结论；情形IV(θ=0)全6点不稳；情形IV(θ=π/2)中平移点L₆与两对非共线点稳，余不稳。

第四节“数值研究—吸引域”：Newton–Raphson法以颜色区隔收敛到各静止点的初始条件域，各情形下稳定点对应有限色块、不稳点域延伸至无穷，耦合情形θ=π/2出现八色域分割。

讨论与结论翻译：本文研究了主天体为伸长线段、次天体具修正牛顿势、无穷小体按Jeans定律变质量的受限三体问题中的微扰效应。研究人员解析导出了运动方程、拟雅可比积分及其他重要动力学量。数值研究检视了势面投影、静止点位置、稳定性及吸引域。结果表明运动方程与微扰平均运动依赖于线段长度与倾角、修正势参数及变质量参数；无微扰时退化为经典R3BP。势面投影随情形而异；静止点数在变质量与线段情形为5，修正势情形6，耦合情形可达8；稳定性在特定参数组合下可稳可不稳；吸引域分析揭示与各静止点相连且延拓至无穷的区域。该成果有助于天体力学与空间动力学研究者理解非标准引力构型与微扰对天体运动的影响。