三体的节拍器 II：一个证明
我证明了：孙博华的量子公式是一条严格定理——它是所有 N 体周期轨道的 \(\tau=T|E|^{3/2}\) 的一条下界，且这条下界是尖锐的： 二体开普勒椭圆（任意偏心率）与三体 Lagrange 等边轨道恰好取到等号。 证明只需四件经典工具：位力定理、Lagrange 恒等式、Hölder 不等式、Wirtinger 不等式。
同时我用精确解析反例否定了孙博华的经典公式作为「定律」的资格： Lagrange 等边轨道的 \(\tau\) 恰好是它的 3 倍，Euler 共线轨道恰好是它的 \(5\sqrt2/2\) 倍。它只是一条更弱的下界。
全部结论都由 125 条精确周期轨道（\(N=2\ldots6\)，质量比跨 0.05–4.4，闭合残差中位数 \(1.8\times10^{-12}\)） 逐步复核，33 项证明链检查零失败。
N = 2…6
位力残差 1.2e-14
0 失败
对孙氏经典公式的精确偏离
ACT I · SETUP要证的到底是什么
设 \(N\) 个质点 \(m_1,\dots,m_N\) 在牛顿引力下作周期运动，周期 \(T\)，总能量 \(E<0\)。 牛顿势是 \(-1\) 次齐次的，因此 \(\mathbf q\to\lambda\mathbf q,\ t\to\lambda^{3/2}t\) 把解映为解， 且 \(E\to\lambda^{-1}E\)、\(T\to\lambda^{3/2}T\)。于是
\(\tau\) 的标度不变性是精确的，这才是「三体也服从开普勒第三定律」这个经验现象的真正内核： 它说明 \(\tau\) 只依赖轨道的形状（与质量），不依赖轨道的大小。 问题于是变成一个纯几何问题：\(\tau\) 在周期解空间上有没有下界？下界是什么？何时取到？
孙博华给出了两个候选答案：
两者在 \(N=2\) 时完全相同（都退化为精确的开普勒第三定律），这是它们都能通过二体检验的原因。 \(N\ge3\) 时由 \(\sum p^3\le(\sum p)^3\) 立得 \(\tau_{\rm C}\le\tau_{\rm Q}\)。 下面证明：\(\tau_{\rm Q}\) 是定理，\(\tau_{\rm C}\) 不是。
ACT II · TWO IDENTITIES两条恒等式
证明。换元给出 \(A_K=T^2\langle K\rangle\)、\(A_U=\langle U\rangle\)（\(\langle\cdot\rangle\) 为时间平均）。 对 \(-1\) 次齐次势的周期解，位力定理给 \(\langle K\rangle=|E|\)、\(\langle U\rangle=2|E|\)。 代入：\(\frac12A_U\sqrt{A_K}=\frac12(2|E|)(T\sqrt{|E|})=T|E|^{3/2}=\tau\)。∎
数值复核：33 组轨道上 \(\tau\) 与 \(\frac12A_U\sqrt{A_K}\) 的相对偏差 \(4\times10^{-14}\)–\(2\times10^{-10}\)（偏差随偏心率升高而增大，纯属积分精度）。
数值复核：最大相对偏差 \(3\times10^{-16}\)–\(8\times10^{-16}\)，即机器精度。
ACT III · THE THEOREM主定理与证明
- S1（重述 T1）。\(\tau=\frac12A_U\sqrt{A_K}\)，且 \(A_K=T^2\langle K\rangle\)、\(A_U=\langle U\rangle\)。
- S2（逐点 Hölder）。令 \(p_{ij}=m_im_j\)、\(d_{ij}=|\mathbf q_{ij}(s)|\)。对每个固定的 \(s\)，用指数 \((3/2,\,3)\) 的 Hölder 不等式：
\[
P=\sum_k p_k=\sum_k\Big(\frac{p_k}{d_k}\Big)^{2/3}\big(p_kd_k^2\big)^{1/3}
\le\Big(\sum_k\frac{p_k}{d_k}\Big)^{2/3}\Big(\sum_kp_kd_k^2\Big)^{1/3}.
\]
两边立方，再用 T2 把 \(\sum_kp_kd_k^2\) 换成 \(M I(s)\)：
\[
\sum_{i
- S3（Jensen）。\(x\mapsto x^{-1/2}\) 在 \((0,\infty)\) 上凸，故对 S2 积分得
\[
\frac{A_U}{G}=\int_0^1\sum_{i
- S4（Wirtinger / Poincaré）。\(\mathbf q_i-\bar{\mathbf q}_i\) 是零均值的 1-周期函数，故 \[ A_K=\frac12\int_0^1\sum_i m_i|\mathbf q_i'|^2ds =\frac12\int_0^1\sum_i m_i|(\mathbf q_i-\bar{\mathbf q}_i)'|^2ds \ \ge\ \frac12\cdot4\pi^2 V=2\pi^2V . \] 等号成立当且仅当每个 \(\mathbf q_i\) 只含单一 Fourier 模（匀速圆周）。
- S5（合并）。 \[ A_U\sqrt{A_K}\ \ge\ \frac{GP^{3/2}}{\sqrt{M\bar I}}\cdot\pi\sqrt2\sqrt V =\pi\sqrt2\,G\,\frac{P^{3/2}}{\sqrt M}\sqrt{\frac{V}{\bar I}}, \] 除以 2 并用 S1 即得主定理。∎
尖锐性：等号在哪里取到
这条下界不是松的。\(N=2\) 时它被每一条开普勒椭圆精确取到，与偏心率无关。 这看起来矛盾——偏心率 \(e>0\) 时 S3 的 Jensen 步与 S4 的 Wirtinger 步都是严格不等式—— 但两个亏格恰好互相抵消：
| \(e\) | Jensen 亏格 | Wirtinger 亏格 | \(V/\bar I\) | \(\tau/\tau_{\rm Q}\)（实测） |
|---|---|---|---|---|
| 0.0 | 1.000000000 | 1.000000000 | 1.000000000 | 1.000000000000 |
| 0.2 | 1.029563014 | 1.030927835 | 0.915094340 | 1.000000000000 |
| 0.5 | 1.172603940 | 1.230769231 | 0.590909091 | 1.000000000000 |
| 0.8 | 1.399999999 | 1.923076923 | 0.265306122 | 1.000000000000 |
| 0.95 | 1.534193600 | 3.094777563 | 0.137280935 | 1.000000000000 |
解析地：二体时 Jensen 亏格为 \(\sqrt{1+\tfrac32e^2}\)，Wirtinger 亏格为 \((1-\tfrac34e^2)^{-1}\)， 而 \(V/\bar I=(1-\tfrac34e^2)/(1+\tfrac32e^2)\)；三者相乘恒等于 1。
\(N=3\) 时等号还被 Lagrange 等边轨道取到（见 ACT IV），而且对任意质量都成立。 \(N\ge4\) 时等号取不到——原因很干净：平面内不存在四点两两等距的位形。
ACT IV · EXACT VALUEShomographic 轨道的精确闭式，以及孙氏量子公式的真身
证明。中心位形满足 \(\nabla_iU=\lambda m_i\mathbf c_i\)，\(\lambda=-U/I\)。 圆轨道角速度 \(\omega^2=U/I\)，故 \(K=U/2\)、\(E=-U/2\)、\(T=2\pi/\omega\)，代入 \(\tau=T|E|^{3/2}\) 得 \(\frac{\pi}{\sqrt2}\sqrt I\,U\)。椭圆 homographic 轨道是同一形状的 Kepler 标度族，而 \(\tau\) 标度不变，故值不变。∎
数值复核：Lagrange 等边、Euler 共线、正 \(N\) 边形（\(N=2\ldots8\)）在 \(e=0,0.4,0.6,0.8\) 下， \(\tau\) 与闭式符合到 \(10^{-15}\)–\(10^{-16}\)，且完全不依赖偏心率。
证明。等边位形 \(d_{ij}\equiv a\)，Lagrange 恒等式给 \(I=a^2P/M\)，而 \(U=GP/a\)， 代入 T3 即得 \(\tau_{\rm hom}=\frac{\pi}{\sqrt2}G\,P^{3/2}/\sqrt M=\tau_{\rm Q}\)。 一般情形正是主定理证明中的 S2 步（Hölder），等号条件即所有 \(d_{ij}\) 相等。∎
这解释了一件此前无人指出的事：Semay 用「包络理论」为量子 N 体推出的公式， 与孙博华用量纲分析猜出的公式，与三体 Lagrange 等边轨道的经典 \(\tau\)， 三者在 \(N=3\) 时是同一个数。包络理论的平均场基态本质上就是一个「尽可能等距」的刚性位形， 所以它必然落在 Hölder 的等号上。 数值：等边 (1,2,3) 质量下 \(\tau_{\rm hom}=33.086305648825\)，\(\tau_{\rm Q}=33.086305649\)，比值 1.000000000
正 \(N\) 边形（等质量、边长参数 \(s\)）另有闭式 \(\tau=\frac{\pi}{4\sqrt2}N^{3/2}S_N\,Gm^{5/2}\)， \(S_N=\sum_{k=1}^{N-1}\csc(\pi k/N)\)。它与两个孙氏公式的比值如下——这正是 \(N\ge4\) 时严格不等号的量化：
| \(N\) | \(\tau_{\rm hom}\)（闭式=数值） | \(\tau_{\rm hom}/\tau_{\rm Q}\) | \(\tau_{\rm hom}/\tau_{\rm C}\) |
|---|---|---|---|
| 2 | 1.570796326795 | 1.000000 | 1.000000 |
| 3 | 6.664324407238 | 1.000000 | 3.000000 |
| 4 | 17.009253552518 | 1.041966 | 6.251795 |
| 5 | 34.184469190917 | 1.088125 | 10.881254 |
| 6 | 59.660041616809 | 1.132368 | 16.985513 |
| 7 | 94.821879793332 | 1.173529 | 24.644105 |
| 8 | 140.987934988088 | 1.211588 | 33.924469 |
ACT V · COUNTEREXAMPLES经典公式为什么不是定律
| 轨道 | \(\tau\)（精确值） | \(\tau_{\rm C}\) | \(\tau/\tau_{\rm C}\) |
|---|---|---|---|
| 二体（任意 \(m_1,m_2\)，任意 \(e\)） | \(\frac{\pi}{\sqrt2}G\frac{(m_1m_2)^{3/2}}{\sqrt{m_1+m_2}}\) | 同左 | 1 |
| Lagrange 等边，\(m_i=m\) | \(\frac{3\pi}{\sqrt2}Gm^{5/2}\) | \(\frac{\pi}{\sqrt2}Gm^{5/2}\) | 3 |
| Euler 共线，\(m_i=m\) | \(\frac{5\pi}{2}Gm^{5/2}\) | \(\frac{\pi}{\sqrt2}Gm^{5/2}\) | \(5\sqrt2/2\) |
| Lagrange 等边，\((1,2,3)\) | 33.086305648825 | 14.367993088 | 2.302779 |
三个反例都是牛顿方程的精确解（中心位形条件残差 \(10^{-16}\)），比值是闭式常数， 不存在数值误差的余地。因子 3 的来源很清楚：它正是 \(N=3\) 的质点对数 \(\binom32\)； Euler 共线的 \(5\sqrt2/2\) 来自其距离比 \(1:1:2\)。
不过 \(\tau_{\rm C}\) 并非毫无意义：由 \(\tau_{\rm C}\le\tau_{\rm Q}\le\tau\)， 经典公式是一条成立但远不尖锐的下界。孙博华的量纲分析抓对了唯一能被量纲分析抓住的东西—— \(\tau\) 的质量量纲与 \(G\) 的幂次；抓不住的正是那个 \(O(1)\) 的对称函数选择， 而「立方和」与「和的立方」的差别恰恰就在这里。数值上 \(2.433\pm0.075\)（廖世俊组） 落在 \(\tau_{\rm C}=2.221\) 与 \(\tau_{\rm Q}=6.664\) 之间，也印证了这一点。
ACT VI · NUMERICS125 条精确周期轨道上的复核
方法：Fourier 参数化的最小作用量搜索（解析梯度，经有限差分校验到 \(10^{-9}\)）， 必要时用变分方程（monodromy）Newton 抛光把轨道精化到真正的 ODE 解。 Newton 里必须消除四个零方向：时间平移、平面旋转、质心平移、Galilei boost （后两者由质心投影消除）。每条轨道都检查闭合残差、位力残差、最小间距、 以及与周期相容的尺度比 \(R_{\rm typ}/R_K\)（后者专门用来剔除 「作用量趋于 0 的逃逸伪解」——在拓扑平凡类中最小化作用量会跑到无穷远，这是搜索中最大的陷阱）。
| 族 | \(N\) | 轨道数 | 质量范围 | \(\tau\) 范围 | \(\min\tau/\tau_{\rm Q}\) | \(\min\) C2 |
|---|---|---|---|---|---|---|
| two_body（\(e=0.6\)，质量延拓） | 2 | 54 | 0.050–4.405 | 0.0242–8.8349 | 0.99999992 | 1.000000 |
| equil（Lagrange 等边，质量延拓） | 3 | 37 | 0.364–1.000 | 3.2823–6.6643 | 0.99999980 | 1.000000 |
| euler_e07（Euler 共线，\(e=0.7\)） | 3 | 23 | 0.552–1.000 | 4.2435–6.6643 | 1.00000217 | 1.000003 |
| euler（Euler 共线，圆） | 3 | 3 | 0.907–1.000 | 7.2366–7.8540 | 1.17631359 | 1.000000 |
| fig8（8 字轨道类，质量延拓） | 3 | 2 | 0.743–0.747 | 5.7817–5.8053 | 1.09981680 | 1.009483 |
| ngon4 | 4 | 2 | 等质量 | 17.009254 | 1.04196589 | 1.000000 |
| ngon5 | 5 | 2 | 等质量 | 34.184469 | 1.08812545 | 1.000000 |
| ngon6 | 6 | 2 | 等质量 | 59.660042 | 1.13236755 | 1.000000 |
其中「C2」指等价命题——质量度量引理 \(\big(\int_0^1 ds/|\mathbf Q|\big)^2\int_0^1|\mathbf Q'|^2ds\ge4\pi^2\)， \(\mathbf Q=(\sqrt{m_1}\mathbf q_1,\dots)\)；它与主定理在 \(V=\bar I\) 时等价。
三点值得注意
(1) 没有反例，只有取等。 所有比值的最小值是 \(1-2\times10^{-7}\)，出现在闭合残差最差（\(1.9\times10^{-6}\)）的那条轨道上； 闭合残差达机器精度（\(\sim10^{-12}\)）的轨道，比值全部 \(\ge1\) 到全精度。 更重要的是，逼近 1 的恰好是理论预言取等号的两族： 全部 54 个二体质量比，以及全部 37 条 Lagrange 等边轨道（质量低至 0.364）—— 后者正是 Hölder 等号条件「所有 \(d_{ij}\) 相等」的体现。 这不是数值巧合，是定理的等号集被数值独立复现。
(2) \(\tau\) 只依赖（同伦类，质量）。 对 8 字轨道分别在 \(T=6.3059\)、\(6.3259\)、\(6.3359\) 三个不同周期下做 Newton 抛光，得到 \[ \tau = 9.237681250724,\quad 9.237682750851,\quad 9.237682750851 , \] 即三条不同大小的轨道给出同一个 \(\tau\) 到 12 位。 这是标度不变性的直接体现，也解释了为什么最小作用量解的 \(\tau\) 是一个良定义的 「类函数」——廖世俊组观测到的 \(\tau/L_f\approx\) 常数，其根源就在这里。
(3) 8 字轨道的精确数据。 从公开初值出发 Newton 抛光到闭合残差 \(7.7\times10^{-14}\)： \(T=6.325913976757529\)，\(E=-1.2871419918\)，\(\tau=9.2376812001\)， 角动量 \(3\times10^{-15}\)，位力残差 \(2\times10^{-14}\)，\(\tau/\tau_{\rm Q}=1.3861391847\)， \(\tau/\tau_{\rm C}=4.158417554\)，\(|\bar{\mathbf q}|=10^{-8}\)（故主定理的无条件形式适用）。
ACT VII · OPEN还剩下什么没被证明
| 命题 | 状态 | 说明 |
|---|---|---|
| \(\tau\) 严格标度不变 | 定理 | 势的齐次性，教科书级 |
| \(\tau=\frac12A_U\sqrt{A_K}\) | 定理 T1 | 位力定理的重新包装 |
| homographic 轨道闭式 \(\frac\pi{\sqrt2}\sqrt I\,U\)，与偏心率无关 | 定理 T3 | 覆盖任意 \(N\)、任意质量 |
| \(\tau_{\rm hom}\ge\tau_{\rm Q}\)，等号 ⟺ 等距 | 定理 T4 | Hölder + Lagrange 恒等式 |
| \(\tau\ge\tau_{\rm Q}\sqrt{V/\bar I}\)（一般）；\(\tau\ge\tau_{\rm Q}\)（当 \(\bar{\mathbf q}=0\)） | 主定理 | 等号集含全部二体椭圆与三体等边轨道 |
| 孙氏经典公式作为普适定律 | 已否定 | 精确反例：因子 3 与 \(5\sqrt2/2\) |
| \(\tau\ge\tau_{\rm Q}\) 对所有周期解（去掉 \(\bar{\mathbf q}=0\)） | 开放 | 125 条轨道全部支持；缺口是 \(\sqrt{V/\bar I}\)，二体情形显示该因子是证明技术的产物而非真实障碍 |
| 质量度量引理 \((\int ds/|\mathbf Q|)^2\int|\mathbf Q'|^2\ge4\pi^2\) | 开放 | 等价于上一条；二体全偏心率取等，几何上无拓扑障碍（\(N\ge3\) 时 \(\mathbb R^{2N-2}\setminus\{0\}\) 单连通），故须用动力学 |
| \(\tau/L_f\approx\) 常数（廖世俊组的经验律） | 开放 | \(\tau\) 是「类+质量」的函数已由数值确认；为何近似正比于自由群字长仍无解释 |
| 量子版本的严格化 | 开放 | 包络理论只给变分上限；且 \(N\ge4\) 时与 homographic 值差 4%–21% |
把主定理推到无条件形式，等价于证明质量度量引理。二体情形提供了强有力的线索： 它对所有偏心率精确取等（数值上 \((\int ds/d)^2\int|q'|^2ds/4\pi^2=1-10^{-13}\)）， 说明这条引理若有证明，等号集必然包含全部开普勒椭圆——这是一个非常强的刚性要求， 也提示正确的证明应当把「径向 Wirtinger」与「角向」两部分以某种恰好互补的方式配对， 正如主定理中 Jensen 亏格与 Wirtinger 亏格在二体情形精确互消那样。
REPRODUCE复现
全部代码在 math/kepler3/：
analytic.py（闭式公式、中心位形求解、Hölder 亏格）、
core.py（N 体流、自适应 Gauss–Legendre 诊断、Fourier 最小作用量、monodromy Newton、形状球拓扑）、
validate.py（V1–V5 精确解校验）、
prove.py（S1–S7 证明链逐步数值复核，33 项检查）、
search.py / continue.py（搜索与质量延拓）、
report.py（汇总与反例筛查）。
计算在 192 核节点上完成。
REFERENCES文献
- 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
- B. H. Sun, Classical and Quantum Kepler's Third Law of N-Body System, Results in Physics 13, 102144 (2019). arXiv:1904.05148
- C. Semay, Quantum support to BoHua Sun's conjecture, Results in Physics 13, 102167 (2019). arXiv:1812.00807
- 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
- C.-Y. Zhao and M.-J. Zhang, A conjecture on Kepler's third law of n-body periodic orbits (2018). arXiv:1811.00735
- 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
- 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
- X. Li and S. J. Liao, Collisionless periodic orbits in the free-fall three-body problem, New Astronomy 70, 22–26 (2019). arXiv:1805.07980
- 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
- V. Dmitrašinović and M. Šuvakov, Topological dependence of Kepler's third law for collisionless periodic three-body orbits, Phys. Lett. A 379, 1939–1945 (2015).
- 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).
- R. Montgomery, A new solution of the three-body problem, Notices AMS 48, 471–481 (2001).（形状球与自由群分类）
- S. J. Liao and X. M. Li, On the periodic solutions of the three-body problem, National Science Review 6, 1070 (2019). PMC8291409
上一篇的结论需要修改一处：孙博华的量子公式不只是「未证明的猜想」，
在 \(\bar{\mathbf q}=0\) 的一大类轨道（含全部等质量 choreography）上它是已证的尖锐下界，
并且它恰好等于三体 Lagrange 等边轨道的精确 \(\tau\)——这是一个可以写成一行恒等式的事实。
经典公式则相反：它被精确反例否定为定律，只能作为一条较弱的下界成立。
仍然开放的，是把下界从 \(\bar{\mathbf q}=0\) 推广到全部周期解——等价于证明一条关于闭曲线的
纯几何-动力学不等式，而二体开普勒椭圆在任意偏心率下精确取等，说明它若成立就必然非常刚性。