2.1 N体从直接求和到树码

2.1 N 体：从直接求和到树码与粒子网格

本节摘要：N 体把引力系统写成点质量的常微分方程。直接求和忠实却按粒子数平方涨；星团用个别时间步把近距遭遇从平方墙里局部挖出来；树码与快速多极用远处成团的多极展开换掉大部分配对；粒子网格把长程力交给傅里叶变换。软化长度、时间积分和周期边界，是这条主线真正的分辨率声明。

先说结论

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

1. 说明直接求和的平方墙为何会先锁死宇宙学，而不是先锁死太阳系
2. 解释 Barnes-Hut 开角准则如何把“远处一团”合法化
3. 区分 PM、P3M、TreePM：短程与长程各交给谁
4. 说出软化长度既防止奇点又规定你能声称的最小结构
5. 判断一篇 N 体论文的粒子数声明在夸耀什么、在隐瞒什么

平方墙是怎么挡住提问的

两个质点，力沿连线，与距离平方成反比。N 个质点，每个都要看其余所有人。配对数目随 N 的平方涨，这不是实现细节，这是提问权。太阳系里 N 很小，直接求和可以配合高阶积分器把轨道送到工程精度。球状星团把 N 抬到十万量级，平方已经痛苦，但近距双星才是真正的杀手：时间步被最硬的一对绑死。宇宙学把 N 抬到十亿以上，平方墙变成政治：你不是算得慢，你是根本问不了“盒子里能否长出星系晕的层次结构”。

Aarseth 一脉的直接 N 体没有否认平方，他们改的是时间。个别时间步让硬双星走小步，外围星走大步，用预测把不同时钟对齐。近距遭遇用正则化把奇点从坐标系里赶出去，而不是靠把步长缩到零。这条路线养活了星团蒸发、双星形成、引力波驱动并合的前身研究。它不适合均匀铺满的宇宙学盒子：那里几乎没有“少数硬对、多数闲人”的结构，每个粒子都在长程势里掉落。

对象                 平方墙首先打在哪         先被改掉的不是力，是时钟或空间
  太阳系             几乎打不到               积分器的长期结构
  星团               近距遭遇                 个别时间步与正则化
  星系               粒子数与盘加热           树码
  宇宙学盒子         均匀长程力               网格势与混合
​

💡 关键直觉：N 体算法史不是“越来越准”，而是“把平方账记到你真正在乎的那一对上，其余用多极或网格赊账”。

树、多极、网格：三种把远处变便宜的办法

Barnes 与 Hut 在 1980 年代把空间做成八叉树。一个节点若张角足够小，其内部粒子对远处某粒子的作用可以用质心甚至更高阶多极代替。开角准则是这套方法的良心：张角越大，便宜越大，力越脏。脏不是道德问题，是你必须用同一开角做收敛测试。快速多极方法把多极之间的平移做得更系统，渐近复杂度可以到接近线性，实现与并行的代价也更高。天文共同体长期更爱树码，不是因为它在渐近符号上更美，而是它与自适应、个别时间步、局部加密更合得来。

粒子网格走另一条河。把质量云到规则格子上，用傅里叶变换解泊松方程，再把力插回粒子。长程力变成一次变换的价格，与粒子数和格子数近线性。短程力被网格尺度抹掉，于是有 P3M：网格管远，直接或修正项管近。TreePM 把短程交给树，长程交给网格，成为后来大型宇宙学模拟的主流骨架之一。Gadget 一脉把这套骨架传播成一种共同体语言：你读到的许多宇宙网图，力并不是按平方配对算出来的。

软化长度常被写成“为了数值稳定”。更老实的说法是：你拒绝解析 1/r 在重合时的无穷大，并用一个核把力在某尺度下压平。这个尺度就是你能声称的最小引力结构。小于软化的“亚晕”或“星团核”，不是发现，是核函数的形状。碰撞 N 体与无碰撞 N 体在这里分道：星团要解析密近，软化必须很小并配合正则化；暗物质晕若物理上就是无碰撞的，过大软化会抑制小晕形成，过小软化会让两点弛豫加热假结构。

图：力计算的复杂度阶梯

每一级台阶都在买粒子数，卖的是短程或长程的某种忠实。没有免费的十亿粒子。

积分器、周期盒子，以及你不能用粒子数撒谎的地方

力算得便宜之后，时间积分仍能把故事写坏。蛙跳积分器便宜且辛，长期能量误差抖动而不漂移，成为无碰撞宇宙学的默认齿轮。高阶积分器在光滑力场里更准，遇到突变的近距力可能更脆。星团路线用埃尔米特方案和正则化，因为那里的力不是光滑场。把宇宙学蛙跳直接搬进密近双星，或把星团正则化搬进周期盒子，都是谱系错位。

周期边界是宇宙学 N 体的另一门方言。盒子代表均匀宇宙里的一块，泊松方程在周期条件下用傅里叶求解最自然。大尺度功率被盒子尺寸截断，小尺度被软化和网格截断。所谓“分辨率”，至少要同时报：粒子数、盒子边长、软化、网格。只报粒子数，等于只报铸币厂的产量，不报面额。

两点弛豫是无碰撞声称的隐形税。有限粒子让近距散射累积，系统会在弛豫时标里变热。物理上无碰撞的晕，数值上可能被加热成另一类物体。加密粒子、加大软化、缩短物理时间，是三种不同的税负转移。论文若只展示终态美丽的密度轮廓，不谈弛豫时标与哈勃时标之比，读者有权怀疑轮廓是被煮过的。

⚠️ 常见坑：把软化长度设成“看起来光滑就行”。光滑往往意味着你已经把要研究的小晕或核区判决出局。

概念账本可以写成：

装粒子：质量、位置、速度
选力法：直接 / 树 / 网格 / 混合
声明：软化、开角、网格、盒子
每个子步：
    算加速度
    蛙跳或埃尔米特更新
    个别时间步则只更新到期的粒子
诊断：能量抖动、动量、功率谱、弛豫时标
收敛：开角减半、软化减半、粒子加倍，结构是否还在
​

GRAPE 和后来的 GPU 把直接求和在中等 N 上救活过一段时间，第 4 章会讲硬件如何反过来挑选算法。对本章而言，只需记住：硬件从未取消阶梯，它只是让你在某一级台阶上多站一会儿。第 2.2 节将看到，流体面对的不是平方墙，而是间断与密度对比——那是另一套阶梯。

从校验粒子到宇宙学骨架的人事分流

直接 N 体没有退出舞台，它改当裁判。新的树码或 GPU 实现，仍要用少量粒子的直接求和对照力向量和能量抖动。星团组把双星、三体交换、引力波驱动的并合当成必须用直接或正则化才能信的过程；宇宙学组把同样的近距物理压进软化，换取盒子里的晕统计。两条职业路径共用牛顿力，却不再共用时间步哲学。

粒子网格的周期势还带来一种特殊的伪力：网格各向异性。同一物理上的球对称晕，在格子对角与轴向可以长得略不同。更高阶的云赋值和力插值能减轻，不能消灭。读晕形状的论文，应当问赋值方案；只问粒子数，会把格子的晶体对称误当成暗物质的自相互作用。

TreePM 的交接半径是另一个被低估的参数。短于它的用树，长于它的用网格。交接处若匹配不好，功率谱会在对应尺度上出现小坎。这些坎有时被写成新物理或重子效应，其实是两套力在握手时较了一下劲。收敛测试必须包括移动交接半径，正如必须包括开角。

个别时间步的书本同步也值得单独记账。不同进程上的粒子处于不同时刻，力是用预测位置算的。预测器的阶数不够，密近处会积累相位误差。星团码为此发展出复杂的块步与外推；宇宙学码通常更粗，因为软化已经把最硬的近距切掉。把星团的正则化期待搬到宇宙学 TreePM 上，会失望；反向亦然。

无碰撞声称的审计：弛豫、各向异性与假亚晕

把暗物质做成粒子，等于用蒙特卡罗抽样一个无碰撞分布函数。抽样不足时，两点散射会在远短于哈勃时标的时间里给晕加热，形状变圆，子结构蒸发。这不是冷暗物质的预言，是有限 N 的税。审计税的办法包括：比较物理时标与弛豫时标、检查速度各向异性、用同一初值跑两档粒子数看小晕是否稳定存在。假亚晕在软化尺度附近成群出现时，应先怀疑核函数，而不是急着写暗物质自相互作用。

力计算的并行图还会引入另一种假力：域边界上的树被截断或重复。正确的实现用本地树加远程多极，错误的实现在边界制造缝。功率谱若在对应进程尺度上出现尖刺，优先查并行而不是查宇宙学。第 4.2 节会把这写成方言问题；在 N 体主线上，它已经是力的一部分。

本节速览

• 平方墙锁的是提问权：太阳系、星团、宇宙学撞墙的部位不同
• 个别时间步改时钟：星团路线的核心不是更快的力，是不让硬双星绑架全场
• 树与网格改远处的价格：开角和格子尺度都是必须收敛的参数
• TreePM 是当代宇宙学骨架：短程树、长程网格，不是一种更神秘的物理
• 软化就是分辨率声明：小于它的结构不是发现
• 粒子数必须和盒子、软化、弛豫时标一起读：单独的 N 是广告

下一节离开点质量，进入连续介质：SPH 与网格如何把同一组欧拉方程切成互相长期看不顺眼的两种流体。