Consistency of standard cosmologies using Bayesian model comparison and tension quantification

技术摘要：使用贝叶斯模型比较与张力量化评估标准宇宙学的一致性

问题陈述
标准宇宙学模型 $\Lambda$CDM 长期以来一直受到来自宇宙微波背景（CMB）、重子声学振荡（BAO）和 Ia 型超新星（SN）的独立约束收敛性的支持。然而，近期的文献报告了不同探测器之间持续存在的差异（“张力”），其中最显著的是哈勃常数（$H_0$）张力、Planck 内部一致性问题（例如低 $\ell$ 与高 $\ell$ 的分裂、CMB 与引力透镜效应之间的差异），以及“曲率张力”（即 Planck 数据倾向于封闭宇宙，而 BAO 数据则倾向于平坦宇宙）。此外，关于这些张力是否需要对 $\Lambda$CDM 进行扩展（如演化的暗能量 $w_0w_a$CDM 或非零中微子质量）也存在持续争论。

一个关键的挑战在于，关于这些张力的定性结论往往对特定的数据产品、似然函数实现以及数据集组合高度敏感。不同的统计指标（证据比、后验偏移、基于拟合的度量）对先验体积和共享子空间的权重分配不同，这使得很难判断某种张力是稳健的物理信号，还是分析选择造成的伪影。本文旨在提供一个统一且一致的贝叶斯评估，以理清这些因素。

方法论
作者采用统一的贝叶斯框架，对五种宇宙学模型进行模型比较和数据集一致性评估，包括基准 $\Lambda$CDM 模型及四种极简扩展模型：

1. $\Lambda$CDM + $\sum m_\nu$（中微子质量之和）。
2. $\Lambda$CDM + $\Omega_K$（空间曲率）。
3. $w$CDM（常数暗能量状态方程）。
4. $w_0w_a$CDM（演化的暗能量）。

数据与似然函数：
本分析使用了一套全面的数据：

• CMB： Planck 公开版本 3 (PR3/2018) 和 Planck 公开版本 4 (PR4/NPIPE)。作者比较了多种高 $\ell$ 似然实现（Plik, CamSpec, HiLLiPoP）以及低 $\ell$ 似限函数（SimAll, LoLLiPoP），同时结合了 CMB 引力透镜似然函数（L3, L4）。
• BAO： 来自 SDSS (DR7, DR12, DR16) 和暗能量光谱仪 (DESI DR2) 的数据。
• SN： 来自暗能量调查 (DES y5)、Pantheon+ 和 Union3 的 Ia 型超新星数据。
• 距离阶梯： 包括局部 $H_0$ 测量值（Riess 2020, Freedman 2020）以作为背景参考。

统计框架：
研究利用了贝叶斯推断的三大支柱：

1. 参数估计： 推导后验分布。
2. 模型比较： 计算贝叶斯证据 ($Z$) 和对数证据差 ($\Delta \ln Z$)，以确定模型偏好，并纳入针对额外参数的 Occam 惩罚项。
3. 张力量化：
   • 证据比 ($R$)： 将两个数据集的联合证据与它们各自证据的乘积进行比较。若 $\ln R \ll 0$，则表明存在张力。
   • 可疑度 ($S$)： 一种源自“Occam 方程”的、独立于先验的基于似然函数的度量。它将张力分解为拟合优度成分和 Occam 惩罚成分。$S$ 通过贝叶斯模型维度的差异 ($\tilde{d}$) 转换为 $p$ 值和 “$\sigma$” 置信度。

作者使用 CosmoPower 模拟器进行快速似然评估，并使用嵌套采样（PolyChord）进行证据计算。

关键结果

• 内部 CMB 一致性： 更新后的 Planck 处理过程 (PR4) 系统性地改善了内部一致性。当从 PR3 (Plik) 转向 PR4 (CamSpec, HiLLiPoP) 时，低 $\ell$ 与高 $\ell$ 多极矩之间，以及原初 CMB 与 CMB 引力透镜效应之间的张力显著降低。例如，在 PR3 中，曲率扩展模型对封闭宇宙的偏好（这曾导致与 BAO 的张力）在采用 PR4 似然函数和更新后的 BAO 数据 (DESI) 后基本消失。
• 曲率张力： “曲率张力”（即 CMB 倾向于 $\Omega_K < 0$ 与 BAO 倾向于平坦的矛盾）被确定很大程度上是早期 Planck 似然实现（PR3/Plik）的一个特征。通过使用 PR4 和现代 BAO 数据 (DESI)，该张力降至 $1\sigma$ 以下，且对曲率的偏好消失了。
• 暗能量扩展：
  • $w$CDM： 用常数 $w$ 取代 $\Lambda$ 并不能稳健地改善拟合或解决张力；它通常仅仅是降低了 CMB 对 $H_0$ 的约束能力，而未能解决潜在的差异。
  • $w_0w_a$CDM： 在 CMB+BAO+SN 组合中观察到的对于演化暗能量的表观证据，对特定似然函数的配对高度敏感。仅在特定的组合下（例如 Planck PR3 + DES y5）才会观察到对 $w_0w_a$CDM 的强烈偏好。切换到 Planck PR4 和/或 Pantheon+ SN 数据会逆转这种偏好，恢复对 $\Lambda$CDM 的适度偏好。
  • DES Dovekie 更新： 文中指出，随后对 DES 超新星似然函数的更新（“DES Dovekie”）移动了 $\Omega_m$ 后验分布，使其与 Pantheon+ 达成一致，并进一步减少了张力，这证实了此前基于 DES y5 声称需要转向 $w_0w_a$CDM 的结论过于草率。
• 中微子质量： 允许自由的中微子质量之和 ($\sum m_\nu$) 倾向于吸收残余的偏移，而不会引入新的不一致性。作者发现，一旦考虑了完整的似然上下文，不存在稳健的、由 $\tau$（光学厚度）驱动的独立差异。
• 度量敏感性： 研究强调，“张力”并非单一的不变量。不同的似然选择（例如 CamSpec 与 HiLLiPoP，两者均基于 PR4）会产生可测量的不同张力统计量（差异约为 $\sim 0.5\sigma$），这强调了即使是 “PR4” 也无法唯一确定高 $\ell$ 似然的选择。

意义与结论
本文得出结论：要求将标准宇宙学模型从 $\Lambda$CDM 更新为动力学暗能量模型（$w_0w_a$CDM）的说法是过早的。作者认为：

1. 许多报告的张力对分析选择（似然实现、数据发布）非常敏感，而非数据的固有属性。
2. 数据处理（Planck PR3 $\to$ PR4）和数据汇编（SDSS $\to$ DESI）的更新系统性地提高了一致性。
3. 扩展模型中表观的“张力缓解”往往是通过降低数据集的约束能力（增加不确定性）实现的，而非通过获得真正改进的同步拟合。
4. “曲率张力”以及关于演化暗能量的最强信号，在合理的各种数据产品选择下并不稳健。

作者主张维持多条独立的分析流水线，以量化残余系统误差，并警告不要过度解读那些可能随着分析成熟和系统误差处理修订而发生变化的差异。他们强调，目前的统计证据并不强烈要求背离 $\Lambda$CDM 范式。