捕捉宇宙重生的余晖：SPT 如何首次"听到"kSZ 信号¶
论文: An Improved Measurement of the Secondary Cosmic Microwave Background Anisotropies from the SPT-SZ + SPTpol Surveys arXiv: 2002.06197 ｜ 作者: Reichardt et al. ｜ 年份: 2020 阅读日期: 2026-03-16
一、一张 138 亿年前的照片，以及它没有告诉我们的事¶
宇宙微波背景（Cosmic Microwave Background, CMB）是大爆炸留给我们的一张"婴儿照"——大约在宇宙诞生后 38 万年，温度降到足以让电子和质子结合成中性氢的那一刻，光子第一次获得自由，从此直线飞行，穿越 138 亿年到达我们的望远镜。[补充]
但 CMB 并不只记录了那个瞬间。在随后的漫长旅途中，光子穿越星系团、穿过正在电离的氢气泡、被高温等离子体散射——这些"路上的遭遇"在 CMB 上留下了微弱的印记，被称为次级各向异性（secondary anisotropies）。[原文] 它们的角度尺度很小——大约几角分，对应多极矩 \(\ell \gtrsim 2000\)——远小于大尺度上那些著名的原初涟漪。[原文]
这些次级信号虽然微弱，却编码了晚期宇宙最重要的信息：结构增长的振幅和再电离（reionization）的持续时间。[原文]
二、两种散射，两个故事¶
在几角分尺度上，最重要的次级效应是两种 Sunyaev-Zel'dovich（SZ）效应，它们的共同点是：CMB 光子被自由电子散射。但散射的物理原因截然不同。[原文]
运动学 SZ 效应（kinematic SZ, kSZ）¶
想象一团自由电子以速度 \(v\) 整体运动。CMB 光子撞上这团电子后，会被多普勒频移——如果电子朝我们运动，光子变蓝；反之变红。信号强度正比于 \((v/c) \, n_e\)，其中 \(n_e\) 是电子数密度。[原文]
kSZ 有一个对测量者来说极为棘手的特征：它不改变 CMB 的黑体频谱。散射后的光子看起来还是"原来的 CMB"，只是温度微微变了。这意味着我们无法通过频率信息把它和原初 CMB 区分开——它像一个完美的伪装者，隐身在原初信号之中。[原文]
kSZ 功率谱有两个物理来源截然不同的贡献：[原文]
-
Patchy kSZ（不均匀再电离 kSZ）：宇宙再电离时期，紫外辐射从第一代恒星和星系出发，逐渐"吹出"电离泡（ionized bubbles）。泡内是电离气体，泡外还是中性氢。电离分数的巨大对比——从 0 跳到 1——结合电离泡的整体运动，产生强烈的 kSZ 信号。patchy kSZ 的振幅直接编码了再电离持续了多久。[原文 + 补充]
-
Homogeneous kSZ（均匀/晚期 kSZ）：再电离完成后，宇宙中所有气体都被电离了，电离分数的对比消失。但大尺度结构的密度涨落依然驱动着物质的整体流动（bulk flows），产生持续的 kSZ 信号。[原文 + 补充]
热 SZ 效应（thermal SZ, tSZ）¶
星系团内部充满了几千万到上亿度的高温气体。当 CMB 光子穿过这些热等离子体时，能量从热电子传递给冷光子——这是逆康普顿散射（inverse Compton scattering）。信号振幅正比于 \((k_B T_e / m_e c^2) \, n_e\)。[原文]
tSZ 效应有一个独特的频谱指纹：净能量传递意味着低频的光子被"搬"到高频。在 217 GHz 处信号恰好为零；低于 217 GHz 光子数减少（负信号），高于 217 GHz 光子数增加（正信号）。[原文] 这个频谱特征是一把"钥匙"——只要在多个频率上观测，就可以利用 tSZ 独特的频率缩放关系把它从其他信号中分离出来。[补充]
tSZ 功率对物质功率谱的归一化参数 \(\sigma_8\) 高度敏感——\(\sigma_8\) 描述的是今天（\(z=0\)）线性物质密度场在 \(8 \, h^{-1}\,\text{Mpc}\) 尺度上的均方根涨落。[原文]
三、拥挤的毫米波天空¶
如果只有 tSZ 和 kSZ，事情会简单很多。但在毫米波段的几角分尺度上，天空远比这拥挤。[原文]
-
尘埃星形成星系（Dusty Star-Forming Galaxies, DSFGs）：它们内部的尘埃被恒星加热后发出热辐射，在 150 和 220 GHz 尤其明亮。最亮的可以逐个探测和遮挡，但每平方角分内有大量暗弱的 DSFG 无法去除。[原文] DSFG 的信号可以分为随机分布的泊松项（Poisson）和空间聚集的成团项（clustering）。[原文]
-
射电星系（Radio galaxies / AGN）：活动星系核的同步辐射，在 95 GHz 更重要。[原文]
-
宇宙红外背景（Cosmic Infrared Background, CIB）：所有 DSFG 的叠加。[补充]
-
tSZ-CIB 反相关：暗物质过密区域既孕育了更多热气体（tSZ 为负信号，因为低于 217 GHz），也拥有更多尘埃星系（CIB 为正信号）。这导致 tSZ 和 CIB 之间存在反相关。[原文]
要从这团"混音"中提取出 kSZ 信号——一个频谱和原初 CMB 完全相同、振幅只有几个 \(\mu\text{K}^2\) 的成分——无异于在管弦乐队的演奏中辨认出一把小提琴的泛音。[补充]
四、以前的困境：kSZ 被 tSZ 挡住了¶
2015 年，George et al.（G15）用 SPT-SZ 巡天的 2540 平方度数据（2540 deg²）进行了当时最精确的小尺度 CMB 功率谱分析。他们在 95、150、220 GHz 三个频率上测量了 6 个频率组合的功率谱，用包含多种前景成分的模型去拟合。[原文]
G15 在不使用外部先验的情况下，kSZ 的约束非常弱。原因在于一个关键的简并（degeneracy）：tSZ 和 kSZ 都在 \(\ell \sim 3000\) 附近贡献几个 \(\mu\text{K}^2\) 的功率。虽然 tSZ 有频率依赖而 kSZ 没有，但 SPT-SZ 在 95 GHz 的噪声太大，无法精确地"钉住"tSZ 的振幅。于是 tSZ 的不确定性直接传递给了 kSZ。[原文 + 补充]
为了打破这个简并，G15 不得不引入外部信息——tSZ 三点函数（bispectrum）的先验。bispectrum 主要由 tSZ 贡献（因为 kSZ 近似高斯分布，对 bispectrum 贡献很小），所以它可以独立约束 tSZ 振幅。加上这个先验后，G15 得到 \(D_{3000}^{\text{kSZ}} = 2.9 \pm 1.3 \, \mu\text{K}^2\)——只有约 \(2\sigma\) 的显著性，不算真正的探测。[原文]
ACT 合作组的结果更弱：\(D_{3000}^{\text{kSZ}} < 8.6 \, \mu\text{K}^2\)（95% CL），只给出了上限。[原文]
kSZ 仍然是一个未被确凿探测到的信号。
五、杠杆点：95 GHz 的灵敏度革命¶
打破困局的关键洞察是：95 GHz 是分离 tSZ 的最佳杠杆点。[补充]
tSZ 效应的频率依赖使得它在 95 GHz 的功率比 150 GHz 高出一倍多——准确地说，SPT 的 95 GHz 波段与 150 GHz 波段的 tSZ 功率比为 2.77。[原文] 而 220 GHz 波段恰好对准了 tSZ 频谱的零点（约 217 GHz），tSZ 功率几乎为零。[原文]
这意味着：如果能在 95 GHz 获得足够低的噪声，tSZ 的振幅就会被 95 GHz 数据"钉住"。一旦 tSZ 确定了，剩余的、在所有频率上都一样的那部分小尺度功率，就自然归因于 kSZ。[补充]
2013–2016 年间，SPT 换装了新一代相机 SPTpol，对 500 deg² 的天区进行了深度观测。SPTpol 在 95 和 150 GHz 的噪声显著低于 SPT-SZ。[原文] 虽然 SPTpol 没有 220 GHz 通道，但这并不影响 tSZ 的约束——因为 tSZ 在 220 GHz 几乎没有信号。[原文]
Reichardt et al. 的关键一步是：把 SPTpol 500 deg² 的低噪声数据与 SPT-SZ 的 2540 deg² 数据合并。在两个巡天重叠的 500 deg² 天区内，他们将时间序列数据（time-ordered data, TOD）按逆方差加权 coadd 成一张图。[原文]
合并后，SPTpol 新数据在全部 2540 deg² 中贡献的统计权重为： - 95 GHz: 83% - 150 GHz: 44% - 220 GHz: 0%（SPTpol 没有 220 GHz）[原文]
换句话说，在 95 GHz，新数据几乎完全主宰了信噪比。这带来了 95×95 GHz bandpower 误差降低约 3 倍的巨大改进，150×150 GHz 也改善了约 30%。[原文]
六、从地图到功率谱：数据分析流水线¶
功率谱估计¶
作者使用 pseudo-\(C_\ell\) 交叉谱方法。[原文] 其核心思想是：
- 给每张地图加一个边缘平滑衰减到零的窗函数，同时遮挡 150 GHz 通量 \(> 6.4\) mJy 的亮点源。[原文]
- 对窗函数化的地图做傅里叶变换，得到 \(\tilde{m}^{\nu_i, A}_\ell\)。[原文]
- 在 \(\ell\)-bin \(b\) 内取加权平均，计算交叉谱：
其中 \(A, B\) 是不同的观测子集，\(\nu_i, \nu_j\) 是频率。[原文]
使用交叉谱（\(A \neq B\)）而非自谱的好处是：噪声在不同观测之间不相关，交叉谱自动消除噪声偏差（noise bias）。[原文]
校准¶
通过将 SPT 地图与 Planck 2018 的 CMB 地图在 \(600 \leq \ell \leq 1500\) 范围内比较来确定绝对校准。最终功率的校准不确定度分别为 0.33%（95 GHz）、0.18%（150 GHz）、0.42%（220 GHz）。[原文]
Beam 和 Transfer function¶
三个频率的 beam 主瓣近似高斯，FWHM 分别为 1.7′、1.2′、1.0′。[原文] 通过 200 组信号注入模拟计算 map-making 和滤波的 transfer function，同时估计样本方差（sample variance）。[原文]
系统误差检验¶
通过扫描方向分裂（scan-direction split）和时间分裂（time split）做 null tests。150 GHz 在 \(\ell < 2500\) 有极小的余量功率（\(<0.1\%\)），统计显著但物理上可忽略，对 \(\ell > 2000\) 的分析无影响。[原文]
最终数据：6 个频率组合（95×95, 95×150, 95×220, 150×150, 150×220, 220×220），在 \(2000 \leq \ell \leq 11000\) 范围内共 88 个 bandpowers。[原文]
七、十个参数的拼图：前景建模¶
固定 \(\Lambda\)CDM 宇宙学参数（取 Planck best-fit 值），用 10 个自由参数描述前景和次级各向异性：[原文]
| 分量 | 参数 |
|---|---|
| tSZ | 振幅 \(D_{3000}^{\text{tSZ}}\)（Shaw 模板） |
| kSZ | 总振幅 \(D_{3000}^{\text{kSZ}}\)（CSF homogeneous + patchy 模板之和） |
| CIB Poisson | 振幅 + 谱指数 \(\beta_p\) |
| CIB 1-halo clustered | 振幅 |
| CIB 2-halo clustered | 振幅 + 谱指数 \(\beta_c\)（1-halo 和 2-halo 共享） |
| Radio galaxies | 泊松振幅 + 谱指数 \(\alpha_{rg}\) |
| tSZ-CIB 相关 | 相关系数 \(\xi\)（自由浮动于 \([-1, 1]\)） |
每个分量只拟合振幅（和频率缩放参数），功率谱形状由理论模板给定。[原文] 作者验证了更换不同 tSZ 模板（Shaw、Battaglia、Bhattacharya、Sehgal）和不同 kSZ 模板（CSF homogeneous、patchy、两者之和）对结果影响很小。[原文]
八、打开信封：首次 ≥3σ 探测 kSZ¶
用 CosmoMC 对 88 个 bandpowers 进行 MCMC 拟合，基线模型的拟合质量为 \(\chi^2 = 99.7\)，自由度 78，PTE = 5.0%——可以接受。[原文]
核心结果：
[原文]
这是一个里程碑：kSZ 从"上限"进入了"测量"时代。[原文 + 补充]
与 G15 的比较：G15 加上 bispectrum 先验后得到 \(D_{3000}^{\text{kSZ}} = 2.9 \pm 1.3 \, \mu\text{K}^2\)，中心值极为一致，但本文的误差小了 30%——而且不需要任何外部先验。[原文] 如果在本文结果上再加上 bispectrum 先验，kSZ 变为 \(2.8 \pm 0.9 \, \mu\text{K}^2\)，即 \(3.1\sigma\) 探测。[原文]
tSZ 的测量与 G15（\(4.38^{+0.83}_{-1.04} \, \mu\text{K}^2\)）、R12（\(4.20 \pm 1.37 \, \mu\text{K}^2\)）和 ACT（\(3.9 \pm 1.7 \, \mu\text{K}^2\)）都在 \(1\sigma\) 内一致，但误差显著缩小。[原文]
tSZ-CIB 相关¶
\(\xi = 0.076 \pm 0.040\)，数据以 98.3% 的置信度排除 \(\xi < 0\)——确认了 DSFG 在星系团中是过密的（positive correlation）。[原文] 增大 \(\xi\) 会缓慢降低推断的 tSZ 功率、同时快速增加推断的 kSZ 功率，因此 tSZ、kSZ 和 \(\xi\) 三者之间存在非平凡的简并。[原文]
为什么 kSZ 终于被探测到了？¶
因果链条很清晰：[原文 + 补充]
- SPTpol 在 95 GHz 贡献了 83% 的统计权重 → 95×95 GHz bandpower 误差降低约 3 倍
- 95 GHz 是 tSZ 信号最强的频率（tSZ 功率比 = 2.77）→ tSZ 振幅被更精确地确定
- tSZ 确定后 → tSZ-kSZ 简并被打破
- kSZ 是唯一剩余的、在所有频率上振幅相同的小尺度功率 → kSZ 被分离出来
九、从功率到物理：再电离持续了多久？¶
kSZ 功率是 homogeneous 和 patchy 两部分之和：
当前数据无法独立分离这两个分量。[原文] 因此作者采用理论模板估计 homogeneous kSZ 的贡献：基于 Shaw et al. (2012) 和 Calabrese et al. (2014) 的公式，在本文宇宙学参数下，\(D_{3000}^{\text{h-kSZ}} \approx 1.65 \, \mu\text{K}^2\)。[原文] 同时还考虑了将其乘以 0.75 和 1.25 的高低估计。[原文]
减去 homogeneous 部分后，剩余的 patchy kSZ 功率 \(D_{3000}^{\text{p-kSZ}} < 2.9 \, \mu\text{K}^2\)（95% CL，取 best estimate 的 homogeneous kSZ）。[原文]
patchy kSZ 功率与再电离的持续时间 \(\Delta z_{\text{re}}\) 之间有一个半解析公式（Calabrese et al. 2014, Eq. 6）：
其中 \(z_{re}\) 是电离分数达到 50% 的红移，\(\Delta z_{re}\) 是从 25% 电离到 75% 电离所经历的红移间隔。[原文]
利用这个公式，可以把 patchy kSZ 功率的约束转化为对 \(\Delta z_{re}\) 的约束：
| 条件 | \(\Delta z_{re}\) 上限（95% CL） |
|---|---|
| 仅功率谱，best-estimate h-kSZ | \(< 5.4\) |
| 加 bispectrum 先验 | \(< 4.1\) |
| 68% 置信区间（含 bispectrum） | \(1.1^{+1.6}_{-0.7}\) |
[原文]
作为对比，G15 在使用相同的定义和模型后（作者做了细致的重新分析以确保可比性），相应的上限为 \(\Delta z_{re} < 8.5\)——本文的约束改善了近一倍。[原文]
\(\Delta z_{re} \sim 1\) 意味着什么？它意味着宇宙从 25% 电离到 75% 电离只用了大约 \(\Delta z \sim 1\) 的时间，大致对应宇宙年龄从约 5 亿年到约 7 亿年之间——再电离发生得相当快（fairly quickly）。[原文] 这与高红移 Lyman-\(\alpha\) 森林和星系紫外光度函数的独立推断一致。[原文 + 补充]
十、其他前景的"体检报告"¶
虽然 kSZ 是这篇论文的明星，但 10 参数拟合同时给出了所有前景的精确约束：
- CIB Poisson 功率在 150 GHz 为 \(7.24 \pm 0.63 \, \mu\text{K}^2\)，在 220 GHz 为 \(61.4 \pm 9.0 \, \mu\text{K}^2\)，与 G15 在 \(0.3\sigma\) 内一致。[原文]
- CIB 聚集项（1-halo + 2-halo）的功率和谱指数也与 G15 一致。[原文] 作者注意到 G15 的 CIB 聚集功率略低于更早的测量，而本文的值向那些更早的结果靠拢。[原文]
- 射电星系功率 \(D_{3000}^r = 1.01 \pm 0.17 \, \mu\text{K}^2\)（150 GHz），比 De Zotti et al. (2005) 模型预测的 \(1.28 \, \mu\text{K}^2\) 低约 25%。谱指数 \(\alpha_{rg} = -0.76 \pm 0.15\)。[原文]
这些结果证明 10 参数模型能自洽地描述所有 6 个频率组合的数据。[原文]
十一、故事的尾声与下一章¶
这篇论文的核心贡献可以用一句话概括：通过加入 SPTpol 的低噪声 95 GHz 数据，首次以 \(\geq 3\sigma\) 的显著性探测到 kSZ 功率谱，并由此约束了宇宙再电离的持续时间。[原文]
但故事远未结束。当前约束仍然相当宽——\(\Delta z_{re} = 1.1^{+1.6}_{-0.7}\) 的 68% 置信区间意味着 kSZ 还不能精确地区分不同的再电离模型。Homogeneous 和 patchy kSZ 目前无法在数据中分离，只能依赖理论模板的估计来间接推断 patchy 部分。[原文]
下一步已在路上。论文发表时，SPT-3G 相机正在进行 1500 deg² 的深度巡天，预期最终温度噪声水平为 3、2、9 \(\mu\text{K-arcmin}\)（95、150、220 GHz）[原文]——比本文的数据再低一个数量级。更远的将来，Simons Observatory 和 CMB-S4 将把测量扩展到更大天区、更低噪声和更多频段。[原文]
当这些下一代实验的数据到来时，kSZ 将不再只是一个"\(3\sigma\) 的探测"，而是一扇透视再电离历史的明亮窗户——我们终将能回答：宇宙是如何在黑暗时代结束时重新点亮的？[补充]
附：关键数字速查¶
| 量 | 值 | 来源 |
|---|---|---|
| \(D_{3000}^{\text{tSZ}}\)（143 GHz） | \(3.42 \pm 0.54 \, \mu\text{K}^2\) | §6 baseline |
| \(D_{3000}^{\text{kSZ}}\) | \(3.0 \pm 1.0 \, \mu\text{K}^2\) | §6 baseline |
| kSZ 探测显著性 | \(\geq 3\sigma\)（首次） | §6 |
| tSZ-CIB 相关 \(\xi\) | \(0.076 \pm 0.040\) | §6 |
| \(\xi < 0\) 被排除 | 98.3% CL | §6 |
| \(D_{3000}^{\text{h-kSZ}}\)（模板值） | \(\approx 1.65 \, \mu\text{K}^2\) | §7 |
| \(D_{3000}^{\text{p-kSZ}}\) 上限 | \(< 2.9 \, \mu\text{K}^2\)（95% CL） | §7 |
| \(\Delta z_{re}\) 上限（含 bispectrum） | \(< 4.1\)（95% CL） | §7 |
| \(\Delta z_{re}\) 68% CI（含 bispectrum） | \(1.1^{+1.6}_{-0.7}\) | §7 |
| Baseline \(\chi^2\) / dof | 99.7 / 78, PTE = 5.0% | §6 |
| 95 GHz tSZ/150 GHz tSZ 功率比 | 2.77 | §5.1 |
| SPTpol 在 95 GHz 的统计权重 | 83% | §2.1 |
| Beam FWHM（95/150/220 GHz） | 1.7′ / 1.2′ / 1.0′ | §2.2 |
| 校准不确定度（95/150/220 GHz） | 0.33% / 0.18% / 0.42% | §2.2 |