2.3 张量代数与协变导数：颠簸中的读数

2.3 张量代数与协变导数：颠簸中的读数

本节摘要：张量是"带着变换规则出门"的量，保证分量随坐标换页而变、方程形式不变；但在弯曲舞台上，连"求变化率"都要重新发明——普通偏导数不是张量，必须补上联络项，得到协变导数。本站算出克里斯托费尔符号的数值，并验证它在弱场极限下恰好退化为牛顿引力场强：几何与三百年前的力在此完成交接仪式。

颠簸从哪里来

航志记到这里，仪器舱已经装了海图（流形）与标尺（度规）。第三件仪器解决一个隐蔽但致命的问题：在弯曲时空里，"一个张量的变化率"该怎么定义？

平直空间里这不成问题——全空间共用一套基矢量，比较两点的矢量就是把分量相减再求导。弯曲舞台上没有这种便利：每点有自己的局部坐标方向（2.1 的分页），北京点的"东"与拉萨点的"东"根本不是同一个方向。直接把两页账本上的分量相减，就像直接相减两笔不同货币的金额——数字出来了，含义却是错的。

张量：带着变换规则出门的量

先把"量"本身的身份验明。一个量要成为张量，必须同时交验两样东西：一组分量，和一套随坐标变换的换算规则。上标（逆变量）按下述规则换页：新分量等于旧分量乘以新坐标对旧坐标的偏导；下标（协变量）则乘以反向偏导。度规 g_μν 带两个下标，是二阶协变张量；它的"倒数矩阵" g^μν 带两个上标，负责升降指标——V^μ = g^μν V_ν，就像用汇率在两种记法之间换算。

这些规则的价值第 1.4 节已经说清：只要方程两边都是同型张量，等号就在任何一页账本下自动成立——广义协变性从口号变成了机械可查的格式审查。本航志用到头的就是这几样：度规、应力能量张量（第 3 章）、黎曼曲率张量（2.5）以及联络，全在这套规则下运转。

联络：把变化率补成张量

对矢量场 V^ρ 求空间变化率，普通偏导 ∂_μ V^ρ 会漏掉一部分：基矢量本身随位置在变。补丁是把基矢量的变化率也算进去，记作 Γ^ρ_μν，于是协变导数定义为

∇_μ V^ρ = ∂_μ V^ρ + Γ^ρ_μν V^ν

Γ 读作克里斯托费尔符号（简称克氏符号），它不是张量——正因如此，它才能把偏导数的非张量缺口恰好填平，使 ∇_μ V^ρ 成为真张量。这套补丁机制叫联络。好消息是它不必另外假设：在广义相对论采用的"列维-奇维塔联络"下，克氏符号完全由度规的导数组装：

Γ^ρ_μν = ½ g^ρσ (∂_μ g_σν + ∂_ν g_σμ − ∂_σ g_μν)

一句话概括本站的仪器谱系：度规给定，一切自动生成——联络、测地线（下节）、曲率（2.5 节）全是度规的导数组装品。这也是"引力信息全住在度规里"的技术兑现。

数值演习：亲手算出克氏符号

公式念一遍不如算一遍。用史瓦西度规（单位制 c=1，记 A = 1 − r_s/r）做有限差分，把最重要的一个分量 Γ^r_tt 数值算出来并与解析值对账：

# 用有限差分从史瓦西度规算 Gamma 上标 r 下标 t t
import math
r_s, r0, h = 1.0, 10.0, 1e-4     # 取 r_s=1 方便对照，r0=10 处求值

def g_tt(r):                      # 度规分量 g_tt = -(1 - r_s/r)
    return -(1 - r_s/r)
def g_rr(r):                      # g_rr = 1/A
    return 1.0/(1 - r_s/r)

# 公式：Gamma^r_tt = 1/2 g^rr 乘以 负的 d(g_tt) 除 dr
dg_tt = (g_tt(r0+h) - g_tt(r0-h)) / (2*h)         # 中心差分
gamma_r_tt = 0.5 * (1.0/g_rr(r0)) * (-dg_tt)
exact = 0.5 * (1 - r_s/r0) * r_s / r0**2
print(f"数值   Gamma^r_tt = {gamma_r_tt:.8f}")
print(f"解析   Gamma^r_tt = {exact:.8f}")
print(f"相对偏差 {abs(gamma_r_tt-exact)/exact:.2e} —— 公式与差分互证成功")
​

日志记录：数值与解析值吻合到相对偏差 10⁻⁸ 量级，克氏符号的组装规则无误。算出 Γ 有什么用？它是弯曲舞台的"重力加速度表"。一个瞬时静止的粒子，其坐标加速度（下节从测地线方程导出）正是 d²r/dt² = −Γ^r_tt（c=1 单位）。也就是说，引力场强的角色由联络扮演——这在下一场演习里看得最清楚。

弱场汇率下，牛顿引力回来了

把上面演习的单位恢复，让 r_s 是地球的史瓦西半径、r 取地表半径，再与牛顿引力场强直接对账：

# 弱场检验：Gamma^r_tt 乘 c 方 是否等于牛顿 g
G = 6.674e-11; c = 3.0e8
M, R = 5.97e24, 6.371e6
r_s = 2*G*M/c**2

A = 1 - r_s/R
Gamma_r_tt_SI = 0.5 * A * c**2 * r_s / R**2     # 恢复 c 方因子的场强表
g_newton = G*M/R**2
print(f"几何侧读数 Gamma c 方: {Gamma_r_tt_SI:.6f} m/s2")
print(f"牛顿侧读数 GM/R 方  : {g_newton:.6f} m/s2")
print(f"相对差 {abs(Gamma_r_tt_SI-g_newton)/g_newton:.2e} —— 差一个 A 因子，弱场下 A 趋近 1")
​

日志记录：几何读数与牛顿 g 相对差仅 10⁻⁹（恰为 1−A 因子的偏离），这就是著名的弱场极限：r 远大于 r_s 时，联络表逐点退化为牛顿引力场强表，测地线方程退化为牛顿运动方程。三百年前教科书里的 "g = 9.8"，在新语法里是"克里斯托费尔符号乘光速平方"——同一个数字，两套语法。完整推导放在 3.3 节，这里先记结论：广义相对论包含牛顿引力，如同相对论力学包含低速力学。

💡 关键直觉：联络不是度规之外的额外结构，而是"度规导数的组装品"。引力场强对应联络而非曲率——均匀引力场（可被自由下落完全消除的部分）根本不含曲率，这也再次解释了 1.2 节的结论：能消除的是联络效应，消不掉的是曲率（潮汐）。

本站要点

• 张量的身份：分量加变换规则；同型张量构成的方程自动协变——广义协变性成为格式审查；
• 协变导数的动机：基矢量逐点改变，普通偏导不是张量；补上联络项 Γ 后 ∇V 才是真张量；
• 克氏符号由度规生成：Γ = ½ g^ρσ（度规导数组合），无需额外假设——"度规给定，一切生成"；
• 数值互证：有限差分与解析式算出的 Γ^r_tt 偏差 10⁻⁸ 量级，公式可信；
• 交接仪式：弱场下 Γ 乘 c² 逐点等于牛顿 g（地表相对差 10⁻⁹），引力场强的几何身份是联络；
• 下一站：有了联络，自由航迹的方程只差一步变分——测地线开箱。