怎样读懂HERA给出的21厘米上限

请打开《First Results from HERA Phase II》的期刊版：2026年2月10日发表于《Astrophysical Journal》998卷33号，共44页。我们读第32页图16，再翻到第41页附录F的表5。机构保存的开放全文可直接查看。旧的arXiv:2511.21289v1于2025年11月26日提交，图16身份相同，篇幅与页码不同；阅读时以这里指明的期刊页码为准。

先认出图里的四种东西

图16有八个小面板，各代表一段观测频率，也对应一段红移。由左上角的56 MHz到右下角的216 MHz，观测的宇宙时期逐渐变晚。这份结果用了约两周的Phase II观测子样本；图中的不同面板来自不同频段，并不是对同一颗气泡拍下的八张照片。

先选下排标着176 MHz、红移约七的面板。横轴是波数 $k$ ，越往右，所对应的空间变化尺度越小。纵轴是

\Delta^2(k)=\frac{k^3P(k)}{2\pi^2},

其中 $P(k)$ 是亮温起伏的三维功率谱。 $\Delta^2$ 的单位是 $\mathrm{mK^2}$ ，可以理解为按对数波数区间衡量起伏强度的一种写法。它不是某处氢气的温度，也不是电离体积分数。

粉色标记是从数据估出的功率。纵向误差棒标的是 $2\sigma$ 误差， $\sigma$ 表示作者估计的热噪声误差。向下箭头表示该估计在这份 $2\sigma$ 误差内与零相容。黑色粗虚线标的是 $1\sigma$ 噪声水平，和粉色点、两倍误差棒各有含义。紫色阶梯则是较早Phase I工作给出的上限，只出现在其中两个面板。要比较新旧结果，应比较同类统计量以及相近的红移和尺度，不能只选图上最低的一条线。期刊版图16与§5.2

误差棒顶端横帽的宽度还带着另一条信息：它是窗口函数的半高全宽。一个功率估计会混合一段邻近尺度的信号，窗口函数描述这份混合权重。横帽标出了这段响应范围的宽度。带着这条信息再看图，点的位置和它实际接收的尺度范围便能分开。

为什么功率估计会是负的

从一张无噪声亮温图计算模平方，结果不会小于零。望远镜里的功率估计多了一步：尽量去掉接收机噪声自身带来的正偏置。

一个简单类比是把同一个天空信号 $s$ 分两次接收，得到 $s+n_1$ 与 $s+n_2$ 。 $n_1$ 、 $n_2$ 是两次接收的噪声。把两个结果交叉相乘，噪声彼此独立、平均为零时，它们的交叉项在多次平均中消失，天空信号的共同部分留下来。有限数据仍会在平均值两边波动，交叉功率估计因而可以为负。HERA这篇工作的交叉谱使用相邻时间交错子集，细节见§3.9.4。

图16用空心粉色标记表示负估计。纵轴取对数，不能把负数直接放在对应的负纵坐标上；应结合空心标记、误差棒和表中的带符号数值来读。负估计说明这次带噪声估计落在零以下；理论天空功率仍是非负量。

翻到表5，自己算两行

表5把每个频段的结果分成四列：波数、带符号的功率估计、一次热噪声误差、两倍误差上限。下面摘出两行，功率三列的单位均为 $\mathrm{mK^2}$ 。

**表 7.1** HERA发表表格中的两行，保留印刷数值。
红移	波数 $[h\,\mathrm{Mpc^{-1}}]$	功率估计	一次误差	发表上限
$16.78$	$0.55$	$252112$	$439236$	$1130584$
$7.05$	$0.70$	$-494$	$887$	$1775$

先把正估计保留，把负估计在计算上限时取零，再加上两倍误差。这两行体现的表列关系是

U=\max\!\left(\widehat{\Delta^2},0\right)+2\sigma.

帽子表示从数据得到的估计， $U$ 表示这里的上限。第一行可以完整回算：

252112+2\times439236=1130584\;\mathrm{mK^2}.

第二行按印刷列计算得到 $2\times887=1774\,\mathrm{mK^2}$ ，发表列是 $1775\,\mathrm{mK^2}$ 。这一末位差与各列分别取整相容；应保留原表数值。若把负估计直接加进去，会得到 $1280\,\mathrm{mK^2}$ ，那就没有采用这张表对负估计的处理方式。

这行上限并没有测出气体温度为1775 mK，也没有给出电离气泡直径。它给出这一频段、尺度窗口和分析方法下的起伏功率上限。把 $\mathrm{mK^2}$ 开平方会得到约 $42.1\,\mathrm{mK}$ 的强度尺度，仍对应这个窗口所描述的空间起伏。

横轴里的小写字母也有数值

这里的 $h$ 是无量纲哈勃参数，定义为 $H_0/(100\,\mathrm{km\,s^{-1}\,Mpc^{-1}})$ 。它不是普朗克常数。若采用示例值 $h=0.674$ ，表中

0.70\,h\,\mathrm{Mpc^{-1}} =0.4718\,\mathrm{Mpc^{-1}},\qquad \lambda=\frac{2\pi}{k}\simeq13.3\,\mathrm{cMpc}.

波长用共动长度表示，已经把宇宙整体膨胀的比例从距离里除去。换算前后的物理尺度相同。

在边长20 cMpc的周期立方体里，允许的波数是离散的，最小非零值为 $2\pi/L\simeq0.314\,\mathrm{Mpc^{-1}}$ 。四组有限模式各取一个平均，和图16中有仪器响应与窗口宽度的估计仍有区别。比较理论与观测，需要先用观测窗口处理理论功率，并统一红移、单位和统计定义。两条曲线在横轴上挨得近，只是比较的起点。

低波数的隆起为什么要先看仪器

再回到图16，观察各面板左侧偏高的点。较小波数对应较大尺度，也更容易受到明亮前景泄漏的影响。论文§5.2把这份Phase II数据的主要困难指向天线之间的互耦：一副天线会把部分接收到的电磁波重新辐射出去，附近天线又接收到它。这些附加路径给频谱带来结构，让原本很亮的前景进入想寻找微弱信号的区域。

这也解释了一个乍看奇怪的现象：仪器升级后可覆盖更宽频率，某些尺度上的上限却可能比旧工作宽松。热噪声灵敏度与前景污染要分别考察。这篇论文同时给出了更宽红移范围的上限，也把互耦列为继续改进分析的关键问题。这里采用的是它报告的两周子样本结果；全部Phase II观测最终能达到怎样的上限，还取决于后续数据与前景处理。

一个上限怎样进入早期宇宙的解释

若某种模型在一个观测窗口中预测很强的起伏，而数据给出的上限很低，模型便会受到限制。要把这句话落实到恒星或气体参数，还需计算该模型的加热、电离和自旋温度历史，再用观测响应和统计方法比较。论文第六节展开的正是这一步；两行表格手算只完成了进入它之前的读数练习。

前一章用高自旋温度假设把中性体积分数直接转成亮温。早期气体若仍然很冷，吸收项会改变起伏幅度。同一份功率强弱不能单独决定电离体积分数，更不能给每个气泡定位。Furlanetto、Oh与Briggs，21厘米信号的温度依赖

即使日后功率谱得到清楚探测，画空间图还需要保留更多信息。前一章的两组条纹已有相同功率、不同相对相位。真实成像还要让有限角分辨率、频率分辨率与前景处理允许这些空间信息留下来。功率谱先回答哪种尺度上有多少起伏，成像再回答它们分布在哪里。

留给下一张图的三个问题

某行印刷值是估计 $-20\,\mathrm{mK^2}$ 、误差 $30\,\mathrm{mK^2}$ 。按上面的表列关系，上限是多少？若直接把负数相加，会得到什么？
已知 $h=0.674$ ，把 $k=0.5\,h\,\mathrm{Mpc^{-1}}$ 换成不含 $h$ 的数值，并求对应共动波长。
两种亮温图在所测窗口中有相同功率，一种含连成片的暗区，另一种形状不同。只用这些功率，能够选出空间图吗？还应保存什么信息？

第一题得到 $60\,\mathrm{mK^2}$ ；直接相加会得到40，采用了另一种处理。第二题得到 $0.337\,\mathrm{Mpc^{-1}}$ ，波长约18.6 cMpc。第三题仅凭给定功率无法选图，空间相位以及未被测量或平均掉的模式也有作用。这个结论并不要求任意两幅图都能有相同功率，只需要承认已有统计没有保存完整空间信息。

再读一篇21厘米结果时，先记下它的版本与观测子样本；看清点、误差和上限各是哪一个量；检查尺度单位与窗口；最后再追踪作者怎样把这些统计量接到物理模型。沿着原图和原表走，这些很大的数字就会逐渐变成可以自己判断的证据。