读懂HERA第一批上限怎样避开前景

阅读前需要了解时延谱和楔形，见一条基线怎样把平滑的前景和氢分开；还需要了解带功率、窗函数和注入检验，见把时延功率变成带功率。至于最终上限怎样约束再电离，请读配套导读读懂HERA的上限。

读的是哪个版本，结论怎样标注

论文是HERA合作组的“First Results from HERA Phase I: Upper Limits on the Epoch of Reionization 21 cm Power Spectrum”，2022年发表于ApJ 925, 221。本导读通读的是arXiv第一版2108.02263v1，没有和期刊版逐行对照，所以下文的节、公式和图的编号都按v1。HERA Collaboration 2022, arXiv v1

下面每一条都带四种标注之一。

• 已发表结果：论文本身测得或写明的内容。
• 准备中的工作：v1发布时论文以“in prep.”引用的工作。依赖其细节之前，先查它后来是否已经发表。
• 模型预测：实验室的模拟，不是HERA数据。
• 未解问题：论文自己也没有解释清楚的现象。

第4.3节与图14、图15：真实数据里的楔形

先读那句“没有做任何前景扣除”。后面的一切都是在“避开”。HERA Collaboration 2022, §4.3

已发表结果。HERA的天线固定朝天，没法转开去躲亮源。团队改为在本地恒星时里挑三段天区：1.25到2.7小时、4.5到6.5小时、8.5到10.75小时，避开射电星系Fornax A和银河系反银心方向。第一天区前景功率最小，给出的上限也最好。

已发表结果。图14和图15画出每个频带、每个天区的柱坐标功率 $P(k_\perp,k_\parallel)$。虚线是地平线，实线是地平线加 $200\,\mathrm{ns}$，在 $z=7.9$ 对应 $\Delta k_\parallel=0.11\,h\,\mathrm{Mpc^{-1}}$，在 $z=10.4$ 对应 $0.10$。下排各图除以预期噪声；$k_\parallel>0.2\,h\,\mathrm{Mpc^{-1}}$ 处，这个比值的均方根在约3%以内等于一，说明高 $k_\parallel$ 处的数据主要是噪声。负的带功率显示为白色格子。

已发表结果，附带一个未解问题。地平线之外也出现了功率，主要在低 $k_\perp$ 处。论文认为其中一部分来自锥化窗的痕迹，但说这解释不了全部，并引用了一种推测：可能是天线之间残余的交叉耦合。第二频带、第二天区在缓冲区之外还有轻微超出，论文说原因不明，未被标记的弱射频干扰是一种可能。

看这两张图时检查三件事：虚线和实线各在哪里；在每个 $k_\perp$ 上，亮区是否留在实线下方；下排噪声比值从哪里开始偏离一。

图16：达到了多大的动态范围

已发表结果。图16不做缓冲区截断，直接按 $|k|$ 分壳平均，再除以 $k=0$ 处的功率。整条流程在前景峰值和高 $k$ 噪声底之间达到约 $10^{9}$ 的动态范围。HERA Collaboration 2022, §4.3, Fig. 16 论文称这是探测的必要条件，但还不充分。

波束色散藏在论文的哪里

v1正文从未提到“波束色散”这个词，得从别的说法里认出这个机制。

一般机制在方法文献里。不论源本身还是仪器响应，频率上的任何结构都会展宽时延核，把前景功率推到地平线以上。Parsons et al. 2012, §2.1 天线旁瓣随频率变化，所以远离指向中心的亮源危害最大。Furlanetto, Oh and Briggs 2006, §9.3

已发表结果。在HERA论文里，同样的物理以“泄漏”的面目出现。第4.3节指出，射频干扰和残余的仪器增益会把前景模式漏到更高的 $k_\parallel$，这也是要挑较暗天区的原因之一。第3.8.1节描述了脉冲星PSR J0742-2822的偏振辐射，法拉第旋转让它随频率振荡。它在线偏振里看得见，在赝斯托克斯I里不明显；论文指出两路馈源波束的不对称可能把它以被压低的程度漏进斯托克斯I。HERA Collaboration 2022, §3.8.1, Fig. 9

模型预测。实验室“随频率变化的天线会漏光”预设，在 $43.8\,\mathrm m$ 基线、$200\,\mathrm{ns}$ 缓冲下，使用宽度按 $1/\nu$ 变化的理想 $14\,\mathrm m$ 艾里斑波束。在分析使用的模式里，漏出的前景是氢功率的 $1.5\times10^{3}$ 倍；把波束固定为频带中心的形状，这个比值降到 $0.20$。这只说明在示意天空里这一机制有多大，并不是对HERA波束的判断。

未解问题。论文没有区分HERA越过地平线的泄漏里，有多少来自波束色散，有多少来自交叉耦合或增益。

第3.8节：压下仪器自身的回波

已发表结果。HERA一期的 $150\,\mathrm m$ 同轴电缆里有反射，子反射分布在 $100$ 到 $1200\,\mathrm{ns}$；还有一种交叉耦合系统误差，位于高时延，几乎不随时间变化，在正负时延上大致对称。团队在 $150$ 到 $1500\,\mathrm{ns}$ 之间拟合了 $25$ 个反射项，并对东西向投影不短于 $14\,\mathrm m$ 的基线滤除交叉耦合。两步合起来，在 $0.1<k<0.8\,h\,\mathrm{Mpc^{-1}}$ 范围内把仪器系统误差压低了约两个数量级。HERA Collaboration 2022, §3.8

读这一节要明白两点。第一，正是这一步让窗口区域能用；没有它，这些回波恰好落在测量氢的地方。第二，它也是信号损失的来源：滤掉高时延处不随时间变化的成分，也可能顺带滤掉一点天空信号，论文在第3.11节对此做了修正。

准备中的工作。论文说这一步在Aguirre等人的验证流程中做过检验，当时该工作尚在准备中。

第3.10、3.11节：流程削掉了多少信号

已发表结果。对近似冗余的基线做相干平均，信号可能部分相消。HERA用指标 $\Delta\chi$ 比较零时延处前景功率的相干平均与非相干平均，平均损失约1%。HERA Collaboration 2022, §3.10, Fig. 10

已发表结果。表2列出四项损失，先写第一频带，括号内为第二频带：绝对定标9%（8%），交叉耦合滤波3%（1%），时间平均1%（1%），冗余平均1%（1%）。这些损失都被当作与尺度和再电离模型无关，有的靠论证，有的经过检验，并且都已在最终功率谱中修正。HERA Collaboration 2022, §3.11, Table 2

准备中的工作。第5节概述了Aguirre等人近乎端到端的模拟，在 $k\ge0.192\,h\,\mathrm{Mpc^{-1}}$ 处恢复出一个真实感较强的模型；还有一次半盲测试，用一条独立流程分析同样的模拟数据，在 $k$ 大于约 $0.2$ 处结果一致。HERA Collaboration 2022, §5

模型预测。实验室说明了方法选择为什么重要：避开楔形能拿回注入信号的 $99.9\,\%$，减去三阶多项式只能拿回 $82\,\%$。

第4.4节与图17、图18：从避开前景到给出上限

已发表结果。把柱坐标功率按 $|k|$ 分壳平均时，$k_\parallel$ 低于地平线加缓冲的格点权重全部为零。图17给出结果，误差棒为 $2\sigma$，负的带功率置零。最好的上限来自第一天区：$z=7.9$、$k=0.192\,h\,\mathrm{Mpc^{-1}}$ 处为 $(30.76\,\mathrm{mK})^2$，$z=10.4$、$k=0.256$ 处为 $(95.74\,\mathrm{mK})^2$。HERA Collaboration 2022, §4.4

已发表结果。图18画出第二频带、第一天区的窗函数。第16到第84百分位的半宽为 $0.031\,h\,\mathrm{Mpc^{-1}}$；由于缓冲区，最低的窗函数在 $0.128$ 以下没有响应。HERA Collaboration 2022, Fig. 18

一个需要细读的地方。正文说带功率的采样间隔是“$\Delta k=0.64\,h\,\mathrm{Mpc^{-1}}$，原生间隔的两倍”。而公开文件和论文自己的表3、表4用的网格是 $0.064$，公开笔记本写的也是 $0.064$。本导读照录印出的数值，不做改动；本系列的方法指南直接使用文件。

另一个细节。第12式把归一化写成 $X^2Y/(\Omega_{pp}B)$，这是常见写法，其中 $X$ 把角度换成距离，$Y$ 把频率换成距离。可第13、14式又用 $c(1+z)^2/(\nu_{21}H)$ 定义 $X$，这是视线方向的换算，而把 $Y$ 定义为横向共动距离。看单位就能判断：只有角度的平方乘频率，才能把 $\mathrm{sr\,Hz}$ 变成 $h^{-3}\mathrm{Mpc^3}$；hera_pspec也正是用横向换算的平方乘以视线方向换算来计算这个因子。hera_pspec conversions.py

自己读论文时的检查清单

• 确认上限靠的是避开而不是扣除，并在图14、图15里找到 $200\,\mathrm{ns}$ 缓冲。
• 对每张图问一句：它展示的是测量、模型还是检验？
• 把论文测到的内容和它以“准备中”引用的内容分开。
• 即使没写“波束色散”，也要从泄漏、增益、反射和交叉耦合这些说法里找出仪器带来的谱结构。
• 把表2读作已经做过的修正，确认你引用的是修正后的上限。
• 比较模型时，用图18或公开文件里的窗函数，并用同一条规则 $U=\max(\hat\Delta^2,0)+2\sigma$。