---
title: "先说清什么叫晕，再读反馈怎样改变聚类"
description: "从晕的密度定义和TNG单位出发，读懂分辨率、重子响应与2026反馈拟合论文中的残差和适用范围。"
canonical_url: "https://sandboxphysics.com/zh/stories/cosmic-web/choose-what-you-call-a-halo/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/cosmic-web/choose-what-you-call-a-halo/"
last_modified: "2026-11-12"
---

[从小起伏到宇宙蛛网](https://sandboxphysics.com/zh/stories/cosmic-web/)

从小起伏到宇宙蛛网 · L4

# 先说清什么叫晕，再读反馈怎样改变聚类

Sandbox Physics · 2026-11-12

密集的结点不自动等于一个定义明确的暗物质晕。要比较模拟或读反馈论文，先固定对象、单位和分母，再问一条曲线究竟支持多强的结论。阅读时打开TNG聚类论文1707.03397v2的图6至7，以及Zhang预印本2609.00807v1的图1、图4和表1。

## 开始前需要什么

读过功率篇，知道密度相对平均背景的含义。本文是带问题的原文阅读，不复现TNG模拟，也不下载其快照。

## 一团物质，怎样划出边界

在图上圈住一个密集结点，圈多大才算晕？球形过密度法先选中心，再从它向外量包围质量，找到球内平均密度达到指定阈值的半径。质量名称必须带着这个阈值；这条边界由测量规则划出，物质仍可穿过它。

$M_{200\mathrm c}$以临界密度为参照，$M_{200\mathrm m}$以同一时刻的平均物质密度为参照。前者的临界密度由当时膨胀率定义，后者还乘上当时的物质密度参数。

$$
\rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G},\qquad \bar\rho_m(z)=\Omega_m(z)\rho_{\rm crit}(z)
$$

两种球形边界分别要求球内平均密度为$200\rho_{\rm crit}$或$200\bar\rho_m$。

$$
M_{200\mathrm c}=\frac{4\pi}{3}R_{200\mathrm c}^3\,200\rho_{\rm crit}
$$

质量与半径需对同一个中心、同一个密度剖面一起求解。

对应阅读：IllustrisTNG official data specifications · Springel et al., First results from IllustrisTNG: matter and galaxy clustering (2018), arXiv:1707.03397v2

## 换参照密度，也会换圈进来的质量

若某时刻$\Omega_m=0.3$，两个阈值分别是临界密度的200倍与60倍。对球内平均密度向外降低的同一剖面，后一个阈值会在更大的半径达到，圈入更多质量。不能把质量简单乘$1/0.3$：球内平均密度变了，半径和包围质量也同时变。

本Lab采用的EdS背景恰有$\Omega_m=1$，所以在相同中心与剖面下，这两种阈值相等。真实模拟若含宇宙学常数，通常不能沿用这个巧合。

朋友链接法按粒子间距连成组，SUBFIND再寻找子结构，球形过密度质量则按选定中心和密度阈值计量。它们回答不同的问题。单凭过密度边界也没有证明所有成员都受束缚，或整团已经达到维里平衡；要给教学场景里的一团亮点命名，也应先确定这些判定规则。

对应阅读：IllustrisTNG official data specifications · Original finite-mesh calculations and independent small-wave Python replay

## 拿到数据以前，先换对单位

TNG快照的Coordinates用$\mathrm{ckpc}/h$：先除以$h$得到共动kpc，再乘当时尺度因子得到物理kpc。固定质量暗物质粒子还要检查header的MassTable；质量单位为$10^{10}M_\odot/h$。不要把不同类型粒子的可变质量字段一概套到暗物质上。 这里的小写$h$是今天的Hubble常数除以$100\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$所得的数。

例如一个说明单位的假想坐标值$100$，取$h=0.6774$、$a=0.5$，对应约$147.6$共动kpc、$73.8$物理kpc。这是算术例题，不是我们读取过的一颗TNG粒子。

快照Velocities存储的是本动速度除以尺度因子平方根，恢复本动速度要乘$\sqrt a$。它与本Lab的正则动量不同，GroupVel等目录字段也需各查自己的定义。看起来都是三列数，并不意味着可以直接互换。

$$
v_{\rm pec}=\sqrt a\,v_{\rm stored},\qquad u_{\rm Lab}=\frac{a\,v_{\rm pec}}{H_\star L}
$$

左式只指上述TNG快照速度字段；右式是本课程的无量纲动量定义。

对应阅读：IllustrisTNG official data specifications · Original finite-mesh calculations and independent small-wave Python replay

## 把分辨率比较，和反馈比较拆开读

打开Springel等人的TNG聚类论文第2节，再看图6：先在同一个盒系列里比较不同分辨率。TNG100与TNG300既改变盒长，也改变质量分辨率。大盒补入更长的波；更细粒子与力分辨率则帮助刻画小结构，两种改进不能互相顶替。

图7转而比较全物理模拟与对应的仅暗物质计算。这里的DMO伴随模拟把总平均物质都作为无碰撞成分表示，不等于只取全物理模拟里的暗物质粒子。比较总物质功率时，分子和分母要对应同一总体。

两图有用的一处细节是对数纵轴。若显示的是自然对数功率比$-0.2$，真正的比值为$e^{-0.2}\simeq0.8187$，约低18.1%。这个数由读图算例指定。

对应阅读：Springel et al., First results from IllustrisTNG: matter and galaxy clustering (2018), arXiv:1707.03397v2 · Original finite-mesh calculations and independent small-wave Python replay

## 2026预印本：残差的百分之几

再读Zhang于2026年9月1日提交的预印本2609.00807v1，重点看图1、图4与表1。它用若干流体模拟及其反馈、宇宙学参数变体研究功率响应。对着表1阅读，先找出各次计算究竟改变了什么。核读截至2026年11月12日，所用版本标注拟投JCAP。

先读轴名：横轴是模拟响应，纵轴是拟合减模拟的绝对响应差。自拟例题$S_{\rm sim}=0.80$、$S_{\rm fit}=0.82$给出残差$0.02$；除以模拟值才得到相对误差$2.5\%$。两者不是同一个“百分之二”。 图1和图4的左右面板分别固定在$k=1h/\mathrm{Mpc}$与$3h/\mathrm{Mpc}$，同一符号系列还会在多个红移出现。

然后比较模型复杂度。固定成单参数时，强反馈例子在$k=3h/\mathrm{Mpc}$处的绝对残差可到约$0.13$。三参数形式在主要测试的$k\leq3h/\mathrm{Mpc}$、$0\leq z\leq2$范围，把最大绝对残差控制在约$0.03$以内，但$\Omega_m=0.1$压力测试仍是例外。它支持在这些模拟与尺度上研究低维近似，不能据此证明所有重子物理只有三个自由度。

$$
S(k,z)=\frac{P_{\rm hydro}(k,z)}{P_{\rm DMO}(k,z)},\qquad \Delta S=S_{\rm fit}-S_{\rm sim}
$$

区分响应本身、绝对响应差与相对误差$\Delta S/S_{\rm sim}$。核读图的颜色区分红移，符号区分模拟系列。

对应阅读：Zhang, Baryonic feedback suppression of the matter power spectrum (2026), arXiv:2609.00807v1 · Original finite-mesh calculations and independent small-wave Python replay

## 最后把结论带回自己的模型

试着写一段自己的结果说明：测量哪一种物质的功率，使用什么晕定义，波数单位是什么，与哪份伴随模拟相除，误差按什么分母报告。然后再补测试范围，以及尚未通过的例子。让这些选择决定一句话能说多远。

本课程的无碰撞PM提供一个可复算的引力基线，没有气体、恒星形成或黑洞反馈。改大初始振幅不能模拟重子排出的作用，也不能从这里的密集结点直接预测TNG的晕目录。要迈向那一步，需要新增物理、合适的空间分辨率及对应的独立验证。

同样，拟合功率响应准确还不等于弱引力透镜的宇宙学参数已经无偏。观测核、噪声、参数退化和校准误差都要进入具体推断；预印本中关于内部缓解的设想，留给这一步继续检验。

对应阅读：Zhang, Baryonic feedback suppression of the matter power spectrum (2026), arXiv:2609.00807v1 · Springel (2005), The cosmological simulation code GADGET-2, astro-ph/0505010v1 · Original finite-mesh calculations and independent small-wave Python replay

## 练习与解释性答案

1. 已知某时刻$\Omega_m=0.3$，能否只用这个数把$M_{200\mathrm c}$换成$M_{200\mathrm m}$？

   展开解释

   不能。还需要中心和密度剖面，求出另一个阈值对应的半径与包围质量。平均密度降低不意味着质量按它的倒数线性增加。

2. 图上给出$\ln(P/P_{\rm ref})=0.1$。这是功率增加10%吗？

   展开解释

   精确比值为$e^{0.1}\simeq1.1052$，增加约10.52%。小对数差可近似百分比，但需要说明近似。

3. 某拟合在$S_{\rm sim}=0.5$处有$\Delta S=0.02$，相对误差是多少？

   展开解释

   是$0.02/0.5=4\%$。绝对响应差仍是$0.02$；若拿它与别处的百分比比较，须先统一分母。

## 动手检验预测

### [写出一个结点尚未告诉你的事](https://sandboxphysics.com/zh/labs/cosmic-web-shell-crossing/)

1. 选“三向会合的结点”，暂停在晚时刻，观察亮点的尺度与共动网格间距。
2. 列出把它称作某个$M_{200\mathrm c}$晕还缺少的定义和分辨率检查。
3. 读原文图6、图7，再写一句能由本Lab支持的引力结论，以及一个需要新增重子物理的问题。

[上一篇：把导出的初值，交给另一个积分器](https://sandboxphysics.com/zh/stories/cosmic-web/rerun-an-export-and-measure-growth/)

## 来源与计算范围

1. [Springel (2005), The cosmological simulation code GADGET-2, astro-ph/0505010v1](https://arxiv.org/abs/astro-ph/0505010v1)

   定向核读第2.1节式1至3、第3.2及3.2.1节式20至23、第4.1节式26至31，查看原PDF第3和11页。用于宇宙学动量、PM窗口与积分方法背景。教学模型另行规定截断和有限梯度，没有短程树力。

2. [Springel et al., First results from IllustrisTNG: matter and galaxy clustering (2018), arXiv:1707.03397v2](https://arxiv.org/abs/1707.03397v2)

   固定v2，核读方法第2.1至2.2节与原图5至7。图6比较分辨率，图7比较全物理与DMO功率，响应坐标使用自然对数。第2.2节讨论有序粒子的散粒噪声限制。这里只读原方法与图，没有重跑TNG或下载快照。

3. [IllustrisTNG official data specifications](https://www.tng-project.org/data/docs/specifications/)

   2026年11月12日核对header的MassTable、暗物质Coordinates/Velocities以及群目录球形过密度定义。快照与群速度约定不同，字段需分别换算。算例使用假想数值，没有提取粒子，也不再分发TNG数据或原素材。

4. [Zhang, Baryonic feedback suppression of the matter power spectrum (2026), arXiv:2609.00807v1](https://arxiv.org/abs/2609.00807v1)

   2026年9月1日提交，截至11月12日核读，版本标为拟投JCAP。定向核读第2至5节、实际PDF与原HTML的图1至4及表1。保留绝对响应残差、单/三参数差别和低物质密度压力测试例外，没有声称独立复现或已验证透镜后验。

5. [Original finite-mesh calculations and independent small-wave Python replay](https://sandboxphysics.com/stories/cosmic-web/materials/cosmic-web-replay-v1.zip)

   自有已验收核/线性、非线性及五预设双精度计算用于课程图与校准。公开小ZIP仅含一份精确小波当前状态、独立NumPy/SciPy轨迹复算及预期结果。示例由Node执行冻结模型和规范序列化函数产生，不冒充浏览器下载。大型原始预设另行留存，包内没有Engine或原论文。