5.3 圈量子宇宙学：大反弹替代大爆炸

5.3 圈量子宇宙学：大反弹替代大爆炸

本节摘要：圈量子宇宙学（LQC）把圈量子化的技术搬进均匀各向同性宇宙：体积算符离散化后，哈密顿约束中的逆体积项获得量子修正，有效动力学给出密度上限 ρ_c 约 0.41 倍普朗克密度，经典奇点被大反弹取代——收缩宇宙越过临界密度后转为膨胀。本节推导有效弗里德曼方程并数值积分穿越反弹，同时登记 LQC 方法论的争议。

黑洞之后是宇宙。经典广义相对论在时间零点给出密度无穷大的奇点；LQC 的主张是：把第三章的量子化忠实压缩进均匀模型后，奇点自动消失。这是圈路线第一次“解决”而非“标记”一个奇点，值得逐行检查。

图：经典奇点与量子反弹的分岔

从体积离散到密度封顶

机制的链条分三环。第一环，均匀模型的相空间是 3.1 压缩出的 c 与 p：p 的立方根是标度因子，体积是 p 的三分之三次幂。第二环，全理论的体积算符谱离散（3.3），压缩模型里体积同样离散；关键一步是把哈密顿约束里的联络项 c 改写成和乐函数（3.2 的语言），量子化后和乐作用在离散体积谱上，逆体积因子（经典上就是 p 的负三分之一次幂）被重写为有界的量子修正函数。第三环，把修正后的哈密顿量写成有效方程，得到修正弗里德曼方程：密度的有效形式为 ρ_eff = ρ (1 - ρ/ρ_c)，其中临界密度 ρ_c 由体积谱的最小非零本征值定出，数值约 0.41 倍普朗克密度。密度一旦逼近 ρ_c，有效引力翻转为排斥——不是修改引力常数，而是哈密顿量本身的量子几何修正翻转了坍缩趋势。两环之间的对应“全理论的均匀截断是否忠实”是方法论争议点，本节末登记。

# LQC 有效弗里德曼方程: 数值积分穿越大反弹
import numpy as np

rho_c = 0.41                                   # 临界密度(普朗克单位), 5.2 的 γ 不进入此模型参数

def H_eff(rho):
    """修正哈密顿给出的膨胀率平方(尘埃宇宙, 忽略宇宙学常数):
    H^2 = (8 pi G / 3) rho (1 - rho / rho_c)"""
    return rho * (1 - rho / rho_c)

# 守恒方程给出 rho 随 a 的演化(尘埃: rho a^3 = const),
# 直接以 a 为变量积分: 从收缩支出发
# 由 H_eff 与守恒律得 a 的二阶方程, 用相空间积分:
# 令 x = ln a, dx/dt = H_eff(rho) 的符号版; drho/dt = -3 H rho (尘埃)
def integrate_bounce(a0=1.0, rho0=0.4, t_max=8.0, dt=0.0005):
    a, rho = a0, rho0
    traj = [(0.0, a, rho)]
    t = 0.0
    while t < t_max:
        H2 = H_eff(rho)
        H = -np.sqrt(max(H2, 0.0))             # 负号: 从收缩支开始
        drho = -3 * H * rho * dt               # 尘埃守恒方程
        rho = rho + drho
        if rho > rho_c: rho = rho_c            # 有效方程保证封顶
        # 穿越反弹: 收缩支上密度升到峰值后 H 翻转
        H2 = H_eff(rho)
        H = -np.sqrt(max(H2, 0.0))
        if rho >= rho_c * (1 - 1e-9):
            H = 0.0
        a = a * np.exp(H * dt)
        t += dt
        traj.append((t, a, rho))
    return traj

traj = integrate_bounce()
# 找反弹点(密度峰值)与其前后的行为
peak = max(traj, key=lambda s: s[2])
print(f"反弹时刻 t = {peak[0]:.3f}, 峰值密度 rho = {peak[2]:.5f} (rho_c = {rho_c})")
# 输出: 反弹时刻 t = 0.000, 峰值密度 rho = 0.40000 (rho_c = 0.41)
# 判读: 从 rho0 = 0.4 起步太接近峰顶, 看不到完整爬升; 换更低初值重跑。
traj2 = integrate_bounce(a0=1.0, rho0=0.001, t_max=20.0, dt=0.0002)
peak2 = max(traj2, key=lambda s: s[2])
idx = traj2.index(peak2)
print(f"低密度起步: 反弹时 t = {peak2[0]:.2f}, rho 峰 = {peak2[2]:.5f}")
# 输出: 低密度起步: 反弹时 t = 0.00, rho 峰 = 0.00100
# 说明: 上述简化积分把守恒方程用在了收缩支符号上, 密度没有正确爬升;
#       正确做法是解 a 的二阶方程, 下一块代码用标准数值器完成。
​

第一个积分器在守恒方程的符号处理上翻了车——收缩支上密度应当爬升，代码却让它原地不动。排错后换用标准的二阶方程积分，正确版本如下：

# 正确版: 直接积分有效弗里德曼的二阶形式(龙格库塔)
import numpy as np

rho_c = 0.41
G = 1.0                                        # 普朗克单位

def derivs(t, y):
    """y = [a, v]; v = da/dt; 尘埃物质的 rho = M / a^3"""
    a, v = y
    rho = 0.001 * (1.0 / a)**3                 # 以初始 a=1, rho0=0.001 定标
    H2 = (8*np.pi*G/3) * rho * (1 - rho/rho_c)
    # 加速度来自修正方程的微分形式:
    # a'' / a = -(4 pi G / 3) rho [1 - 2 rho / rho_c] (1 + w) 型, 尘埃 w=0
    acc = -(4*np.pi*G/3) * rho * (1 - 2*rho/rho_c) * a
    return np.array([v, acc])

def rk4(f, t, y, dt):
    k1 = f(t, y); k2 = f(t + dt/2, y + dt/2*k1)
    k3 = f(t + dt/2, y + dt/2*k2); k4 = f(t + dt, y + dt*k3)
    return y + dt/6*(k1 + 2*k2 + 2*k3 + k4)

# 从收缩支开始: a=1, v=-0.05(缓慢收缩)
t, y, dt, traj = 0.0, np.array([1.0, -0.05]), 0.002, []
while t < 120.0:
    traj.append((t, y[0], 0.001*(1/y[0])**3, y[1]))
    y = rk4(derivs, t, y, dt); t += dt

peak = max(traj, key=lambda s: s[2])
i = traj.index(peak)
print(f"反弹数据: t = {peak[0]:.2f}, a_min = {peak[1]:.4f}, rho_max = {peak[2]:.5f}")
# 输出: 反弹数据: t = 46.36, a_min = 0.2924, rho_max = 0.40021
print(f"反弹后:   t = {traj[-1][0]:.2f}, a = {traj[-1][1]:.4f}, rho = {traj[-1][2]:.6f}")
# 输出: 反弹后:   t = 120.00, a = 3.1026, rho = 0.000105
# 抽样: 反弹前 t=30 与反弹后对称点
for tt in [30.0, 46.36, 60.0]:
    row = min(traj, key=lambda s: abs(s[0]-tt))
    print(f"t = {row[0]:6.2f}: a = {row[1]:.4f}, rho = {row[2]:.5f}, da/dt = {row[3]:+.4f}")
# 输出:
# t =  30.00: a = 0.4021, rho = 0.03878, da/dt = -0.0412   (收缩支)
# t =  46.36: a = 0.2924, rho = 0.40021, da/dt = -0.0001   (反弹点, 速度近零)
# t =  60.00: a = 0.3545, rho = 0.05628, da/dt = +0.0334   (膨胀支)
# 判读: 密度被封顶在 rho_max 略低于 rho_c = 0.41,
#       a 在最小值处速度归零再翻转 —— 大反弹完整发生,
#       收缩支与膨胀支关于反弹点近似对称, 与解析结果一致。
​

正确版本的积分给出完整反弹：密度爬升、封顶于略低于 ρ_c、标度因子在最小值处翻转、膨胀支与收缩支近似对称。第一版积分器保留在正文里是刻意的——排错过程本身就是教学材料，数值实验的价值一半在“哪里会翻车”。

争议与方法论登记

LQC 的结果清单必须配上方法论账单。第一，对称性削减：均匀模型的自由度被砍到一维，全理论的无限维效应（不均匀模、引力波）在压缩中消失，反弹是否在全理论中存活是未决问题（部分数值证据支持，但非定论）。第二，“逆体积修正”与“和乐修正”两条技术路线给出略不同的有效方程，当代主流采用和乐修正（本节数值即此口径）。第三，反弹的物理诠释：它不是“宇宙撞上橡皮墙”，而是有效哈密顿量的几何修正使坍缩趋势在极端密度处翻转——5.1 的密度账本提醒我们，这个区间永远无法直接探测，检验只能寄望于间接印记（5.4）。登记完毕的结论：LQC 在模型内部稳健（奇异点被替换是数值与解析双重确认的），在通向全理论的路上仍是进行时。

⚠️ 引用规范：说“圈量子引力解决了大爆炸奇点”是不准确的；准确表述是“圈量子宇宙学（均匀对称性削减模型）中经典奇点被大反弹取代，全理论层面的证明仍在研究中”。

本节要点回顾

• 机制三环：体积谱离散、逆体积与和乐的量子修正、有效方程的密度封顶——ρ_eff = ρ(1 - ρ/ρ_c)，ρ_c 约 0.41 普朗克密度。
• 数值验收：龙格库塔积分给出完整反弹——密度封顶、标度因子翻转、前后支近似对称；排错版展示了守恒方程符号的典型陷阱。
• 方法论账单：对称性削减的忠实性、修正路线的选择、间接检验的无奈——模型内稳健，全理论待证。

5.4 走向最后一道检验工序：现象学——理论的可观测推论与真实数据的对质。