Constraining viscous fluid models in $f(Q)$ gravity with data

该研究利用 DESI、宇宙学时钟、超新星及结构增长等多组观测数据，结合 Akaike 和贝叶斯信息准则，在 $f(Q)$ 引力框架下约束了包含体粘性的宇宙模型，发现仅非粘性的幂律模型获得了观测支持，而指数和对数模型均被统计检验排除。

要点

这篇论文就像是一场宇宙侦探小说。侦探们（作者）试图解开宇宙加速膨胀和星系形成的谜题，他们不满足于现有的“标准答案”（$\Lambda$CDM 模型），而是去测试几种新的“替代理论”（$f(Q)$ 引力理论），并特别考察了一个叫“粘滞性”（Bulk Viscosity）的因素是否在其中扮演了关键角色。

为了让你轻松理解，我们可以把宇宙想象成一个正在膨胀的巨型面团，把引力理论想象成揉面团的配方。

1. 背景：为什么我们要找新配方？

目前的宇宙学“标准配方”叫 $\Lambda$CDM 模型。它假设宇宙里充满了看不见的“暗能量”（像是一种神秘的推手，把宇宙越推越快）和“暗物质”。这个配方虽然好用，但有些问题解释不通，比如宇宙膨胀的速度（哈勃常数 $H_0$）在不同测量方法下对不上号，就像用不同的尺子量同一个东西，结果却不一样。

于是，物理学家们开始尝试新配方。这篇论文关注的是 $f(Q)$ 引力理论。

• 旧配方（广义相对论）： 认为引力是时空的“弯曲”（像保龄球压在蹦床上）。
• 新配方（$f(Q)$ 理论）： 认为引力其实是时空的“非度量性”（Non-metricity）。想象一下，如果你拿着尺子在宇宙里走，尺子的长度会随着位置变化，这种“尺子变来变去”的特性就是 $Q$。$f(Q)$ 理论就是研究这种特性如何影响宇宙。

2. 核心变量：宇宙面团的“粘滞性”

论文引入了一个关键概念：体粘滞性（Bulk Viscosity）。

• 比喻： 想象你在搅拌蜂蜜（高粘滞性）和水（低粘滞性）。蜂蜜流动时内部有摩擦，会消耗能量并产生热量。
• 在宇宙中： 宇宙中的流体（物质）如果具有粘滞性，它在膨胀时内部会有摩擦。这种摩擦会产生一种“有效压力”，可能会改变宇宙膨胀的速度，甚至不需要“暗能量”就能解释宇宙为什么在加速膨胀。
• 研究目的： 作者想知道，如果在 $f(Q)$ 的新配方里加入“粘滞性”这个调料，能不能让理论更完美地符合观测数据？

3. 实验过程：三种新配方 vs. 真实数据

作者挑选了三种具体的 $f(Q)$ 理论模型（就像三种不同的新面团配方）：

1. 幂律模型 ($f_1$CDM)： 像是一个简单的数学幂次方关系。
2. 指数模型 ($f_2$CDM)： 像是一个指数增长的关系。
3. 对数模型 ($f_3$CDM)： 像是一个对数关系。

然后，他们收集了宇宙中所有的“目击证词”（观测数据）：

• 宇宙时钟 (CC)： 通过老恒星的年龄来推算宇宙膨胀速度。
• 声波指纹 (BAO)： 宇宙早期留下的声波印记，像尺子一样测量距离。
• 超新星 (PantheonP + SH0ES)： 宇宙中的“标准烛光”，用来测量距离。
• 星系生长率 (f 和 $f\sigma_8$)： 观察星系团是如何在引力作用下聚集长大的。

作者利用超级计算机（MCMC 模拟），让这三种新配方在“有无粘滞性”两种情况下，去拟合这些真实数据。

4. 调查结果：谁赢了？

经过严格的统计打分（使用 AIC 和 BIC 标准，类似于给模型打分，分数越低越好，且要惩罚过于复杂的模型），结果如下：

• 关于“粘滞性”：

  • 结论： 加不加“粘滞性”调料，并没有让理论变得更好。
  • 比喻： 就像你试图在蛋糕里加一种神秘的香料，结果发现加了之后，蛋糕不仅没变好吃，反而因为增加了不必要的成分（参数），让评委（统计数据）觉得这个配方太复杂且不必要。
  • 数据表现： 加入粘滞性后，模型的统计分数（$\Delta$AIC 和 $\Delta$BIC）变差了，说明数据并不支持宇宙流体有显著的粘滞性。

• 关于三种模型：

  • 指数模型 ($f_2$) 和 对数模型 ($f_3$)： 被“淘汰”了。 无论有没有粘滞性，它们都无法很好地解释观测数据，被统计标准直接拒绝。
  • 幂律模型 ($f_1$CDM)： 唯一的“幸存者”。
    • 在没有粘滞性的情况下，这个模型表现最好，甚至能和一些观测数据达成“实质性支持”（Substantial observational support）。
    • 一旦加入粘滞性，它的表现就稍微变差了，虽然没被完全淘汰，但不再是最优解。

5. 最终结论：侦探的总结

1. 新理论有潜力： $f(Q)$ 引力理论（特别是幂律模型）是一个很有希望的替代方案，它不需要暗能量就能解释宇宙加速膨胀，且在没有粘滞性的情况下，能很好地拟合数据。
2. 粘滞性不是关键： 在这项研究中，宇宙流体的“粘滞性”并没有被证实是解释宇宙加速膨胀的关键因素。加入它反而让模型变得臃肿且不符合数据。
3. 未来的路： 虽然幂律模型表现不错，但作者也指出，要完全取代 $\Lambda$CDM 模型，还需要更多测试（比如结合宇宙微波背景辐射 CMB 的数据）。

一句话总结：
这篇论文就像是在测试几种新的宇宙“引擎”设计。结果发现，其中一种设计（幂律 $f(Q)$）在没有额外“摩擦”（粘滞性）的情况下，跑得最稳、最符合观测；而试图给引擎加“摩擦”不仅没帮助，反而让设计变得多余。其他两种设计则因为跑得太偏，直接被刷掉了。

技术摘要

这是一份关于在 $f(Q)$ 引力框架下利用观测数据约束粘性流体模型的详细技术总结。

1. 研究背景与问题 (Problem)

• $\Lambda$CDM 模型的局限性：尽管 $\Lambda$CDM 模型是目前解释宇宙现象的主流模型，但其面临的挑战（如哈勃常数 $H_0$ 张力和 $S_8$ 张力）促使物理学家寻求替代理论。
• $f(Q)$ 引力理论：作为广义相对论（GR）的修正理论之一，$f(Q)$ 引力基于非度量性（non-metricity, $Q$）而非曲率或挠率来描述引力相互作用。它旨在在不引入暗能量成分的情况下解释宇宙加速膨胀。
• 粘性流体的作用：在宇宙流体中引入体粘性（bulk viscosity）已被证明可以改变背景动力学，解释加速膨胀和熵产生，并影响结构形成。
• 核心问题：目前尚不清楚在 $f(Q)$ 引力框架下，引入体粘性是否能改善模型对观测数据的拟合度，以及不同的 $f(Q)$ 模型（幂律、指数、对数）在考虑粘性后的宇宙学可行性。

2. 方法论 (Methodology)

• 理论框架：
  • 基于 $f(Q)$ 引力理论，假设物质流体具有体粘性，有效压力修正为 $\bar{p} = p_m + \zeta H$（其中 $\zeta$ 为粘性系数，$H$ 为哈勃参数）。
  • 推导了修正后的弗里德曼方程（Friedmann equations）和标量微扰演化方程。
  • 利用 1+3 协变形式（1+3 covariant formalism）和准静态近似（quasi-static approximation），推导了密度对比度 $\delta(z)$、增长率 $f(z)$ 和红移空间畸变 $f\sigma_8(z)$ 的演化方程。
• 模型选择：研究了三种典型的 $f(Q)$ 模型：
  1. $f_1$CDM：幂律模型 ($F(Q) = \alpha Q^n$)。
  2. $f_2$CDM：指数模型 ($F(Q) = \beta Q_0 (1 - e^{-p\sqrt{Q/Q_0}})$)。
  3. $f_3$CDM：对数模型 ($F(Q) = \epsilon \ln(\Gamma Q/Q_0)$)。
• 观测数据：
  • 宇宙时计 (CC)：31 个 $H(z)$ 数据点。
  • 重子声学振荡 (BAO)：来自 DESI 巡天的距离和关联测量。
  • Ia 型超新星 (SNIa)：PantheonP + SH0ES 样本（1701 条光变曲线）。
  • 结构增长数据：增长率 $f(z)$ 和红移空间畸变 $f\sigma_8(z)$ 数据（来自 VIPERS 和 SDSS）。
  • 进行了多种数据集的联合分析（如 CC+BAO, PantheonP+SH0ES+f, PantheonP+SH0ES+f$\sigma_8$ 等）。
• 统计方法：
  • 使用 MCMC 模拟工具 Kosmulator 进行参数约束。
  • 拟合参数包括：$\Omega_m, H_0, r_d, M_{abs}, \sigma_8, \gamma$ 以及模型参数 $n, p, \Gamma$ 和粘性系数 $\zeta$。
  • 使用 AIC (Akaike Information Criterion) 和 BIC (Bayesian Information Criterion) 进行模型选择，以 $\Lambda$CDM 为基准比较各模型的优劣（$\Delta IC$ 判据）。

3. 关键贡献 (Key Contributions)

1. 首次系统性研究：这是文献中首次同时研究体粘性对 $f(Q)$ 引力框架下宇宙晚期加速膨胀和物质密度涨落（大尺度结构形成）的影响。
2. 理论推导完善：在 $f(Q)$ 引力中引入了粘性流体，推导了包含粘性项的修正弗里德曼方程和密度微扰演化方程，并给出了三种典型模型的具体解析形式。
3. 全面的数据约束：利用最新的 DESI BAO 数据、PantheonP+SH0ES 超新星数据以及最新的增长率数据，对模型进行了严格的参数约束。
4. 模型选择与粘性效应评估：通过 AIC/BIC 统计标准，量化了引入粘性系数 $\zeta$ 是否增加了模型的统计显著性，并对比了三种 $f(Q)$ 模型的可行性。

4. 主要结果 (Results)

• 参数约束：
  • $f_1$CDM (幂律模型)：在无粘性情况下，该模型表现出稳健的参数估计，且未遭受到任何数据集的明确拒绝。引入粘性后，$\Omega_m$ 倾向于更高值，$S_8$ 略微降低，但统计性能略有下降。
  • $f_2$CDM (指数模型) 和 $f_3$CDM (对数模型)：无论是否引入粘性，这两个模型在多个数据集（特别是包含 PantheonP+SH0ES 的数据集）上均被 AIC/BIC 标准拒绝（$\Delta IC \ge 10$）。
  • 粘性系数 $\zeta$：在所有模型中，引入粘性系数 $\zeta$ 作为一个自由参数，并没有显著改善 $\chi^2$ 拟合优度，反而增加了 $\Delta AIC$ 和 $\Delta BIC$ 值。这表明从统计角度看，引入体粘性缺乏正当性（即模型过于复杂且未带来更好的拟合）。
• 张力缓解情况：
  • 粘性效应在某些数据集中略微增加了 $\Omega_m$ 的估计值，可能有助于缩小局部距离阶梯测量与 CMB 测量之间的差异，但在其他数据集（如 $f$ 和 $f\sigma_8$）中，这种效应反而加剧了张力。
  • 粘性导致 $S_8$ 值降低，这在一定程度上缓解了 $S_8$ 张力，但 $f_3$CDM 模型在粘性下得到的 $S_8$ 值过低，超出了合理范围。
• 结构形成：
  • 无粘性的 $f(Q)$ 模型（特别是 $f_1$CDM）在红移空间畸变 $f\sigma_8(z)$ 的预测上与观测数据吻合良好，甚至在某些红移范围内优于 $\Lambda$CDM。
  • 粘性对增长率 $f(z)$ 的影响较小，主要表现为对密度微扰演化的额外阻尼。

5. 结论与意义 (Significance)

• 模型优选：在研究的三种 $f(Q)$ 模型中，无粘性的幂律模型 ($f_1$CDM) 是唯一一个在统计上稳健且未遭明确拒绝的替代模型。它在多个数据集上获得了“实质性观测支持”（Substantial observational support）。
• 粘性效应的否定：研究结果表明，在当前观测数据的精度下，在 $f(Q)$ 引力框架下引入体粘性不能显著改善模型对宇宙的拟合。相反，它增加了模型的复杂性而未带来统计收益，导致模型被 BIC 标准拒绝。
• 理论启示：虽然粘性流体在早期宇宙或特定理论场景中可能重要，但在解释晚期宇宙加速膨胀和大尺度结构形成的 $f(Q)$ 模型中，简单的无粘性修正可能更为有效。
• 未来方向：研究建议未来需要结合 CMB 物理等更多观测手段，对 $f(Q)$ 模型进行更严格的约束，以进一步验证粘性在修正引力理论中的潜在作用。

总结：该论文通过严谨的统计分析和多源观测数据，确立了 $f(Q)$ 幂律模型作为 $\Lambda$CDM 的有力替代者，同时指出在当前数据水平下，引入体粘性并不能提升这些修正引力模型的宇宙学解释能力。