科学网—万有引力和三体（2）
精选
已有 275 次阅读
2026-11-7 07:54
| 系统分类: 科研笔记
万有引力和三体（2）
张武昌 2026年11月7日星期六
第五章 瑞典 - 挪威联合王国的悬赏
在整个 19 世纪中后期的科学界、思想界与精英社会中，人类已经形成一种近乎信仰的绝对共识：宇宙是一台结构精密、秩序恒定、规律铁律、完全可被人类数学穷尽的巨型钟表。只要掌握足够的算力、足够精准的初始条件，人类便可以 绝对精确地推演宇宙过去一切状态、预知宇宙未来一切演化 。
上帝彻底退场，误差彻底消失，扰动全部可解，未知全部可测。在维多利亚时代科学乐观主义的顶峰， 人类坚信自己已经站在了洞悉寰宇终极真理的终点 。
然而，历史最宏大的戏剧性永远蕴藏在极致巅峰之后的骤然坠落。就在经典力学权威万古独尊、确定性宇宙观统治人类认知顶点的时代，一场由王室主导、举国瞩目、面向全球学界的 世纪顶级数学悬赏竞赛 ，即将亲手击碎人类两百年的科学信仰，掀开现代科学最震撼、最颠覆、最悲壮的认知革命。
瑞典 - 挪威联合王国的国王奥斯卡二世，年轻的时候正经学过数学，不是附庸风雅的门外汉。他在大学修过数学、天文，看得懂天体力学的大意，知道牛顿、拉普拉斯那条线索，明白 “ 太阳系稳不稳定 ” 是整个时代悬着的大难题。 他长期资助瑞典当时最核心的数学家米塔格 ‑ 列夫勒创办的顶尖数学期刊《 Acta Mathematica 》，本身很愿意以 “ 支持科学的君主 ” 的形象示人。在当时欧洲，君王赞助数学、天文竞赛是一种传统。
1887 年，为纪念国王 奥斯卡二世的 六十寿诞 （1889 年） ，米塔格 ‑ 列夫勒提议和操办了一场轰动欧洲的国际数学悬赏，以抬高自己作为主编的期刊《 Acta Mathematica 》《数学学报》的地位。 他拉上当时最顶尖的两位大牛：魏尔斯特拉斯、埃尔米特，一起拟定四道难题， 第一道，就是 N 体问题：严格证明太阳系是否稳定，能不能写出收敛的级数解 ，另外三道是函数论的题目，只是这第一道名气最大，后世只记住三体这个了。
王室斥巨资设立 人类近代史上规格最高、关注度最广、目标最终极的科学悬赏 。国王以一国君主之尊，向 全世界所有数学家、力学家、天文学家 公开征解那个自牛顿时代起、困扰人类两百年、贯穿莱布尼茨论战、牵动宇宙秩序根本的终极难题： N 体问题通解与太阳系绝对稳定性的 严格数学证明 。
这场皇家悬赏的规格与时代分量，远超普通学术竞赛。它不是小众学界的内部切磋，而是 国家王权背书、全欧洲精英围观、代表人类文明叩问宇宙终极秩序的巅峰答辩 。奥斯卡二世本身极具科学素养，崇尚数理思辨，深知这道问题的终极重量：若能彻底给出通用解析解，就意味着人类 彻底、永久、绝对地锁死宇宙秩序 ，为持续两百年的宇宙稳定性之争画上终极句号，将牛顿 – 拉普拉斯的确定性宇宙观推向绝对真理的神坛，让人类掌握统治星辰运转的终极密钥。
大赛规则庄严而决绝：面向全球学者开放征解，不限国籍、不限资历，以严密的数学证明为唯一标尺，优胜者将获得巨额王室奖金与欧洲学界至高荣誉。消息一经发布，瞬间席卷整个西方学术圈层，成为 19 世纪末最盛大、最万众期待的科学盛事。
彼时的学界万众笃定、信心爆棚。所有人都默认：三体、多体问题只是 人类算力不足、技巧未臻极致 留下的遗留难题，绝非自然规律的本质局限。既然连隐匿深空的未知行星都能被笔尖推演捕获，那区区多体摄动的通用解，必然只是时间与技巧问题。整个科学界已然做好准备，迎接 经典力学最后一块拼图的圆满归位 ，迎接人类确定性宇宙体系的最终完满。
在全球无数顶尖学者的角逐之中，彼时崭露锋芒、天赋冠绝欧洲的法国天才数学家 亨利 ・ 庞加莱 入局攻坚。作为横跨数学、物理、天体力学的旷世奇才，庞加莱耗费数年心血，穷尽当时人类全部顶级数学工具，对限制性三体问题展开极致推演。在漫长的孤绝演算之后，庞加莱交出了一份逻辑严密、体系恢弘、推演精致、看似完美无缺的重磅论文，一度被评审团认定为 完美解决世纪难题的终极答卷 。
举国欢庆、学界沸腾，欧洲学界一致认定：两百年悬案终被攻克，人类宇宙秩序彻底尘埃落定。王室奖金即将颁发，经典力学的终极加冕仪式即将圆满落幕，人类即将迎来 完全确定、完全可控、完全可预测的终极宇宙图景 。
可就在论文付梓印刷、即将公之于众、载入史册的最后关头，一个细微到极致、却足以颠覆整个人类科学根基的 致命漏洞 ，被校对学者意外发现。
庞加莱在复盘推演时骤然洞悉了一个碾压时代、颠覆认知的恐怖真相： 三体及多体系统之中，存在一类对初始条件的极致敏感依赖性 。在多体引力耦合的复杂系统里，星体初始位置、速度的 无穷微小偏差 ，并不会随时间推移平均消解，反而会持续放大、指数级叠加、层层裂变，最终导致整个天体系统的演化轨迹彻底偏移、全然不可预测。
这不是计算误差，不是方法缺陷，不是人力局限 —— 这是宇宙系统与生俱来、根植于自然法则深处的本质属性 。
这一刻，庞加莱亲手推开了人类认知宇宙最幽暗、最真实、最颠覆的一扇大门。
他推翻了自己所有的既有结论，撕碎了即将封神的完美答卷，以极致的学术坦诚，向全世界公布了这个震碎时代的终极真相： 不存在通用的、永久收敛的三体解析通解；多体系统不具备长期可精确预测性；牛顿 – 拉普拉斯构建的绝对确定性宇宙，彻底崩塌。
这场奥斯卡二世皇家悬赏，就此迎来人类科学史上最悲壮、最戏剧性的惊天反转。
原本为 终结难题、圆满经典力学、固化确定性宇宙信仰 而生的王室大赛，最终非但没有为旧体系加冕，反而亲手 掘开了经典决定论的坟墓 。
1988 Wisdom & Sussman 算出冥王星轨道是混沌的，算是庞加莱混沌预言验证的序曲。
1989 年，拉斯卡（ Jacques Laskar ） ：在《 Nature 》发表，用平均久期方程积分 2 亿年，第一次给出内太阳系混沌的数值证据，证明行星轨道初值敏感、长期不可预报，印证庞加莱的判断。
1992 年，萨斯曼、威斯多姆（ Sussman & Wisdom ） ：直接完整 N 体辛积分，不做平均近似，在《 Science 》确认整个太阳系混沌，是大家最常引用的 “ 计算机模拟证实太阳系混沌 ” 标志性工作。
Laskar, J. (1989). A numerical experiment on the chaotic behaviour of the Solar System. Nature , 338(6212), 237–238.
Sussman, G. J., & Wisdom, J. (1992). Chaotic Evolution of the Solar System. Science , 257(5066), 56–62.
引力多体系统长期轨道无法精确预报。 计算机数值模拟不是数学证明，只是独立的数值实验，直观复现了庞加莱预言的定性行为。
第六章 一般三体
庞加莱在奥斯卡二世悬赏的震撼结论，为 19 世纪的天体力学画上了一道决绝的休止符。他已经证明，人们梦寐以求的一般三体通用解析解并不存在，对初始条件的敏感依赖，注定了多体引力系统长期轨道的不可精确预报。但这并不意味着人类对三体问题的探索就此止步；恰恰相反，庞加莱的定性分析，为后来者开辟了一条全新的道路。
此前一个多世纪里，学界的目光大多集中在 限制性三体问题 —— 也就是其中一个天体质量微乎其微，几乎不干扰另外两个大天体的运动，这是一种刻意简化、用来降低数学难度的近似模型。而悬在所有力学家头顶的终极难题，依旧是 一般三体问题 ：三个天体都拥有可观质量，彼此之间互相施加引力，没有谁可以被当成无足轻重的尘埃。
进入二十世纪，随着数学工具的持续拓展，加上电子计算机这一划时代工具的诞生，人类终于拥有了直面这一难题的全新武器。研究者开始有意识地把研究重心，从 “ 小质量天体在两大天体引力场中的运动 ” ，转向三个质量相当、彼此地位对等的三体构型。这是一次研究范式上的重大转向：限制性三体是宇宙的局部近似，而三个质量相近的天体共舞，才是三体问题最原始、最完整、最残酷的本来面貌。
在解析的道路上，数学家依旧在艰难搜寻零星的特殊周期轨道。 1912 年，桑德曼给出了一般三体问题的级数解，这个结果在理论层面意义重大，但它有着与生俱来的致命短板：级数收敛速度慢到难以想象，没有任何实际的计算价值，无法用来追踪天体漫长的轨道演化。它在数学上宣告解的存在，却无法满足天文学家想要看见轨道、预判走向的诉求。桑德曼的成果像是一个象征性的界碑，它告诉世人，解析道路的潜力已经走到了边界，想要窥见一般三体纷繁的演化图景，必须另寻出路。
真正的变革，伴随着计算机的诞生徐徐拉开序幕。在二十世纪上半叶，数值计算依靠人力手摇计算机，计算成本极高，只能做很短时间尺度的轨道推演，稍微长一点的演化，海量的四则运算足以耗尽一名计算员数年的光阴。二十世纪中叶之后，电子数字计算机登上舞台，算力的爆发式增长，彻底改写了三体研究的格局。人们不再执着于写出一个可以一劳永逸描述轨道的数学公式，转而拥抱数值模拟这一条全新的求索路径：给定一组精确的初始位置、初始速度，依靠微分方程的分步迭代，一步一步追随三个天体在引力拉扯下的运动，把抽象的微分方程，转化为屏幕上可以动态观看的轨道轨迹。
等质量三体，成为这一轮数值探索里最核心的试验场。相较于质量悬殊的构型，三个质量近乎相等的天体，不存在谁主导、谁依附，引力耦合强度旗鼓相当，混沌效应表现得最为剧烈、最为纯粹。在模拟之中，研究者亲眼见证了庞加莱理论预言的图景：仅仅是初始速度、位置上毫厘之差，在一段时间之后，整个系统的演化走向就会彻底分道扬镳。一部分三体系统会维持一段时间的束缚状态，在反复的靠近、绕转、引力弹弓式的相互抛掷之后，最终走向散射解体：其中一个天体获得足够动能，被另外两个天体甩出去，剩下两者形成稳定的双星系统，这就是三体散射最经典的结局。还有极少数特殊的初始条件，能够落入稳定的周期轨道，三个天体周而复始，沿着固定的闭环轨迹循环运动，这类解极为稀少，如同混沌汪洋之中零星的孤岛。
二十世纪下半叶，一大批统计性的模拟工作陆续铺开，研究者不再纠结于单个轨道的长期预报，转而转向统计规律。他们批量生成成千上万组随机初始条件，对大量等质量三体系统开展大规模数值试验，统计系统解体的概率、轨道寿命、逃逸天体的速度分布。
这是一次深刻的认知转变：既然我们无法精准预言某一个三体系统千万年之后的命运，那我们可以去寻找混沌背后的统计秩序。 Hut 与 Bahcall 在 1983 年的经典工作，正是这一方向的里程碑，他们把三体散射当成一个统计物理问题来处理，建立起三体相互作用的统计框架，把三体问题从纯粹的天体力学，和恒星动力学、星团演化紧密联结在了一起。
二十一世纪，刘慈欣的《三体》走到台前。
作家拾起了这个在天体力学里沉寂了三百年的古老难题，将抽象的微分方程、混沌效应、不可预报的轨道演化，转化成一幅可以被所有人感知的宏大图景：一颗行星，被三颗质量可观、互相牵引的恒星裹挟，在毫无规律的恒星引力场里浮沉。没有稳定的年、固定的季节，恒纪元与乱纪元交替降临，文明一次次在炙烤与冰封之中覆灭、重启、挣扎求生。小说里三体人所有的苦难、焦虑、向外寻找新家园的执念，根源正是一般三体问题与生俱来的混沌本性 —— 不存在一套可以永久预知的节律，没有能够一劳永逸掌握的星辰时刻表。
这并非简单的借用名词，而是一次跨越学科边界的深情呼应。限制性三体是学者简化出来的模型，桑德曼的级数只是理论上的存在，计算机模拟里成千上万次的散射与逃逸，是冰冷的数字结果；而《三体》把这一切动力学的本质，赋予了文明的重量、生命的悲欢。庞加莱笔下对初始条件敏感的混沌，不再只是数学上的一句结论，它化作整个种族悬在头顶永恒的生存危机。
作品一经问世，浪潮席卷而来。它先是在国内掀起阅读热潮，继而跨越语言的壁垒，走向世界文坛，斩获国际科幻领域至高的荣誉。数以千万计的读者，第一次因为一部小说，听见了那个从牛顿时代就开始回响的古老疑问：当三颗星辰彼此拉扯，宇宙的秩序究竟还能不能被把握？很多读者并不能读懂摄动理论、看不懂桑德曼级数，也无法复现三体数值模拟的迭代过程，但他们能读懂混沌带来的绝望，读懂人类数百年来想要读懂星辰、预测命运的执着。
科学与文学，在这里完成了一次动人的合鸣。
第七章 太阳系的质心
两个天体围绕公共质心旋转；两个天体的质量相近时质心悬浮在两天体之间；质量差距逐步拉大，质心向质量更大天体偏移；当质量悬殊足够大，质心落入大天体内部，表现为小天体环绕大天体运行。双星系统、恒星 – 行星系统，没有天然的物理分割，只是同一套动力学连续谱上不同的表象。
太阳系中由于太阳的质量比行星大很多，所以太阳绕着质心旋转，太阳的旋转半径（太阳质心和太阳系质心之间的距离）很小，有时候太阳系质心落在太阳内部，有时候在外面。太阳系八大行星看上去就象小天体沿着椭圆轨道环绕主星太阳旋转。
在《自然哲学的数学原理》第三卷 “ 宇宙的系统 ” 之中，牛顿已经明确论证 N 体引力系统的公共质心是惯性中心： 太阳并不静止在太阳系的几何中心，太阳与所有行星共同绕系统公共质心运动；该质心要么静止，要么匀速直线运动
牛顿已经掌握这个力学原理，但是他并没有做定量数值计算，没有画出太阳摆动的轨迹、没有分解周期。在牛顿的时代，天文学家的计算重心，全部放在行星轨道的摄动上面；太阳本身的摆动，只是一条推论，没有单独作为研究对象。 在整个 18 世纪、 19 世纪，天体力学的核心目标是行星历表、行星摄动、彗星轨道、水星近日点这类问题。太阳绕质心的摆动，仅仅是编制行星历表时一个附带的坐标修正项，不是独立课题。勒威耶、纽康姆在计算水星轨道摄动的时候，他们的历表底层已经隐含这套质心变换，但他们不会专门去研究太阳本身的这条回旋曲线。
太阳质心摆动，在地球上直接光学观测极难测出 。它的视线方向速度幅度只有大约 12 m/s ，远远小于当年光谱仪、天体测量仪器的误差阈值；和行星摄动不一样，行星轨道偏差可以靠凌日、冲日、子午环位置长年累积读出，太阳自身的这一晃动，在地面传统天体测量里，被太阳本身的视差、自转、仪器系统误差掩盖，没有直接观测证据，长期属于 纯理论计算的运动 ，不是直接观测得到的现象。
1965 年， Paul D. Jose 在 The Astronomical Journal 发表《 Sun’s motion and sunspots 》，是这个方向的里程碑 。 Jose 使用 Eckert 等人 1951 的巨行星积分结果，向前回溯到 1610 年、向后推演至 2010 年，完整绘制太阳相对于太阳系质心的运动轨迹；他识别出由木星（ 11.86 a ）、土星（ 29.46 a ）会合周期叠加形成的 约 178.7 年准周期 （后世称为 Jose 周期），描述轨迹会在三叶形卷曲形态和无序回旋之间切换，同时讨论太阳绕质心运动带来的角动量变化，并且试探性把它和太阳黑子周期联系起来。 从这篇文章开始，太阳绕质心摆动，才从历表里面一个不起眼的坐标修正项，变成一个独立、专门的动力学研究对象。
Fairbridge & Shirley （ 1987 ）、 Charvátová （ 1988 ）跟进 Jose 的工作，进一步确认 179 年准周期的轨迹交替特征：一段时间太阳轨迹是对称三叶环，一段时间是缠绕、没有闭合趋势的无序段； Charvátová 后续长期研究这两种模态的切换。
同时在这一段时期，这套恒星摆动的力学原理，成为 视向速度法搜寻系外行星 的底层思想模板：遥远恒星被行星拉扯而 wobble ，和太阳被木星土星拉扯而 wobble 是完全一样的物理。
地面传统光学天体测量，始终无法直接测出太阳的质心摆动。
真正的直接测量，来自甚长基线干涉测量 VLBI 、盖亚卫星 Gaia 的天球参考架。 VLBI 建立国际天球参考架 ICRF （ 1997 ，后续 ICRF2 、 ICRF3 ），在构建惯性参考框架的时候，必须把太阳系质心作为基准原点，并且把太阳相对于质心的运动作为模型项纳入，从射电源的位置残差中间接约束太阳质心运动。
Gaia （ 2013 升空， 2016 之后持续释放数据），凭借微角秒级的天体测量精度，能够在惯性框架下定量约束太阳系质心的长期运动、太阳相对于质心的微小位移，把从前纯粹数值积分的轨道，放到直接观测的约束之下。
太阳质心环绕太阳系质心的运动没有一个恒定的轨道半径，不是固定圆周运动。太阳相对太阳系质心的晃动，是巨行星引力持续调制出来的复杂摆线。当木星、土星等大质量行星分布在太阳两边、引力矩互相抵消时，太阳系质心可以落回太阳中心，二者距离降到接近 0 。当四颗巨行星挤在同一侧，太阳质心和太阳系质心的最大距离可达约 1 530 540 千米，而太阳半径为 695 700 千米，所以这时太阳光球位于太阳系质心之外。
日地距离 1.496 亿千米，太阳 wobble 带来的 太阳在天球上的角位移 很小，只有大约 ±0.01 角秒级别。而我们人眼的分辨极限大概是 60 角秒（ 1 角分），差了几千倍，肉眼完全察觉不到太阳位置的晃动。
时至今日，我们对太阳摆动的轨迹，主要仍然依靠 JPL DE 系列、 IMCCE VSOP 系列行星历表做数值积分；观测的作用，更多是校准历表、修正外海王星天体带来的微小偏差，而不是像行星轨道那样，靠肉眼和望远镜直接看见这个晃动。
终 章
从牛顿提笔演算万有引力的那一刻起，人类就踏上了追问三体的旅途。这条路，有神学的牵绊，有王权赞助的竞赛，有笔尖预言行星的辉煌，有混沌理论突如其来的颠覆，有计算机无数次模拟的漫长积累。而在二十一世纪，《三体》为这段史诗画上了一个面向公众的句点。
科学的探索不会就此停下，新的三体轨道还在持续被发现，新的数值模拟不断揭示多体系统更多隐秘的统计规律，天体力学的远征依旧向前。但属于 “三体” 这个名字的宏大叙事，至此完成了一次完整的闭环：它始于一个困扰近代科学的数学难题，在三百年一代代学者的推演与思辨中层层展开，最终借着文学的光芒，抵达无数人的心底。
星辰的纠缠依旧在宇宙深处无声上演，而人类对秩序、对未知、对浩瀚苍穹的好奇，将和这个三体的故事一道，长久流传下去。
画蛇添足：混沌
混沌不是世界的混沌，（宏观）系统在某一时刻都有一个确定的状态，薛定谔的猫是微观系统，所以可以出现叠加态（两个或多个状态的概率之和为 1 ）。
混沌是数学里的状态，是系统随时间的延长，其不确定性增加，超出科学家的忍受范围。
数学里混沌的出现是现实和理论上都无法避免的，任一参数（参量）的初始值都是一个确定的值，这个值叫做真值，真值的确定按照小数点后有效数字的长短（多少）可以无穷取值，在数学工具里的数值永远达不到真值，都是在不同容忍度下取一定的有效数字，随着时间的延伸，数学计算的结果会出现在越来越粗的管道内，视为混沌。这就是天气预报短期预报可靠性高，长期预报可参考，再长的时间就是混沌、无法预报的原因。
举一个简单的例子，军工厂按照 50 克装药的标准制造子弹，所有子弹打出的轨迹在 200 米内几乎完全相同，但是此后会发散开来，部分原因是 50 克的装药是相同的，但是每个子弹的装药在 50.01 和 50.99 之间有发散。这个精度上的问题永远存在，在足够长的时间内轨道结果就会发散，所以特定子弹的一定时间后的轨道混沌无解（是 数学上 无法预测精确轨迹，不是其轨迹变为薛定谔的猫）。
对于天气预报而言，解决办法就是每天更新初始参数。
上述轨道问题类似预报精度的问题，是提高不同时间精度的轨道所在范围的精度，（例如）从 1 年内的粗轨道发展到一天之内的细轨道位置范围。
蝴蝶翅膀效应是指在（气候模拟的）初始状态参量取值中不能体现出的极小的被舍弃的数位中的由蝴蝶扇动翅膀造成的微小的变化，经过足够长的数值模拟都能造成结果的巨大（不可忽视、科学家不能接受的）和实际值的偏离。其实这个偏离也是在系统振动的大管道内。不是蝴蝶翅膀的扇动通过一个必定的机制掀起了风暴。
转载本文请联系原作者获取授权，同时请注明本文来自张武昌科学网博客。 链接地址： https://blog.sciencenet.cn/blog-39359-1560073.html
上一篇： 万有引力和三体（1） 欢迎参加科学网十佳博文评选活动！ 主办单位： 
支持单位： 