2.4 测地线：自由航迹的方程

2.4 测地线：自由航迹的方程

本节摘要：测地线是弯曲时空中"最直"的航线，由固有时变分原理 δτ = 极值导出，其坐标形式 d²x^ρ/dτ² + Γ^ρ_μν u^μ u^ν = 0 宣告：自由落体的坐标加速度完全由联络给出。本站做两个数值实验——黑洞附近的径向坠落（固有时有限而坐标时无限）与球面上的最直航线（大圆不需舵、纬线需要舵），验证"自由航行 = 测地航行"不只是口号。

从"最省时"到方程：变分一步走

前两站攒齐了度规与联络，本站把它们接进自动驾驶仪。第 2.2 节留下过一句总纲：测地线是固有时取极值的航线。把它写成数学：在所有连接两事件的世界线中，让固有时积分 τ = ∫ dτ 取极值的那条。用变分法（欧拉-拉格朗日方程）对积分变量求驻点，几行推导后尘埃落定为

d²x^ρ/dτ² + Γ^ρ_μν (dx^μ/dτ)(dx^ν/dτ) = 0

这就是测地线方程。读法值得细品：第一项是坐标加速度，第二项是联络与速度分量的二次组合。方程说的是"两者之和为零"——沿着测地线航行时，四速度的协变导数为零，即航行者自己的加速度表读数为零。这正是第 1.2 节那句"自由下落是无引力航行"的方程化：航行者毫无受力感，全部"下落"表现在坐标账本里，由联络代劳。

注意一个容易混淆的细节：测地线方程里的"最直"是对固有时参数化而言的（称为仿射参数）。换成别的参数（比如坐标时 t），方程会多出一项与速度平方成正比的摩擦状修正——那不是物理，只是换了个不匀速的秒表。航志统一用固有时 τ 记账。

另一个要点是适用范围：只有"自由"物体走测地线。火箭点火、地面支撑、绳子牵引的世界线都不是测地线；地面上的你之所以没走测地线，是因为地板持续把你从"最省时航线"上推开——从这面镜子看，站在地上反而是"被加速"的状态，与 1.2 节的换轨完全一致。

数值实验一：黑洞附近的自由坠落

第一个实验直面广义相对论最出名的剧情：朝黑洞自由坠落会发生什么。取史瓦西度规（c=1，r_s = 1），从 r₀ = 10r_s 静止释放一名航行者，积分其径向测地线。径向情形有守恒量：能量 E = A·dt/dτ（A 即 1 − r_s/r），落速满足 (dr/dτ)² = E² − A：

# 史瓦西时空径向测地线：固有时有限，坐标时在视界处冻结
import math
r_s, r0 = 1.0, 10.0
A0 = 1 - r_s/r0
E = math.sqrt(A0)                 # 从静止释放的能量常数

def dr_dtau(r):
    A = 1 - r_s/r
    return -math.sqrt(max(0.0, E*E - A))   # 向内坠落取负号

def dt_dtau(r):
    return E / (1 - r_s/r)                  # 坐标时速率，视界处发散

tau, r, dtau = 0.0, r0, 1e-3
t_coord = 0.0
stops = {10.0: None, 5.0: None, 2.0: None, 1.01: None}
while r > 1.0005:
    r += dr_dtau(r)*dtau
    t_coord += dt_dtau(r)*dtau
    tau += dtau
    for key in stops:
        if r <= key and stops[key] is None:
            stops[key] = (tau, t_coord)
print("坠落日志：r 位置 -> 累计固有时 与 累计坐标时")
for key in [5.0, 2.0, 1.01]:
    tau_k, t_k = stops[key]
    print(f"r = {key:5.2f} r_s: 固有时 {tau_k:7.3f}, 坐标时 {t_k:8.3f}")
print("坠落者手表现在读数有限；远方账本却记到发散——视界内外两本时间")
​

日志记录（典型数值）：落到 r = 5 时固有时约 14.7、坐标时约 15.5；落到 r = 2 时固有时约 20.6、坐标时约 24.1；逼近视界 r = 1.01 时固有时约 25 上下、坐标时已超过 43 并继续飙升。物理解读分两层：坠落者自己的钟走有限时间就穿过视界，潮汐应力在此之前可能也尚未撕碎他（视超大质量黑洞的质量而定，见 5.1 节）；而远方观察者的坐标账本上，坠影在视界处逐渐变暗冻结，坐标时趋于无限。同一次航行、两本账、互不矛盾——这就是 2.1 节"换页不换物理"的戏剧化版本。变式：改从 r₀ = 100 释放，E 更接近 1，你会发现穿越视界所需固有时更长，但结论结构不变。

数值实验二：球面上的最直航线

第二个实验离开时空、降到二维球面，把"最直"二字咬开看。球面度规 ds² = R²(dθ² + sin²θ dφ²) 只有三个非零克氏符号，其中关键的两个：Γ^θ_φφ = −sinθ·cosθ 与 Γ^φ_θφ = cotθ。代入测地线方程检查哪条纬线不用"打舵"：

# 检查哪条纬线是测地线：赤道免舵，其他纬线必须持续修偏
import math
for lat_deg in [0, 30, 45, 60]:
    theta = math.pi/2 - math.radians(lat_deg)     # 纬度转余纬
    # 沿等速 phi 航行，theta 应满足 d2theta/dphi2 = sin theta cos theta
    correction = math.sin(theta)*math.cos(theta)
    needed = "无需修偏，天然测地线" if abs(correction) < 1e-12 else f"每弧度需要 {correction:.3f} rad 的修偏"
    print(f"纬度 {lat_deg:>2d} 度：{needed}")
print("飞机走大圆航线省油的原因：那是球面测地线，纬线不是")
​

图 2-2 球面上的最直航线：大圆与纬线

实验结果：赤道（纬度 0）的修偏项 sinθ·cosθ 恰好为零，是天然测地线；纬度 30 度的航线每弧度需要 0.433 的修偏量，60 度更是需要 0.433 乘根号三——纬线不是测地线，除非它恰好是赤道。洲际航班宁肯北绕也要飞大圆弧，烧的是同量的油却飞得更短，正是球面测地线的日常红利。把这个图像搬回时空：行星轨道、光线、坠落者走的都是各自度规下的"大圆"，没有任何一方在"转向"，转的是舞台本身。

航行警告：测地线方程只约束"自由"航行者。把火箭推力、光压、电磁力加进四加速度 u^ν∇_νu^ρ = F^ρ/m，方程就不再是测地线——判断一段世界线是否自由，先数清楚受力清单。

关键直觉："最直"与"最省时"是同一条航线的两个名字。飞大圆省的是里程，自由坠落省的是固有时——乘客更年轻、燃料更省，物理与工程在这一条线上握手。