摘要:广义相对论把引力从力的名单除名,改记为时空曲率:物质能量通过爱因斯坦场方程决定洛伦兹度量的曲率,自由粒子沿测地线滑行。本节复盘从等效原理到场方程的推导主线,用施瓦西解核算三大经典检验与 GPS 修正量级,并以黑洞与引力波收束,完成本册仪器在物理工地的验收。
封闭电梯里的观察者无法区分「静止于引力场」与「无引力空间中匀加速」——这是等效原理的通俗版本,爱因斯坦称之为一生最幸福的想法。几何转译:引力效应可以通过选择适当的局部参照系完全消除,恰似弯曲流形上每点都有局部平直的坐标卡(第一章)。由此定调:引力若要几何化,几何必须是洛伦兹的(第五章第三节),而「受力运动」应改写为「无力的测地线运动」(第三章)。
场方程的推导主线三步。第一步定左端:场方程应是关于度量的二阶方程且自动满足守恒——第 3.2 节的比安基恒等式恰好保证爱因斯坦张量 G_{ab} = Ric_{ab} - \tfrac12 Scal\, g_{ab} 无散,别无分店。第二步定右端:物质能量动量张量 T_{ab} 是广义协变的守恒量,两端成比例即自动守恒。第三步定常数:牛顿极限(弱场低速)对表万有引力公式,解出比例常数 \kappa = 8\pi(几何单位制)。合并成一行:
左边是曲率,右边是物质,左边随右边流动:第 3.2 节预装配的部件在此正式上岗。惠勒的格言式总结——「物质告诉时空如何弯曲,时空告诉物质如何运动」——的后半句正是第三章的测地线方程。
场方程最简单的非平庸解由施瓦西在战争年代给出:球对称真空(T_{ab} = 0)时度量为
先按本册流程验货:r 大时退化为球坐标下的闵可夫斯基度量(前两项号差一负一正、后两项正是第 1.4 节手推的球面度量)——量具交接无断裂;真空场方程 G_{ab} = 0 可由它直接验算(代入克氏符号与曲率公式,机器核算见下方代码);r = 2M 处分量发散,经第 4.1 节的临床分辨,那是坐标账本到期(换自由下落坐标即痊愈),而 r = 0 是真奇点。
三大经典检验的核算(只报结论与领头项,推导见各检验对应章的作业):一、行星进动——在施瓦西时空中解测地线方程,近日点每圈进动 \tfrac{6\pi M}{a(1-e^2)}(a 半长轴、e 偏心率),水星数值为每百年四十三角秒,与观测天衣无缝,这是广义相对论的首胜。二、光线偏折——光子沿类光测地线经过太阳边缘,偏折角 \tfrac{4M}{R_\odot} 约一点七五角秒,一九一九年日食观测证实,倍的因子(等效原理只给一半)恰来自空间曲率部分——第 3.2 节「里奇之外还有截面」的实测注脚。三、引力红移——静止于引力场中的钟走得慢,r 处与无穷远的钟速比为 \sqrt{1 - \tfrac{2M}{r}},第 5.3 节「类时测地线原时最长」的静站版本。
GPS 核算:导航卫星钟每天比地面钟快约三十八微秒(狭义速度效应慢七微秒 + 引力势效应快四十五微秒,后者正是红移公式在地球质量下的读数)。若不修正,光速乘三十八微秒即每天约十一公里的定位漂移——引力时间效应是 GPS 的生命线。这个量级核算也提醒我们:黎曼几何不是奢侈品,是把卫星导航误差控制在米级以内的日常必需品。

from sympy import symbols, sin, Matrix, simplify, Rational r, th, ph, t, M = symbols('r theta varphi t M', positive=True) f = 1 - Rational(2)*M/r coords = [t, r, th, ph] g = Matrix([[-f, 0, 0, 0], [0, 1/f, 0, 0], [0, 0, r**2, 0], [0, 0, 0, r**2*sin(th)**2]]) ginv = g.inv() # 克氏符号(第 2.2 节公式) G = {(k, i, j): simplify(sum(ginv[k, l] * (g[l, i].diff(coords[j]) + g[l, j].diff(coords[i]) - g[i, j].diff(coords[l])) for l in range(4)) / 2) for k in range(4) for i in range(4) for j in range(4)} # 里奇张量 R_kj(第 3.2 节缩并公式,逐项按定义累加) def Ricci(k, j): expr = 0 for m in range(4): # - d_m Gamma^m_{kj} expr -= G[(m, k, j)].diff(coords[m]) for l in range(4): # + d_l Gamma^l_{jk} expr += G[(l, j, k)].diff(coords[l]) for m in range(4): # 乘积差两项 expr += G[(l, m, k)] * G[(m, j, l)] expr -= G[(l, m, j)] * G[(m, k, l)] return simplify(expr) Ric_m = Matrix(4, 4, lambda k, j: Ricci(k, j)) Scal = simplify(sum(ginv[i, i] * Ric_m[i, i] for i in range(4))) G_ein = Ric_m - Scal * g / 2 # 爱因斯坦张量 print([simplify(G_ein[i, j]) for i in range(4) for j in range(4) if i == j]) # 输出: [0, 0, 0, 0] —— 真空场方程 G = 0 成立
符号计算确认爱因斯坦张量的四个对角分量全部为零(非对角分量由对称性同样为零),施瓦西度量确实解开真空场方程。手写循环的目的是让读者看清汤是如何煮出来的——每个分量都是第 2.2、3.1、3.2 节公式的机械代入;正式工程应改用专门的张量计算系统,第七章工具链会给出选型。
物理工地的验收结论:第一章(坐标卡与图册)支撑了任意坐标系下的协变表述;第二章(克氏符号)是行星轨道方程的全部「力项」来源——它取代了牛顿引力;第三章(曲率与测地线)承担三大检验的计算主体;第四章(奇点与完备性)给出黑洞视界与奇点的分辨标准;第五章(洛伦兹号差)是因果结构与红移的母体。没有一件仪器缺席,也没有一件仪器需要临时发明——这正是黎曼当年演讲的预言在物理侧的兑现:几何先于物理准备好了一切语法。