6.5 数值相对论与量子引力：未完成的海图
6.5 数值相对论与量子引力：未完成的海图
本节摘要：场方程的非线性让解析解只覆盖极少数对称情形，两支工程队伍分别从两端逼近其余部分。数值相对论把时空切成空间薄片逐层演化（3 加 1 分解与 BSSN 变量），三十年稳定性战争后于 2005 年攻克双黑洞并合，产物就是 6.2 节的模板库；量子引力则在普朗克地界（长度 10⁻³⁵ 米、时间 10⁻⁴⁴ 秒）处理几何的量子涨落，弦论与圈量子引力是两条主路线，黑洞信息疑难成为共同的试金石。航志在此收笔，附开放海域清单。
把场方程交给机器：3 加 1 分解
第 5、6 章的精确解全都依赖对称性（球对称、均匀各向同性），真实宇宙很少这么客气——双黑洞并合这类"最值钱"的过程毫无解析解可言，只能数值求解。主流方案是 3 加 1 分解（ADM 形式）：把四维时空切成一叠三维空间薄片，每片配一个时刻标签；度规拆成"片上的空间度规 γ_ij、片间如何粘接的移位 β_i、片面倾斜的 lapse 函数 α"三组变量。场方程随之拆成两类：
- 演化方程：告诉下一片的空间度规与外曲率长什么样（含时间导数，可逐步推进）；
- 约束方程：每一片都必须满足的相容条件（哈密顿约束与动量约束），像每步都要复查的预算平衡表。
配上坐标规范（如何选下一片、如何定移位——著名的 BSSN 变量组与规范条件是三十年经验的结晶），演化就能逐片推进：给定初始数据（两颗黑洞的位置动量），算出并合、提取 6.2 节波形模板。
数值稳定性的三十年战争
这事难在哪？1970 年代就写出了演化方程，此后三十年计算机翻了几个数量级，双黑洞轨道却总是算几十圈就崩——数值不稳定性像瀑布一样吞掉解。转机出现在世纪之交：BSSN 重新组合变量让方程更"驯服"、广义谐和坐标提供更稳的规范，2005 至 2006 年三个小组（普雷托里斯；贝克与同事；卡米内利与路易斯组）先后独立完成稳定演化。此后行业进入工业化：模板库填满探测器的匹配滤波器（6.2 节），并合波形的最后几圈不再靠近似拼凑。理解"为什么稳定格式至关重要"最好的办法是自己摔一跤：
# 演化格式的分野：显式欧拉 vs 龙格库塔（黑洞旁圆轨道） import math GM = 1.0; c = 1.0 # 几何单位 r0 = 8.0 # 稳定圆轨道半径，远在 ISCO 外 v0 = math.sqrt(GM/r0) def deriv(s): x, y, vx, vy = s rr = math.hypot(x, y) f = -GM/rr**3 return (vx, vy, f*x, f*y) def energy(s): x, y, vx, vy = s return 0.5*(vx*vx+vy*vy) - GM/math.hypot(x, y) def euler(s, dt): d = deriv(s) return [s[i] + dt*d[i] for i in range(4)] def rk4(s, dt): k1 = deriv(s) k2 = deriv([s[i]+dt/2*k1[i] for i in range(4)]) k3 = deriv([s[i]+dt/2*k2[i] for i in range(4)]) k4 = deriv([s[i]+dt*k3[i] for i in range(4)]) return [s[i]+dt/6*(k1[i]+2*k2[i]+2*k3[i]+k4[i]) for i in range(4)] T_orbit = 2*math.pi*r0/v0 dt = T_orbit/200 s_e = s_r = [r0, 0.0, 0.0, v0] E0 = energy(s_r) for n in range(int(20*T_orbit/dt)): s_e = euler(s_e, dt) s_r = rk4(s_r, dt) print(f"20 圈后能量相对漂移：") print(f" 显式欧拉: {(energy(s_e)-E0)/abs(E0):+.2e} （轨道肉眼可见地外扩）") print(f" 四阶龙格库塔: {(energy(s_r)-E0)/abs(E0):+.2e} （守恒到机器精度附近）") print("数值相对论的三十年，就是把这种漂移按进约束方程里的三十年")
日志读数：同样的步长，显式欧拉二十圈后能量漂移千分之一量级（轨道外扩），四阶龙格库塔守恒到近机器精度——这不是"精度"问题，是"约束是否被维持"的问题。数值相对论的全部艺术就在于：选变量、选规范、离散化，让哈密顿与动量约束在数万圈演化后仍被压在误差棒内。
量子引力：两条主路线
奇点定理（6.4 节）与霍金辐射（5.3 节）从两个方向指认同一处边界：广义相对论在极强场与极小尺度交权。交权给谁？两条主路线竞争了四十年。
弦论：出发点是把粒子换成振动的弦，统一全部相互作用；代价是额外维度（十维时空）与海量的可能真空。它对航志主线最大的回礼是全息原理与 ADS/CFT 对偶——某个引力理论可以完全等价于其边界上低一维的量子场论。5.3 节的黑洞信息疑难在这条路上拿到了至今最像答案的进展：岛公式把黑洞辐射的熵算出了符合幺正性的曲线（2019 至 2020 年的工作），"信息没有丢"从直觉升格为半严格计算。
圈量子引力：出发点不另加粒子，直接把几何本身量子化——面积与体积有离散谱，时空由自旋网络织成。它对主线的回礼是大反弹宇宙学：把 6.3 节的大爆炸奇点替换为前一收缩相的反弹，经典奇点在量子几何中消解。预言的可观测印记（原初引力波谱的修改）尚在检验能力边缘。
两条路线风格迥异（统一优先对几何优先），却共享同一个困境：量子引力的特征尺度远离实验台。算一下这笔账就明白：
# 普朗克地界：量子引力的"门牌号" import math hbar = 1.0546e-34; G = 6.674e-11; c = 2.99792458e8 l_p = math.sqrt(hbar*G/c**3) t_p = math.sqrt(hbar*G/c**5) m_p = math.sqrt(hbar*c/G) rho_p = m_p/l_p**3 print(f"普朗克长度 {l_p:.2e} m（质子半径的 1e-20 倍）") print(f"普朗克时间 {t_p:.2e} s") print(f"普朗克质量 {m_p:.2e} kg（约 22 微克——宏观可称量！）") print(f"普朗克密度 {rho_p:.2e} kg/m3") print("把质子放大到可观测宇宙，也才放大到普朗克长度的量级")
读数值得玩味：普朗克长度小到荒谬（质子的一亿亿亿分之一），但普朗克质量只有 22 微克——一粒盐的质量。量子引力难的不是质量，而是把这么小的质量压缩到这么小的尺度（需要的质量密度 10⁹⁶ 千克每立方米）。这解释了为什么"实验窗口"如此稀缺：LHC 的能标差了十五个数量级，自然实验室只剩早期宇宙与黑洞视界附近。
航志收笔：开放海域清单
广义相对论的主航线至此闭合：从自由落体的电梯出发，途经仪器舱、场方程、四大检验、黑洞与涟漪，抵达宇宙学与因果边界。收笔处把未绘制的海域开列如下，作为下一册航志的候选航线：
- 暗能量的本质：宇宙学常数还是动力学场？状态方程测量是首要任务（6.3 节的预算表悬而未决）；
- 黑洞信息疑难：岛公式之后，幺正性的完整证明仍差最后一步（5.3 与本节弦论路线）；
- 宇宙监督与裸奇点：数值实验继续搜证（6.4 节悬案）；
- 桌面量子引力：让两枚微小质量引力纠缠的实验方案，试图把 22 微克的门槛搬进实验室；
- 引力波色散检验：引力子质量上限、偏振模式清点——6.1 节的"两种偏振"若被推翻，修改引力翻身；
- 多信使标准警笛：引力波测哈勃常数与距离阶梯的分歧仲裁（6.2 与 6.3 节的接口）。
这册航志没有结束语式的总结——广义相对论是一座仍在施工的理论，任何"总结"都会在下一场并合事件或下一次巡天数据发布后过期。合上日志时请记住第 1 章的那次换轨：你已经知道如何用自由落体的眼睛看世界，剩下的海图，留给下一次出航。
本站要点
- 3 加 1 分解：时空切片、度规拆三组变量、方程拆演化与约束两类——约束维持是数值稳定的生命线；
- 三十年战争：BSSN 变量与规范条件驯服不稳定性，2005 至 2006 年双黑洞并合被攻克，产物即 6.2 节的模板库；
- 格式实验：同一步长下欧拉漂移千分之一、龙格库塔近机器精度——数值相对论的艺术就是把这种漂移按进约束里；
- 两条量子路线：弦论凭全息与岛公式回击信息疑难，圈量子引力以自旋网络给出反弹宇宙；普朗克地界（10⁻³⁵ 米、22 微克、10⁹⁶ 密度）远离实验台是共同困境；
- 开放清单：暗能量、信息疑难、宇宙监督、桌面量子引力、偏振与色散检验、标准警笛——六片未绘制的海域，各配了本航志的对应锚点。