SPT-3G D1: A Measurement of Secondary Cosmic Microwave Background Anisotropy Power — 五问总纲¶
- arXiv: 2601.20551 (Chaubal, Huang, Reichardt et al., SPT-3G Collaboration, 2026 投稿 JCAP)
- 数据: SPT-3G D1（2019–2020 冬季）95 / 150 / 220 GHz, 1646 deg² 南天
- 测量: 6 个频率自/互谱 \(1700 \le \ell \le 11{,}000\)
- 目标: 在 arcmin 尺度限制次级 CMB（tSZ、kSZ）、CIB、射电星系功率，并由 kSZ 推 reionization 时长
1. 问题是什么（What problem）¶
一句话版：用 SPT-3G 两年数据把 \(\ell\sim 3000\) 处毫米波前景拆得更细，重新校准 tSZ/kSZ 功率并据此约束再电离时长——结果发现"答案高度依赖你假设的模板"。
领域原来有什么问题¶
毫米波小角尺度功率谱里同时叠着多个分量（透镜化主 CMB、tSZ、kSZ、CIB Poisson + clustering、tSZ–CIB 相关、射电星系、银河 cirrus）。前人在 \(\ell=3000\) 处通过单一振幅参数加 simulation-based 模板做联合拟合，给出 \(D_{3000}^{\rm tSZ}\)、\(D_{3000}^{\rm kSZ}\)，再用 kSZ 推再电离时长 \(\Delta z_{\rm re}\)。
为什么重要¶
- kSZ at \(\ell=3000\) 是除 CMB optical depth \(\tau\) 之外几乎唯一能直接限制 reionization "持续多久"的观测量。
- tSZ amplitude 对 \(\sigma_8\) 极敏感（\(\propto \sigma_8^{7\sim8}\)），是大尺度结构生长的"放大器"。
- tSZ–CIB correlation 既是天体物理信号本身，也是把 kSZ/tSZ 分开来的关键 nuisance。
以前的方法为什么不够¶
- SPT-SZ + SPTpol（R21）和 ACT DR6（Louis25, Beringue25）单参数模板拟合 \(D_{3000}\)，结果差异在误差棒内尚可接受。
- SPT-3G 把 \(\ell\sim3000\) 处 95×95、150×150、220×220 GHz 的不确定性分别压缩了 2.0/1.4/2.9 倍，到这一精度时不同模拟模板（Agora vs G15/Shaw+Zahn）给出的 \(D_{3000}^{\rm tSZ}\)、\(D_{3000}^{\rm kSZ}\) 出现 \(\sim 5\sigma\) 差异——模板本身成了系统误差主导。
2. 核心想法（Why it works）¶
这个方法之所以有效，是因为 数据精度足以反过来约束 SZ 和 CIB 在 \(\ell\) 上的形状：把以前固定的"模板形状"开放成 power-law 修正 (\(\ell^\alpha\)) 或 monotonic cubic spline，让数据自己告诉我们 tSZ-CIB correlation 怎么随 \(\ell\) 变。
关键洞察¶
- tSZ 与 kSZ 在不同 \(\ell\) 上的"形状"差异，比 \(\ell=3000\) 处的振幅更难造假——所以引入 \(\ell^{\alpha_{\rm tSZ}}\)、\(\ell^{\alpha_{\rm kSZ}}\) 两个自由度可以吸掉绝大部分模板差异。
- 60% 的 SZ 模板依赖性来自 tSZ–CIB 相关函数 \(\xi(\ell)\) 形状的不同（Z12 形态随 \(\ell\) 上升 vs Agora 在 \(\ell\sim 2000\) 见峰后下降）。把 \(\xi(\ell)\) 也自由化，反而约束变更稳。
- 三档模型梯子（Agora 模板 → free CIB + \(\ell^\alpha\) SZ → free CIB+SZ）：从最受先验束缚到最 agnostic。三者结果在 \(\ell\sim3000\) 处一致 (\(\lesssim 1\sigma\))，差异主要出现在 \(\ell\) 边缘。
与现有方法的本质区别¶
| 旧方法（R21、ACT 旧分析） | 本文 | |
|---|---|---|
| tSZ/kSZ 形状 | 固定模拟模板 | 模板 + \(\ell^\alpha\) 修正 / 7 节点 spline |
| tSZ–CIB 相关 | 单参数缩放 Z12 模板 | 8 节点 spline，\(\xi(\ell)\ge 0\) |
| CIB clustering | 1-halo + 2-halo 模板 | 5 节点 spline，\(\beta_C(\ell)\) 也 spline |
| 模型选择敏感度 | \(\le 1\sigma\)（在 SPT-SZ 精度下） | \(\sim 5\sigma\)（必须开放形状） |
3. 技术实现（How it is done）¶
方法流水线¶
problem → observation → idea → method → result：
- Problem：SPT-SZ 时代的单模板单振幅 SZ 拟合在 SPT-3G 精度下崩了。
- Observation：模板间差异主要在 \(\ell\) 形状，不在 \(\ell=3000\) 处的总功率。
- Idea：让数据约束 \(\ell\) 形状。
- Method：pseudo-\(C_\ell\) 估计 → 200 bundle 交叉谱 → mode-coupling + transfer function → 多模板 / power-law / spline 三档前景模型 → cobaya MCMC 联合拟合 → 后处理出 \(\Delta z_{\rm re}\)。
- Result：\(\dtsz=4.91\pm0.37\)、\(\dksz=1.75\pm0.86\,\mu{\rm K}^2\)（free CIB + \(\ell^\alpha\) G15 SZ）；推得 \(\Delta^{50}z_{\rm re}<3.8\) (Calabrese14) 或 \(\Delta^{90}z_{\rm re}<6.1\) (AMBER, Kramer25)，95% C.L.
解析部分¶
数据在 200 个时间 bundle 之间做 pseudo-\(C_\ell\) 互相关 (Hivon02 + Tristram05)，掩模 W 由 150 GHz 权重图加 6 mJy 点源掩蔽构成；mode-coupling 矩阵用 NaMaster 算。Bandpower 反卷核 \(K_{bb'} = P_{b\ell}(M_{\ell\ell'}[\mathbf{W}]\,F_{\ell'}\,B_\ell^i B_\ell^j) Q_{\ell'b'}\) 同时吸收 mask、滤波传递函数 \(F_\ell\)、波束 \(B_\ell\)、像素 window。波束由场内 AGN (\(r\lesssim 1.\!'5\)) 和土星 (\(r\gtrsim 3'\)) 拼接，按四种 SED 分别构造。
数值部分¶
- 协方差：100 套高斯 + Poisson 模拟提取 sample variance；用 \(A+B/(\ell-500)\) 拟合 \(\nu_b^{-1}\)；噪声项从 bundle cross-spectra 取，\(\ell<6000\) 允许频率间相关，\(\ell>6000\) 强制对角。
- MCMC：cobaya 跑各模型链，Planck 2018 \(\Lambda\)CDM 高斯先验 (
plikHM_TTTEEE_lowl_lowE)；备选 CMB-SPA / no-tau 先验做稳健性检查。 - 再电离后处理：从 \(\dksz\) 减去 \(\dhksz=1.65(\sigma_8/0.8)^{4.46}\)（Calabrese14 Eq.5）得 patchy kSZ；再用 patchy kSZ 公式（Calabrese14 Eq.6 或 Kramer25 AMBER fit）反推 \(\Delta z_{\rm re}\)；用 uniform \(\Delta z\) 先验重采样去除 prior volume effect。
关键近似和假设¶
- 非相对论 tSZ 频谱（数据不限制电子温度）。
- \(z_{\rm re}=7.68\)（Planck 中位）固定后才解 \(\Delta z_{\rm re}\)。
- 95 GHz CIB 由 150/220 + 单 \(\beta\) SED 外推，不是独立测量。
- \(\xi(\ell)\ge 0\)（spline 节点取非负）。
- 5%–95% 与 25%–75% 两种 \(\Delta z_{\rm re}\) 定义不可直接比对。
- 拟合质量 borderline：模板 fit PTE 极低；free CIB / free SZ PTE 仅 6%——可能是协方差/波束被低估，或模型仍漏分量。
4. 证据（What evidence）¶
最重要的 3 张图¶
- Fig.~3 (
tsz-ell-lots/ksz-ell-lots)：tSZ、kSZ 功率谱形状在三种模型下的中位线 + 68% CL 带；\(\ell\sim3000\) 处稳定，但小角尺度上模型差异显现。 - Fig.~4 (
sztriangleb)：\(\dksz\)、\(\dtsz\)、\(\xi_{3000}\) 三角形 posterior；展示 tSZ–kSZ–tSZ-CIB 强简并 + 三模型一致。 - Fig.~7 (
dz1d_50/dz1d_90)：再电离时长 1D posterior；柔性模型下偏好"快速再电离"，模板模型偏好慢。
关键数字¶
| 量 | 值 | 模型 |
|---|---|---|
| \(\dtsz\) at 143 GHz | \(4.91\pm 0.37\,\mu{\rm K}^2\) | free CIB + \(\ell^\alpha\) G15 SZ（baseline 报告值） |
| \(\dksz\) | \(1.75\pm 0.86\,\mu{\rm K}^2\) | 同上 |
| \(\alpha_{\rm tSZ}\), \(\alpha_{\rm kSZ}\) | \(0.37\pm 0.15\), \(0.0\pm 1.2\) | 同上 |
| \(\dtsz\) (Agora 模板) | \(4.28\pm 0.37\) | Agora templates |
| \(\dksz\) (Agora 模板) | \(3.96\pm 0.82\) | Agora templates |
| \(\dtsz\) (R21 G15 模板) | \(5.47\pm 0.20\) | G15 templates（baseline 对照） |
| \(\dksz\) (R21 G15 模板) | \(<1.4\) (95% C.L.) | G15 templates |
| \(\xi_{3000}\) | \(0.091\pm 0.025\) (Agora) → \(0.036\pm 0.021\) (\(\ell^\alpha\)) | 模板依赖 |
| \(\Delta^{50}z_{\rm re}\) | \(0.67^{+1.34}_{-0.51}\)，95% 上限 \(<3.8\) | Calabrese14, \(\sigma_8=0.812\) |
| \(\Delta^{90}z_{\rm re}\) | \(1.5^{+2.2}_{-1.1}\)，95% 上限 \(<6.1\) | AMBER (Kramer25), \(\sigma_8=0.812\) |
| 不确定性压缩 vs R21 (\(\ell\sim3000\)) | 95×95: 2.0×, 150×150: 1.4×, 220×220: 2.9× | — |
| 总 CIB power at \(\ell=3000\) | \(\sim 10\,\mu{\rm K}^2\) (150 GHz), \(\sim 108\,\mu{\rm K}^2\) (220 GHz) | 各模型 \(<1\sigma\) 一致 |
5. 结论与影响¶
核心结论（≤3 条）¶
- SPT-3G D1 给出迄今最精的 \(\ell\gtrsim 3000\) 毫米波温度功率谱（6 个 cross-spectra），在 \(\ell\sim 3000\) 处不确定性较 SPT-SZ/SPTpol 减小 1.4–2.9×。
- 次级 CMB 拆分严重依赖 tSZ–CIB 相关性的形状假设：单参数模板 fit 在不同模拟下给出 \(\sim 5\sigma\) 不同的 \(\dksz\)、\(\dtsz\)；引入 \(\ell^\alpha\) 或 spline 后回到一致。Agora-flavor 的"\(\xi(\ell)\) 随 \(\ell\) 下降"形状被数据偏好。
- 再电离倾向短：在 \(\sigma_8=0.812\) 下，free CIB + \(\ell^\alpha\) G15 SZ 给 \(\Delta^{50}z_{\rm re}<3.8\)；与 Raghunathan24 的 kSZ trispectrum 结果 (\(\Delta^{50}z_{\rm re}<4.5\)) 自洽。Agora 模板 fit 偏好更长的 \(\Delta z_{\rm re}\approx 7\)。
对领域的影响¶
- 把"前景模板的 \(\ell\)-dependence"从 nuisance 升格为数据可约束的物理对象——下一代 CMB（SPT-3G+, SO, CMB-S4, CMB-HD）必须在这个层面建模。
- 给 reionization 提供一条独立于 \(\tau\)、独立于 high-\(z\) Lyman-α 的时间尺度限制；与 21 cm/Lyman 联合分析 (e.g., arXiv:2504.09725) 互补。
- tSZ–CIB 相关性的"\(\ell\) 上升 vs \(\ell\) 下降"问题成为 CIB 半解析模型（Z12, Shang12 vs Agora/Omori24）的判据。
局限性和开放问题¶
- 拟合质量 borderline（PTE 低至 ~6%）：可能是波束/协方差被低估，或模型漏了 CO 线、辐射星系-tSZ 相关等成分。
- 再电离时长仅用 \(\ell=3000\) 处单点 kSZ，未充分利用 spline kSZ 在不同 \(\ell\) 上的信息——文中明言"超出本文范围"。
- 95 GHz CIB 不是独立约束（SED 推导）。
- \(\sigma_8\) 取值在 0.77–0.812 之间会让 \(\Delta z_{\rm re}\) 上限改变 ~25%。
- \(z_{\rm re}=7.68\) 固定值用于 \(\Delta z_{\rm re}\) 后处理，未联合放开。
6.（进阶）如果我是作者，我怎么想到这个 idea？¶
problem → observation → insight → method 的思维链：
- problem：SPT-3G D1 误差棒比 SPT-SZ 小一半以上，但单参数模板拟合仍给"漂亮的"\(\dtsz/\dksz\)——一旦换 Agora 模板，结果就跳 \(\sim 5\sigma\)。
- observation：模板差异在不同 \(\ell\) 上不是均匀的——Z12 的 tSZ–CIB 上升 vs Agora 下降；Shaw vs Agora tSZ 形状也不一样。
- insight：所以"模板假设"实际上是把 \(\ell\)-shape 锁死成先验，把所有 \(\ell\) 上的 mismatch 强行塞回 \(\ell=3000\) 振幅。如果让 \(\ell\)-shape 自由，振幅自然会回到正确值。
- method：先用最便宜的两参数 power-law 修正吃掉 80% 的差异 → 再用 monotonic cubic spline 提供 model-agnostic 上限 → 最后比较 PTE 与 posterior 一致性证明 spline 不是过拟合。
副产物：spline tSZ–CIB correlation 反而成了 Agora 模拟的判据，不只是 nuisance。