1.1 从开普勒到数字实验

1.1 从开普勒到数字实验

本节摘要：计算天体物理学的起点不是某台机器的开机仪式，而是解析天体力学把能积的积分积完之后，剩下那些不可积的运动。开普勒给出几何轨道，牛顿给出力的语言，拉普拉斯用摄动把太阳系维持成可预报的钟；庞加莱则证明一般三体没有全局解析积分。数字实验从此成为合法手段：在守恒律看守下，让方程自己往前走。

学习目标

阅读完本节，你应当能够：

1. 用两体可积与三体不可积，说明“有方程”和“有封闭解”不是一回事
2. 说出摄动论在太阳系预报里完成了什么，以及它为何无法推广到星团与星系
3. 解释数字实验与演示动画的差别：前者受守恒与初值约束，后者只求看起来像
4. 指出二十世纪中叶恒星结构与早期 N 体各自接过了哪一段解析欠账
5. 为后文的离散化做好动机：我们不是喜欢网格，是解析已经停在河边

钟表宇宙是怎么被上紧发条的

把计算天体物理的第一幕写成“计算机发明之后”，会漏掉真正决定方向的那两百年。开普勒用第谷的火星数据逼出椭圆，几何已经够用：太阳在一个焦点，面积速度守恒，周期平方与半长轴立方成正比。牛顿接着做了一件更危险的事——他把几何改成力。危险之处在于：力是瞬时的、成对的、可叠加的。两体问题因此可以化成一个等效的单体制心运动，轨道再次闭合。太阳系看起来像一座被校准过的航海钟。

航海钟这个类比不是修辞。经度问题逼着十八世纪把天体位置算到角分，月球理论、木星卫星掩食都是预报工具。拉普拉斯把行星相互摄动写成级数，勒威耶与亚当斯后来能指出海王星该在哪一带，靠的就是这座钟还走得动。摄动论的工作方式很像给已经上好油的齿轮加一丁点侧向力：主运动仍是椭圆，修正项按质量比和小偏心率展开。只要扰动小、时间不太长，纸和笔就赢了。

赢是有边界的。边界不在“公式写得不够美”，而在叠加假设本身。星团里没有一颗占绝对优势的中心天体；星系盘里恒星数以千亿计；宇宙学里每个质点都在其余所有质点的势里掉落。摄动级数在这些地方不是收敛得慢，而是展开点根本不存在。

解析传统能做的                    解析传统做不到的
  两体 → 圆锥曲线                    一般三体的全局积分
  小摄动 → 级数修正                  密近遭遇、交换、逃逸
  周期轨道的稳定性估计               星团蒸发、星系并合
  太阳系数十年预报                   十亿年的结构形成
​

💡 关键直觉：解析天体力学不是失败了，它是把能积的部分积完了。计算登场，是去接管积分表空白的那些页。

图：从封闭轨道到受纪律的数字实验

解析传统把能积的积分积完之后，提问方式才从“寻找封闭形式”换成“设计可重复演化”。中间不是机器发明，是不可积把门关上了。

庞加莱切断的不是希望，是封闭形式

十九世纪末，瑞典国王悬赏三体问题。庞加莱交上的答案让悬赏的含义反过来：一般三体不存在足够多的解析积分，相空间里密布着对初值敏感的轨道。钟表隐喻在这里断裂。断裂之后，物理学家并没有停止写方程，他们停止了“写完方程就等于理解运动”这种自信。

不可积不是不可算。差分可以一步一步走，只是每一步都在借未来的精度还现在的债。二十世纪上半叶，人们用机械积分仪和人工计算组做限制性三体、周期轨道搜索。电子计算机到来之后，同一笔债换了一种偿还速度。冯·诺依曼一代人清楚：机器擅长的是重复的算术，不是重新发现一个守恒量。

天体物理比纯天体力学更早感到饥渴。爱丁顿的恒星结构已经是一组微分方程：流体静力学平衡、质量连续、能量输运、物态。主序还可以靠相似性论证和手算积分；到了壳层燃烧、对流边界、热脉冲，手算积分变成了一种礼貌的拖延。亨耶伊方法把恒星结构变成边值问题的松弛求解，这是计算天体物理第一次以“内部结构”而不是“点质量轨道”的身份出现。与此同时，奥尔特、霍姆伯格、图姆雷用极少的粒子去看星系盘会不会自己长出旋臂。粒子少得可笑，问题却已经不是摄动论能提问的类型。

所谓数字实验，指的是：给定方程、给定初边值、给定守恒检查，让系统演化，再问结果是否依赖步长、粒子数、边界。它像实验室，因为你可以改一个参数再跑；它又不像实验室，因为你改的是模型而不是自然。混淆这两点，会把模拟图当成照片。

数字实验的纪律：什么可以变，什么不能装看不见

早期计算的名声并不好。有人觉得那是“不会积分的人在蛮干”。这种批评有一半是对的：没有收敛检查的演化，只是昂贵的故事。另一半则过时了。当解析已经证明不存在全局积分，蛮干之外的选择只剩下两种——放弃提问，或接受受纪律约束的推演。

纪律有三条，后来整个领域都在重复它们，只是换了名字。

第一条，守恒是账本。两体问题里能量和角动量可以精确守；N 体里它们变成诊断量。能量漂移过大，不是宇宙在告诉你新物理，是你的积分器或时间步在撒谎。星系模拟里角动量泄漏会让盘变厚，看起来像“加热”，其实是算法在磨损齿轮。

第二条，初值不是装饰。宇宙学模拟的初始条件来自早期宇宙的涨落谱；恒星碰撞的初值来自轨道参数；超新星的初值来自一维演化码交来的球对称模型。换一套初值，故事可以完全不同。把某次著名模拟的终态当成“宇宙必然如此”，等于把一次实验的样品当成物种。

第三条，分辨率是提问权。你分辨不了的尺度，你就没有资格断言那上面发生了什么，除非你显式引入亚网格模型，并承认那是模型。第 5 章会回到这里；现在只需记住：数字实验的资格来自分辨，不来自配色。

⚠️ 常见坑：把可视化里好看的旋臂、喷流、细丝直接写成“已被模拟证实的自然”。先问：同样的结构在两倍粒子数或两倍网格下还在不在。

二十世纪六十到七十年代，纪律逐渐变成共同体习惯。Aarseth 的直接 N 体把近距遭遇当作必须认真积分的事件，而不是麻烦。Miller 和 Prendergast、Hohl 的盘模拟开始讨论热稳定性。White 和 Rees 稍后给出星系形成的层次图景，计算与理论第一次在宇宙学尺度上互相需要。这些工作今天看起来粒子数寒酸，方法史的意义却很大：它们证明，提问方式已经从“寻找积分”换成“设计可重复的演化”。

概念上，一次最小的数字实验可以写成下面这样的账本，而不是一段能跑的工程脚本：

给定：粒子或网格上的状态
每个时间步：
    算力或通量
    更新速度与位置，或更新守恒量
    记下能量、角动量、质量
结束时：
    对照解析极限或线性理论
    对照更粗和更细的分辨率
    只有三份对照都还活着，才讨论物理
​

这本账看起来朴素，却把计算从“画图”里救了出来。第 2 章的树码、粒子网格、SPH、自适应网格，都是在这本账越来越厚、对象越来越多时，去改“算力或通量”那一行的成本。没有这本账，复杂度再低也只是更快地画出不可信的故事。

数字实验如何进入天文台的日常

二十世纪中叶之后，数字实验不再只活在数学系的三体搜索里。射电巡天给出的中性氢速度场，逼着动力学家用自洽的盘模型去对；球状星团的光度函数和蓝离散星，逼着碰撞 N 体去解释质量分层；超新星遗迹的 X 射线温度，逼着流体码去对塞多夫解以外的非均匀介质。观测仪器每提高一档角分辨率，解析模型能调的自由参数就显得更不够用。计算因此不是从理论单向输出到观测，而是被观测的细节拖进必须逐步推进的地带。

这一拖，改变了论文的写法。早期天体力学论文以定理和级数为主；计算论文必须交代初值、步长、守恒诊断、对照的解析极限。读者若跳过这些段落只看图，等于跳过实验记录只看照片。共同体后来把“代码比对项目”当成一种集体体检：同一激波管、同一静水球、同一宇宙学初值，多家方言分别提交。差异被写进方法，而不是被修辞抹平。第 1 章后面的离散化与稳定性，正是为了让你能读懂这些体检表。

航海图的隐喻在这里仍然有用。解析积分是沿岸的灯塔，数字实验是离岸之后必须打开的航位推算。推算要记海流（耗散）、记罗经差（约束误差）、记海图比例尺（分辨率）。不记这三项，船也可以走，只是你不知道自己在哪。

要点速记

• 解析先完成能完成的：两体与小摄动把太阳系维持成可预报的钟，这是计算的前提，不是计算的反面
• 不可积切断封闭形式：一般三体没有全局解析积分，星团与星系更没有展开点
• 数字实验是受纪律的推演：守恒、初值、分辨率三条，缺一就退回演示
• 恒星结构与早期 N 体同时登场：计算天体物理一开始就不是单一的轨道问题
• 可视化不是证据：结构必须在分辨率对照下存活
• 后面所有算法都在还同一笔债：让不可积系统在可检查的账本里往前走

下一节将把“必须算”落实到方程本身：牛顿运动、流体守恒、辐射转移，为何写得出却解不出，离散化究竟改写了对象的哪一层。