跳到主要内容
Sandbox Physics
从小起伏到宇宙蛛网

从小起伏到宇宙蛛网 · L4

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

Sandbox Physics ·

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

一团物质，怎样划出边界

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

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

ρcrit(z)=3H(z)28πG,ρˉm(z)=Ωm(z)ρcrit(z)\rho_{\rm crit}(z)=\frac{3H(z)^2}{8\pi G},\qquad \bar\rho_m(z)=\Omega_m(z)\rho_{\rm crit}(z)
两种球形边界分别要求球内平均密度为200ρcrit200\rho_{\rm crit}或200ρˉm200\bar\rho_m。
M200c=4π3R200c3 200ρcritM_{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

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

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

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

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

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

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

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

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

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

vpec=a vstored,uLab=a vpecH⋆Lv_{\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-0.2，真正的比值为e−0.2≃0.8187e^{-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。

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

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

S(k,z)=Phydro(k,z)PDMO(k,z),ΔS=Sfit−SsimS(k,z)=\frac{P_{\rm hydro}(k,z)}{P_{\rm DMO}(k,z)},\qquad \Delta S=S_{\rm fit}-S_{\rm sim}
区分响应本身、绝对响应差与相对误差ΔS/Ssim\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. 已知某时刻Ωm=0.3\Omega_m=0.3，能否只用这个数把M200cM_{200\mathrm c}换成M200mM_{200\mathrm m}？

    展开解释

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

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

    展开解释

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

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

    展开解释

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

动手检验预测

写出一个结点尚未告诉你的事

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

来源与计算范围

  1. Springel (2005), The cosmological simulation code GADGET-2, 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

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

  3. IllustrisTNG official data specifications

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

  4. Zhang, Baryonic feedback suppression of the matter power spectrum (2026), arXiv: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

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

反馈这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策