三体问题

三体问题（英语：Three-body problem）是天体力学中的基本力学模型。它是指三个质量、初始位置和初始速度都是任意的可视为质点的天体，在相互之间万有引力的作用下的运动规律问题。

初始位置在斜三角形顶点，且初始速度均为零的三个相同物体的近似轨迹。按照动量守恒定律，质心仍然存在。

内容

三体问题是多体问题的一个特例。对于一般多体问题而言，不存在一般的解析解，即难以预测所有三体问题的数学情景，但存在较为复杂的广义解析解：松德曼-汪秋栋级数解。现在更多是研究其特殊情况或使用数值方法求解。

例如太阳系中，考虑太阳、地球和月球的运动，它们彼此以万有引力相吸引，若假设三个星球都可设为质点，并且忽略其他星球的引力，太阳、地球和月球的运动即可以视为三体问题。

三体问题也被用于模拟经典力学或量子力学中三个粒子的运动状态。

数学概述

三体问题可以用三个质量为$m_{i}$的相互作用的物体的矢量位置$\mathbf {r_{i}} =(x_{i},y_{i},z_{i})$的牛顿运动方程数学表示： $${\begin{cases}{\ddot {\mathbf {r} }}_{\mathbf {1} }&=-Gm_{2}{\frac {\mathbf {r_{1}} -\mathbf {r_{2}} }{|\mathbf {r_{1}} -\mathbf {r_{2}} |^{3}}}-Gm_{3}{\frac {\mathbf {r_{1}} -\mathbf {r_{3}} }{|\mathbf {r_{1}} -\mathbf {r_{3}} |^{3}}},\\{\ddot {\mathbf {r} }}_{\mathbf {2} }&=-Gm_{3}{\frac {\mathbf {r_{2}} -\mathbf {r_{3}} }{|\mathbf {r_{2}} -\mathbf {r_{3}} |^{3}}}-Gm_{1}{\frac {\mathbf {r_{2}} -\mathbf {r_{1}} }{|\mathbf {r_{2}} -\mathbf {r_{1}} |^{3}}},\\{\ddot {\mathbf {r} }}_{\mathbf {3} }&=-Gm_{1}{\frac {\mathbf {r_{3}} -\mathbf {r_{1}} }{|\mathbf {r_{3}} -\mathbf {r_{1}} |^{3}}}-Gm_{2}{\frac {\mathbf {r_{3}} -\mathbf {r_{2}} }{|\mathbf {r_{3}} -\mathbf {r_{2}} |^{3}}}.\end{cases}}$$ 其中$G$为万有引力常数。这是一组9个二阶微分方程构成的方程组。

这个问题也可以用哈密顿形式等价表示，此时可以用一组18个一阶微分方程表示，这些方程分别对应于$\mathbf {r_{i}}$位置和动量$\mathbf {p_{i}}$的一个分量： $${\frac {d\mathbf {r_{i}} }{dt}}={\frac {\partial {\mathcal {H}}}{\partial \mathbf {p_{i}} }},\qquad {\frac {d\mathbf {p_{i}} }{dt}}=-{\frac {\partial {\mathcal {H}}}{\partial \mathbf {r_{i}} }},$$ 此处${\mathcal {H}}$为哈密顿量： $${\mathcal {H}}=-{\frac {Gm_{1}m_{2}}{|\mathbf {r_{1}} -\mathbf {r_{2}} |}}-{\frac {Gm_{2}m_{3}}{|\mathbf {r_{3}} -\mathbf {r_{2}} |}}-{\frac {Gm_{3}m_{1}}{|\mathbf {r_{3}} -\mathbf {r_{1}} |}}+{\frac {\mathbf {p_{1}} ^{2}}{2m_{1}}}+{\frac {\mathbf {p_{2}} ^{2}}{2m_{2}}}+{\frac {\mathbf {p_{3}} ^{2}}{2m_{3}}}.$$

限制性三体问题

圆型限制性三体问题是太阳系中椭圆轨道的有效近似，这里可以用两个主天体的引力以及它们旋转的离心效应而产生的电势组合实现可视化（科里奥利效应是动态的，此处未显示）。然后，可以将拉格朗日点视为合成表面上梯度为零的五个位置（用蓝线显示），表明此处处于平衡状态。

当所讨论的三个天体中﹐三体中其中两体的质量极大，以至于第三体的质量几乎不能对其造成任何扰动；或是有一个天体的质量与其他两体的质量相比小到可以忽略时，可将三体问题简化为二体问题的变型。这样的三体问题称为限制性三体问题。一般地把小质量的天体称为无限小质量体﹐简称小天体；把另外两个大质量的天体称为有限质量体。

设$m_{1,2}$为两个有限质量体的质量，其平面坐标分别为$(x_{1},y_{1})$和$(x_{2},y_{2})$，另设小天体坐标为$(x,y)$。所选单位应该要确保两有限质量体之间距离和万有引力常数均为1，这样数学处理会比较简单。则小天体的运动可以以下公式数学表示： $${\begin{cases}{\frac {d^{2}x}{dt^{2}}}=-m_{1}{\frac {x-x_{1}}{r_{1}^{3}}}-m_{2}{\frac {x-x_{2}}{r_{2}^{3}}},\\{\frac {d^{2}y}{dt^{2}}}=-m_{1}{\frac {y-y_{1}}{r_{1}^{3}}}-m_{2}{\frac {y-y_{2}}{r_{2}^{3}}}.\end{cases}}$$ 其中$r_{i}={\sqrt {(x-x_{i})^{2}+(y-y_{i})^{2}}}.$

在此数学表示中，运动方程通过坐标$x_{i}(t),y_{i}(t)$具有明确的时间相关性。但可以通过转换为旋转参考系来消除这一时间相关性，从而简化了后续分析。

历史

1499年，韦斯普奇最早提出三体问题。他利用月球位置讯息来确定自己在巴西的位置。这一方法适用于海上确定经度导航，因此在17世纪20年代，此法变得尤为实用。实际上，约翰·哈里森发明的航海经线仪解决了确定经度的问题。但是，由于太阳和行星对月球绕地公转的扰动效应，月球理论的准确性很低。

亚美利哥·韦斯普奇和伽利略·伽利莱随后也提出了三体问题。伽利略确定所有天体的下落速度都是相同均匀的，但他并没有将其应用到行星运动。

1687年，当时艾萨克·牛顿发表了《自然哲学的数学原理》。在《原理》第一卷的第66号命题及其22个推论中，牛顿首次定义和研究了三个受相互扰动的重力吸引影响的巨大物体的运动问题。在第三册的第25至35条命题中，牛顿将他的66号命题结果用于月球在地球和太阳的引力影响下的运动。

让·勒朗·达朗贝尔和亚历克西·克洛德·克莱罗两个长期竞争者都试图找到该问题的某种通性。他们竞相在1747年向法国皇家科学研究院提交了他们的第一批三体问题分析。“三体问题”（法语：Problème des trois Corps）这个名字至此开始广泛使用。达朗贝尔于1761年发布的文章表明，该名称最早使用于1747年。

1887年，为了祝贺自己的60岁寿诞，瑞典国王奥斯卡二世赞助了一项现金奖励的竞赛，征求太阳系的稳定性问题的解答，这是三体问题的一个变型。

法国数学家亨利·庞加莱简化了问题，提出了限制性三体问题：即三体中其中两体的质量极大，以至于第三体的质量完全不能对其造成任何扰动。面对这个问题，庞加莱运用了他发明的相图理论，并且最终发现了混沌理论。虽然庞加莱没有成功给出一个完整的解答，不过他的工作令人印象深刻，因此他于1888年赢得了奖金。

庞加莱发现这个系统的演变经常是混沌的，即如果初始状态有一个小的扰动，例如个体的初始位置有一个小的变动，则后来的状态可能会有极大的不同。如果这一小变动不能被测量仪器所探测，那么我们将不能预测最终状态为何。

著名数学家兼竞赛裁判卡尔·魏尔施特拉斯说：“这个工作不能真正视为对所求的问题的完善解答，但是它的重要性使得它的出版将标志着天体力学的一个新时代的诞生。”

魏尔斯特拉斯并不知道他自己的预测有多准确。在庞加莱的论文中，他描述了例如同宿点（homoclinic points）之类的新思想。这些概念在1890年的瑞典《数学学报》备忘录中出版，后来该书在编辑途中被发现一个理论上的错误，然而该错误实际上导致了庞加莱一些进一步的发现，这些发现在现代视为混沌理论的开端。

19世纪末到20世纪初，科学家们研究使用短程二体吸力来解决三体问题的方法。19世纪后期，乔治·威廉·希尔在用金星和水星的运动来研究限制性三体问题，不久进一步延伸到四体问题以计算月球绕地球轨道和行星环绕恒星轨道。

1907年和1909年时，芬兰数学家卡尔·弗里肖夫·松德曼分别发表两篇论文，证明三体问题存在一个t^1/3的幂级数解，除初始状态是角动量为零的情况（这种情况极其罕见，只有当物体勒贝格测度为0时才会出现）外，对所有实数t始终收敛。松德曼的研究成果收录于1912年的瑞典《数学学报》上 。

1930年，数学家大卫·贝洛里奇（David Beloriszky）指出，松德曼级数收敛极慢，如果将松德曼级数解用于天文观测，则计算有效项数将至少达到10^{7006800000000000000♠8000000}项。

1970年，前苏联理论物理学家维塔利·尼古拉耶维奇·叶菲莫夫预测当三个相同的玻色子之间的对相互作用接近共振时，三体谱表现出无限的束缚态序列。这一预测于2006年证实，因此也称为叶菲莫夫效应。

20世纪70年代，米歇尔·赫农和罗杰A.布鲁克分别找到了一套特殊解族，这些特殊解族构成了同系列特殊解族: 布鲁克-赫农-哈德吉德梅特里奥解族（Broucke–Henon–Hadjidemetriou family）。

1979年，前苏联数学家L·K·巴巴德赞詹兹（LK Babadzanjanz）发表论文《Existence of the continuations in the N-body problem（多体问题连续性的存在性）》， 后又于1993年发表论文《On the global solution of the N-body problem（论多体问题全局解）》。

1985年，中国数学家汪秋栋发表论文《多体问题全局解的存在性》 ， 后又于1991年发表论文《The global solution of the n-body problem（多体问题全局解）》。

三体问题在单周期T≃6.3259时的“8”字型解动画。

1993年，圣塔菲研究所的美国物理学家克里斯·摩尔构想了一种零角动量解，该解适用于三个相等质量物体做“8”字型运动。 此情形在2000年由数学家阿兰·契纳和理查德·蒙哥马利（Richard Montgomery）证明。 该解在数值上证明了对于质量和轨道参数小扰动的稳定性，这增加了在物理宇宙中可以观察到这种轨道的可能。但有人认为不太可能发生这种情况，因为其稳定的范围较小，因此不太可能发生这种情况。此外，双星散射事件会导致“8”字形轨道的出现概率估计不足1%。

2013年，贝尔格莱德物理研究所的物理学家米洛万·乌瓦科夫（Milovan Šuvakov）与维利科·德米特拉·伊诺维（Veljko Dmitrašinović）发现了等质量零角动量三体问题的13种新的解族。

2015年，物理学家安娜·胡多马尔（Ana Hudomal）发现了14种等质量零角动量三体问题的新解族。

2017年，中国计算机科学家李晓明和中国数学家廖世俊发现了669组等质量零角动量三体问题的新周期轨道。 2018年，不等质量的零动量系统又新增了1223组新解。

2018年，李晓明和廖世俊发布了234个不等质量“自由落体”三体问题的解。三体问题的自由落体公式从所有三个静止的物体开始。正因为如此，质量在一个自由落体配置不在一个闭合的“循环”轨道上运行，而是沿着一个开放的“轨道”向前和向后运行。

2019年，爱丁堡大学菲利普·布林（Philip Breen）等人训练了一种用于三体问题的快速神经网络。

2022年，廖世俊、杨宇和李晓明在《新天文学》上提出了一个获得“三体问题”周期解的路线图。

2023年，伊万·赫里斯托夫（Ivan Hristov）、拉多斯拉娃·赫里斯托娃（Radoslava Hristova）、维利科·德米特拉希诺维奇（Veljko Dmitrašinović）和谷川清隆（Kiyotaka Tanikawa）发表了一篇关于等质量三体系统的周期轨道研究，发现了12409个不同轨道。

解法

一般性解法

由三个物体相互作用构成的引力系统是混沌系统，不过由三个物体弹性相互作用构成的系统则不然。

三体序列的收敛半径由到最近奇点的距离决定。因此，必须要考虑三体问题的可能奇点。三体问题中的唯二奇点是二元碰撞（两个天体间的瞬时碰撞）和三元碰撞（三个天体间的瞬时碰撞），下面对此类状态进行简单讨论。

实际上多体间不太可能发生碰撞（包括二体和三体），因为此前已证明碰撞状态对应勒贝格测度为0的一组初始状态。然而，由于没有对这一初始状态进行研究，以避免相应解出现碰撞奇点。因此，松德曼的求解思路包含以下几大要点:

1. 通过适当的变量替换，用正则化去分析二元碰撞以外的解。
2. 用保罗·潘勒韦方程的思路证明只有在角动量L=0时才会发生三元碰撞。将初始状态限制为L ≠ 0，从三体问题的变换方程中去除所有实奇点。
3. 证明如果L≠0，则不仅不存在三元碰撞，而且系统严格有界远离三元碰撞。通过对微分方程使用柯西存在性定理，证明在以实轴为中心的复平面（符合柯西—柯瓦列夫斯卡娅定理的阴影）中，证明存在一带状区域（与L取值有关），其中不包含复奇点。
4. 构造保角变换h(z)=$e^{\frac {\pi \omega }{2\Omega }}$，把这一带状区域映射到单位圆盘。他证明了如果ω = t^1/3（正则化后的新变量），且| ω- ω*| ≤ Ω，则映射可由下面两式得出：

$${\begin{cases}\omega ={\frac {2\Omega }{\pi }}log{1+\tau \over 1-\tau },\\\tau ={\frac {e^{\frac {\pi \omega }{2\Omega }}-1}{e^{\frac {\pi \omega }{2\Omega }}+1}}.\end{cases}}$$ 由此，松德曼证明，l和Ω一旦确定，就可以通过对τ的连续微分计算各项发展，从而证明了三体问题存在收敛级数解。

1985年，中国数学家汪秋栋发表论文《多体问题全局解的存在性》，后又于1991年发表论文《The global solution of the n-body problem（多体问题全局解）》， 将解法推广至多体问题中，并给出了角动量L为0的初始状态下的解法。

特殊解法

传统理论

20组三体问题周期性特解示例。

1767年，莱昂哈德·欧拉提出了三个周期解序列，其中三个物体在任意时刻均共线。具体内容参见欧拉三体问题。

1772年，拉格朗日找到了一系列由三个质点组成的等边三角形解。这些解与欧拉的共线解一起构成了三体问题的中心构形。该系列解对于任何质量比的天体均有效，并且都沿开普勒椭圆运动。这四个解族是唯一有明确解析式的已知解族。在圆型限制性三体问题的特殊情况下，这些解在随原点变换的旋转坐标轴中称为拉格朗日点L₁, L₂, L₃, L₄和L₅，其中L₃, L₄是拉格朗日对称解的实例。

1892年至1899年，庞加莱为限制性三体问题找到了一组无穷多的周期解，并将这组解推广至一般三体问题。

1893年，恩斯特·迈塞尔提出了三体勾股问题：将质量之比为3：4：5的三个质点分别置于3：4：5直角三角形的顶点处，研究其轨迹。卡尔·伯劳在1913年进一步研究了这个问题。

1967年，维克多·塞贝赫利和弗雷德里克·彼得斯（C. Frederick Peters）利用数值积分理论建立了这个问题的最终逃逸模型，同时找到了附近的周期解。

20世纪70年代，米歇尔·赫农和罗杰A.布鲁克分别找到了一套特殊解族，这些特殊解族构成了同系列特殊解族: 布鲁克-赫农-哈德吉德梅特里奥解族。在这一解族中，这三个物体都具有相同的质量，可以表现出逆行和直行两种形式。在布鲁克的一些解中，两个物体延相同的轨道运行。

数值分析方法

尽管高精度需要大量的CPU时间，但是可以用数值积分在计算机上得到多体问题的任意高精度解。现已有人尝试创建计算机程序，结合狭义相对论等现代物理学理论，以数字方式解决涉及电磁相互作用和引力相互作用的三体问题（以及扩展的多体问题）。此外，使用随机游走理论，可以计算三体处于不同位置的概率。

2017年，廖世俊和李晓明发明了“精准数值模拟（the clean numerical simulation，CNS）”的混沌系统数值模拟新策略，用超级计算机成功获得了695个等质量三体系统的周期轨道。

2019年，布林等人训练了一种用于三体问题的快速神经网络，该神经网络使用数字积分器进行训练。

2023年，伊万·赫里斯托夫、拉多斯拉娃·赫里斯托娃、维利科·德米特拉希诺维奇和谷川清隆发表了一篇关于等质量三体系统的周期轨道研究。其中，他们发现了12409个不同轨道。

多体问题

三体问题是多体问题的一个特例，它描述了多个物体如何在一种相互作用（如引力）下移动。这些问题具有收敛幂级数形式的全局解析解，正如松德曼（n = 3）和汪秋栋（n >3）所证明的那样（详见N体问题）。然而，松德曼-汪秋栋级数收敛得如此之慢，以至于它们没有什么实用性。因此，目前还需要用数值积分形式的数值分析方法求取近似解，或在某些特殊情况下，使用经典三角级数近似（参见多体模拟）。原子系统，例如原子、离子和分子，可以用量子多体问题来处理。在经典力学系统中，多体系统通常是指一个星系或星系团；行星系统，如恒星、行星及其卫星等也可视为多体系统。一部分可以用扰动理论简化处理，其中系统被视为二体系统加上一个假定导致偏离无扰动双体轨迹的额外力。

文化领域

中国作家刘慈欣围绕三体问题展开创作了《三体》三部曲中的第一部作品《三体》，向公众科普了三体问题。

另见

• N体问题
• 多体问题
• 混沌理论
• 三体 (小说)
• 多体模拟

延伸阅读

• Aarseth, S. J. Gravitational n-Body Simulations. New York: Cambridge University Press. 2003. ISBN 978-0-521-43272-6. 
• Bagla, J. S. Cosmological N-body simulation: Techniques, scope and status. Current Science. 2005, 88: 1088–1100. Bibcode:2005CSci...88.1088B. arXiv:astro-ph/0411043 . 
• Chambers, J. E.; Wetherill, G. W. Making the Terrestrial Planets: N-Body Integrations of Planetary Embryos in Three Dimensions. Icarus. 1998, 136 (2): 304–327. Bibcode:1998Icar..136..304C. CiteSeerX 10.1.1.64.7797 . doi:10.1006/icar.1998.6007. 
• Efstathiou, G.; Davis, M.; White, S. D. M.; Frenk, C. S. Numerical techniques for large cosmological N-body simulations. Astrophysical Journal. 1985, 57: 241–260. Bibcode:1985ApJS...57..241E. doi:10.1086/191003. 
• Hulkower, Neal D. The Zero Energy Three Body Problem. Indiana University Mathematics Journal. 1978, 27 (3): 409–447. Bibcode:1978IUMJ...27..409H. doi:10.1512/iumj.1978.27.27030 . 
• Hulkower, Neal D. Central Configurations and Hyperbolic-Elliptic Motion in the Three-Body Problem. Celestial Mechanics. 1980, 21 (1): 37–41. Bibcode:1980CeMec..21...37H. S2CID 123404551. doi:10.1007/BF01230244. 
• Li, Xiaoming; Liao, Shijun. On the stability of the three classes of Newtonian three-body planar periodic orbits. Science China Physics, Mechanics & Astronomy. 2014, 57 (11): 2121–2126. Bibcode:2014SCPMA..57.2121L. S2CID 73682020. arXiv:1312.6796 . doi:10.1007/s11433-014-5563-5. 
• Moore, Cristopher. Braids in Classical Dynamics (PDF). Physical Review Letters. 1993, 70 (24): 3675–3679 [2016-01-01]. Bibcode:1993PhRvL..70.3675M. PMID 10053934. doi:10.1103/PhysRevLett.70.3675. （原始内容 (PDF)存档于2018-10-08）. 
• Poincaré, H. New Methods of Celestial Mechanics 3 vol. English translated. American Institute of Physics. 1967. ISBN 978-1-56396-117-5.  含有内容需登入查看的页面 (link)
• Šuvakov, Milovan; Dmitrašinović, V. Three Classes of Newtonian Three-Body Planar Periodic Orbits. Physical Review Letters. 2013, 110 (10): 114301. Bibcode:2013PhRvL.110k4301S. PMID 25166541. S2CID 118554305. arXiv:1303.0181 . doi:10.1103/PhysRevLett.110.114301. 

外链

• Chenciner, Alain. Three body problem. Scholarpedia. 2007, 2 (10): 2111. Bibcode:2007SchpJ...2.2111C. doi:10.4249/scholarpedia.2111 . 
• Physicists Discover a Whopping 13 New Solutions to Three-Body Problem （页面存档备份，存于互联网档案馆） (Science)
• 3body simulator （页面存档备份，存于互联网档案馆） - an example of a computer program that solves the three-body problem numerically