CMBComp: A Simple and Accurate Compressed CMB Likelihood for Dark Energy, Curvature, and Massive Neutrinos

技术摘要：“CMBComp：一种用于暗能量、曲率及质量中微子的简单且精确的压缩 CMB 似然函数”

问题陈述
完整的宇宙微波背景（CMB）似然函数在计算上非常昂贵，且结构复杂，这使得它难以部署在探索性分析中，特别是对于晚期现象学宇宙模型。这种局限性阻碍了在大型模型空间中的快速扫描，以及与低红移几何探测器（如重子声学振荡 BAO）的联合分析。虽然对于晚期推断而言，主导性的 CMB 信息在概念上被理解为编码在少量的几何量中（声波视界角尺度、最后散射面的共动距离以及物理密度参数），但目前仍需要一种实用的、针对模型优化的且准确的压缩似密函数，以在不牺牲精度的前提下替代完整的似然函数。

方法论
作者提出了 CMBComp，这是一套紧凑的、高斯压缩的 CMB 似然函数，其源自组合的 SPA 数据集（SPT-3G D1 + ACT DR6 + Planck PR3 原初各向异性 + Planck PR4/NPIPE CMB 引力透镜）。该压缩策略依赖于这样一个事实：对于标准的平滑暗能量和几何驱动扩展模型，CMB 主要作为一个高红移几何锚点。

该框架构建了五种不同的压缩似然函数，每一种都针对特定的宇宙学模型空间。核心几何观测量是移动参数 $R$ 和声学角尺度 $\ell_a$。根据模型的不同，这些量会通过增加一个或两个物理密度参数来构成数据向量 $\mathbf{v}$。似然函数定义为：
$$-2 \ln \mathcal{L}_{\text{comp}} = (\mathbf{v}_{\text{th}} - \bar{\mathbf{v}})^\top \mathbf{C}^{-1} (\mathbf{v}_{\text{th}} - \bar{\mathbf{v}})$$
其中 $\bar{\mathbf{v}}$ 和 $\mathbf{C}$ 是从特定模型的 SPA 链中提取的平均向量和协方差矩阵，而 $\mathbf{v}_{\text{th}}$ 是通过 CAMB 计算的理论预测值。

五种具体的压缩方式如下：

1. CMB-3： 一个三参数压缩 $(R, \ell_a, \Omega_b h^2)$，适用于空间平坦的 $\Lambda$CDM 模型。
2. CMB-3w： 将相同的三参数向量应用于动力学暗能量 $w_0 w_a$CDM 模型，其均值和协方差是从 $w_0 w_a$CDM 链中重新提取的，以适应扩展后的参数空间。
3. CMB-4$\nu$： 一个四参数压缩 $(R, \ell_a, \Omega_b h^2, \Omega_c h^2)$，适用于质量中微子 $\nu\Lambda$CDM 模型。引入 $\Omega_c h^2$ 对于正确追踪物质预算至关重要，因为当中微子质量总和 ($\sum m_\nu$) 变化时，这可以避免将质量中微子错误地视为复合时期的冷物质。
4. CMB-4k： 一个四参数压缩 $(R, \ell_a, \Omega_b h^2, \Omega_\Lambda)$，适用于非平坦的 $o\Lambda$CDM 模型。至关重要的是，对于该模型，最后散射红移 $z_*$ 是在每个 MCMC 步骤中从 CAMB 自洽评估的，而不是固定的，以解释由曲率引起的声波视界映射偏移。
5. CMB-4$\nu w$： 一个四参数压缩 $(R, \ell_a, \Omega_b h^2, \Omega_c h^2)$，适用于联合质量中微子与动力学暗能量的 $\nu w_0 w_a$CDM 模型，其统计量从相应的联合链中重新提取。

主要贡献

• 模型适配的似然函数： 本文提供了五种特定的、可直接使用的压缩似然函数（均值向量和协方差矩阵），这些函数源自最新一代的 CMB 数据（SPA）。
• 实现独立性： 作者提供了一个独立的、轻量级的 Python 实现（CMBComp），它不需要完整的玻尔兹曼求解器（如 CAMB）或复杂的推断框架（如 Cobaya 或 CosmoMC）即可运行。它包含一个自定义的背景演化计算器和一个自适应 Metropolis–Hastings 采样器。
• 验证： 通过将压缩似然函数与 DESI DR2 BAO 似然函数结合，并将所得的边缘化后验与使用完整 SPA CMB 似然函数生成的参考链进行对比，对压缩似然函数进行了严格验证。

结果
验证表明，在所有五种模型类别中，压缩似然函数都能高精度地重现完整 CMB 分析的边缘化约束和二维轮廓形状：

• $\Lambda$CDM： 与完整 CMB 相比，压缩后验在 $H_0$ 和 $\Omega_m$ 上显示出微小的残余偏移（$\sim 0.6\sigma$）。作者将其归因于三参数压缩无法完全模拟 CMB 引力透镜效应，而引力透镜效应倾向于在完整似然函数中将哈勃常数向下移动。
• $w_0 w_a$CDM： 一致性极佳。额外的暗能量参数吸收了在 $\Lambda$CDM 中观察到的偏差，使得所有参数的偏移都变得微乎其微（$< 0.3\sigma$）。
• $\nu\Lambda$CDM： 几何参数（$\Omega_m, H_0$）的重现误差 $< 0.12\sigma$。中微子质量上限比完整 CMB 参考值略弱（$\sim 13\%$），这反映了通过间接物质预算关系对 $\sum m_\nu$ 的约束力有所减弱，但定性结构得到了保留。
• $o\Lambda$CDM： 一致性几乎是精确的，所有均值偏移均低于 $0.13\sigma$。自洽计算 $z_*$ 以及引入 $\Omega_\Lambda$ 消除了在 $\Lambda$CDM 案例中观察到的几何泄漏。
• $\nu w_0 w_a$CDM： 暗能量和几何参数的一致性在 $0.09\sigma$ 以内。由于暗能量参数开启了额外的简并方向，中微子质量约束比参考值稍宽（$\sim 14\%$），但后验结构得到了忠实重现。

意义与主张
本文声称，CMBComp 特别适用于修改晚期膨胀历史的模型进行宇宙学推断。它能够以极低的计算开销将准确的 CMB 约束纳入新的分析中，且无需依赖基于完整玻尔兹曼求解器的流水线。

作者强调，这些压缩函数对于那些效应局限于晚期膨胀历史和背景演化的模型是有效的，这些模型保留了标准的早期宇宙物理学。他们明确警告，对于引入了新早期物理、修改了复合过程或显著改变了原初扰动部门的模型，仍应使用完整的 CMB 似然函数。这项工作为研究人员、学生和教育工作者提供了一个实用的、易于获取的工具，用于执行高精度、简化、透明且计算高效的宇宙学推断。