广义相对论基础
广义相对论基础(General Relativity Basics)
本章目标:从「等效原理」出发,建立广义相对论的最小数学工具箱(张量、度规、协变导数、曲率),给出爱因斯坦场方程的常见推导动机与弱场极限,并以 Schwarzschild 时空为主线串起经典检验(引力红移、光线偏折、近日点进动、Shapiro 延迟)与典型考题.
为方便计算,本文会在不同小节切换单位:
- 几何单位:𝑐 =𝐺 =1
(常用于黑洞/测地线推导) - SI 单位:保留 𝑐,𝐺
(适合数值估算与实验量级)
符号约定(请先看)
| - 指标:希腊字母 $\mu,\nu,\dots=0,1,2,3$,其中 $x^0=ct$.
- 度规号:本文默认 **$(-,+,+,+)$**(时间分量为负).
- 爱因斯坦求和约定:重复上下指标自动求和.
- 4-速度:$u^\mu=\dfrac{dx^\mu}{d\tau}$,其中 $\tau$ 为固有时.
|
0. 从「引力是什么」到「时空几何」
牛顿引力把引力看作一种作用于物体的力:
𝐅=𝑚𝐠,∇⋅𝐠=−4𝜋𝐺𝜌.
然而两个关键事实迫使我们升级理论:
- 惯性质量 = 引力质量(实验精度极高)
- 光也会被引力影响(偏折、延迟、红移),而光「没有静止质量」
广义相对论的核心观点是:
引力不是一种「额外的力」,而是 时空几何(度规)的体现;自由落体沿着时空中的测地线 运动.
1. 等效原理:广义相对论的出发点
1.1 弱等效原理(WEP)
所有自由落体在同一重力场中具有相同的加速度,与其内部组成无关.
这等价于「惯性质量 𝑚𝑖
与引力质量 𝑚𝑔
相等」,从而
𝑚𝑖𝐚=𝑚𝑔𝐠 ⇒ 𝐚=𝐠.
1.2 爱因斯坦等效原理(EEP)
在任意一点附近,总能选择一个局部惯性系,使得该点附近的物理定律与狭义相对论相同.
直观表述:
- 封闭电梯里做实验:无法区分「处于均匀重力场静止」与「在无重力的空间做匀加速运动」.
1.3 用等效原理推导引力红移(最经典的 1 页推导)
考虑一部在深空中以加速度 𝑎
匀加速上升的电梯.电梯底部发出频率 𝜈emit
的光子向上,顶部接收.
顶部在光子飞行时间 Δ𝑡 ≈ℎ/𝑐
内获得速度增量 Δ𝑣 ≈𝑎Δ𝑡 =𝑎ℎ/𝑐
.
对顶部而言,接收到的光发生相对论多普勒效应(弱速):
𝜈rec𝜈emit≈1−Δ𝑣𝑐=1−𝑎ℎ𝑐2.
由等效原理,把匀加速与匀重力 𝑔 =𝑎
等价,得到引力场中「向上爬升」会 红移:
Δ𝜈𝜈≈−𝑔ℎ𝑐2 (𝑔ℎ≪𝑐2)
同理可得时间膨胀:上方的钟走得更快.
2. 最小数学工具箱:度规、张量与协变导数
2.1 间隔与度规(metric)
狭义相对论中 Minkowski 间隔:
𝑑𝑠2=𝜂𝜇𝜈𝑑𝑥𝜇𝑑𝑥𝜈=−𝑐2𝑑𝑡2+𝑑𝑥2+𝑑𝑦2+𝑑𝑧2.
广义相对论把 𝜂𝜇𝜈
推广为随位置变化的 𝑔𝜇𝜈(𝑥)
:
𝑑𝑠2=𝑔𝜇𝜈(𝑥)𝑑𝑥𝜇𝑑𝑥𝜈 
物理意义:度规决定「距离/时间」的测量方式,因此也决定自由粒子的运动.
2.2 固有时与 4 - 速度
对时类轨迹(可作为粒子世界线),定义固有时:
𝑑𝜏=1𝑐√−𝑑𝑠2.
因此
𝑢𝜇=𝑑𝑥𝜇𝑑𝜏,𝑢𝜇𝑢𝜇=−𝑐2.
2.3 指标升降
𝑣𝜇=𝑔𝜇𝜈𝑣𝜈,𝑣𝜇=𝑔𝜇𝜈𝑣𝜈,
其中 𝑔𝜇𝜈
是 𝑔𝜇𝜈
的矩阵逆:𝑔𝜇𝛼𝑔𝛼𝜈 =𝛿𝜇 𝜈
.
2.4 为什么需要协变导数?
在曲面/曲时空中,普通偏导 𝜕𝜇
作用在张量上不再保持「张量变换律」.需要引入带「连接」的导数:
∇𝜇𝑣𝜈=𝜕𝜇𝑣𝜈+Γ𝜈 𝜇𝜆𝑣𝜆.
其中 Γ𝜈 𝜇𝜆
为 Christoffel 符号(Levi-Civita 连接).
2.5 Christoffel 符号的公式
在无挠、度规相容(∇𝛼𝑔𝜇𝜈 =0
)的假设下:
Γ𝜌 𝜇𝜈=12𝑔𝜌𝜎(𝜕𝜇𝑔𝜈𝜎+𝜕𝜈𝑔𝜇𝜎−𝜕𝜎𝑔𝜇𝜈) 
小提醒:Christoffel 不是张量
| 它依赖坐标系选择(可以在某点选取局部惯性系使其为 0),但它组合成的曲率张量是坐标不变量.
|
3. 测地线:自由落体的方程
3.1 从「极值固有时」推导测地线方程
自由粒子在两事件间走「固有时最大」(或作用量极值)的路径:
𝑆=−𝑚𝑐∫𝑑𝜏=−𝑚𝑐∫√−𝑔𝜇𝜈˙𝑥𝜇˙𝑥𝜈𝑑𝜆,
其中 ˙𝑥𝜇 =𝑑𝑥𝜇𝑑𝜆
,𝜆
是参数.
对该拉格朗日量做变分(略去常数并选取仿射参数),可得欧拉 - 拉格朗日方程最终化为:
𝑑2𝑥𝜌𝑑𝜏2+Γ𝜌 𝜇𝜈𝑑𝑥𝜇𝑑𝜏𝑑𝑥𝜈𝑑𝜏=0 
这就是测地线方程:自由落体「没有外力」,但坐标加速度可不为零.
3.2 牛顿极限:测地线如何还原 𝐚 =𝐠
在弱场、低速下,度规可写成
𝑔00≈−(1+2Φ/𝑐2),|Φ|≪𝑐2.
测地线的空间分量(只保留一阶小量)给出
𝑑2𝑥𝑖𝑑𝑡2≈−𝜕𝑖Φ,
即恢复牛顿引力势 Φ
的运动方程.
4. 曲率:引力的「强度」如何编码
4.1 曲率来自「绕一圈不回原位」
在平直空间,把向量平移绕一圈会回到原向量;在曲空间则会发生偏转.数学上表现为协变导数不对易:
(∇𝜇∇𝜈−∇𝜈∇𝜇)𝑣𝜌=𝑅𝜌 𝜎𝜇𝜈𝑣𝜎.
这里 𝑅𝜌 𝜎𝜇𝜈
是 Riemann 曲率张量.
4.2 Riemann、Ricci 与标量曲率
从 Christoffel 可写出:
𝑅𝜌 𝜎𝜇𝜈=𝜕𝜇Γ𝜌 𝜈𝜎−𝜕𝜈Γ𝜌 𝜇𝜎+Γ𝜌 𝜇𝜆Γ𝜆 𝜈𝜎−Γ𝜌 𝜈𝜆Γ𝜆 𝜇𝜎.
收缩得到 Ricci 张量与标量曲率:
𝑅𝜇𝜈=𝑅𝜌 𝜇𝜌𝜈,𝑅=𝑔𝜇𝜈𝑅𝜇𝜈.
4.3 爱因斯坦张量与 Bianchi 恒等式
定义
𝐺𝜇𝜈=𝑅𝜇𝜈−12𝑔𝜇𝜈𝑅.
它满足重要恒等式(几何上「守恒」):
∇𝜇𝐺𝜇𝜈=0 
这将迫使右边的物质源也满足能动量守恒:∇𝜇𝑇𝜇𝜈 =0
.
5. 爱因斯坦场方程:几何 = 物质
5.1 形式与常数
广义相对论最核心的方程:
𝐺𝜇𝜈=8𝜋𝐺𝑐4𝑇𝜇𝜈 
其中 𝑇𝜇𝜈
为能动量张量,描述能量密度、动量流、压强/剪应力等.
5.2 为什么是 𝐺𝜇𝜈
?(最常用的「推导动机」)
我们希望左边是:
- 由度规及其一、二阶导数组成(对应「局域」引力)
- 对称二阶张量(与 𝑇𝜇𝜈
匹配) - 满足协变守恒 ∇𝜇(左边) =0

- 弱场极限还原泊松方程 ∇2Φ =4𝜋𝐺𝜌

在这些条件下,最自然且最简单的选择就是 𝐺𝜇𝜈
.
更系统的推导:Hilbert–Einstein 作用量
| 取总作用量
$$
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 常见物质:理想流体的 𝑇𝜇𝜈
宇宙学与恒星结构常用「理想流体」:
𝑇𝜇𝜈=(𝜌𝑐2+𝑝)𝑢𝜇𝑢𝜈+𝑝𝑔𝜇𝜈,
其中 𝜌
为质量密度(或能量密度除以 𝑐2
),𝑝
为各向同性压强.
6. 弱场近似与线性化引力(选学但很实用)
设
𝑔𝜇𝜈=𝜂𝜇𝜈+ℎ𝜇𝜈,|ℎ𝜇𝜈|≪1.
在线性近似与合适规范(如 Lorenz 规范)下,场方程会变成「波动方程」形式,真空中存在引力波解.
6.1 线性化方程的标准形(了解即可)
定义迹反转(trace-reversed)扰动
¯ℎ𝜇𝜈≡ℎ𝜇𝜈−12𝜂𝜇𝜈ℎ,ℎ≡𝜂𝛼𝛽ℎ𝛼𝛽.
在 Lorenz 规范
𝜕𝜇¯ℎ𝜇𝜈=0
下,Einstein 方程线性化为
◻¯ℎ𝜇𝜈=−16𝜋𝐺𝑐4𝑇𝜇𝜈 (◻≡−1𝑐2𝜕2𝑡+∇2)
真空区 𝑇𝜇𝜈 =0
给出波动方程:◻¯ℎ𝜇𝜈 =0
.
6.2 引力波最小图景:TT 规范与应变
对沿 +𝑧
传播的平面引力波,可在横向 - 无迹(TT)规范下写成
𝑑𝑠2≈−𝑐2𝑑𝑡2+(1+ℎ+)𝑑𝑥2+(1−ℎ+)𝑑𝑦2+2ℎ×𝑑𝑥𝑑𝑦+𝑑𝑧2,
其中 ℎ+,ℎ×
为两种偏振.它们对应干涉仪(LIGO 类)测到的 应变 ℎ ∼Δ𝐿/𝐿
.
对于入门最关键的是:在静态弱场、低速下
𝑔00≈−(1+2Φ/𝑐2),∇2Φ=4𝜋𝐺𝜌.
因此广义相对论与牛顿引力在日常尺度高度一致,只在精密实验/强引力处显著.
7. Schwarzschild 时空:最重要的精确解
7.1 度规(真空、球对称、静态)
对质量 𝑀
的球对称真空外部,Einstein 方程给出 Schwarzschild 解:
𝑑𝑠2=−(1−2𝐺𝑀𝑟𝑐2)𝑐2𝑑𝑡2+(1−2𝐺𝑀𝑟𝑐2)−1𝑑𝑟2+𝑟2(𝑑𝜃2+sin2𝜃𝑑𝜙2)
定义 Schwarzschild 半径
𝑟𝑠=2𝐺𝑀𝑐2.
当 𝑟 =𝑟𝑠
时,𝑔𝑡𝑡
变号、𝑔𝑟𝑟
发散(在该坐标系下).这对应事件视界;更好的坐标(Eddington–Finkelstein、Kruskal)能消除坐标奇点.
7.2 静止观测者的引力时间膨胀
对固定 𝑟,𝜃,𝜙
的静止钟,𝑑𝑟 =𝑑𝜃 =𝑑𝜙 =0
,有
𝑑𝜏=√1−𝑟𝑠𝑟𝑑𝑡.
因此远处观察者看到近处的钟变慢.

7.3 能量与角动量守恒(测地线的第一积分)
Schwarzschild 度规对 𝑡
、𝜙
不显含,导致守恒量:
𝐸≡(1−𝑟𝑠𝑟)𝑐2𝑑𝑡𝑑𝜏,𝐿≡𝑟2𝑑𝜙𝑑𝜏.
在赤道面 𝜃 =𝜋/2
可将测地线化为「一维有效势」问题.
7.4 有效势(时类测地线)
在几何单位 𝐺 =𝑐 =1
下,𝑟𝑠 =2𝑀
,并可写成
(𝑑𝑟𝑑𝜏)2+𝑉eff(𝑟)=𝐸2,
其中(典型形式)
𝑉eff(𝑟)=(1−2𝑀𝑟)(1+𝐿2𝑟2).

7.5 光子球与 ISCO(竞赛很常考的两个半径)
在几何单位 𝐺 =𝑐 =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)
两静止观测者分别在 𝑟1
(发射)与 𝑟2
(接收,通常更远)处:
𝜈2𝜈1=√1−𝑟𝑠/𝑟11−𝑟𝑠/𝑟2
当 𝑟2 →∞
:
𝜈∞𝜈𝑟=√1−𝑟𝑠𝑟.
8.2 光线偏折(弱场近似)
对冲量参数(impact parameter)𝑏
,弱场一阶:
𝛼≈4𝐺𝑀𝑏𝑐2 

8.3 近日点进动(弱场、近开普勒轨道)
椭圆轨道半长轴 𝑎
、偏心率 𝑒
,每一周额外进动角
Δ𝜙≈6𝜋𝐺𝑀𝑎(1−𝑒2)𝑐2 
8.4 Shapiro 时间延迟(选学)
雷达信号掠过质量 𝑀
的天体会产生额外往返时间延迟.对「从 𝑟1
发出、掠过最近距离 𝑏
、到 𝑟2
接收」的几何,弱场近似常见写法为
Δ𝑡Shapiro≈2𝐺𝑀𝑐3ln(4𝑟1𝑟2𝑏2) 
(不同教材对几何定义略有差异,但共同点是:它随 𝐺𝑀/𝑐3
线性增长,并含对数项.)
8.5 宇宙学一瞥:FRW 度规与 Friedmann 方程(选读)
在「大尺度均匀各向同性」的宇宙假设下,度规可写为 FRW 形式:
𝑑𝑠2=−𝑐2𝑑𝑡2+𝑎(𝑡)2[𝑑𝑟21−𝑘𝑟2+𝑟2(𝑑𝜃2+sin2𝜃𝑑𝜙2)],![ds^2=-c^2dt^2+a(t)^2\left[\dfrac{dr^2}{1-kr^2}+r^2(d\theta^2+\sin^2\theta\,d\phi^2)\right],](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
其中 𝑎(𝑡)
为尺度因子,𝑘 =0, ±1
表示空间曲率符号.
代入 Einstein 方程并取物质为理想流体,可得 Friedmann 方程(最常用版本):
(˙𝑎𝑎)2=8𝜋𝐺3𝜌−𝑘𝑐2𝑎2+Λ𝑐23 
以及加速度方程
¨𝑎𝑎=−4𝜋𝐺3(𝜌+3𝑝𝑐2)+Λ𝑐23 
并由 ∇𝜇𝑇𝜇𝜈 =0
得到连续性方程
˙𝜌+3˙𝑎𝑎(𝜌+𝑝𝑐2)=0.
9. 例题(含全过程推导与计算)
例题 1:用弱场近似估算地球表面与高空的引力红移
已知地球半径 𝑅 ≈6.37 ×106 m
,质量 𝑀 ≈5.97 ×1024 kg
.在高度 ℎ ≪𝑅
处的钟与地面钟相比,频率相对变化 Δ𝜈/𝜈
约为多少?
解:
弱场中
Δ𝜈𝜈≈−ΔΦ𝑐2.
地球引力势(取无穷远为 0)为 Φ(𝑟) = −𝐺𝑀/𝑟
,因此
ΔΦ=Φ(𝑅+ℎ)−Φ(𝑅)=−𝐺𝑀𝑅+ℎ+𝐺𝑀𝑅.
当 ℎ ≪𝑅
,用展开
1𝑅+ℎ≈1𝑅(1−ℎ𝑅),
得
ΔΦ≈−𝐺𝑀𝑅(1−ℎ𝑅)+𝐺𝑀𝑅=𝐺𝑀𝑅2ℎ=𝑔ℎ.
所以
Δ𝜈𝜈≈−𝑔ℎ𝑐2 
若取 ℎ =1000 m
、𝑔 ≈9.8 m/s2
,则
∣Δ𝜈𝜈∣≈9.8×1000(3.0×108)2∼1.1×10−13.
** 结论:** 千米量级高度差带来 10−13
量级的频移(原子钟可测).
例题 2:Schwarzschild 时空中的静止钟差(精确公式)
在 Schwarzschild 外部时空中,两只静止钟分别位于 𝑟1
与 𝑟2
(𝑟2 >𝑟1 >𝑟𝑠
).若远处坐标时间过了 Δ𝑡
,两钟的固有时分别为多少?两者比值?
解:
静止钟有 𝑑𝑟 =𝑑𝜃 =𝑑𝜙 =0
,于是
𝑑𝜏=√1−𝑟𝑠𝑟𝑑𝑡.
积分得到
Δ𝜏(𝑟)=√1−𝑟𝑠𝑟Δ𝑡.
所以
Δ𝜏1Δ𝜏2=√1−𝑟𝑠/𝑟11−𝑟𝑠/𝑟2 
与引力红移公式一致(频率与时间尺度互为倒数).
例题 3:弱场光线偏折(给出标准推导骨架)
证明在弱场近似下,掠过质量 𝑀
的光线偏折角为 𝛼 ≈4𝐺𝑀𝑏𝑐2
.
解题思路骨架:
- 弱场近似度规(各向同性形式)可写为
𝑑𝑠2≈−(1+2Φ/𝑐2)𝑐2𝑑𝑡2+(1−2Φ/𝑐2)(𝑑𝑥2+𝑑𝑦2+𝑑𝑧2),Φ=−𝐺𝑀𝑟.
- 对光线 𝑑𝑠2 =0
,可得到等效「折射率」
𝑛(𝑟)≈1−2Φ𝑐2=1+2𝐺𝑀𝑟𝑐2.
- 用几何光学近似(费马原理)
𝛿∫𝑛𝑑𝑙=0
把光线当作在介质中传播的射线. 4. 在小偏折近似中,偏折角约为
𝛼≈∫+∞−∞𝜕𝜕𝑏(2𝐺𝑀𝑐2√𝑥2+𝑏2)𝑑𝑥=4𝐺𝑀𝑏𝑐2.
例题 4:近日点进动公式(从有效势到近似解)
给出椭圆轨道的相对论修正导致的每周进动
Δ𝜙≈6𝜋𝐺𝑀𝑎(1−𝑒2)𝑐2.
推导要点:
- 对 Schwarzschild 测地线,可得到轨道方程的近似形式(令 𝑢 =1/𝑟
)
𝑑2𝑢𝑑𝜙2+𝑢=𝐺𝑀𝐿2+3𝐺𝑀𝑐2𝑢2.
- 把 3𝐺𝑀𝑐2𝑢2
视为小扰动,设开普勒解
𝑢0=𝐺𝑀𝐿2(1+𝑒cos𝜙).
- 扰动会把角频率从 1 改成 1 −𝛿
,使得一周后相位多走
Δ𝜙≈2𝜋𝛿=6𝜋𝐺𝑀𝑎(1−𝑒2)𝑐2.
(竞赛层面通常要求记结论并能解释「为何是小扰动导致的频率偏移」.)
例题 5:求 Schwarzschild 的光子球半径与 ISCO 半径
在几何单位 𝐺 =𝑐 =1
下,证明 Schwarzschild 时空中光子的圆轨道半径为 𝑟 =3𝑀
,时类最内稳定圆轨道半径为 𝑟 =6𝑀
.
解(思路 + 关键步骤):
- 对赤道面运动,利用守恒量可把径向运动写成
(𝑑𝑟𝑑𝜆)2+𝑉eff(𝑟)=𝐸2.
- 圆轨道条件:𝑟 =常数
意味着
𝑑𝑟𝑑𝜆=0⇒𝑉eff(𝑟)=𝐸2,
且要「保持在圆轨道上」还需
𝑑𝑉eff𝑑𝑟=0.
- ** 光子(零质量)** 的有效势(形式上)为
𝑉(null)eff(𝑟)=(1−2𝑀𝑟)𝐿2𝑟2.
求导并令 𝑑𝑉/𝑑𝑟 =0
得
𝑟=3𝑀.
- ** 时类(有质量)** 圆轨道的稳定性由二阶导数决定:稳定要求 𝑑2𝑉eff𝑑𝑟2 >0
. 对 Schwarzschild,稳定圆轨道存在于 𝑟 >6𝑀
;临界点为
𝑟ISCO=6𝑀.
例题 6:径向自由落体到达视界需要多少固有时?坐标时间会怎样?
考虑从 𝑟 =𝑟0
(静止释放)径向下落到 𝑟 =𝑟𝑠
.定性说明:
- 自由落体者的固有时 𝜏
是否有限? - 远处静止观察者的坐标时间 𝑡
是否有限?
解(定性 + 关键结论):
- 对下落者,𝜏
是沿世界线积分得到的真实物理时间. - 在 Schwarzschild 坐标下,𝑡
在靠近 𝑟𝑠
时会出现对数发散(这是坐标效应). 结论:
固有时有限:自由落体者在有限 𝜏
内穿过视界.
- 坐标时间发散:远处观察者用 𝑡
描述会看到物体「永远接近但不穿过」视界.
(若需要具体积分表达式,可在后续补充「径向测地线第一积分」并做积分计算.)
例题 7:掠日光线偏折角的数量级估算
取 𝑏 ≈𝑅⊙ ≈7.0 ×108 m
,𝑀 ≈𝑀⊙ ≈2.0 ×1030 kg
.
解:
𝛼≈4𝐺𝑀𝑏𝑐2.
代入 𝐺 ≈6.67 ×10−11
、𝑐 ≈3.0 ×108
,得
𝛼∼4×6.67×10−11×2.0×10307.0×108×(3.0×108)2∼8.5×10−6rad.
换算为角秒(1 rad ≈206265 arcsec
):
𝛼∼1.75arcsec.
这就是经典的「掠日偏折约 1.75 角秒」.
例题 8:水星近日点进动的数量级估算
已知水星轨道参数(可查表):半长轴 𝑎 ≈5.79 ×1010 m
,偏心率 𝑒 ≈0.206
,太阳质量 𝑀 ≈𝑀⊙
.估算每周进动角.
解:
Δ𝜙≈6𝜋𝐺𝑀𝑎(1−𝑒2)𝑐2.
代入数值可得每一周(每绕太阳一圈)约 ∼5 ×10−7 rad
量级,对应约 0.1 arcsec
量级.累积到每世纪(约 415 圈)即得到著名的 ∼43 arcsec/century
.
10. 习题(含提示/答案要点)
10.1 基础题
- 证明:在局部惯性系中,某点可令 Γ𝜌 𝜇𝜈 =0
,但一般不能令 𝜕𝜎𝑔𝜇𝜈 =0
同时在邻域处处为 0. - 设弱场度规 𝑔00 = −(1 +2Φ/𝑐2)
,证明静止钟的频移满足 Δ𝜈/𝜈 ≈ −ΔΦ/𝑐2
. - 写出理想流体 𝑇𝜇𝜈
并解释各项物理意义.
10.2 计算题(更贴近竞赛/考题)
黑洞视界尺度估算:太阳质量 𝑀⊙ ≈2.0 ×1030 kg
,求 𝑟𝑠
的数量级.
答案要点:𝑟𝑠 =2𝐺𝑀/𝑐2 ≈3 km
.
引力时间膨胀:若某天体 𝑟𝑠 =10 km
,在 𝑟 =40 km
静止的钟相对于无穷远的钟慢多少?
提示:𝑑𝜏/𝑑𝑡 =√1−𝑟𝑠/𝑟
.
光偏折角数值:太阳半径 𝑅⊙ ≈7.0 ×108 m
,质量 𝑀⊙
,估算掠日光线偏折角.
提示:𝛼 ≈4𝐺𝑀/(𝑏𝑐2)
,取 𝑏 ≈𝑅⊙
.
光子球数值:对某黑洞 𝑟𝑠 =30 km
,求光子球半径 𝑟ph
与 ISCO 半径 𝑟ISCO
.
答案要点:𝑟ph =32𝑟𝑠
,𝑟ISCO =3𝑟𝑠
.
引力波应变量级:若某引力波通过地球时的应变 ℎ ∼10−21
,LIGO 级别的臂长 𝐿 ∼4 km
,估算长度变化 Δ𝐿
.
答案要点:Δ𝐿 ∼ℎ𝐿 ∼4 ×10−18 m
.
11. 用 Python 复现/绘制示意图(可选)
本文插图由 Python 生成并保存在 docs/modern/images/:
schwarzschild_time_dilation_factor.pngschwarzschild_effective_potential.pnglight_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(现代写法,推导清晰)
本页面最近更新:2026/7/3 22:14:58,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, 匿名同学
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用