2.3 可积性障碍：布伦斯、庞加莱与桑德曼的三重变奏
2.3 可积性障碍：布伦斯、庞加莱与桑德曼的三重变奏
本节摘要："可积"在天体力学里是一张执照，持有它的系统才能被守恒量彻底驯服。本节先给执照下严格定义，再依次过堂三份历史文件：布伦斯证明再无新的代数积分，庞加莱证明一般情形再无新的解析积分并顺手发现混沌的雏形，桑德曼的收敛级数则是一次"账面和解"——存在性赢了，实用性输了。这三份文件合起来，就是"三体无通解"这句流行语的全部内涵。
行话定义：「可积」是一张什么执照
上一节结尾留下的悬念，用一个行话来正式收束：可积性。它的严格说法是：一个自由度为 n 的哈密顿系统若存在 n 个彼此独立的、两两对合（泊松括号为零）的守恒量，则称其为可积系统。这句话的含金量在于它的兑换表：满足条件后，运动可以被"作用量-角度变量"完全拆开——每个守恒量钉住一层不变的环面，轨迹贴着环面匀速打转，永不出圈。预测变成查表：给初值，查它落在哪张环面，算出打转的频率，未来任意时刻的位置一式写出。
二体问题持照上岗：化约后自由度只剩极少几个，动量、角动量、能量恰好凑够，于是每条轨道都是一张标准环面——椭圆就是环面在坐标面上的影子。第 1 章推过的"圆锥曲线通解"，从可积性语言看就是这张执照的实物。
三体问题申请同一张执照时，账本危机发作：九个自由度需要九笔互相兼容的守恒账，人类全部库存凑不出这个数。于是问题变成一个精确的数学追问——**到底是不容易找，还是根本不存在？**十九世纪末，两份否定性文件先后给出答案。
布伦斯障碍：代数积分先出局
1887 年，海因里希·布伦斯证明：三体问题不存在经典的十个（动量、角动量、能量、质心类）之外的代数首次积分——即不能由坐标与动量通过加、乘、开方等代数运算构造出来的新守恒量。这份文件的力度在于它封的是一整类：任何指望"再找一个多项式形状的守恒量补账本"的方案，从 1887 年起就在数学上关门了。
代数积分出局，非代数的积分还开着门——比如无穷级数、隐函数、极限过程。庞加莱的文件处理的就是这个更宽的口子。
庞加莱障碍：解析积分也出局，顺手捡到混沌
1890 年那篇因国王悬赏而重写的论文里，庞加莱考虑的是限制性三体问题的近可积情形：两主星之外的小天体质量趋于零时，系统近似可积；把小质量作为小参数逐级往上加，人类两千年来对付这种问题的看家本领是摄动级数——把解展开成小参数的幂级数，一项项修正。开普勒轨道之所以算得准，靠的就是这套手艺对太阳系各行星逐级修正。
庞加莱证明的事情用现代语言说就是：这套级数在关键处发散。不存在除经典账本之外、对小参数一致收敛的解析首次积分。更深刻的发现在证明路上：发散的机制不是"误差慢慢变大"，而是一种几何灾难——他追踪的某族轨道的同宿流形（稳定与不稳定轨迹的接缝）发生了横截相交。一旦横截相交，接缝附近就会无穷次折叠，织成今天叫"同宿缠绕"的精细网络：任何解析积分想穿过这张网，都会被折叠得支离破碎。混沌不是从庞加莱的论文里"被预言"，而是直接在纸面上被画了出来——这是它第一次以数学对象的身份出现。
两个障碍合起来的准确措辞值得原样记住：一般三体问题不存在足够多的、足够正则的（代数的或解析的）首次积分。它没有否定：数值可算、特殊解存在、统计规律可用——这正是第 3、4、5 章各自接手的活。
桑德曼级数：一次账面和解
判决之后，1912 年芬兰数学家桑德曼给出了迟到二十年的"和解方案"：他证明三体问题的解可以写成时间的幂级数——技术上需要先把时间变量开立方再展开——且级数对所有时间收敛。存在性问题就此终结：三体的解不但存在、唯一，还能写下来。
欢呼很快冷场。收敛慢到什么程度？后人估算，用桑德曼级数把轨道算到实用天文精度，需要的项数在 10^{8\,000\,000} 量级——把已知宇宙的原子全变成纸也写不下这些项。级数收敛"在数学上"，与"在你有生之年算出来"，中间隔着一条欲望的鸿沟。1991 年王琪东（Q. Wang）把存在性推广到一般 N 体情形，补齐了纯数学的版图；但计算实践早就集体转向了另一条路——数值积分，也就是第 4 章的主题。
动手复盘：近距交会为什么是级数与积分器的共同天敌
实验背景：桑德曼级数之所以要把时间开立方，根源在碰撞奇点——两体距离趋于零时引力势发散，解虽可延拓但导数不再光滑。这不是纸面上的洁癖，它直接决定了数值积分器在近距交会时的表现。操作：让两体做一条近距交会的轨道（近日距设得很小），分别用固定大步长与固定小步长积分穿过交会点，监测能量误差。代码如下：
import math mu = 4 * math.pi ** 2 def run(dt, r_peri): """二体椭圆轨道，近日距 r_peri，固定步长 dt，蛙跳积分一个周期""" a = 1.0 e = 1 - r_peri / a pos = [r_peri, 0.0] vel = [0.0, math.sqrt(mu * (1 + e) / (a * (1 - e)))] def acc(p): r = math.hypot(p[0], p[1]) f = -mu / r ** 3 return [f * p[0], f * p[1]] E0 = 0.5 * (vel[0] ** 2 + vel[1] ** 2) - mu / math.hypot(*pos) acc_now = acc(pos) t, period = 0.0, 2 * math.pi * math.sqrt(a ** 3 / mu) while t < period: vel = [vel[i] + 0.5 * dt * acc_now[i] for i in range(2)] pos = [pos[i] + dt * vel[i] for i in range(2)] acc_now = acc(pos) vel = [vel[i] + 0.5 * dt * acc_now[i] for i in range(2)] t += dt E1 = 0.5 * (vel[0] ** 2 + vel[1] ** 2) - mu / math.hypot(*pos) return abs((E1 - E0) / E0) print("近日距 0.5，大步长:", run(0.01, 0.5)) print("近日距 0.01，大步长:", run(0.01, 0.01)) print("近日距 0.01，小步长:", run(0.0005, 0.01))
结果：近日距 0.5 时大步长的能量误差在 10^{-6} 量级；近日距缩到 0.01，同一大步长的误差暴涨几个数量级；步长缩小二十倍后误差才压回来。解读：误差不是随轨道变扁线性变差，而是被近距点的力梯度陡增引爆——力在交会瞬间变化太快，固定步长"踩不准"。桑德曼在纸面上用重标时间对付的正是这个奇性；数值实践的对策则有两派：自适应步长（交会时自动加密）与碰撞正则化（KS 变换等专用工具，把奇点在坐标变换里抹平），第 4 章会再遇它们。变式：把近日距逐步减半重跑，观察误差爆炸的临界步长如何随距离缩小——你会亲眼看到"奇点附近一切常规手段失灵"这句警告的真实分量，也就理解了为什么三体散射的数值研究都要郑重声明处理了碰撞。
三重变奏的总谱
把三份文件叠起来看，十九世纪末到二十世纪初的这组"否定性成果"其实分工明确：布伦斯封掉代数积分，庞加莱封掉解析积分并交付混沌的几何机制，桑德曼保住存在性却放弃实用性。三者共同的效果是把问题从"求解"改写成了三个更可行的小问题——找特殊解、做数值实验、问统计规律。第 3 章沿第一条路走，第 4 章沿第二条，第 5 章的混沌理论就是第三条的入场券。
⚠️ 复述定理时别越界：庞加莱证明的是"不存在某些正则类型的首次积分"，不是"三体系统不可预测"，更不是"三体问题被解决或被证死"。每一版流行转述多越一次界，读者就多一分对混沌的误解。
本节要点回顾
- 可积执照：足够多、彼此对合的守恒量换来环面上的规则运动，二体持照、三体缺账；
- 布伦斯文件：代数首次积分只有经典十笔，1887 年后"再找一个多项式守恒量"这条路关门；
- 庞加莱文件：解析首次积分同告出局，同宿缠绕是发散的几何机制，也是混沌的首次亮相；
- 桑德曼变奏：收敛级数终结存在性争议，天文数字的项数宣告它与实践无缘；
- 近距交会：碰撞奇点是级数发散与积分器失效的共同根源，重标时间与正则化是对策；
- 路线图：特殊解、数值实验、统计规律三条侧路就此正式通车。
障碍立在眼前，缝隙也从这里开始。下一章进入猎物图谱：在混沌的海面上，先去看那五座人尽皆知的灯塔。