---
title: "怎样读懂HERA给出的21厘米上限"
description: "从HERA的功率谱图和两行原表入手，分清带符号的估计、误差、窗口与上限，再看这些数字怎样限制早期宇宙。"
canonical_url: "https://sandboxphysics.com/zh/stories/reionization/reading-hera-upper-limits/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/reionization/reading-hera-upper-limits/"
last_modified: "2026-11-12"
---

[星光怎样改变星系间的氢](https://sandboxphysics.com/zh/stories/reionization/)

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

Sandbox Physics · 2026-11-12

HERA给出的上限，限制的是某些尺度上21厘米亮温起伏的强度。要读懂它，先看每个点究竟估计了什么，再看误差与仪器怎样进入这个估计。即使把强度量得很准，要画出气泡的位置还需要空间相位。

请打开《First Results from HERA Phase II》的[期刊版](https://doi.org/10.3847/1538-4357/ae2d54)：2026年2月10日发表于《Astrophysical Journal》998卷33号，共44页。我们读第32页图16，再翻到第41页附录F的表5。[机构保存的开放全文](https://api.repository.cam.ac.uk/server/api/core/bitstreams/d98f6b1c-b4da-41a6-ad94-e4343baacde5/content)可直接查看。旧的[arXiv:2511.21289v1](https://arxiv.org/abs/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](https://doi.org/10.3847/1538-4357/ae2d54)

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

## 为什么功率估计会是负的

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

一个简单类比是把同一个天空信号 $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厘米信号的温度依赖](https://arxiv.org/abs/astro-ph/0608032v2)

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

## 留给下一张图的三个问题

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

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

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

## 来源

1. [期刊版](https://doi.org/10.3847/1538-4357/ae2d54)
2. [机构保存的开放全文](https://api.repository.cam.ac.uk/server/api/core/bitstreams/d98f6b1c-b4da-41a6-ad94-e4343baacde5/content)
3. [arXiv:2511.21289v1](https://arxiv.org/abs/2511.21289v1)
4. [Furlanetto、Oh与Briggs，21厘米信号的温度依赖](https://arxiv.org/abs/astro-ph/0608032v2)

## 继续阅读

星光怎样改变星系间的氢 · 阅读路径

### 从光子、电离体积到无线电信号

1. [文章 **星系的光怎样改变宇宙中的氢**](https://sandboxphysics.com/zh/stories/reionization/starlight-and-hydrogen/)
2. [文章 **电离一半究竟是什么意思**](https://sandboxphysics.com/zh/stories/reionization/ionized-volume-and-local-fraction/)
3. [文章 **一片膨胀的气体需要多少光**](https://sandboxphysics.com/zh/stories/reionization/photons-in-an-expanding-box/)
4. [文章 **两个电离泡相遇以后，光去了哪里**](https://sandboxphysics.com/zh/stories/reionization/when-ionized-bubbles-overlap/)
5. [文章 **从中性氢的信号读出电离泡**](https://sandboxphysics.com/zh/stories/reionization/reading-a-neutral-hydrogen-map/)
6. [文章 **亲手复算一段电离历史**](https://sandboxphysics.com/zh/stories/reionization/recompute-a-fixed-source-history/)
7. [当前页面 **怎样读懂HERA给出的21厘米上限**](https://sandboxphysics.com/zh/stories/reionization/reading-hera-upper-limits/)