3.3 宇宙学 N 体：从粒子网格到流体宇宙

3.3 宇宙学 N 体：从粒子网格到流体宇宙

本节摘要：宇宙学模拟的问题很明确：线性涨落如何长成网、晕、空洞。粒子网格先让周期盒子里的长程力可算；更高分辨率的树与 TreePM 让晕的内部结构进入议程；千禧年一代把半解析星系挂上晕树，证明骨架已经能支撑巡天式的统计；Illustris、EAGLE、TNG 再把流体和配方装进同一只盒子。世代更替改写的是“结构”一词包含多少重子，而不是把前一代作废。

本节目标

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

1. 说明宇宙学盒子的周期边界与共动坐标解决了什么时钟问题
2. 按世代排列：PM、高分辨纯 N 体、半解析挂载、流体宇宙
3. 解释功率谱、晕质量函数、相关函数为何是公开路标
4. 区分线性理论可检验的大尺度与配方主导的小尺度
5. 判断一张宇宙网可视化图里，细丝的颜色来自暗物质还是来自气体冷却

盒子、共动坐标，以及骨架为何先被算出来

牛顿 N 体加上膨胀背景，粒子的 peculiar 运动才是结构形成。共动坐标把哈勃流拆走，周期边界让傅里叶泊松自然落地。盒子边长截断最大模式，粒子数和软化截断最小模式。宇宙学模拟因此天生是带窗函数的实验：你问的任何统计量，都要声明窗。

1970–80 年代的 PM 与 P3M 工作，包括 Davis、Efstathiou、Frenk、White 一脉，证明冷暗物质可以长出网状结构，而不是纯粹的大团块。这是对模型的定性胜利。定量上，功率谱的形状、晕的密度轮廓、子结构的丰度，要等粒子数涨上去才吵得起来。Navarro、Frenk、White 的密度轮廓成为无碰撞晕的共同体语言，随后的争议是：轮廓在很小半径是否依赖软化、粒子数和暖暗物质一类的扩展。

纯 N 体不包含星系。要把巡天比上，需要一种把晕翻译成星系的规则。半解析模型挂在并合树上：冷却、星形成、反馈写成常微分方程，沿树积分。千禧年模拟把这套做法做到巡天尺度：骨架是 N 体，星系是树上的账本。流体尚未进入盒子，统计已经可以和星系聚类对话。这是方法史上极重要的一课——不是所有对比都值得立刻上流体。

世代                 盒子里实际在演化的           星系如何出现
  PM/P3M             暗物质粒子                   几乎不出现，或事后指认
  高分辨纯 N 体      晕与子晕                     子晕或半解析
  千禧年式           暗物质加并合树               半解析账本
  流体宇宙           气体加星加黑洞               粒子或网格上的重子
​

💡 关键直觉：半解析不是流体失败后的退路，它是骨架已经可信时的正当翻译层。

流体宇宙：同一只盒子开始长颜色

一旦气体进入周期盒子，第 3.2 节的整座栈都要装下，还要跑完从高红移到今天。这是算力与工程的双重墙。Illustris 展示形态多样性可以在宇宙学体积里出现；其反馈也暴露了气体分数等量上的张力。TNG 调整磁场与反馈，EAGLE 把校准写进设计，SIMBA、Astrid 等继续改风与黑洞。流体宇宙的产品不再只是功率谱，而是可与光谱巡天、弱引力透镜、星系群 X 射线对比的合成观测。合成观测本身是一层方法：把模拟粒子投到光锥上、做尘埃衰减、做仪器响应。对比“像不像”之前，先问投影和噪声有没有被同等对待。

大尺度上，流体对物质功率谱的修正是后来的研究主题：反馈把气体推出晕，改变弱透镜信号。这迫使宇宙学与星系形成重新握手——你不能再假装小尺度配方与大尺度参数估计无关。方法史在这里拐弯：宇宙学模拟不再只为“证明冷暗物质长网”，而要为参数推断提供可靠的理论模板。

图：宇宙学模拟世代

每一代都留下仍在使用的遗产。流体宇宙没有把千禧年式的树作废，参数扫描仍常常走半解析。

初值、宇宙学参数，以及不要编造的对照

初值来自早期宇宙的涨落谱，用 Zel'dovich 或更高阶拉格朗日微扰把粒子放下。起始红移太低，小尺度已经非线性；过高则步数浪费。二阶拉格朗日初值减少瞬变，成为后来的标配。宇宙学参数——物质密度、谱指数、涨落振幅——改变骨架。流体宇宙若只跑一组参数，就把配方与宇宙学锁死在同一次实验里。要解开，需要网格状的参数扫描，或第 5 章会提到的仿真器。

对照必须用公开事实。线性功率谱在足够大的尺度上应贴住理论；晕质量函数有拟合形式可对；重子声振荡的尺度是宇宙学标准尺。具体某次巡天的拟合优度数字，若不是你正在分析的那份数据，不要随手编。方法教程只需要说：这些路标存在，模拟必须先过路标，再谈形态学的审美。

⚠️ 常见坑：把流体宇宙里好看的细丝，直接写成冷暗物质被证实。细丝在纯 N 体里已经存在；颜色来自气体时，你在看反馈。

概念账本：

选盒子边长与粒子数 → 声明窗函数
给初值：涨落谱 + 拉格朗日位移
纯 N 体跑到今天 → 晕寻找与并合树
可选：半解析 或 流体栈
统计：功率谱、质量函数、相关函数
合成观测前：先过线性与晕的路标
​

宇宙学把骨架做到最大。下一节走到相反的极端：视界附近，牛顿势不再够用，网格要铺在弯曲时空上。

从证明模型到服务推断：角色转换

早期盒子的任务是定性：冷暗物质能否长网，热暗物质是否把小结构抹掉。任务完成后，角色转为给观测推断提供模板：弱透镜的物质谱、重子声振荡的标准尺、星系偏置、21 厘米功率谱。模板必须覆盖参数空间，单点流体宇宙不够用。于是出现模拟套件：同一条款下改宇宙学，或同一宇宙学下改条款。套件的存在，说明领域承认单次著名模拟不再等于宇宙。

重子对物质功率谱的抑制或增强，尺度大约在兆秒差距以下到晕尺度，精确形状随风与星形成条款而变。参数估计若忽略这一层，会把条款误差吞进宇宙学参数。当代做法是边缘化条款，或用仿真器在条款空间插值。第 5 章的机器学习岗位，有很大一块是为这个角色转换服务的。

初始条件的二阶拉格朗日位移、起始红移、转移函数的辐射成分，都会在小尺度留下瞬变。纯 N 体的小尺度收敛因此有起点效应，不只是粒子数效应。把起点红移写进清单，和写软化一样基本。流体宇宙的起点还要决定何时打开冷却与紫外背景，打开得太早或太晚会改小晕的气体分数。这些是世代成熟后才被认真对待的细节，早期定性胜利并不需要它们。

寻找器、定义质量与“同一晕”的幻觉

不同晕寻找器对桥接物质的归属不同，球形超密与椭圆友元给出的质量可以对不上。子晕在潮汐里被剥，寻找器可能提前宣判死亡。并合树跨快照匹配时，过疏的全家福会把飞掠当成并合。半解析吃的是树，流体宇宙的卫星统计吃的是寻找器。对比观测的卫星计数，必须先对齐定义。所谓“同一晕在不同模拟里的命运”，常常是不同寻找器的命运。

质量定义还有孔径问题：维里半径、200 倍临界密度、绑定粒子，给出不同数字。弱透镜测的是投影质量，X 射线测的是热气体质量。合成观测必须把定义推到与观测同一孔径。只在模拟内部用一种定义画质量函数，再与另一种定义的观测比，是窗函数错误，不是宇宙学错误。

问题：流体宇宙是不是让千禧年式的半解析过时了？

没有。参数扫描、并合树统计、许多聚类问题，半解析仍然是正当翻译层。流体宇宙把配方变成动力学，也把条款推进宇宙学误差预算。两者分工：骨架与树回答“晕如何长”，流体回答“气体如何被按住”。用一张彩色细丝废除半解析，会把定性骨架胜利和条款不确定度混为一谈。细丝颜色若来自气体，你在审冷却与反馈；网的存在本身，纯 N 体已经给出。

寻找器与孔径定义能让“同一晕”变成幻觉。维里质量、两百倍临界密度、绑定粒子、投影透镜质量，数字可以对不上。半解析吃树，卫星计数吃寻找器，弱透镜吃投影。对比前先对齐定义，否则窗函数错误会冒充宇宙学错误。质量函数图上的每一根柱，都已经过翻译，不是粒子数的直接读出。骨架世代回答网如何长，流体世代回答气体如何被按住，两代试卷不能互相废除。

世代史还可以压成一条河床：粒子网格长出网，树码解析晕，并合树挂上半解析星系，流体把配方推进误差预算，光锥把盒子翻译成巡天。后一代不废除前一代，只是把提问权往更贵的层推进。读一张彩色宇宙网，先问颜色来自暗物质还是气体，再问寻找器叫什么名字。大尺度对线性理论，小尺度对配方，两张卷混审，就会把窗函数错误写成新物理。

核心回顾

• 周期盒子加共动坐标：先解决时钟与长程力
• 骨架先于星系：纯 N 体与半解析是正当世代，不是过时玩具
• 千禧年证明树能支撑巡天统计：流体不是对比的唯一门票
• 流体宇宙让配方进入宇宙学误差预算：透镜时代尤其如此
• 大尺度对线性理论，小尺度对配方：不要混审
• 细丝颜色要问来源：暗物质骨架与气体冷却是两代产物

下一节进入致密天体：后牛顿、数值相对论、GRMHD 如何把视界从禁区改写成实验室。