三体的节拍器

广义开普勒第三定律猜想（Sun's conjecture）现状综述 · 2026 年 9 月

ACT I · BACKGROUND从开普勒到混沌：三百年只找到三族解

1609 年，开普勒发表行星运动定律；1687 年，牛顿用万有引力定律证明它们可以从二体问题中严格推出。 对二体束缚轨道，第三定律有一个用能量书写的等价形式：

但加一颗星，天就变了。三体问题的 18 个变量无法像二体那样约化，庞加莱 1889 年证明了它的混沌本质。 此后三百年间，人类只找到过 3 族周期解（Euler、Lagrange 的共线/等边解与 Broucke–Hénon 族）。 1993 年 Moore 用最小作用量原理数值发现了著名的「8 字轨道」， 2000 年 Chenciner 与 Montgomery 在 Annals of Mathematics 上给出了它的严格存在性证明 ——这是数学家在周期解方向上的里程碑，但证明的是「存在」，而非任何关于周期与能量的定律。

真正的转折发生在 2013 年：Šuvakov 与 Dmitrašinović 用数值方法一举找到等质量平面三体问题的 13 条新周期轨道（11 个新族）。2017 年起，上海交通大学廖世俊团队用「干净数值模拟」 （Clean Numerical Simulation, CNS）配合超级计算机，把周期解的发现推进到了工业化规模。 把三百年的积累与最近十年的井喷放在同一条时间轴上，才能看出这场爆炸有多剧烈：

三百年攒下 3 族，十年井喷到两千余族。样本量从个位数变成千位数， 统计层面的规律——开普勒不变量近似为常数——才第一次有可能被看见。 正是这批海量周期解，让一条隐藏的规律浮出了水面。

ACT II · THE PATTERN数值中浮现的「准开普勒第三定律」

2015 年，Dmitrašinović 与 Šuvakov 首先注意到：等质量、零角动量三体周期轨道的标度不变周期 \(T|E|^{3/2}\) 近似是常数，但带有拓扑依赖——它与轨道绕三个两体碰撞点环绕方式的 「自由群元素」长度 \(L_f\) 成正比。也就是说，纯组合 \(T|E|^{3/2}\) 并非常数， 要除以拓扑因子 \(L_f\) 之后才趋于普适。

2017 年，李晓明与廖世俊在 695 族周期解上验证了这个规律：定义平均周期 \(\bar{T}=T/L_f\)，则标度不变平均周期

他们称之为「准开普勒第三定律」（quasi Kepler's third law）。 2018 年，对 \(m_1=m_2=1\)、\(m_3\) 变化的 1349 族轨道，他们又给出线性经验律 \(T|E|^{3/2}\approx 3.074\,m_3-0.617\)。2019 年，连初始条件为「三体从静止自由释放」的 313 族任意质量比周期解也近似服从同样的结构。廖世俊在 National Science Review 的观点文章中总结：两千多个周期解全部近似满足这条广义定律——三体问题 混沌的表象之下，存在某种至今无法解释的秩序。

ACT III · THE CONJECTURE孙博华猜想：量纲分析给出的显式公式

孙博华，1989 年兰州大学博士，曾任南非开普半岛理工大学机械工程系长聘教授， 2010 年当选南非科学院院士，2018 年回国后任西安建筑科技大学力学技术研究院创院院长、 中科院北京纳米能源所资深研究员。他的研究风格鲜明——用量纲分析寻找普适标度律 （湍流、毛细上升、多米诺传播皆属此列），还写过专著《量纲分析与 Lie 群》。

2018 年，他在 SCIENCE CHINA Physics, Mechanics & Astronomy 发表 《Kepler's third law of n-body periodic orbits in a Newtonian gravitation field》， 思路分三步：

第一步：Buckingham π 定理

对 N 体周期轨道，选取基本参数：引力场参数 \(\alpha_n\sim Gm_im_j\)、约化质量 \(\mu_n\)、 周期轨道围成的面积 \(A_n\)。量纲分析给出两条独立关系： \(|E_n|A_n^{1/2}/\alpha_n=\text{const}\) 与 \(T_n A_n^{-3/4}\alpha_n^{1/2}\mu_n^{-1/2}=\text{const}\)。 消去面积 \(A_n\)（这正是开普勒第三定律「3/4 次方面积律」的多体影子），得 \(T_n|E_n|^{3/2}=\text{const}\times\alpha_n\sqrt{\mu_n}\)。

第二步：二体匹配定常数

要求 N=2 时退化为经典开普勒第三定律，定出 \(\text{const}=\pi/\sqrt{2}\)。

第三步：质量乘积对称性假设

这是整个推导中唯一不从方程出发的一步。牛顿引力中，力正比于所有两两质量乘积 \(m_im_j\) 的线性组合；孙博华据此假设 N 体推广中的质量因子应取所有 \(\binom{n}{2}\) 个 质量乘积的立方和。最终猜想：

三体情形的非平凡预言为 \(T_3|E_3|^{3/2}=\frac{\pi}{\sqrt{2}}G\big[\frac{(m_1m_2)^3+(m_1m_3)^3+(m_2m_3)^3}{m_1+m_2+m_3}\big]^{1/2}\)。 当 \(m_3\to 0\) 时自动退化回二体定律；等质量时质量因子为常数，回到廖世俊组观察到的经验现象。

与数值的对照

等质量三体：公式预言 \(\pi/\sqrt{2}=2.2214\)，而 695 族数值拟合为 \(2.433\pm0.075\) ——偏差约 8.7%（孙博华原文自报 5.8%）。不等质量情形：公式给出关于 \(m_3\) 的非线性曲线， 与李-景-廖的线性经验律 \(3.074m_3-0.617\) 趋势一致但不重合。此外，一个重质量 \(M\) 加 \(n-1\) 个单位质量的极限下，公式给出 \(T_n|E_n|^{3/2}\approx\sqrt{(n-1)/2}\;\pi GM\)， 是关于 \(M\) 的线性律。

注意一个概念上的错位：数值上的「普适常数」2.433 是在除以拓扑因子 \(L_f\) 之后才出现的； 孙博华的公式里没有 \(L_f\)。两者的对照因此只能是近似层面的——这正是「准」字的含义， 也是后来量子版本与经典版本对不上账的伏笔。

ACT IV · QUANTUM量子多体版本：另一条更微妙的猜想

量子系统没有轨道，周期无从谈起，但这条定律仍然被延伸到了量子世界，靠的是两次构造：

2018 年，Semay（比利时蒙斯大学）在 《Quantum support to BoHua Sun's conjecture》中首次指出： 对 \(N\) 个全同自引力粒子，若给「周期」一个合适的量子定义，则存在类似的「开普勒不变量」， 为孙博华猜想提供了量子支持。随后他与 Willemyns 在 2019 年给出了定态「等效周期」的正式定义 \(T_q=\pi I/\langle K\rangle\)（作用量比平均动能）。

2019 年，孙博华在 《Classical and Quantum Kepler's Third Law of N-Body System》 中顺势提出量子 N 体猜想——注意质量因子的形式变了，从「立方和」变成「和的立方」：

2021 年，Semay 与 Willemyns 在 《Quasi Kepler's third law for quantum many-body systems》 中用「包络理论」（envelope theory，一种给多体能谱提供解析上限的近似方法）为一般质量体系 推出了与孙博华完全相同的公式，并用自引力玻色子基态的高精度数值能量做了检验，符合良好。 但论文同时揭示了一个深层裂缝：

对 \(N\ge 3\)，量子开普勒不变量并不等于经典值——全同粒子情形下恰好是经典值的 \(N(N-1)/2\) 倍。作者的解释：经典不变量的数值检验依赖拓扑因子 \(L_f\)， 而 \(L_f\) 在量子力学中没有对应物。两条「定律」连彼此都对不上。 Semay & Willemyns, Eur. Phys. J. Plus 136, 342 (2021)

此外，包络理论给出的只是「上限的上限」，量子周期本身是半经典构造。所以量子版本的证据强度， 比经典数值版本还要弱一档。

ACT V · STATUS现状：证据、争议与为什么还没有证明

三个版本的法律地位

为什么至今没有证明

其一，对象本身不可解析。三体问题混沌、不可积；连周期解的完整分类都远未完成 （2022 年廖世俊组的路线图仍需「超算 + 机器学习」）。要对整个周期解空间断言一条普适定律， 先得刻画这个空间——这是前提性的缺失。

其二，数值上它就不是严格常数。695 族轨道的 \(\bar{T}^{\,*}=2.433\pm0.075\) 自带约 3% 散布；Dmitrašinović 等人 2018 年的专门数值检验表明，开普勒不变量对轨道拓扑 存在残余依赖。一条「近似成立」的规律，可能根本没有精确成立的版本可供证明—— 这决定了它的归宿更可能是「标度律」而非「定理」。

其三，推导本身不唯一。量纲分析只能定出函数形式到「质量的对称函数」为止； 「立方和」与「和的立方」在 N=2 时完全等价，都满足所有约束。孙博华经典版选了前者， 量子包络理论倾向后者——这正说明猜想的形式是构造性选择，而非方程的必然推论。 赵长印与张明江 2018 年的评论文章 （arXiv:1811.00735）从理论与应用两方面 梳理了该猜想，同样以「conjecture」定性。

其四，数学界的火力在别处。天体力学中严格的现代成果集中在周期解的存在性 （Chenciner–Montgomery 2000 之于 8 字轨道）、变分构造、稳定性与拓扑分类上； 尚无人把这条周期-能量关系作为定理目标。一个可对照的事实是：证明一条 8 字轨道的 「存在」已值得一篇 Annals of Mathematics——为无穷多族轨道证明一条统一标度律， 难度只会更高。

它有什么用

尽管没有证明，这条猜想并非空转。它把周期解搜索的未知量从两个（轨道形状与周期）压缩到 一个：给定质量与能量，公式直接预言周期量级，可作为数值搜索的初值筛选器—— 孙博华原文称之为寻找周期解的「捷径」（shortcut）。对航天任务设计中的周期轨道族普查， 这类先验估计有实际价值。量子版本则揭示了一件理论上耐人寻味的事： 反平方力在多体量子体系中似乎也保留了某种经典不变量的影子——Semay 与 Willemyns 的原话是， 「something special seems to happen with an inverse-square law force」。

开放问题

若要推动它从猜想走向定理（或明确的近似界），至少有三条路都还没有人走通： (1) 从 N 体哈密顿量的标度对称性（\(r\to\lambda r,\ t\to\lambda^{3/2}t\)）出发， 解释为什么拓扑修正后的不变量近似普适——同伦群平均作用量或许是桥梁； (2) 给出不变量对 \(L_f\) 残余依赖的严格上界，把「准」定量化； (3) 在量子侧，寻找 \(L_f\) 的真正对应物（或证明它不存在）， 从而解释 \(N(N-1)/2\) 因子，并判断经典与量子两条猜想能否统一。

ACT VI · THE TWIST · 2026-09-09 补记反转：精确等式被证伪

本文上线当天，我们收到一份由 GPT 给出的严格推导；本节已逐步独立复核 （区间运算全部用有理数精确核验，另做了独立能量交叉验证），确认其成立。 上一节问「为什么没有证明」，答案是：因为它不成立。 精确等式（Sun 2018 猜想）在两种解读下分别被严格反驳； 但同一推导也给出了这个方向上第一个可以完整证明的定理。

反例一：Lagrange 等边三角形旋转解（纯解析）

三个任意正质量置于边长 \(\ell\) 的等边三角形顶点、绕质心匀速旋转。 由质心条件 \(\sum_i m_i\mathbf q_i=0\)，牛顿方程化为 \(\ddot{\mathbf q}_i=-(GM/\ell^3)\,\mathbf q_i\)，取 \(\omega^2=GM/\ell^3\) 即得精确无碰撞周期解。 其能量 \(E=-GS/(2\ell)\)（记 \(S=\sum_{i<j}m_im_j\)），消去 \(\ell\) 得

精确解出现的是「和的立方」\(S^3\)，而猜想是「立方的和」； 正质量下 \(S^3>\sum_{i<j}(m_im_j)^3\)，二者严格矛盾。 孙原文只检验了 \(m_3\to0\) 与「一重多轻」两个极限——恰好都是两个质量分支无法区分的退化情形， 漏掉了三体问题里最经典的精确解。

反例二：8 字轨道 + 计算机辅助证明（拓扑归一化也救不回来）

等质量 8 字轨道是无碰撞、零角动量周期解，拓扑字 \(abAB\)，故 \(L_f=4\)。 Kapela–Simó 的严格区间证书（\(G=m=1\)、\(T=2\pi\) 归一化）给出真实解初值的包围盒； 即使把区间大幅放宽，仍可推出 \(K<1.374\)、\(r_{12},r_{13}<1.651\)、\(r_{23}<0.688\)， 从而 \(|E|>1.29\)。于是

把 \(T\) 解读成 \(T/L_f\)，猜想依然不成立——且这个反例恰好落在经验拟合所在的零角动量类内。 整串比较只涉及有理数四则与开方，可机器核验；区间中点给出的能量 \(-1.29298\) 与 8 字轨道已知值 \(-1.29297\) 四位有效数字一致，佐证了区间证书的有效性。

漏洞定位：量纲分析丢掉了轨道形状

孙推导的关键步骤 \(|E|\sqrt A/\alpha=\mathrm{const}\) 连二体都不成立： 由 \(A=\pi a^2\sqrt{1-e^2}\)、\(|E|=\alpha/(2a)\)，它等于 \(\frac{\sqrt\pi}{2}(1-e^2)^{1/4}\)， 依赖偏心率。量纲分析只允许 \(T|E|^{3/2}=GM^{5/2}\,\Phi(\text{质量比};\text{无量纲轨道数据})\)，排除不了形状与拓扑依赖； 二体的偏心率恰好消去是特殊性质，不能推广。更具讽刺意味的是：取正则的广义约化质量 \(\mu=S/M\)（\(n=2\) 时回到真正的约化质量），孙自己的公式骨架 \(\alpha\sqrt\mu\) 给出的恰是「和的立方」分支——他在量纲分析的两条合法分支里选错了那一条。

可证明的修正版：周期–能量不等式

定理（附条件）。设无碰撞 \(T\)-周期解在质心系中满足 每个质点的时间平均位置为零：\(\frac1T\int_0^T\mathbf q_i\,dt=0\) （该条件涵盖全部追逐型轨道——含 8 字轨道——与 Lagrange 旋转解）。则

证明骨架：Lagrange–Jacobi 恒等式 \(I''=4K-2U\) 周期积分给出位力关系 \(2\langle K\rangle=\langle U\rangle\)，故作用量 \(\mathcal A=\int_0^T(K+U)\,dt=3|E|T\)； 质心系动能按质点对分解后，对每对（零均值、周期的）相对矢量联用 Wirtinger 不等式与 \(y^{-1/2}\) 的 Jensen 不等式，再对唯一自由参数取极小， 每对贡献 \(\ge\frac32(GM)^{2/3}(2\pi/T)^{2/3}\)，装配即毕。

限定与余波

条件充分而非必要。二体偏心轨道的时间平均位置为 \(-\frac32ae\,\hat e\)（指向远拱点），违反定理条件，但不等式对一切偏心率取等号。 放松适用条件是明确可做的问题。

被证伪的是精确等式，经验规律还活着。 \(\bar T|E|^{3/2}\approx2.433\) 的近似恒常性（百分之几散布）依然成立、依然无解释—— 它现在成为唯一真正待证的目标：由 \(\mathcal A=3|E|T\) 有 \(\bar T|E|^{3/2}=\mathcal A\sqrt{|E|}/(3L_f)\)，即把「每单位拓扑字长的作用量近似恒定」 变成带误差界的命题。

一个耐人寻味的巧合。可证下界的常数恰好等于孙博华 2019 年的 量子版公式；等质量时「下界 ÷ 经典猜想值」精确等于 \(n(n-1)/2\) ——正是 Semay–Willemyns 发现的量子-经典偏差因子。这个因子由此获得纯经典解释， 也再次指向「和的立方」才是物理上正确的分支。

反思：如此明显的错误，为什么没有被数值证伪

反例一只需小学代数，反例二的轨道是全领域最出名的数值解——为什么八年里无人察觉？ 这不是数值不够多，而是错误被五层结构恰好藏住了：

其一，被测样本里没有反例的位置。两千多个「验证」轨道 （695、1349、313 族）全部是零总角动量情形——CNS 搜索就是按此设计的； 而决定性反例 Lagrange 旋转解是给定能量下角动量最大的轨道，位于轨道空间的另一头。 因子三的分歧只活在旋转轨道上；在零角动量类内，8 字轨道只差 4%，695 族拟合只差 9.5%。 反例不是被忽略了，是它从来不在数据集里。

其二，两个结构因素把错误压缩成了「小散布」。 等质量把错误的质量因子退化成一个纯数字——选错分支只是常数差一个 O(1) 因子； 拓扑归一化又把量程压进窄带：实测 \(2.433\)、孙预言 \(2.2214\)、可证下界除以 \(L_f\) 后 \(1.666\)，全是 O(2) 的数。常数错 10% 在物理学家眼里是「近似吻合」，不是「结构错误」。 数值拟合的对象 \(\bar T\) 与公式里的 \(T\) 本不是同一个量，这个错位反而成了错误公式的藏身所。

其三，公式把量纲分析真正决定的部分全做对了。 \(T|E|^{3/2}=GM^{5/2}\Phi(\text{质量比})\) 的骨架——标度律、对 \(G\)、总质量与能量的依赖—— 严格正确，任何分支都给得出来。数值对照确认的是「标度对、量级对」； 唯一出错的是量纲分析无权决定、只能由动力学裁决的那个 O(1) 因子。

其四，质量依赖的检验区间太窄。1349 族的线性经验律 \(3.074m_3-0.617\) 只在 \(m_3\approx1\) 附近拟合（外推到 \(m_3\to0\) 甚至给负值）。 孙的非线性曲线在这个窄窗里也近似线性，差约 10%——「趋势一致但不重合」， 再次被读成吻合。要判分支错误，需要跨越量级质量比的数据，而这样的数据不存在。

其五，最便宜的检验没人做，而线索早已公开。 Lagrange 解的周期–能量关系人人「知道」，恰恰因此没人花五分钟去算它的质量因子是 \(S^3/M\) 还是 \(\sum(m_im_j)^3/M\)。孙自己 2019 年的量子公式用的就是「和的立方」， Semay–Willemyns 也发现量子与经典常数差 \(N(N-1)/2\)——他们把差异归因于 \(L_f\) 没有量子对应物，而没有怀疑经典公式选错了分支。 证伪所需的全部信息早已公开，缺的是把两件事实放在一起的那个念头。

一句话心得：被检验的命题版本（零角动量、等质量、除以 \(L_f\)、允许 10% 散布） 恰好是错误公式的最佳藏身所；而致命检验不需要超算——只需要一个 250 岁的精确解， 和对质量分支的一点怀疑。证伪从来不缺数据，缺的是视角。

REFERENCES主要文献

1. B. H. Sun, Kepler's third law of n-body periodic orbits in a Newtonian gravitation field, Sci. China-Phys. Mech. Astron. 61, 054721 (2018). arXiv:1807.10685
2. B. H. Sun, Classical and Quantum Kepler's Third Law of N-Body System, Results in Physics 13, 102144 (2019). arXiv:1904.05148
3. X. M. Li and S. J. Liao, More than six hundred new families of Newtonian periodic planar collisionless three-body orbits, Sci. China-Phys. Mech. Astron. 60, 129511 (2017). arXiv:1705.00527
4. X. Li, Y. P. Jing and S. J. Liao, Over a thousand new periodic orbits of a planar three-body system with unequal masses, Publ. Astron. Soc. Japan 70, 64 (2018). arXiv:1709.04775
5. X. Li and S. J. Liao, Collisionless periodic orbits in the free-fall three-body problem, New Astronomy 70, 22 (2019). arXiv:1805.07980
6. S. Liao, X. Li and Y. Yang, Three-body problem — from Newton to supercomputer plus machine learning, New Astronomy 96, 101850 (2022). arXiv:2106.11010
7. M. Šuvakov and V. Dmitrašinović, Three classes of Newtonian three-body planar periodic orbits, Phys. Rev. Lett. 110, 114301 (2013).
8. V. Dmitrašinović and M. Šuvakov, Topological dependence of Kepler's third law for collisionless periodic three-body orbits with vanishing angular momentum and equal masses, Phys. Lett. A 379, 1939–1945 (2015).
9. V. Dmitrašinović, A. Hudomal, M. Shibayama and A. Sugita, Linear stability of periodic three-body orbits with zero angular momentum and topological dependence of Kepler's third law: a numerical test. arXiv:1705.03728
10. C. Semay, Quantum support to BoHua Sun's conjecture, Results in Physics 13, 102167 (2019). arXiv:1812.00807
11. C. Semay and C. T. Willemyns, Quasi Kepler's third law for quantum many-body systems, Eur. Phys. J. Plus 136, 342 (2021). arXiv:2005.06973
12. C.-Y. Zhao and M.-J. Zhang, A conjecture on Kepler's third law of n-body periodic orbits (2018). arXiv:1811.00735
13. A. Chenciner and R. Montgomery, A remarkable periodic solution of the three-body problem in the case of equal masses, Annals of Mathematics 152, 881–901 (2000).
14. S. J. Liao and X. M. Li, On the periodic solutions of the three-body problem, National Science Review 6, 1070 (2019). PMC8291409
15. T. Kapela and C. Simó, Computer assisted proofs for nonsymmetric planar choreographies and for stability of the Eight, Nonlinearity 20, 1241–1255 (2007).（ACT VI 反例二的严格区间证书来源）

本文基于公开文献调研撰写，公式与数值均核对自原始论文。 孙博华 2018 年论文被引约 18 次，量子扩展方向被引均在个位数——这个小而活跃的领域， 等待的可能正是一位愿意把它当成定理来证的数学家。