Lambda CDM模型中H/H0随a(t)变化绘图异常的技术问询

嘿，这个问题我之前也纠结过，核心其实是混淆了哈勃参数作为尺度因子a的函数和作为时间t的函数的行为，或者是对论文里的图的轴/度量理解错了。咱们一步步理清楚：

1. Lambda CDM中H(a)的正确行为 #

标准Lambda CDM模型里，哈勃参数作为尺度因子a的表达式是完全正确的：

H(a)/H₀ = √(Ωᵣ a⁻⁴ + Ωₘ a⁻³ + Ωₖ a⁻² + Ω_Λ)

用当前Planck 2018的观测参数（Ωᵣ≈8.5×10⁻⁵，Ωₘ≈0.31，Ω_Λ≈0.69，Ωₖ≈0）代入计算，这个函数随a增大是单调递减的：

• 早期宇宙（a→0）：辐射主导，H(a)∝a⁻²，快速下降；
• 物质主导阶段：H(a)∝a⁻¹.⁵，下降速率变慢；
• 暗能量主导阶段：H(a)趋近于常数H₀×√(Ω_Λ)≈0.83H₀，下降速率无限趋近于0，但绝不会回升。

2. 为什么论文里会出现“先降后升”？ #

最可能的几个原因：

• 看错了图的度量或坐标轴：论文里的图大概率展示的不是H(a)，而是其他相关量：
  • 比如宇宙减速参数q(t)：它在a≈0.7时从正变负，对应宇宙从减速膨胀切换到加速膨胀，视觉上会有一个明显的“转折”；
  • 或者是H(t)随时间t的变化：虽然H(t)整体还是下降的，但在a≈0.7之后，下降速率变得极慢，视觉上可能被误判为“回升”；
• 轴的方向搞反了：如果图的横轴是从a=1到a=0（从现在回溯到过去），那H(a)会随a减小而增大，看起来像是“先升后降”，和你描述的“先降后升”刚好相反。

3. 验证你的绘图结果 #

你可以手动计算几个关键节点的H(a)/H₀值来确认：

• 当a=0.7时：√(8.5e-5×0.7⁻⁴ + 0.31×0.7⁻³ + 0.69) ≈ 1.263
• 当a=1时：√(Ωᵣ+Ωₘ+Ω_Λ)=1
• 当a=2时：√(8.5e-5×2⁻⁴ + 0.31×2⁻³ + 0.69) ≈ 0.854

显然，从a=0.7到a=1，H(a)/H₀是从1.263降到1，完全是持续下降的，你的绘图结果是正确的。

总结 #

Lambda CDM模型中，H(a)/H₀随a增大确实是单调递减的，你的绘图没有问题。论文里的“先降后升”应该是对图的内容或坐标轴的误解，建议再仔细核对论文的图例和标注。

内容的提问来源于stack exchange，提问作者Vishal Verma