3.2 八字轨道：数值发现与变分法确权
3.2 八字轨道：数值发现与变分法确权
本节摘要：等质量三体可以沿着同一条 8 字形曲线互相追逐、周而复始——这就是八字轨道（figure-eight orbit）。它 1993 年由摩尔在计算机上撞见，2000 年由尚晨与蒙哥马利用变分法证明存在并确认线性稳定。本节给出可直接运行的初值数据，亲手复现这个周期解，并用量化实验看清它"稳定得有多勉强"。
一个打破预期的猎物
提起三体轨道，多数人的想象是一团甩开的乱麻。1993 年，克里斯托弗·摩尔在研究粒子系统的数值实验里撞见了一个反例般干净的图形：三个等质量天体排成一列，沿同一条 8 字形闭曲线鱼贯而行——你追我、我追他、他追你，一圈之后严丝合缝地回到起点。这类"多个天体沿同一条曲线排队追逐"的解后来得到一个专用名字：编排解（choreography）。八字轨道是编排解家族中最著名的一员。
它的稀奇之处值得摆上桌面。其一，初值毫无花样：三个质量全相等，总角动量恰好为零——而零角动量在直觉里近乎"最混沌"的初始条件。其二，它居然稳定：给初值一个微小扰动，轨道并不立刻散架，天体会自行靠拢回 8 字形附近——这在周期解里是稀缺品质，大多数周期解一碰就碎。其三，它的存在性证明用上了当时三体研究里最锋利的数学工具，成了一场"数值发现加数学确权"的标准示范。
从发现到确权：一场跨越七年的接力
摩尔的原始贡献是数值的：他在模拟中观察到 8 字形闭轨，并测出了近似初值。但数值实验有个天然软肋——你看到的"周期轨道"可能是积分误差编造出来的赝品，第 4 章会专门解剖这类风险。要让它成为数学事实，需要一套不依赖积分器的证明。
2000 年，阿兰·尚晨与理查德·蒙哥马利在《数学年鉴》上给出了证明，工具是变分法：把"周期轨道"翻译成"作用量的极小值点"——在所有具有 8 字形对称性的闭合路径中，找出使牛顿力学作用量取极小的那条，再证明这条极小路径不经过碰撞，因而是真实解。整套证明不写一行积分代码，纯靠变分原理压轴。随后西莫（Carles Simó）补上第三块拼图：数值线性化分析确认 8 字轨道是线性稳定的，并沿"质量不再全相等"的方向把它延拓出一大批亲戚。
这场接力之所以值得细看，是因为它后来成了范本：计算机负责"找到嫌疑解"，变分法负责"定罪为真解"，线性化负责"评估稳定成色"。本章下一节的上千族周期解，用的正是这套三人流水线。
亲手复现：把 8 字跑出来
实验背景：验证八字轨道的周期性与"稳定得有多勉强"。操作：采用通行的归一化（质量全为 1、引力常数为 1），初值取精修后的高精度数据；用第 2 章的受力零件加 RK4 积分器推进一个周期，检查三个天体是否回到出发点；然后故意污染初值，观察轨道能撑住几个周期。
import math G = 1.0 m = 1.0 # 八字轨道的高精度初值（归一化单位） q1 = [-0.97000436, 0.24308753] q2 = [ 0.97000436, -0.24308753] q3 = [ 0.0, 0.0] v1 = [ 0.4662036850, 0.4323657300] v2 = [ 0.4662036850, 0.4323657300] v3 = [-0.9324073700, -0.8647314600] state = q1 + v1 + q2 + v2 + q3 + v3 # 12 个数：3 体 × (2 位置 + 2 速度) T = 6.32591398 # 周期（归一化时间单位） def accel(q): """q 为 6 个数（3 体位置），返回 6 个加速度分量""" out = [0.0] * 6 for i in range(3): for j in range(3): if i == j: continue dx = q[2 * j] - q[2 * i] dy = q[2 * j + 1] - q[2 * i + 1] r = math.hypot(dx, dy) + 1e-12 f = G * m / r ** 3 out[2 * i] += f * dx out[2 * i + 1] += f * dy return out def rk4_step(s, dt): q, v = s[:6], s[6:] a1 = accel(q) q2_ = [q[i] + 0.5 * dt * v[i] for i in range(6)] v2_ = [v[i] + 0.5 * dt * a1[i] for i in range(6)] a2 = accel(q2_) q3_ = [q[i] + 0.5 * dt * v2_[i] for i in range(6)] v3_ = [v[i] + 0.5 * dt * a2[i] for i in range(6)] a3 = accel(q3_) q4_ = [q[i] + dt * v3_[i] for i in range(6)] v4_ = [v[i] + dt * a3[i] for i in range(6)] a4 = accel(q4_) qn = [q[i] + dt / 6 * (v[i] + 2 * v2_[i] + 2 * v3_[i] + v4_[i]) for i in range(6)] vn = [v[i] + dt / 6 * (a1[i] + 2 * a2[i] + 2 * a3[i] + a4[i]) for i in range(6)] return qn + vn # 主实验：推进一个周期，测回归误差 dt = 1.0e-4 s = list(state) for _ in range(int(round(T / dt))): s = rk4_step(s, dt) q_err = max(abs(s[i] - state[i]) for i in range(6)) print("一周期后位置回归误差:", q_err) # 总角动量（应为 0）：逐体算位置与速度的叉积分量再累加 Lz = state[0] * state[7] - state[1] * state[6] \ + state[2] * state[9] - state[3] * state[8] \ + state[4] * state[11] - state[5] * state[10] print("总角动量 z 分量:", Lz) # 变式：污染初值 1e-5，看几个周期后开始失真 s = [x + (1e-5 if i == 0 else 0.0) for i, x in enumerate(state)] kept = 0 for cycle in range(1, 200): for _ in range(int(round(T / dt))): s = rk4_step(s, dt) if max(abs(s[i] - state[i]) for i in range(6)) > 0.1: print("扰动轨道在约", cycle, "个周期后偏离 8 字形") break kept = cycle
结果：干净初值下一个周期后的位置回归误差在 10^{-7} 量级（受步长与精度限制，步长折半误差平方缩小），总角动量在浮点意义下为零——8 字确属零角动量解；被污染初值的轨道能在 8 字形附近坚持成百个周期才慢慢漂走。解读：第一层，回归误差证明这是真周期轨道，不是数值赝品——注意这里的前提是步长足够小，把步长放大十倍重跑，"周期"会明显变形，这正是第 4 章误差对账的预告。第二层更有意思：扰动轨道的"存活期"长而有限，说明 8 字轨道的稳定不是超稳定——它处在相空间里一个很小的稳定岛中央，岛内指数不发散，岛外就是混沌海。蒙哥马利有一句流传很广的调侃：三体的稳定解"像针尖上的舞池"。变式：把污染幅度从 10^{-5} 换成 10^{-3} 与 10^{-7} 各跑一次，存活周期数应大致按对数规律变化——这是李雅普诺夫指数概念的感性前奏，第 5 章会给它定量化。
图 3-3 八字轨道与编排
为什么这个解值得单独一节
八字轨道在猎物清单上的地位，超过它作为"一个漂亮解"的本身。对数学，它验证了变分法在三体这类强非线性问题上的杀伤力——作用的极小化天然偏爱对称构型，8 字之后，同一条流水线沿着对称群的不同选择，量产出一大批新编排解。对物理直觉，它修正了一个流行误读："零角动量加等质量"远不是混沌的同义词，规则运动在初值空间里的分布比想象的更细碎、更顽强。对工程想象，它展示了周期解作为"免费基础设施"的潜质：一条稳定的闭轨，就是一个可以反复利用的运输走廊——这个念头在第 6 章的行星际高速路一节会重新出现。
本节要点回顾
- 猎物身份：等质量、零角动量、平面编排解，三体沿同一条 8 字曲线互相追逐，周期约六个归一化时间单位；
- 流水线范本：数值发现（1993）、变分法确权（2000）、线性稳定分析（西莫），三步接力成为此后找解的标准流程；
- 可复现数据：本节初值可直接运行，一周期回归误差随步长平方缩小，是验证积分器的好标尺；
- 稳定的成色：线性稳定但稳定岛极小，扰动轨道可存活成百周期后漂走——存活期随扰动幅度对数变化；
- 观念修正：零角动量不等于必然混沌，规则孤岛在初值空间的分布远比直觉细密。
8 字证明了孤岛存在，但它是"猜出来"的孤岛。下一节看猎具升级后的收获：当发现流程变成参数扫描加精修，孤岛从一座变成上千座。