4.3 量刑推演：密度演化与宇宙命运

4.3 量刑推演：密度演化与宇宙命运

本节摘要：把 w 的口供外推到遥远未来，宇宙的终局只有三种剧本：w 等于负一是永恒加速的冷寂，w 大于负一减速回归或停滞，w 小于负一则密度发散、宇宙在有限时间内被大撕裂——幻影能量 w 取负 1.1 与负 1.5 的撕裂时刻分别在约 1160 亿年与 230 亿年后。本节完成这三条时间线的数值推演，把"暗能量是什么"与"宇宙怎么收场"焊在一起。

审讯的最后一问是量刑：如果口供属实，宇宙会怎么收场？本节把第四章的参数推到未来极限，三条终局时间线全部算出。这一节也顺带回收旧卷宗（原教程"宇宙命运"主题在此合并审理）。

剧本一：恒定负一，永恒的冷寂

宇宙学常数的口供推到未来：物质继续稀释、真空能恒定，两者的比值只增不减。加速永远持续，但膨胀率趋近一个非零的下限——宇宙以指数方式越来越空、越来越冷。星系彼此退行出视界，恒星的燃料账本逐项耗尽，最后的能量流动归于涨落。这不是爆炸而是漫长的散场，业内叫法很多：热寂、大冻结、渐冷终局。

# 演算 1：冷寂剧本的账目推演
import math
Om, OL, H0 = 0.315, 0.685, 67.4
t_H = 14.51                      # 哈勃时间（Gyr）
def t_a(a):                      # 由尺度因子算宇宙年龄（a>1 为未来）
    n = 40000; da = (a-1.0)/n; s = 0.0
    for i in range(n):
        am = 1.0+(i+0.5)*da
        s += 1.0/(am*math.sqrt(Om*am**-3+OL))
    return 13.80 + s*da*t_H
a100 = (100*OL/Om)**(1/3)        # 真空能比物质大百倍的时刻
print("真空能超出物质百倍：尺度因子 a = %.2f" % a100)
print("那时宇宙年龄 t = %.1f Gyr（距今 %.1f Gyr）" % (t_a(a100), t_a(a100)-13.80))
print("渐近膨胀率 = H0×√OL = %.1f km/s/Mpc" % (H0*math.sqrt(OL)))
# 输出:
# 真空能超出物质百倍：尺度因子 a = 6.01
# 那时宇宙年龄 t = 44.1 Gyr（距今 30.3 Gyr）
# 渐近膨胀率 = H0×√OL = 55.8 km/s/Mpc
​

三百亿年后真空能对物质的压制就到百倍量级，加速进入"无人驾驶"——引力束缚系统内部不受影响（星系、星系团靠自身引力锁住），但一切未束缚的天体将永久退场。有趣的是这条剧本的"本地新闻"并不惊悚：太阳系照常，银河系与仙女座照常并合，只是可观测宇宙的内容清单越来越短。

剧本二：w 大于负一，减速回归或悬停

动力学口供若 w 恒大于负一（如解冻类场的晚期值负 0.9），密度的稀释虽慢于物质但仍在进行。两条支线：温和情形下加速逐渐减弱，宇宙滑向"悬停"（膨胀率趋零但不坍缩，视 w 的渐近值而定）；若 w 渐近回到负三分之一以上，减速回归，宇宙进入幂律的缓慢稀释。这条剧本的关键可观测后果是结构生长不会完全停摆——真空能主导的冷寂剧本里生长终将冻结（有效引力为零），w 大于负一的剧本里冻结不完全。遥远未来的星系形成史因此带着 w 的签名，尽管实测它们远超当前能力。

剧本三：w 小于负一，有限时间的大撕裂

幻影口供的推演最戏剧化。密度随膨胀增长（1.3 的正指数），膨胀率越来越大，尺度因子在有限时间内发散——不是渐近而是爆炸式的收场，业内称大撕裂。撕裂的时刻可以用一个收敛积分精确算出。

# 演算 2：大撕裂时刻的积分（数值段 + 解析尾巴）
import math
Om = 0.315
Ode = 1-Om
t_H = 14.51
def t_rip(w):
    # H(a)=H0·√(Om·a^-3 + Ode·a^{-3(1+w)})；t_rip = ∫_1^∞ da/(aH)
    # 数值段到 a=1000（对数步长），幻影项主导的尾巴解析积分
    p = -3*(1+w)/2          # 幻影密度指数的一半，w<-1 时为正
    n = 20000; la1 = math.log(1000.0); dla = la1/n; s = 0.0
    for i in range(n):
        a = math.exp((i+0.5)*dla)
        H = math.sqrt(Om*a**-3 + Ode*a**(-3*(1+w)))
        s += dla/H          # da/(aH) = dln(a)/H
    tail = 1000.0**(-p)/(p*math.sqrt(Ode))
    return (s+tail)*t_H
for w in [-1.05, -1.10, -1.20, -1.50]:
    print("w = %.2f -> 距今 %.1f Gyr 后大撕裂" % (w, t_rip(w)))
# 输出:
# w = -1.05 -> 距今 232.6 Gyr 后大撕裂
# w = -1.10 -> 距今 115.9 Gyr 后大撕裂
# w = -1.20 -> 距今 57.5 Gyr 后大撕裂
# w = -1.50 -> 距今 22.7 Gyr 后大撕裂
​

撕裂时刻对 w 极度敏感：从负 1.5 到负 1.05，时刻从 230 亿年推到 2330 亿年——w 每靠近负一一点，死刑日期指数推迟。撕裂的过程同样有次序：先撕散星系团，再撕散星系（撕裂前约六千万年），然后是行星系统（前三个月）、恒星与行星（前一分钟）、原子（前十的负十九次方秒）。当前观测的 w 误差棒（负一加减正负 0.05）对应的撕裂时刻全部远在千亿年外，但这不妨碍幻影剧本在理论上被严肃对待——它同时是第五章稳定性关卡的重点审查对象（负动能的鬼场问题）。

图 三条终局时间线的对照

量刑的现实意义

有人问：终局在千亿年外，量刑有什么用？有用在反向约束：撕裂时刻对 w 的指数敏感性意味着，只要把 w 钉得更接近负一，就等于给宇宙开出了更长的寿命证明；反过来，任何声称 w 显著小于负一的观测都在缩短宇宙的死刑日期——这是把宇宙学参数测量与"宇宙命运"直接挂钩的少数场景。此外终局剧本还牵动一个哲学问题（7.1 会接手）：如果终局是永续冷寂，可观测宇宙的熵预算、生命与信息的远期前景都有硬约束——科学史上少有哪次"量刑推演"如此直白地通向存在层面的追问。

⚠️ 常见坑：把大撕裂渲染成"马上发生的爆炸"。即便幻影口供属实，现有误差棒下的撕裂时刻也在千亿年量级；"宇宙将死于撕裂"与"宇宙近期死于撕裂"是两个完全不同的命题，科普报道最爱混掉的正是这一条。

💡 关键直觉：三剧本的分岔点全在 w 等于负一这条线上——密度恒定（冷寂）、密度仍降（悬停）、密度反升（撕裂）。一个参数的三种取值区间，对应存在层面的三种终局，这就是暗能量问题分量之所在。

本节要点回顾

• 冷寂剧本：真空能恒定，约 300 亿年后对物质压制达百倍，之后指数膨胀永续
• 悬停剧本：w 大于负一，减速回归或趋停，结构生长不完全冻结
• 撕裂剧本：w 小于负一，密度发散，宇宙死于有限时间
• 敏感性：撕裂时刻随 w 指数变化，负 1.5 对应 230 亿年，负 1.05 对应 2330 亿年
• 反向约束：钉住 w 等于给宇宙续命或定罪，参数测量与终局推演互为表里

审讯全部结束：读数是"接近负一的常数"，动机问题原封未动。下一章清点三桩悬案——那是比终局更让人睡不着的东西。