跳转至

广义相对论基础

广义相对论基础(General Relativity Basics)

本章目标:从「等效原理」出发,建立广义相对论的最小数学工具箱(张量、度规、协变导数、曲率),给出爱因斯坦场方程的常见推导动机与弱场极限,并以 Schwarzschild 时空为主线串起经典检验(引力红移、光线偏折、近日点进动、Shapiro 延迟)与典型考题.

为方便计算,本文会在不同小节切换单位:

  • 几何单位c=G=1(常用于黑洞/测地线推导)
  • SI 单位:保留 c,G(适合数值估算与实验量级)
符号约定(请先看)
1
2
3
4
- 指标:希腊字母 $\mu,\nu,\dots=0,1,2,3$,其中 $x^0=ct$.  
- 度规号:本文默认 **$(-,+,+,+)$**(时间分量为负).  
- 爱因斯坦求和约定:重复上下指标自动求和.  
- 4-速度:$u^\mu=\dfrac{dx^\mu}{d\tau}$,其中 $\tau$ 为固有时.

0. 从「引力是什么」到「时空几何」

牛顿引力把引力看作一种作用于物体的力:

F=mg,g=4πGρ.

然而两个关键事实迫使我们升级理论:

  1. 惯性质量 = 引力质量(实验精度极高)
  2. 光也会被引力影响(偏折、延迟、红移),而光「没有静止质量」

广义相对论的核心观点是:

引力不是一种「额外的力」,而是 时空几何(度规)的体现;自由落体沿着时空中的测地线 运动.


1. 等效原理:广义相对论的出发点

1.1 弱等效原理(WEP)

所有自由落体在同一重力场中具有相同的加速度,与其内部组成无关.

这等价于「惯性质量 mi 与引力质量 mg 相等」,从而

mia=mgg  a=g.

1.2 爱因斯坦等效原理(EEP)

在任意一点附近,总能选择一个局部惯性系,使得该点附近的物理定律与狭义相对论相同.

直观表述:

  • 封闭电梯里做实验:无法区分「处于均匀重力场静止」与「在无重力的空间做匀加速运动」.

1.3 用等效原理推导引力红移(最经典的 1 页推导)

考虑一部在深空中以加速度 a 匀加速上升的电梯.电梯底部发出频率 νemit 的光子向上,顶部接收.

顶部在光子飞行时间 Δth/c 内获得速度增量 ΔvaΔt=ah/c

对顶部而言,接收到的光发生相对论多普勒效应(弱速):

νrecνemit1Δvc=1ahc2.

由等效原理,把匀加速与匀重力 g=a 等价,得到引力场中「向上爬升」会 红移

 Δννghc2 (ghc2)

同理可得时间膨胀:上方的钟走得更快.


2. 最小数学工具箱:度规、张量与协变导数

2.1 间隔与度规(metric)

狭义相对论中 Minkowski 间隔:

ds2=ημνdxμdxν=c2dt2+dx2+dy2+dz2.

广义相对论把 ημν 推广为随位置变化的 gμν(x)

 ds2=gμν(x)dxμdxν 

物理意义:度规决定「距离/时间」的测量方式,因此也决定自由粒子的运动.

2.2 固有时与 4 - 速度

对时类轨迹(可作为粒子世界线),定义固有时:

dτ=1cds2.

因此

uμ=dxμdτ,uμuμ=c2.

2.3 指标升降

vμ=gμνvν,vμ=gμνvν,

其中 gμνgμν 的矩阵逆:gμαgαν=δ νμ

2.4 为什么需要协变导数?

在曲面/曲时空中,普通偏导 μ 作用在张量上不再保持「张量变换律」.需要引入带「连接」的导数:

μvν=μvν+Γ μλνvλ.

其中 Γ μλν 为 Christoffel 符号(Levi-Civita 连接).

2.5 Christoffel 符号的公式

在无挠、度规相容(αgμν=0)的假设下:

 Γ μνρ=12gρσ(μgνσ+νgμσσgμν) 
小提醒:Christoffel 不是张量
1
它依赖坐标系选择(可以在某点选取局部惯性系使其为 0),但它组合成的曲率张量是坐标不变量.

3. 测地线:自由落体的方程

3.1 从「极值固有时」推导测地线方程

自由粒子在两事件间走「固有时最大」(或作用量极值)的路径:

S=mcdτ=mcgμνx˙μx˙νdλ,

其中 x˙μ=dxμdλλ 是参数.

对该拉格朗日量做变分(略去常数并选取仿射参数),可得欧拉 - 拉格朗日方程最终化为:

 d2xρdτ2+Γ μνρdxμdτdxνdτ=0 

这就是测地线方程:自由落体「没有外力」,但坐标加速度可不为零.

3.2 牛顿极限:测地线如何还原 a=g

在弱场、低速下,度规可写成

g00(1+2Φ/c2),|Φ|c2.

测地线的空间分量(只保留一阶小量)给出

d2xidt2iΦ,

即恢复牛顿引力势 Φ 的运动方程.


4. 曲率:引力的「强度」如何编码

4.1 曲率来自「绕一圈不回原位」

在平直空间,把向量平移绕一圈会回到原向量;在曲空间则会发生偏转.数学上表现为协变导数不对易:

(μννμ)vρ=R σμνρvσ.

这里 R σμνρ 是 Riemann 曲率张量.

4.2 Riemann、Ricci 与标量曲率

从 Christoffel 可写出:

R σμνρ=μΓ νσρνΓ μσρ+Γ μλρΓ νσλΓ νλρΓ μσλ.

收缩得到 Ricci 张量与标量曲率:

Rμν=R μρνρ,R=gμνRμν.

4.3 爱因斯坦张量与 Bianchi 恒等式

定义

Gμν=Rμν12gμνR.

它满足重要恒等式(几何上「守恒」):

 μGμν=0 

这将迫使右边的物质源也满足能动量守恒:μTμν=0


5. 爱因斯坦场方程:几何 = 物质

5.1 形式与常数

广义相对论最核心的方程:

 Gμν=8πGc4Tμν 

其中 Tμν 为能动量张量,描述能量密度、动量流、压强/剪应力等.

5.2 为什么是 Gμν?(最常用的「推导动机」)

我们希望左边是:

  1. 由度规及其一、二阶导数组成(对应「局域」引力)
  2. 对称二阶张量(与 Tμν 匹配)
  3. 满足协变守恒 μ(左边)=0
  4. 弱场极限还原泊松方程 2Φ=4πGρ

在这些条件下,最自然且最简单的选择就是 Gμν

更系统的推导:Hilbert–Einstein 作用量
1
2
3
4
5
6
7
取总作用量  
$$  
S=\dfrac{c^3}{16\pi G}\int R\sqrt{-g}\,d^4x+S_{\text{matter}}.  
$$  
对 $g_{\mu\nu}$ 变分可得 $G_{\mu\nu}=\frac{8\pi G}{c^4}T_{\mu\nu}$.

(完整变分细节较长,竞赛/入门通常只要求理解「为何是 $R$、为何要 $\sqrt{-g}$、为何能得到守恒」.)

5.3 常见物质:理想流体的 Tμν

宇宙学与恒星结构常用「理想流体」:

Tμν=(ρc2+p)uμuν+pgμν,

其中 ρ 为质量密度(或能量密度除以 c2),p 为各向同性压强.


6. 弱场近似与线性化引力(选学但很实用)

gμν=ημν+hμν,|hμν|1.

在线性近似与合适规范(如 Lorenz 规范)下,场方程会变成「波动方程」形式,真空中存在引力波解.

6.1 线性化方程的标准形(了解即可)

定义迹反转(trace-reversed)扰动

h¯μνhμν12ημνh,hηαβhαβ.

在 Lorenz 规范

μh¯μν=0

下,Einstein 方程线性化为

 h¯μν=16πGc4Tμν (1c2t2+2)

真空区 Tμν=0 给出波动方程:h¯μν=0

6.2 引力波最小图景:TT 规范与应变

对沿 +z 传播的平面引力波,可在横向 - 无迹(TT)规范下写成

ds2c2dt2+(1+h+)dx2+(1h+)dy2+2h×dxdy+dz2,

其中 h+,h× 为两种偏振.它们对应干涉仪(LIGO 类)测到的 应变 hΔL/L

对于入门最关键的是:在静态弱场、低速下

g00(1+2Φ/c2),2Φ=4πGρ.

因此广义相对论与牛顿引力在日常尺度高度一致,只在精密实验/强引力处显著.


7. Schwarzschild 时空:最重要的精确解

7.1 度规(真空、球对称、静态)

对质量 M 的球对称真空外部,Einstein 方程给出 Schwarzschild 解:

ds2=(12GMrc2)c2dt2+(12GMrc2)1dr2+r2(dθ2+sin2θdϕ2)

定义 Schwarzschild 半径

rs=2GMc2.

r=rs 时,gtt 变号、grr 发散(在该坐标系下).这对应事件视界;更好的坐标(Eddington–Finkelstein、Kruskal)能消除坐标奇点.

7.2 静止观测者的引力时间膨胀

对固定 r,θ,ϕ 的静止钟,dr=dθ=dϕ=0,有

dτ=1rsrdt.

因此远处观察者看到近处的钟变慢.

Schwarzschild 引力时间膨胀因子

7.3 能量与角动量守恒(测地线的第一积分)

Schwarzschild 度规对 tϕ 不显含,导致守恒量:

E(1rsr)c2dtdτ,Lr2dϕdτ.

在赤道面 θ=π/2 可将测地线化为「一维有效势」问题.

7.4 有效势(时类测地线)

在几何单位 G=c=1 下,rs=2M,并可写成

(drdτ)2+Veff(r)=E2,

其中(典型形式)

Veff(r)=(12Mr)(1+L2r2).

Schwarzschild 有效势示意

7.5 光子球与 ISCO(竞赛很常考的两个半径)

在几何单位 G=c=1 的 Schwarzschild 时空:

  • 光子球(photon sphere):不稳定的圆形光子轨道 $$ \boxed{ r_{\text{ph}}=3M=\dfrac{3}{2}r_s } $$
  • 最内稳定圆轨道(ISCO,时类):稳定圆轨道存在的内边界 $$ \boxed{ r_{\text{ISCO}}=6M=3r_s } $$

直观理解:越靠近黑洞,时空曲率越强,有效势的「势阱」结构改变,稳定圆轨道会消失.


8. 经典检验与常用结果(入门必会)

8.1 引力红移(精确形式:Schwarzschild)

两静止观测者分别在 r1(发射)与 r2(接收,通常更远)处:

ν2ν1=1rs/r11rs/r2

r2

ννr=1rsr.

8.2 光线偏折(弱场近似)

对冲量参数(impact parameter)b,弱场一阶:

 α4GMbc2 

弱场光线偏折角

8.3 近日点进动(弱场、近开普勒轨道)

椭圆轨道半长轴 a、偏心率 e,每一周额外进动角

 Δϕ6πGMa(1e2)c2 

8.4 Shapiro 时间延迟(选学)

雷达信号掠过质量 M 的天体会产生额外往返时间延迟.对「从 r1 发出、掠过最近距离 b、到 r2 接收」的几何,弱场近似常见写法为

 ΔtShapiro2GMc3ln(4r1r2b2) 

(不同教材对几何定义略有差异,但共同点是:它随 GM/c3 线性增长,并含对数项.)

8.5 宇宙学一瞥:FRW 度规与 Friedmann 方程(选读)

在「大尺度均匀各向同性」的宇宙假设下,度规可写为 FRW 形式:

ds2=c2dt2+a(t)2[dr21kr2+r2(dθ2+sin2θdϕ2)],

其中 a(t) 为尺度因子,k=0,±1 表示空间曲率符号.

代入 Einstein 方程并取物质为理想流体,可得 Friedmann 方程(最常用版本):

 (a˙a)2=8πG3ρkc2a2+Λc23 

以及加速度方程

 a¨a=4πG3(ρ+3pc2)+Λc23 

并由 μTμν=0 得到连续性方程

ρ˙+3a˙a(ρ+pc2)=0.

9. 例题(含全过程推导与计算)

例题 1:用弱场近似估算地球表面与高空的引力红移

已知地球半径 R6.37×106m,质量 M5.97×1024kg.在高度 hR 处的钟与地面钟相比,频率相对变化 Δν/ν 约为多少?

解:

弱场中

ΔννΔΦc2.

地球引力势(取无穷远为 0)为 Φ(r)=GM/r,因此

ΔΦ=Φ(R+h)Φ(R)=GMR+h+GMR.

hR,用展开

1R+h1R(1hR),

ΔΦGMR(1hR)+GMR=GMR2h=gh.

所以

 Δννghc2 

若取 h=1000mg9.8m/s2,则

|Δνν|9.8×1000(3.0×108)21.1×1013.

** 结论:** 千米量级高度差带来 1013 量级的频移(原子钟可测).


例题 2:Schwarzschild 时空中的静止钟差(精确公式)

在 Schwarzschild 外部时空中,两只静止钟分别位于 r1r2r2>r1>rs).若远处坐标时间过了 Δt,两钟的固有时分别为多少?两者比值?

解:

静止钟有 dr=dθ=dϕ=0,于是

dτ=1rsrdt.

积分得到

Δτ(r)=1rsrΔt.

所以

 Δτ1Δτ2=1rs/r11rs/r2 

与引力红移公式一致(频率与时间尺度互为倒数).


例题 3:弱场光线偏折(给出标准推导骨架)

证明在弱场近似下,掠过质量 M 的光线偏折角为 α4GMbc2

解题思路骨架:

  1. 弱场近似度规(各向同性形式)可写为
ds2(1+2Φ/c2)c2dt2+(12Φ/c2)(dx2+dy2+dz2),Φ=GMr.
  1. 对光线 ds2=0,可得到等效「折射率」
n(r)12Φc2=1+2GMrc2.
  1. 用几何光学近似(费马原理)
δndl=0

把光线当作在介质中传播的射线. 4. 在小偏折近似中,偏折角约为

α+b(2GMc2x2+b2)dx=4GMbc2.

例题 4:近日点进动公式(从有效势到近似解)

给出椭圆轨道的相对论修正导致的每周进动

Δϕ6πGMa(1e2)c2.

推导要点:

  1. 对 Schwarzschild 测地线,可得到轨道方程的近似形式(令 u=1/r
d2udϕ2+u=GML2+3GMc2u2.
  1. 3GMc2u2 视为小扰动,设开普勒解
u0=GML2(1+ecosϕ).
  1. 扰动会把角频率从 1 改成 1δ,使得一周后相位多走
Δϕ2πδ=6πGMa(1e2)c2.

(竞赛层面通常要求记结论并能解释「为何是小扰动导致的频率偏移」.)


例题 5:求 Schwarzschild 的光子球半径与 ISCO 半径

在几何单位 G=c=1 下,证明 Schwarzschild 时空中光子的圆轨道半径为 r=3M,时类最内稳定圆轨道半径为 r=6M

解(思路 + 关键步骤):

  1. 对赤道面运动,利用守恒量可把径向运动写成
(drdλ)2+Veff(r)=E2.
  1. 圆轨道条件r=常数 意味着
drdλ=0Veff(r)=E2,

且要「保持在圆轨道上」还需

dVeffdr=0.
  1. ** 光子(零质量)** 的有效势(形式上)为
Veff(null)(r)=(12Mr)L2r2.

求导并令 dV/dr=0

r=3M.
  1. ** 时类(有质量)** 圆轨道的稳定性由二阶导数决定:稳定要求 d2Veffdr2>0. 对 Schwarzschild,稳定圆轨道存在于 r>6M;临界点为
rISCO=6M.

例题 6:径向自由落体到达视界需要多少固有时?坐标时间会怎样?

考虑从 r=r0(静止释放)径向下落到 r=rs.定性说明:

  • 自由落体者的固有时 τ 是否有限?
  • 远处静止观察者的坐标时间 t 是否有限?

解(定性 + 关键结论):

  1. 对下落者,τ 是沿世界线积分得到的真实物理时间.
  2. 在 Schwarzschild 坐标下,t 在靠近 rs 时会出现对数发散(这是坐标效应).
  3. 结论:

  4. 固有时有限:自由落体者在有限 τ 内穿过视界.

  5. 坐标时间发散:远处观察者用 t 描述会看到物体「永远接近但不穿过」视界.

(若需要具体积分表达式,可在后续补充「径向测地线第一积分」并做积分计算.)


例题 7:掠日光线偏折角的数量级估算

bR7.0×108mMM2.0×1030kg

解:

α4GMbc2.

代入 G6.67×1011c3.0×108,得

α4×6.67×1011×2.0×10307.0×108×(3.0×108)28.5×106rad.

换算为角秒(1rad206265arcsec):

α1.75arcsec.

这就是经典的「掠日偏折约 1.75 角秒」.


例题 8:水星近日点进动的数量级估算

已知水星轨道参数(可查表):半长轴 a5.79×1010m,偏心率 e0.206,太阳质量 MM.估算每周进动角.

解:

Δϕ6πGMa(1e2)c2.

代入数值可得每一周(每绕太阳一圈)约 5×107rad 量级,对应约 0.1arcsec 量级.累积到每世纪(约 415 圈)即得到著名的 43arcsec/century


10. 习题(含提示/答案要点)

10.1 基础题

  1. 证明:在局部惯性系中,某点可令 Γ μνρ=0,但一般不能令 σgμν=0 同时在邻域处处为 0.
  2. 设弱场度规 g00=(1+2Φ/c2),证明静止钟的频移满足 Δν/νΔΦ/c2
  3. 写出理想流体 Tμν 并解释各项物理意义.

10.2 计算题(更贴近竞赛/考题)

  1. 黑洞视界尺度估算:太阳质量 M2.0×1030kg,求 rs 的数量级.

  2. 答案要点:rs=2GM/c23km

  3. 引力时间膨胀:若某天体 rs=10km,在 r=40km 静止的钟相对于无穷远的钟慢多少?

  4. 提示:dτ/dt=1rs/r

  5. 光偏折角数值:太阳半径 R7.0×108m,质量 M,估算掠日光线偏折角.

  6. 提示:α4GM/(bc2),取 bR

  7. 光子球数值:对某黑洞 rs=30km,求光子球半径 rph 与 ISCO 半径 rISCO

  8. 答案要点:rph=32rsrISCO=3rs

  9. 引力波应变量级:若某引力波通过地球时的应变 h1021,LIGO 级别的臂长 L4km,估算长度变化 ΔL

  10. 答案要点:ΔLhL4×1018m


11. 用 Python 复现/绘制示意图(可选)

本文插图由 Python 生成并保存在 docs/modern/images/

  • schwarzschild_time_dilation_factor.png
  • schwarzschild_effective_potential.png
  • light_deflection_weak_field.png
示例代码(生成弱场光偏折图)
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
import numpy as np
import matplotlib.pyplot as plt

M = 1.0
b = np.linspace(4.0, 80.0, 800)
alpha = 4 * M / b
plt.plot(b, alpha)
plt.xlabel(r"$b/M$")
plt.ylabel(r"$\alpha\;\mathrm{(rad)}$")
plt.title(r"弱场近似:$\alpha \approx 4M/b$")
plt.grid(True, alpha=0.25)
plt.show()

参考与延伸阅读

  • Misner, Thorne, Wheeler,Gravitation(体系最全,但很厚)
  • Schutz,A First Course in General Relativity(入门友好)
  • Carroll,Spacetime and Geometry(现代写法,推导清晰)


评论