12月3号，也就是今天，国产科幻动画《三体》本应该要开播。可惜好事多磨，昨天传来了动画暂缓播出的消息。不过，作为《三体》粉丝，我觉得是时候写点什么纪念一下科幻史上这辉煌灿烂的一笔。
遥想四年前我读完《三体》的那个夏天，那时我还是头发茂密的物理系大一本科生。我第一次感受到科幻不仅能对技术做出展望，还能引发对外星文明、宇宙社会、甚至回归到人性自身的探讨。这次我们先不聊太空电梯、超强度纳米线、智子、引力波通讯、费米悖论、降维打击等天马行空的东西（也许播出后结合剧情一起写写），就先聊聊《三体》的基础背景设定，也就是在三星系统。
三体人的故乡 | 图源：《未来漫游指南》第2集38:16
1、 三星系统和混沌
《三体》假设在距离地球4.3光年的半人马座三星系统中，存在一种外星文明。现实中，这个三星系统是真实存在的。目前的观测结果表明，半人马座α星是质量为的黄绿色主序星，表面温度；半人马座β星是的橙色主序星，表面温度；比邻星是的红矮星，表面温度，它距离前两颗恒星有0.2光年，相当于太阳到奥尔特云的距离。所以目前的状态是，比邻星远远地望着α和β缠缠绵绵。
三星系统与我们的位置关系
2016年，人们发现比邻星具有一个位于宜居带的行星——Proxima Centauri b。它距离比邻星0.05个天文单位，质量至少为地球的1.07倍，半径约为地球的1.3倍，公转周期11.2个地球日。据估算，它的平均表面温度为-39℃。如果大气中存在一些温室气体，还是有可能存在碳基生命的。
已经确认的比邻星b、c和尚待确认的比邻星d | 图源：KOSMO
但我得吐槽一句，半人马座三星系统位于南半球天区，赤纬61°~63°S，而红岸基地位于43°N以北的大兴安岭，所以这三颗星全天都在红岸基地的地平线之下。只有27°N以南的华南地区，才有可能直接收到三体文明的信号。
说回正题，因为这三颗恒星目前可以看成一颗小红星远远地绕着双星系统转，所以万幸在接下来的百万年间，它们的轨道都比较稳定。
三体星系的轨迹示意。从该频道的动画来看轨道还算稳定，周期预计54万年 | 图源：KOSMO
但是！别看刚刚的示意图很有序，但这很可能就像天气预报那样不靠谱。而且，模拟时也会出现这种稳定之下暗流汹涌的结果。
稳定了，但不完全稳定 | 图源[6]
原著中，三体的一轮文明设计了人阵列计算机来计算运行轨迹，结果计算出恒纪元后没几天就变成了乱纪元，最终导致了这轮文明的覆灭。我随手跑了个程序，感觉当时的轨迹应该是类似这样子的。
横纵轴单位为A.U.,本次乱纪元持续约9个世纪
这里借用Zzzyt205做的更漂亮的三体模拟视频。虽然计算机做数值计算会存在精度问题，不等于实际发生的情况，但这个模拟结果足够说明三体星系多不适合生存了。
红黄蓝三颗恒星和一颗白色行星 | 特别感谢：bilibili@Zzzyt205
比如50秒左右，行星离三颗恒星都很近，此时三日凌空，行星变成炼狱。1分钟以后，行星不断远离三颗恒星，那就会出现漫长的冰期。又比如在第29秒和38秒，行星从一颗恒星转移到另一颗恒星，如果轨道高度低于洛希极限，那就会出现大撕裂。
未播出海报，目测应该是正剧中的大撕裂场面
那么，有没有可能把它们的轨道精确地写成解析式呢？早在牛顿时代，人们就发现两个物体在引力作用下的运动可以找到解析解，也就是用精确的数学表达式写出这两个物体的轨迹。这件事并不难，我们经常能在高中物理题里见到弱化版的双星模型。如果你有志于做物理研究的话，这类弱化成正圆轨道的题应该是十拿十稳的送分题。
送分题，请高中生自行学习
但是当问题变成三个物体时，别说高中生，就是一流数理学家，都束手无策。1890年，法国数学家庞加莱发表论文，说明对于三体系统而言，只要系统的初始状态发生一点点改变，就会对后续的运动产生极大的影响。1961年，美国气象学家洛伦茨在研究气象预报时，也发现稍微改变一点输入数据，预测结果就会大相径庭。11年后，他在美国科学发展学会第139次会议上提出著名的蝴蝶效应，“巴西的一只蝴蝶扇动翅膀，就可能在得克萨斯州掀起一场龙卷风”。由此，人们正式对混沌学科开展系统性的研究，三体问题也成了一个世界公认的难题。我们再重申一遍问题：三个物体在平方反比引力的作用下能否找到它们轨迹的解析通解？
目前的结论是不能，代数上只能用无限多的级数去逼近，数值上也只能不断地迭代做有限差分求解。可以说，三体作为一个混沌系统，除了质心是稳定的以外，再也不能找出一点秩序了。
三体系统唯一稳定的质心 | 图源：Wikipedia
不过，我们总可以根据最新的观测结果更新模型的初始数据，从而实现在短期内做到一定程度的准确预报。随着算力的提升，迭代精度也可以不断增加，这在另一方面也能减小数值解的误差。
2、 稳定的特解
其实，三体问题也可以有稳定的周期解。最有代表性的就是欧拉-拉格朗日族、8字族和布鲁克-赫农族这三族(family)解。
8字族
布鲁克-赫农族 | 图源：[4]
不过直到十年前，人们也就只找到了这三族解。注意这里的用词，是三族解，不是三个解。
2013年，塞尔维亚贝尔格莱德物理研究所[7]的MilovanŠuvakov和VeljkoDmitrašinović，利用计算机，成功地找出了13族新的三体问题的周期性特殊解。他们使用的方法，是先从一个已知的周期解出发，然后在计算机模拟中不断调整系统内部的参数，改变初始条件，直到发现新的、稳定的三体运动模式。2015年，AnaHudomal用同样的方法新发现了14族解。
双弯曲三角+大圆|图源[6]
2017年，上海交通大学的廖世俊教授[8]利用天河二号超级计算机，新发现了669个周期解。从这里开始解的发现就不再按族算了，而是按个。比如这个扭曲的8字形轨道和糖果轨道。
扭曲的八字型轨道 | 图源[8]
糖果轨道 | 图源[8]
2018年，他们又新发现1349个解。2021年，他们又基于机器学习搞出了135445个解，而且称数值精度比普朗克尺度还小……果然人算还是不如超算啊。感兴趣的读者可以去参考资料[8]看看其他轨道。
3、 拉格朗日点
如果三体系统中，有一个天体的质量小到相对于另外两个天体可以忽略不计，那么对于这个小质量天体，就存在五个位置，可以让系统稳定地做周期性转动。这就是刚刚提到的欧拉-拉格朗日族中的一种情况，这五个点就是拉格朗日点。下面的推导非常优美，推荐有高一基础的读者看看，不感兴趣的读者可以直接跳到#晕轮轨道小节。
五个拉格朗日点

#

在前面的假设中，第三个物体对另外两个物体的引力可以忽略，所以我们先求解相距为的和的两体运动。
解得
然后再在和之间放上，要求它能在万有引力的作用下一起绕质心公转
这是个关于的五次方程，没有解析的求根公式。所以我们需要上传统艺能——小量近似，。再把和代入，就会发现关于的二次项也会消掉，只剩下一个可以轻松解出的三次项

#

和类似，列出向心运动方程
做传统艺能，，解得
在这种小量近似下，非常优美。

#

这次我们选取的参考点不在地球，而是在相对于的对侧
做传统艺能，，，解得

#和

这部分是我觉得拉格朗日最牛的地方。先想象在可以构成等边三角形的地方放置一个，然后我们反过来证明，和在这个地方产生的引力刚好能为提供向心力。
首先是证明合力方向指向公转中心O。由于，又因为，所以力合成的平行四边形和平行四边形CDOE是相似的，这样就证明了合力指向O点。
接着是证明合力大小刚好提供向心力。你当然可以套余弦公式把合力的大小和的长度算出来，来验证它们是符合圆周运动的运动方程。但更优雅的方式是，验证它们所受到的引力加速度总与公转半径成正比常数。第一个等号已经成立，而，其中、分别是、的单位矢量。代入第一个等号，。你会发现，用矢量的语言不仅说明了大小关系，还再一次验证了合力方向指向O点。
再仔细想想，上述利用矢量关系的推导过程是否绝对依赖于和这个假设呢？你不妨试试把系统的公转中心放在三者的质心上，并重新考虑对另外两个物体的引力，会发现上述关系依然成立，之后总可以求出一个角速度，使得三者一起绕着质心公转。所以对于任意三体系统，这种等边三角形的圆周运动模式都是成立的。
木星轨道的特洛伊小行星群
木星的轨道上存在两群小行星，统称为特洛伊群，其中运行在木星“前方”的为“希腊群”，运行在木星“后方”的为“纯特洛伊群”。这两群小行星所处的位置，正是拉格朗日点和。上图中紫色的部分是希尔达小行星群，它们的三角形轨道与木星形成2:3的轨道共振。希尔达小行星的远日点可能是、、这三个拉格朗日点，当它们连续绕太阳旋转三圈，会相继通过这三个点。
最后，恭喜你看完优美推导！再复习一遍五个拉格朗日点和相应的势能面，因为有意思的部分马上要来了。
再复习一遍五个拉格朗日点 | 图源：Wikipedia

#晕轮轨道

拉格朗日点理论上只是一个没有大小的点，但它对于航天探测有重大意义，因为航天器可以稳定地悬停在它周围的晕轮轨道上。1978年，人们将国际彗星探测器发往日地点，开启了对拉格朗日点的第一次应用。尽管它在2014年9月16日失联，但在这36年1个月3天中，它都恪尽职守地保持在点附近的晕轮轨道上。
你可能会好奇，点虽然在切向上看是势能极小点，但在径向上是势能极大点。也就是说，只要在径向有一点点偏差，就应该会万劫不复地坠入深渊。但上节的示意图考虑的是引力势+随日地系统同步公转的等效离心势，实际的运动轨迹并不是严格地一同公转。这就为晕轮轨道的存在提供了空间。比如，1995年发射的太阳和太阳圈探测器(SOHO)，它的轨道俯视图就长这样：
SOHO的晕轮轨道
日地系统的点能在相对固定的位置避开部分太阳辐射，有利于天文探测，所以詹姆斯·韦布望远镜就选择在这点附近的晕轮轨道。
再想象一下，如果我们从日地系统切换到地月系统，那么地月系统的点就会是对月背而言极佳的观测点。这里需要补充一个背景知识，月亮因为潮汐锁定，它的公转周期和自转周期已经同步了，也就是说，它只有一个面固定朝向地球。
潮汐锁定 | 图源：星球研究所
鹊桥号所在的晕轮轨道 | 图源：星球研究所
例如鹊桥号通信中继卫星就放置在地月系统的拉格朗日点附近，为嫦娥四号在月背软着陆提供稳定的地月通讯信号转播服务。具体地说，鹊桥号在绕点的晕轮轨道上，随地月系统一起做周期运动。这样不仅能实现在月背上空的悬停，还能最大程度地节省轨道调整的燃料。
鹊桥号
到这里，我们讲了三体问题作为混沌体系在数学上没有解析通解，但存在一些稳定的特殊解，着重介绍了其中的五个拉格朗日点及其在航天上的应用。希望三体人的智子在扫描这篇文章时，能收获一些延长恒纪元的灵感，愿三体文明和地球文明和谐共存。
+
+
参考资料：
[1] 既然三体运动非常不稳定，那太阳系为什么能稳定存在？- 知乎@UFO
[2] 拉格朗日点（平动点）的近似推导 - 知乎@Dr.Space
[3] 三体运动模拟 含行星修正版 - bilibili@Zzzyt205
[4] 上海交大科学家领衔提出求解三体问题周期解之路线图
[5] Alpha Centauri. What Does the Closest Stellar System Conceal? - youtube@Kosmo
[6] 为什么三体不稳定，我们的太阳系如此稳定？
[7] Three-Body Gallery
[8] Movies of the Periodic Planar Collisionless Three-Body Orbits in Real Space or on Shape Sphere - 廖世俊教授
>End
>>>                        
本文转载自“中科院物理所”，原标题《三体开播前，帮三体人找找稳定解》。
为分享前沿资讯及有价值的观点，太空与网络微信公众号转载此文，并经过编辑。
支持保护知识产权，转载请注明原出处及作者。
部分图片难以找到原始出处，故文中未加以标注，如若侵犯了您的权益，请第一时间联系我们。
HISTORY/往期推荐
王兆耀：国家战略性力量——商业航天定位的再思考（一）
王兆耀：国家战略性力量——商业航天定位的再思考（二）
王兆耀：国家战略性力量——商业航天定位的再思考（三）
杨长风：航天商业化发展和话语体系转变的思考
>>>                  
充满激情的新时代，
充满挑战的新疆域，
与踔厉奋发的引领者，
卓尔不群的企业家，
一起开拓，
一起体验，
一起感悟，
共同打造更真品质，
共同实现更高价值，
共同见证商业航天更大的跨越！
——《卫星与网络》，观察，记录，传播，引领。
>>>                                           
· 《卫星与网络》特别顾问：王兆耀
· 《卫星与网络》编辑委员会
主任：陈玉忠
高级顾问：王国玉、刘程、童旭东、相振华、王志义、杨烈
编辑委员：曹志刚、陈嵩辉、邓中亮、荆继武、景贵飞、郎燕、刘进军、刘天雄、宁振波、秦智、汪春霆、吴季、徐小舒、阎丽娟、赵敏、肇启明、周建华、朱铎先
· 《卫星与网络》创始人：刘雨菲
· 《卫星与网络》副社长：袁鸿翼
·  微信公众号（ID：satnetdy）团队
编辑：艳玲、哈玫
主笔记者：李刚、魏兴、张雪松、霍剑、乐瑜、刀子、赵栋
策划部：杨艳
视觉总监：董宁
专业摄影：冯小京、宋伟
设计部：顾锰、潘希峎、杨小明
行政部：姜河、林紫
业务部：王锦熙、瑾怡
原创文章转载授权、转载文章侵权、投稿等事宜，请加微信：18600881613
商务合作；展览展厅设计、企业VI/CI及室内设计、企业文化建设及品牌推广；企业口碑传播及整体营销传播等，请加微信：13811260603
杂志订阅，请加微信：
wangxiaoyu9960
· 卫星与网络各分部：
成都分部负责人：沈淮
长沙分部负责人：宾鸿浦
西安分部负责人：郭朝晖
青岛分部负责人：江伟
· 卫星与网络总部负责人：农燕
· 会议活动部负责人：乔颢益
· 投融资及战略层面合作：刘雨菲
· 本平台签约设计公司：一画开天（北京）文化创意设计有限公司
· 航天加（深圳）股权投资基金管理负责人：杨艳
继续阅读
阅读原文
关键词
太阳
就是
颗恒星
地球
引力