第八学期 · 第1章 · 广义相对论 General Relativity

等效原理与时空几何 Equivalence Principle & Spacetime Geometry

🎯 学习目标 Learning Objectives

学完本章后，你将能够：

📖 核心概念 Core Concepts

1. 弱等效原理：自由下落普适性 The Weak Equivalence Principle (WEP)

弱等效原理 (weak equivalence principle, WEP) 又称 Universality of Free Fall (自由下落普适性)：在相同引力场中，任何检验粒子（无论成分、质量、内部结构）的加速度都相同。伽利略传说中比萨斜塔的实验、牛顿的单摆实验、Eötvös 1908 年的扭秤实验，直至 2017 年 MICROSCOPE 卫星，都证实了下式定义的 Eötvös 比值极小：

\[ \boxed{\;\eta\;\equiv\;2\,\frac{|a_1-a_2|}{|a_1+a_2|}\;\lesssim\;10^{-15}\quad\text{(MICROSCOPE 2017, Ti vs Pt)}\;} \]

意思是：钛与铂在地球引力场中下落的加速度之差小于百亿亿分之一。这强烈提示——引力质量 \(m_g\)（出现在 \(F_g=m_g g\) 中）与惯性质量 \(m_i\)（出现在 \(F=m_i a\) 中）严格相等：\(m_g=m_i\)。这一点不能由牛顿理论解释，它是广义相对论的实验基石。

等效原理的另一种表述（爱因斯坦称之为「我最快乐的想法」）：一个在均匀引力场中自由下落的密封舱里的观察者，感觉不到引力——舱内所有物体相对舱体都漂浮着，与外太空中无引力的惯性系局部无法区分 (locally indistinguishable)。

weak equivalence principle (弱等效原理), universality of free fall (自由下落普适性), Eötvös ratio \(\eta\) (Eötvös 比值), gravitational mass \(m_g\) vs inertial mass \(m_i\) (引力质量与惯性质量)

2. 爱因斯坦与强等效原理 The Einstein & Strong Equivalence Principles

WEP 只谈力学。爱因斯坦把它推广到全部物理（电磁、核力……），给出：

本章主要讨论 EEP。关键点：等效是局部的——只在一点（严格说是无穷小邻域）成立。延展区域内必出现潮汐效应（第 9 节），它是真正的引力几何信号，无法消除。

Einstein Equivalence Principle EEP (爱因斯坦等效原理), Strong Equivalence Principle SEP (强等效原理), local inertial frame LIF (局部惯性系), local Lorentz invariance LLI (局部洛伦兹不变性), local position invariance LPI (局部位置不变性)

3. 引力不是力，是时空几何 Gravity = Spacetime Geometry, Not a Force

等效原理的深刻推论：既然自由下落可以局部「消除」引力，引力就不能是牛顿意义上的「力」——你总不能通过换参考系就把电磁力或摩擦力变没。引力的本质是时空几何的弯曲：粒子在弯曲时空中走「直」的路径（测地线），看上去却被我们这些坐标观察者解读为「被引力拉弯」。

类比：地球表面是二维球面（弯曲）。两架沿「直线」飞行的飞机（在各自局部看是直的，即沿大圆 great circle 飞行），从赤道同一点出发、向北平行飞行，最终会在北极汇合——它们之间的距离随时间减小。从「平」的地图（墨卡托投影）上看，仿佛两机受了一个相互吸引的「力」；实际上根本没有力，它们各自沿球面几何上的「直线」走，是曲面本身的曲率让它们汇合。引力对物体轨迹的作用与此完全类似。

所以广义相对论的纲领：物质告诉时空如何弯曲；弯曲的时空告诉物质如何运动（Wheeler）。本章只讨论后者——给定时空几何（度规），粒子如何运动。下章再讨论前者（爱因斯坦方程）。

geometry not force (几何而非力), geodesic (测地线), great circle (大圆), embedding diagram (嵌入图), Wheeler's slogan (惠勒纲领)

4. 度规张量与原时 The Metric Tensor \(g_{\mu\nu}\) & Proper Time

时空中一切距离信息都编码在一个对称的 \(4\times4\) 张量——度规张量 \(g_{\mu\nu}=g_{\nu\mu}\)（metric tensor）里。它给出原时 (proper time, 检验粒子自身时钟所记录的时间)：

\[ \boxed{\;c^{2}\,\mathrm{d}\tau^{2}\;=\;g_{\mu\nu}\,\mathrm{d}x^{\mu}\,\mathrm{d}x^{\nu}\;=\;g_{\mu\nu}\,\dv{x^{\mu}}{\lambda}\dv{x^{\nu}}{\lambda}\,\mathrm{d}\lambda^{2}\;} \]

采用 mostly-minus 符号差 \((+,-,-,-)\)：类时世界线（粒子）\(\mathrm{d}\tau^{2}>0\)，类光（光子）\(\mathrm{d}\tau=0\)，类空 \(\mathrm{d}\tau^{2}<0\)。一条世界线 \(x^{\mu}(\lambda)\) 的累计原时是

\[ \tau_{AB}\;=\;\int_{A}^{B}\mathrm{d}\tau\;=\;\frac{1}{c}\int_{\lambda_A}^{\lambda_B}\!\sqrt{\,g_{\mu\nu}\,\dot x^{\mu}\dot x^{\nu}\,}\;\mathrm{d}\lambda. \]

度规就是广义相对论的「引力势」：\(g_{\mu\nu}\) 在每点的取值告诉你该处时空的几何（间隔的尺与钟）。自由粒子选择使原时取极值（通常极大值）的路径——这就是测地线（第 7 节）。

度规分量依赖坐标系：同一个时空，用不同坐标（球面、笛卡尔、各向同性、自由下落……）写出来的 \(g_{\mu\nu}\) 形式不同，但时空几何本身不变（坐标规范自由度）。

metric tensor \(g_{\mu\nu}\) (度规张量), proper time \(\tau\) (原时), line element \(\mathrm{d}s^{2}\) (线元), mostly-minus signature (mostly-minus 符号差), coordinate freedom (坐标自由度)

5. 闵可夫斯基度规（狭义相对论回顾）Review: Minkowski Metric (Special Relativity)

无引力（或局部惯性系下）时空是平直的，度规为闵可夫斯基度规 (Minkowski metric)：在笛卡尔坐标 \((ct,x,y,z)\) 下

\[ \boxed{\;\mathrm{d}s^{2}\;=\;c^{2}\mathrm{d}t^{2}-\mathrm{d}x^{2}-\mathrm{d}y^{2}-\mathrm{d}z^{2}\;=\;\eta_{\mu\nu}\,\mathrm{d}x^{\mu}\mathrm{d}x^{\nu},\qquad \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1)\;} \]

对应洛伦兹变换 (Lorentz transformation) 保持 \(\eta_{\mu\nu}\) 形式不变。沿速度 \(\vec{v}\) 运动的粒子，原时与坐标时关系为时间膨胀：

\[ \mathrm{d}\tau\;=\;\sqrt{1-\dfrac{v^{2}}{c^{2}}}\;\mathrm{d}t\;\equiv\;\frac{\mathrm{d}t}{\gamma},\qquad \gamma=\frac{1}{\sqrt{1-v^{2}/c^{2}}}. \]

闵可夫斯基时空中两点之间的测地线是直线 \(x^{\mu}(\tau)=u^{\mu}\tau+x_0^{\mu}\)（\(u^{\mu}\) 是四速，常矢量），相应于匀速直线运动——正是牛顿第一定律。所以狭义相对论是广义相对论在「无曲率」极限下的特例。

注意：用球坐标写闵可夫斯基度规（描述平直时空）时形式复杂：\(\mathrm{d}s^{2}=c^{2}\mathrm{d}t^{2}-\mathrm{d}r^{2}-r^{2}\mathrm{d}\theta^{2}-r^{2}\sin^{2}\!\theta\,\mathrm{d}\phi^{2}\)。度规分量随 \(r\) 变化，但时空仍平直（曲率张量为零）——这是「坐标非平凡但时空平直」的典型例子（详见第 10 节）。

Minkowski metric \(\eta_{\mu\nu}\) (闵可夫斯基度规), flat spacetime (平直时空), Lorentz transformation (洛伦兹变换), time dilation \(\gamma\) (时间膨胀), special relativity review (狭义相对论回顾)

6. 弯曲时空一例：施瓦氏度规 A First Curved Example: Schwarzschild Metric

球对称、无转动、无电荷的质量 \(M\)（如太阳、地球外）外部的时空由施瓦氏度规 (Schwarzschild metric, 1916) 描述：

\[ \boxed{\;\mathrm{d}s^{2}\;=\;\Bigl(1-\dfrac{r_s}{r}\Bigr)c^{2}\mathrm{d}t^{2}\;-\;\Bigl(1-\dfrac{r_s}{r}\Bigr)^{-1}\!\mathrm{d}r^{2}\;-\;r^{2}\mathrm{d}\Omega^{2},\qquad r_s=\dfrac{2GM}{c^{2}}\;} \]

其中 \(\mathrm{d}\Omega^{2}=\mathrm{d}\theta^{2}+\sin^{2}\!\theta\,\mathrm{d}\phi^{2}\)，而 \(r_s=2GM/c^{2}\) 是施瓦氏半径 (Schwarzschild radius)，对应黑洞的视界。对地球 \(r_s\approx 8.87\,\mathrm{mm}\)，对太阳 \(r_s\approx 2.95\,\mathrm{km}\)。\(\sqrt{1-r_s/r}\) 因子告诉我们：在 \(r\) 处的坐标钟 \(t\) 与原时 \(\tau\) 不同步——靠近 \(M\) 的钟走得慢（引力红移）。本章只把施瓦氏度规作为示例；下章会从爱因斯坦方程推导它。

Schwarzschild metric (施瓦氏度规), Schwarzschild radius \(r_s=2GM/c^{2}\) (施瓦氏半径), gravitational time dilation (引力时间膨胀), event horizon (事件视界)

7. 测地线：弯曲时空中的「直线」Geodesics: Straightest Paths in Curved Spacetime

自由粒子（不受非引力作用）的世界线使原时取极值 \(\delta\!\int\mathrm{d}\tau=0\)。变分给出测地线方程 (geodesic equation)：

\[ \boxed{\;\dv[2]{x^{\mu}}{\tau}+\Gamma^{\mu}_{\;\alpha\beta}\,\dv{x^{\alpha}}{\tau}\dv{x^{\beta}}{\tau}=0\;} \]

其中 克里斯托费尔符号 (Christoffel symbols / 连络系数) 完全由度规及其导数决定：

\[ \Gamma^{\mu}_{\;\alpha\beta}\;=\;\frac{1}{2}\,g^{\mu\nu}\Bigl(\partial_{\alpha}g_{\beta\nu}+\partial_{\beta}g_{\alpha\nu}-\partial_{\nu}g_{\alpha\beta}\Bigr). \]

测地线方程的几何意义：粒子在每个局部惯性系里都瞬时地匀速直线运动，但从一点到另一点把所有这些局部「直」段拼起来，在弯曲时空中就形成了曲线。这正是「最直的可能路径」(straightest possible path) 的含义——它局部看是直的，整体看却是弯的。地球绕太阳公转就是这样的测地线：地球并未被任何「力」拉住，它只是在弯曲时空中沿测地线自由下落，那个测地线恰好是闭合的椭圆。

geodesic equation (测地线方程), Christoffel symbols \(\Gamma^{\mu}_{\;\alpha\beta}\) (克里斯托费尔符号), connection coefficients (连络系数), straightest possible path (最直路径)

8. 测地线 = 牛顿第一定律的推广 Geodesic as Generalization of Newton's 1st Law

牛顿第一定律：「不受力的粒子做匀速直线运动」——平直时空中测地线方程的解。弯曲时空中，「不受力的粒子」沿测地线走，而测地线在曲率存在时表现为曲线。所以测地线方程就是把牛顿第一定律推广到弯曲时空。在弱场、低速极限下（见推导 1.B），测地线方程退化为牛顿第二定律：

\[ \dv[2]{\vec{x}}{t}\;=\;-\vec{\nabla}\Phi\quad\Longleftrightarrow\quad\vec{F}=m\vec{a}\quad(\Phi=\text{牛顿引力势}). \]

这表明：「引力 = 时空曲率」在弱场下重现牛顿引力，而在强场、高速下给出可观测的修正（水星近日点进动、光线偏折、引力波……，下两章详谈）。

Newton's first law (牛顿第一定律), Newtonian limit (牛顿极限), weak-field slow-motion (弱场低速)

9. 潮汐力：曲率的指纹 Tidal Forces — The Fingerprint of Curvature

等效原理告诉我们：在一点可以选自由下落系消除引力。但在延展区域（即使很小），不同位置的「下落方向」并不一致——这就是潮汐效应。考察一束最初静止、彼此相近的自由下落检验粒子。在平直时空中它们永远保持相对静止；在弯曲时空中它们的相对距离随时间变化——相对加速度正比于黎曼曲率张量 (Riemann curvature tensor) \(R^{\mu}_{\;\nu\alpha\beta}\)：

\[ \boxed{\;\Delta a^{i}\;=\;-c^{2}\,R^{i}_{\;0j0}\,\Delta x^{j}\quad(i,j=1,2,3)\;} \]

这是测地偏移方程 (geodesic deviation equation) 的非相对论极限。关键点：

tidal force (潮汐力), Riemann curvature tensor \(R^{\mu}_{\;\nu\alpha\beta}\) (黎曼曲率张量), geodesic deviation (测地偏移), intrinsic curvature (内禀曲率)

10. 坐标选择与规范自由度 Coordinate Choices & Gauge Freedom

广义相对论是广义协变 (general covariance) 的：物理结果与坐标选择无关，但具体计算中合理选坐标能极大简化问题。常见坐标选择：

重要：「度规分量是位置的函数」\(\neq\)「时空弯曲」。判断是否真有曲率要用黎曼张量（度规的二阶导数组合）。球坐标下的平直时空度规有非平凡 \(g_{\mu\nu}(r)\)，但 \(R^{\mu}_{\;\nu\alpha\beta}\equiv0\)。

general covariance (广义协变), Schwarzschild coordinates (施瓦氏坐标), isotropic coordinates (各向同性坐标), local inertial frame LIF (局部惯性系), coordinate singularity vs physical singularity (坐标奇点与物理奇点)

🔭 推导 Derivation

推导 1.A — 从原时极值导出测地线方程 Geodesic Equation from Extremal Proper Time

步骤 1（作用量）：粒子的世界线 \(x^{\mu}(\lambda)\) 的累计原时为

\[ \tau_{AB}=\frac{1}{c}\int_{\lambda_A}^{\lambda_B}\!L\,\mathrm{d}\lambda,\qquad L\equiv\sqrt{g_{\alpha\beta}\,\dot x^{\alpha}\dot x^{\beta}},\qquad \dot x^{\mu}\equiv\dv{x^{\mu}}{\lambda}. \]

步骤 2（欧拉–拉格朗日方程）：原时极值 \(\delta\!\int L\,\mathrm{d}\lambda=0\) 给出欧拉–拉格朗日方程

\[ \dv{}{\lambda}\Bigl(\frac{\partial L}{\partial\dot x^{\mu}}\Bigr)-\frac{\partial L}{\partial x^{\mu}}=0. \]

步骤 3（以原时本身为参数）：取 \(\lambda=\tau\)（对类时测地线总可以），则 \(g_{\alpha\beta}\dot x^{\alpha}\dot x^{\beta}=c^{2}\) 是常数（沿世界线）。计算偏导：

\[ \frac{\partial L}{\partial\dot x^{\mu}}=\frac{1}{2L}\cdot 2\,g_{\mu\beta}\dot x^{\beta}=\frac{g_{\mu\beta}\dot x^{\beta}}{L},\qquad \frac{\partial L}{\partial x^{\mu}}=\frac{1}{2L}\,(\partial_{\mu}g_{\alpha\beta})\,\dot x^{\alpha}\dot x^{\beta}. \]

步骤 4（代入化简）：代入 E–L 方程并两边乘 \(L\)：

\[ \dv{}{\tau}\bigl(g_{\mu\beta}\dot x^{\beta}\bigr)-\frac{1}{2}(\partial_{\mu}g_{\alpha\beta})\dot x^{\alpha}\dot x^{\beta}=0. \]

展开第一项：\(g_{\mu\beta}\ddot x^{\beta}+(\partial_{\gamma}g_{\mu\beta})\dot x^{\gamma}\dot x^{\beta}-\tfrac{1}{2}(\partial_{\mu}g_{\alpha\beta})\dot x^{\alpha}\dot x^{\beta}=0\)。把 \(\dot x^{\beta}\) 写成 \(\dv{x^{\beta}}{\tau}\)，并利用度规对称性合并 \(\dot x^{\alpha}\dot x^{\beta}\) 项（中括号内 \(\alpha,\beta\) 对称）：

\[ g_{\mu\beta}\,\ddot x^{\beta}+\frac{1}{2}\bigl(\partial_{\alpha}g_{\mu\beta}+\partial_{\beta}g_{\mu\alpha}-\partial_{\mu}g_{\alpha\beta}\bigr)\dot x^{\alpha}\dot x^{\beta}=0. \]

步骤 5（升指标得测地线方程）：两边乘逆度规 \(g^{\mu\nu}\) 求和，并利用 \(g^{\mu\nu}g_{\mu\beta}=\delta^{\nu}_{\;\beta}\)：

\[ \boxed{\;\ddot x^{\nu}+\underbrace{\frac{1}{2}g^{\nu\mu}\bigl(\partial_{\alpha}g_{\mu\beta}+\partial_{\beta}g_{\mu\alpha}-\partial_{\mu}g_{\alpha\beta}\bigr)}_{\displaystyle\Gamma^{\nu}_{\;\alpha\beta}}\dot x^{\alpha}\dot x^{\beta}=0\quad\Longleftrightarrow\quad \dv[2]{x^{\nu}}{\tau}+\Gamma^{\nu}_{\;\alpha\beta}\dv{x^{\alpha}}{\tau}\dv{x^{\beta}}{\tau}=0.\;} \]

这正是测地线方程，且同时给出克里斯托费尔符号的表达式。注意 \(\Gamma^{\nu}_{\;\alpha\beta}\) 是度规一阶导数的组合，而曲率（黎曼张量）含二阶导数——故克里斯托费尔符号本身可被局部坐标变换消除（局部惯性系），而曲率不能。

推导 1.B — 牛顿极限：测地线方程 → 牛顿第二定律 Newtonian Limit of the Geodesic Equation

步骤 1（弱场低速假设）：取静态弱场度规 \(g_{00}=-(1+2\Phi/c^{2})\)（符号差 \(-,+,+,+\)，本章正文用 \((+,-,-,-)\)，等价讨论；此处为与传统教材一致采用 mostly-plus）……

为避免符号混淆，采用本章统一的 mostly-minus 符号差 \((+,-,-,-)\)。设粒子低速 \(\dv{x^{i}}{\tau}\ll c\dv{t}{\tau}\)，所以测地线方程中 \(\alpha=\beta=0\) 项主导：

\[ \dv[2]{x^{\mu}}{\tau}\;\approx\;-\Gamma^{\mu}_{\;00}\Bigl(\dv{x^{0}}{\tau}\Bigr)^{2}=-\Gamma^{\mu}_{\;00}\,c^{2}\Bigl(\dv{t}{\tau}\Bigr)^{2}. \]

步骤 2（度规的弱场形式）：取 \(g_{00}=1+2\Phi/c^{2}\)（\(|\Phi|/c^{2}\ll1\)），空间分量近似 \(g_{ij}\approx-\delta_{ij}\)。则 \(\Gamma^{i}_{\;00}=\tfrac{1}{2}g^{ij}(-\partial_{j}g_{00})=-\tfrac{1}{2}\delta^{ij}\partial_{j}(1+2\Phi/c^{2})=-\partial_{i}\Phi/c^{2}\)（用到 \(g^{ij}\approx-\delta^{ij}\)，符号：\(-\delta^{ij}\cdot(-\partial_{j}(2\Phi/c^{2}))\cdot\tfrac{1}{2}=\partial_{i}\Phi/c^{2}\)……

步骤 3（小心符号）：仔细：\(g_{ij}=-\delta_{ij}\Rightarrow g^{ij}=-\delta^{ij}\)。代入公式：\(\Gamma^{i}_{\;00}=\tfrac{1}{2}g^{i\nu}(\partial_{0}g_{0\nu}+\partial_{0}g_{0\nu}-\partial_{\nu}g_{00})\)。静态 \(\partial_{0}=0\)，所以 \(\Gamma^{i}_{\;00}=-\tfrac{1}{2}g^{ij}\partial_{j}g_{00}=-\tfrac{1}{2}(-\delta^{ij})(2\partial_{j}\Phi/c^{2})=+\partial_{i}\Phi/c^{2}\)。

步骤 4（代入）：空间分量方程：

\[ \dv[2]{x^{i}}{\tau}\;\approx\;-\Gamma^{i}_{\;00}\,c^{2}\Bigl(\dv{t}{\tau}\Bigr)^{2}\;=\;-\partial_{i}\Phi\,\Bigl(\dv{t}{\tau}\Bigr)^{2}. \]

低速下 \(\dv{t}{\tau}\to 1\)，左边 \(\dv[2]{x^{i}}{\tau}\to\dv[2]{x^{i}}{t}\)。故

\[ \boxed{\;\dv[2]{\vec{x}}{t}\;=\;-\vec{\nabla}\Phi\quad\Longleftrightarrow\quad m\,\ddot{\vec{x}}=-m\,\vec{\nabla}\Phi=\vec{F}_{\text{Newton}}\;} \]

步骤 5（结论）：测地线方程在弱场低速极限下精确还原牛顿第二定律加牛顿引力。两个要点：(1) 引力势 \(\Phi\) 编码在度规的时间分量 \(g_{00}=1+2\Phi/c^{2}\) 中——「时间弯曲」主导了日常引力（苹果下落、行星公转）；(2) 空间曲率（\(g_{ij}\) 偏离 \(-\delta_{ij}\)）贡献小，但对光线偏折有可观测效应（下章）。

✏️ 典型例题 Worked Examples

例 1.1 Eötvös 比值：MICROSCOPE 卫星检验 The Eötvös Ratio (MICROSCOPE Test)

MICROSCOPE 卫星（2017 年发表）用一对 Ti/Pt 检验质量在地球轨道上做自由下落，测量它们的相对加速度。结果给 \(\eta_{\text{Ti,Pt}}=(1.0\pm 9.1_{\text{stat}}\pm 8.9_{\text{syst}})\times 10^{-15}\)。设在地球轨道高度 \(r=7.0\times 10^{6}\,\mathrm{m}\) 处，地球引力给两质量的共同加速度 \(a_{\text{Earth}}=GM_{\oplus}/r^{2}\)。(a) 求该处 \(a_{\text{Earth}}\)。(b) 把 \(\eta=10^{-15}\) 翻译成「Ti 与 Pt 加速度之差 \(|a_{\text{Ti}}-a_{\text{Pt}}|\)」(c) 这种差异比什么更小？给一个尺度感的对照。

(a) 地球引力加速度：\(GM_{\oplus}=3.986\times 10^{14}\,\mathrm{m^{3}/s^{2}}\)。

\[ a_{\text{Earth}}=\frac{GM_{\oplus}}{r^{2}}=\frac{3.986\times 10^{14}}{(7.0\times 10^{6})^{2}}=\frac{3.986\times 10^{14}}{4.90\times 10^{13}}=\mathbf{8.13\,\mathrm{m/s^{2}}}. \]

(b) 翻译为加速度差：由 \(\eta\) 定义（\(a_1\approx a_2\)，\(a_1+a_2\approx 2a\)）：

\[ |a_{\text{Ti}}-a_{\text{Pt}}|\;\approx\;\eta\cdot a_{\text{Earth}}=10^{-15}\times 8.13=\mathbf{8.1\times 10^{-15}\,\mathrm{m/s^{2}}}. \]

(c) 尺度感：这个加速度差相当于一颗原本静止的检验质量自由下落整整一年（\(T=3.16\times 10^{7}\,\mathrm{s}\)）后，Ti 比 Pt 多走的距离

\[ \Delta x=\tfrac{1}{2}\,\Delta a\,T^{2}=\tfrac{1}{2}(8.1\times 10^{-15})(3.16\times 10^{7})^{2}\,\mathrm{m}=\tfrac{1}{2}(8.1\times 10^{-15})(10^{15})\,\mathrm{m}\approx\mathbf{4\,\mathrm{m}}. \]

即在「最严格的检验」下，Ti 和 Pt 一年内下落距离之差还不到一层楼的高度——而它们各自（按牛顿匀加速估算）一年内下落约 \(\tfrac{1}{2}\cdot 8.13\cdot(3.16\times 10^{7})^{2}\approx 4.1\times 10^{15}\,\mathrm{m}\approx 0.43\,\mathrm{光年}\)。这个比值 \(4\,\mathrm{m}/(0.43\,\mathrm{光年})\sim 10^{-15}\) 正是 \(\eta\) 的大小，也是「引力质量 = 惯性质量」的精度——现代物理学中最精确的实验事实之一。

答案：(a) \(a_{\text{Earth}}=8.13\,\mathrm{m/s^{2}}\)；(b) \(|\Delta a|\approx 8\times 10^{-15}\,\mathrm{m/s^{2}}\)；(c) 一年下落差 \(\sim 4\,\mathrm{m}\)，比 4 光年的总下落距离小 \(10^{-16}\) 倍。

例 1.2 引力红移：Pound–Rebka 实验 Gravitational Redshift (Pound–Rebka)

1960 年 Pound 与 Rebka 在哈佛 Jefferson 实验室塔顶 (\(H=22.5\,\mathrm{m}\)) 用 ⁵⁷Fe 的 \(14.4137\,\mathrm{keV}\) Mössbauer \(\gamma\) 射线做实验：底部发射、顶部吸收（或反之），通过控制源的多普勒速度来补偿引力频移。\(g=9.81\,\mathrm{m/s^{2}}\)，\(c=2.998\times 10^{8}\,\mathrm{m/s}\)。(a) 用等效原理推导光子从底到顶的分数频移 \(\Delta\nu/\nu\)。(b) 计算数值。(c) 求补偿所需的多普勒速度 \(v\) 并与 Mössbauer 自然线宽对应的速度比较。

(a) 等效原理推导：把光子视为能量 \(E=h\nu\)、等效质量 \(m_{\text{eff}}=E/c^{2}=h\nu/c^{2}\) 的粒子。从底升到顶 \(H\) 需克服引力做功 \(m_{\text{eff}}\,gH\)，能量减少：

\[ \Delta E=-m_{\text{eff}}\,gH=-\frac{h\nu}{c^{2}}\,gH\quad\Rightarrow\quad\frac{\Delta\nu}{\nu}=\frac{\Delta E}{E}=-\frac{gH}{c^{2}}. \]

负号表示光子向上传播时频率降低（红移）。也可从施瓦氏度规 \(g_{00}=1-r_s/r\) 严格得到：\(\nu(r)\propto\sqrt{g_{00}}\)，弱场下 \(\Delta\nu/\nu=\Delta\Phi/c^{2}=-gH/c^{2}\)。两者一致——这就是局部位置不变性 (LPI) 的可观测检验。

(b) 数值：

\[ \frac{\Delta\nu}{\nu}=-\frac{gH}{c^{2}}=-\frac{(9.81)(22.5)}{(2.998\times 10^{8})^{2}}=-\frac{220.7}{8.988\times 10^{16}}=\mathbf{-2.46\times 10^{-15}}. \]

(c) 多普勒补偿速度：把放射源以小速度 \(v\) 朝吸收体方向移动，产生一阶多普勒频移 \(\Delta\nu/\nu=v/c\)。令 \(v/c=gH/c^{2}\)：

\[ v=\frac{gH}{c}=\frac{220.7}{2.998\times 10^{8}}=7.36\times 10^{-7}\,\mathrm{m/s}=\mathbf{0.74\,\mu\mathrm{m/s}}. \]

对照 Mössbauer 自然线宽：⁵⁷Fe 激发态寿命 \(\tau\approx 1.4\times 10^{-7}\,\mathrm{s}\)，自然线宽 \(\Gamma=\hbar/\tau\approx 4.7\times 10^{-9}\,\mathrm{eV}\)。对应速度宽度 \(\Gamma_v/c=\Gamma/E=4.7\times 10^{-9}/14413.7=3.3\times 10^{-13}\)，即 \(\Gamma_v\approx 98\,\mathrm{\mu m/s}\)。引力位移对应的 \(v=0.74\,\mathrm{\mu m/s}\) 仅是线宽的 \(\sim 1/130\)。Pound–Rebka 通过对完整共振曲线的扫描与积分，1960 年测得结果符合预期至 \((\pm 10)\%\)；1965 年 Pound–Snider 改进到 \((\pm 1)\%\)；现代 H/maser 实验已达 \(10^{-4}\) 量级。

答案：(a) \(\Delta\nu/\nu=-gH/c^{2}\)；(b) \(-2.46\times 10^{-15}\)；(c) 补偿速度 \(v\approx 0.74\,\mathrm{\mu m/s}\)（线宽的 \(\sim 1/130\)）。

例 1.3 用度规计算原时：双生子「旅行」Proper Time from the Metric (Twin Travel)

用闵可夫斯基度规 \(\mathrm{d}s^{2}=c^{2}\mathrm{d}t^{2}-\mathrm{d}x^{2}\) 直接计算两根世界线的原时。事件 A：\((t,x)=(0,0)\)；事件 B：\((t,x)=(10\,\mathrm{y},0)\)（坐标时 10 年）。观察者 Alice 静止在 \(x=0\)。观察者 Bob 以 \(v=+0.6c\) 走 5 年（坐标时）到 \((5\,\mathrm{y},3\,\mathrm{ly})\)（注意 \(0.6c\times 5\,\mathrm{y}=3\,\mathrm{ly}\)），瞬间反向，再以 \(v=-0.6c\) 走 5 年回到 B。求：(a) Alice 的原时 \(\tau_A\)。(b) Bob 的原时 \(\tau_B\)（沿每段折线用度规积分）。(c) 两人重逢时的年龄差。

(a) Alice：静止 → \(\mathrm{d}x=0\)，所以 \(c^{2}\mathrm{d}\tau^{2}=c^{2}\mathrm{d}t^{2}-0\)，\(\mathrm{d}\tau=\mathrm{d}t\)。

\[ \tau_A=\int_{0}^{10\,\mathrm{y}}\mathrm{d}t=\mathbf{10\,\mathrm{y}}. \]

(b) Bob 第一段（\(v=+0.6c\)）：\(\mathrm{d}x=0.6c\,\mathrm{d}t\)。

\[ c^{2}\mathrm{d}\tau^{2}=c^{2}\mathrm{d}t^{2}-(0.6c\,\mathrm{d}t)^{2}=c^{2}\mathrm{d}t^{2}(1-0.36)=c^{2}\mathrm{d}t^{2}\cdot 0.64. \]
\[ \mathrm{d}\tau=\mathrm{d}t\,\sqrt{0.64}=0.8\,\mathrm{d}t.\quad\Rightarrow\quad \Delta\tau_{B,\text{去}}=0.8\times 5\,\mathrm{y}=\mathbf{4\,\mathrm{y}}. \]

第二段（\(v=-0.6c\)）：\(\mathrm{d}x=-0.6c\,\mathrm{d}t\)，平方后相同：\(\Delta\tau_{B,\text{回}}=0.8\times 5=4\,\mathrm{y}\)。

\[ \tau_B=\Delta\tau_{B,\text{去}}+\Delta\tau_{B,\text{回}}=4\,\mathrm{y}+4\,\mathrm{y}=\mathbf{8\,\mathrm{y}}. \]

(c) 年龄差：两人从同一事件 A 出发、在事件 B 重逢。Alice 老 10 岁，Bob 老 8 岁。Bob 年轻了 2 岁。

关键观察：两根世界线连接同样的两个事件，但原时不同——原时是世界线的函数，不是端点的函数。Alice 的世界线（直线，匀速 = 0）是闵可夫斯基时空的测地线，对应极大原时。Bob 的折线（非测地线）原时较小。这与「欧氏几何里直线是最短」恰好相反——闵可夫斯基时空中直线（测地线）给出原时极大值（这是洛伦兹符号差的结果）。注意：Bob 的瞬间反向需要加速，所以两兄弟的运动学不对称——这就是「双生子佯谬」并不矛盾的原因。

答案：(a) \(\tau_A=10\,\mathrm{y}\)；(b) \(\tau_B=8\,\mathrm{y}\)；(c) Bob 比 Alice 年轻 2 岁。

例 1.4 球面上的测地线：北京—旧金山大圆距离 Geodesic on a Sphere (Great Circle)

把地球视为半径 \(R=6371\,\mathrm{km}\) 的理想球面（暂时忽略赤道隆起）。两城市：北京 \((\lambda_1=116.4^{\circ}\mathrm{E},\,\varphi_1=39.9^{\circ}\mathrm{N})\)，旧金山 \((\lambda_2=122.4^{\circ}\mathrm{W},\,\varphi_2=37.8^{\circ}\mathrm{N})\)。球面度规 \(\mathrm{d}s^{2}=R^{2}\mathrm{d}\varphi^{2}+R^{2}\cos^{2}\!\varphi\,\mathrm{d}\lambda^{2}\)。(a) 写出球面上的测地线方程并说明大圆 (great circle) 是其解。(b) 用球面余弦定理 \(\cos\Delta\sigma=\sin\varphi_1\sin\varphi_2+\cos\varphi_1\cos\varphi_2\cos(\Delta\lambda)\) 计算北京—旧金山大圆的角距 \(\Delta\sigma\) 与飞行距离 \(L=R\,\Delta\sigma\)。(c) 若错误地按「同纬度航线」（沿 \(\varphi=39^{\circ}\mathrm{N}\) 纬线）飞，需要多远的距离？比大圆远多少？

(a) 测地线方程 → 大圆：变分 \(\delta\!\int\mathrm{d}s=0\)（球面度规）给出测地线方程。几何上等价于：测地线是「过球心平面与球面的交线」，即大圆（半径为 \(R\) 的圆）。纬线（除赤道外）半径 \(R\cos\varphi小圆，不是测地线——沿纬线飞需要持续改变航向（被「力」拉离直线），所以不是最短路径。

(b) 大圆角距与距离：\(\Delta\lambda=116.4^{\circ}-(-122.4^{\circ})=238.8^{\circ}\)。

\[ \cos\Delta\sigma=\sin 39.9^{\circ}\sin 37.8^{\circ}+\cos 39.9^{\circ}\cos 37.8^{\circ}\cos 238.8^{\circ}. \]
\[ =(0.6414)(0.6129)+(0.7672)(0.7902)(-0.5210)=0.3932-0.3157=\mathbf{0.0775}. \]
\[ \Delta\sigma=\arccos(0.0775)=1.4932\,\mathrm{rad}=85.56^{\circ}. \]
\[ L=R\,\Delta\sigma=6371\times 1.4932=\mathbf{9512\,\mathrm{km}}. \]

(c) 沿纬度 39°N 的「纬线」距离：纬线半径 \(R_\varphi=R\cos 39^{\circ}=6371\times 0.7771=4950\,\mathrm{km}\)。但两城市不在同一纬度（39.9°N 与 37.8°N），先把旧金山「搬到」39°N 上做粗略估计：纬线上的经度差仍为 238.8°，但绕纬线时只能走 \(\le 360^{\circ}\) 的弧。最短纬线弧角 \(=360^{\circ}-238.8^{\circ}=121.2^{\circ}\)（往东飞经欧洲）。距离：

\[ L_{\text{纬线}}=R_\varphi\times(121.2^{\circ}\times\pi/180)=4950\times 2.115=\mathbf{10470\,\mathrm{km}}. \]

比大圆远 \(10470-9512=\mathbf{958\,\mathrm{km}}\)，即多飞约 10%。这就是洲际航班总沿大圆航线的原因——节省燃料与时间。注意大圆的最高纬度（顶点）约为 \(61^{\circ}\mathrm{N}\)（深入西伯利亚/北极圈附近），尽管两城市都在 \(38^{\circ}\mathrm{N}\) 附近——这印证「球面上两点最短路径向上凸向极区」。

答案：(a) 测地线 = 大圆 ✓；(b) 大圆距离 \(L\approx 9510\,\mathrm{km}\)；(c) 纬线粗估 \(10470\,\mathrm{km}\)，比大圆远约 \(960\,\mathrm{km}\)（\(\sim 10\%\)）。

例 1.5 自由下落实验室中的潮汐加速度 Tidal Acceleration in a Freely-Falling Lab

国际空间站 (ISS) 在地球表面上方 \(h=400\,\mathrm{km}\) 做近似自由下落（圆轨道），半径 \(r=R_{\oplus}+h=6371+400=6771\,\mathrm{km}\)。舱内一个 \(L=2.0\,\mathrm{m}\) 的实验装置，沿径向（指向地心方向）摆放。求：(a) 装置两端的潮汐加速度差 \(\Delta a_r\)（径向拉伸）。(b) 若装置改为沿横向（垂直于径向）摆放，潮汐加速度 \(\Delta a_\perp\)（横向压缩）。(c) 站在地球表面静止（不是自由下落）的相同 2 m 装置上的潮汐加速度（作对照）。已知 \(GM_{\oplus}=3.986\times 10^{14}\,\mathrm{m^{3}\!,s^{-2}}\)。

(a) 径向潮汐（ISS 高度）：由 \(a_r=-GM/r^{2}\)，对 \(r\) 求导：\(\dv{a_r}{r}=+2GM/r^{3}\)。两端距离差 \(\Delta r=L=2\,\mathrm{m}\)，相对加速度

\[ \Delta a_r=\dv{a_r}{r}\,\Delta r=\frac{2GM}{r^{3}}\,L=\frac{2(3.986\times 10^{14})(2.0)}{(6.771\times 10^{6})^{3}}. \]
\[ (6.771)^{3}=310.3\Rightarrow (6.771\times 10^{6})^{3}=3.103\times 10^{20}. \]
\[ \Delta a_r=\frac{1.594\times 10^{15}}{3.103\times 10^{20}}=\mathbf{5.14\times 10^{-6}\,\mathrm{m/s^{2}}}=5.1\,\mathrm{\mu m/s^{2}}. \]

这是拉伸（远离地心的那一端比靠近地心的那一端下落得慢，相对而言装置被「拉长」）。注意 ISS 整体加速度 \(a=GM/r^{2}=8.70\,\mathrm{m/s^{2}}\)，潮汐仅是其 \(6\times 10^{-7}\) 倍——所以 ISS 内的宇航员确实几乎「无重力」，但 2 m 装置两端有 \(\sim 5\,\mathrm{\mu m/s^{2}}\) 的相对拉伸，是真实存在且无法通过自由下落消除的曲率信号。

(b) 横向潮汐（ISS 高度）：横向（切向）方向的引力分量把粒子向中心轴拉回。简单推导（在过中心的平面内取两个横向邻居，离中心轴距离 \(b\)，到 \(M\) 距离 \(\sqrt{r^{2}+b^{2}}\)，方向偏角 \(\sim b/r\)）：

\[ \Delta a_\perp=-\frac{GM}{r^{3}}\,b\quad\Rightarrow\quad\text{对两端 }b=L=2\,\mathrm{m}\text{，}\;\;|\Delta a_\perp|=\frac{GM}{r^{3}}\,L=\frac{\Delta a_r}{2}=\mathbf{2.57\times 10^{-6}\,\mathrm{m/s^{2}}}. \]

横向是压缩（两端向中间轴靠拢），大小是径向拉伸的一半。三个主方向（一径向 + 两横向）合起来，潮汐张量的迹为零 \(\mathrm{tr}\,\mathcal{R}=0\)（这是真空爱因斯坦方程 \(R_{\mu\nu}=0\) 的体现，见下章）。

(c) 地球表面（\(r=R_{\oplus}=6.371\times 10^{6}\,\mathrm{m}\)）：注意——站在地面不是自由下落，地面支撑力把整个人托起。但潮汐张量只依赖 \(r\) 与 \(M\)，与参考系无关：

\[ \Delta a_r^{\text{地面}}=\frac{2GM}{R_{\oplus}^{3}}\,L=\frac{2(3.986\times 10^{14})(2.0)}{(6.371\times 10^{6})^{3}}=\frac{1.594\times 10^{15}}{2.585\times 10^{20}}=\mathbf{6.16\times 10^{-6}\,\mathrm{m/s^{2}}}. \]

仅比 ISS 处大 20%（因为 \(r\) 只差 6%）。\(\Delta a_\perp^{\text{地面}}=-3.08\times 10^{-6}\,\mathrm{m/s^{2}}\)。思考：2 m 高的人站在地面上，头比脚受到的地球引力略小（差 \(6\,\mathrm{\mu m/s^{2}}\)）——这就是为什么「直立时身高略被引力拉长」（实测效应约 \(\sim 10\) 微米量级，被血液循环等其他效应掩盖）。

答案：(a) ISS 处径向 \(\Delta a_r=+5.1\,\mathrm{\mu m/s^{2}}\)（拉伸）；(b) 横向 \(\Delta a_\perp=-2.6\,\mathrm{\mu m/s^{2}}\)（压缩，半大小）；(c) 地面 \(\Delta a_r=+6.2\,\mathrm{\mu m/s^{2}}\)。

⚠️ 常见误区 Common Pitfalls

🧮 练习题 Practice Problems

1. 太阳与地球的施瓦氏半径 Schwarzschild Radii of Sun & Earth
\(r_s=2GM/c^{2}\)。\(M_{\odot}=1.989\times 10^{30}\,\mathrm{kg}\)，\(M_{\oplus}=5.972\times 10^{24}\,\mathrm{kg}\)，\(G=6.674\times 10^{-11}\,\mathrm{N\,m^{2}\!,kg^{-2}}\)，\(c=2.998\times 10^{8}\,\mathrm{m/s}\)。
太阳：\(r_{s,\odot}=\dfrac{2(6.674\times 10^{-11})(1.989\times 10^{30})}{(2.998\times 10^{8})^{2}}=\dfrac{2.654\times 10^{20}}{8.988\times 10^{16}}=\mathbf{2.95\,\mathrm{km}}\)。
地球：\(r_{s,\oplus}=\dfrac{2(6.674\times 10^{-11})(5.972\times 10^{24})}{(2.998\times 10^{8})^{2}}=\dfrac{7.971\times 10^{14}}{8.988\times 10^{16}}=\mathbf{8.87\,\mathrm{mm}}\)。
对照：太阳半径 \(R_{\odot}=6.96\times 10^{5}\,\mathrm{km}\)，\(r_{s,\odot}/R_{\odot}\sim 4\times 10^{-6}\)（弱场）；中子星 \(r_s/R\sim 0.3\)（强场）；黑洞 \(r_s=R\)（视界即表面）。
思考：如果把太阳压缩到 2.95 km 半径以内（不改变其质量），它会变成黑洞——光也无法逃出。 答：\(r_{s,\odot}\approx 2.95\,\mathrm{km}\)；\(r_{s,\oplus}\approx 8.87\,\mathrm{mm}\)。
2. 中子星表面的引力红移 Redshift from a Neutron Star Surface
施瓦氏度规下静态观察者测得的光子频率 \(\nu(r)\propto\sqrt{g_{00}(r)}=\sqrt{1-r_s/r}\)。中子星 \(M=1.4\,M_{\odot}\)、半径 \(R=10\,\mathrm{km}\)。
\(r_s=2.95\,\mathrm{km}\times 1.4=4.13\,\mathrm{km}\)，\(r_s/R=0.413\)（强场）。
红移因子：\(\sqrt{1-r_s/R}=\sqrt{1-0.413}=\sqrt{0.587}=0.766\)。
从表面发出的光子到达无穷远观察者时频率变为 \(\nu_{\infty}=0.766\,\nu_{\text{em}}\)，即红移 \(z=1/0.766-1=0.305\)（波长增加 30%）。
能量减少：\(E_{\infty}/E_{\text{em}}=0.766\)，每光子损失约 23% 能量。
对照：太阳表面 \(r_s/R_{\odot}=4\times 10^{-6}\)，红移 \(z\approx 2\times 10^{-6}\)（极小，需要谱线 Doppler 补偿测）；白矮星 \(z\sim 10^{-4}\)；中子星 \(z\sim 0.1{-}0.3\)；黑洞视界 \(z\to\infty\)。
思考：当发射体位于视界 \(R=r_s\)，\(\sqrt{1-r_s/R}\to 0\)，频率红移到零——光「无法」从视界逃出。 答：红移因子 \(0.766\)，分数红移 \(z\approx 0.305\)。
3. 2D 球面上纬线不是测地线 A Parallel of Latitude is NOT a Geodesic
用球面度规 \(\mathrm{d}s^{2}=R^{2}\mathrm{d}\theta^{2}+R^{2}\sin^{2}\!\theta\,\mathrm{d}\phi^{2}\)。考察固定纬度 \(\theta=\theta_0\)、\(\dv{\theta}{s}=0\)、\(\dv{\phi}{s}=1/(R\sin\theta_0)\) 的曲线，代入测地线方程，验证除 \(\theta_0=\pi/2\)（赤道）外都不满足。
两条测地线方程（拉格朗日量 \(L=\tfrac{1}{2}(R^{2}\dot\theta^{2}+R^{2}\sin^{2}\!\theta\,\dot\phi^{2})\)，仿射参数为弧长 \(s\)）：
(i) \(\ddot\theta-\sin\theta\cos\theta\,\dot\phi^{2}=0\)；
(ii) \(\dv{}{s}(\sin^{2}\!\theta\,\dot\phi)=0\)。
沿纬线 \(\theta=\theta_0\)：\(\dot\theta=0,\;\ddot\theta=0\)，但 \(\dot\phi=\) 常数 \(\ne 0\)。
代入 (i)：\(0-\sin\theta_0\cos\theta_0\,\dot\phi^{2}=0\)。除 \(\sin\theta_0=0\)（极点）或 \(\cos\theta_0=0\)（\(\theta_0=\pi/2\)，赤道），方程不成立。
结论：球面上唯一作为测地线的「纬线」是赤道（大圆）。其他纬线（如北回归线、北京所在的 39°N）都不是测地线——沿纬线匀速航行的飞机必须不断向南（向赤道）调整航向，否则会被「甩出」纬线。这就是上一节例 1.4 中沿纬度航线不是最短路径的数学原因。
思考：这恰好类比弯曲时空中的「自由下落」——沿测地线不需要任何力，偏离测地线才需要力。 答：除赤道外所有纬线都不是测地线（方程 (i) 不满足）。
4. 月球对地球的潮汐力（拉伸海洋）Lunar Tides on Earth
月球质量 \(M_{\text{月}}=7.35\times 10^{22}\,\mathrm{kg}\)，到地球距离 \(d=3.84\times 10^{8}\,\mathrm{m}\)。计算月球引力在地球表面（\(R_{\oplus}=6.37\times 10^{6}\,\mathrm{m}\)）最近月点和与月心连线垂直方向上的潮汐加速度，并对比地球自身引力 \(g=9.81\,\mathrm{m/s^{2}}\)。
径向（朝月一侧，潮汐拉伸）：\(\Delta a_r=2GM_{\text{月}}R_{\oplus}/d^{3}\)。
\(d^{3}=(3.84\times 10^{8})^{3}=5.66\times 10^{25}\,\mathrm{m^{3}}\)。
\(\Delta a_r=\dfrac{2(6.674\times 10^{-11})(7.35\times 10^{22})(6.37\times 10^{6})}{5.66\times 10^{25}}=\dfrac{6.25\times 10^{19}}{5.66\times 10^{25}}=\mathbf{1.10\times 10^{-6}\,\mathrm{m/s^{2}}}\)。
横向（地球两端与月心连线垂直的方向，潮汐压缩）：\(\Delta a_\perp=-\tfrac{1}{2}\Delta a_r=-\mathbf{0.55\times 10^{-6}\,\mathrm{m/s^{2}}}\)。
对比 \(g\)：\(\Delta a_r/g\approx 1.1\times 10^{-7}\)，月球潮汐只有地球自引力的千万分之一。
但海水（自由流体）对这一微小差异极敏感：每秒钟被「拉长」\(\tfrac{1}{2}(1.1\times 10^{-6})(1)^{2}\approx 0.5\,\mathrm{\mu m}\)，几小时下来就是米级的海面隆起——这就是潮汐。一日两次涨潮（地球自转把每个海岸带入两个隆起带：一朝月、一背月）。
注：太阳对地球的潮汐 (\(\sim 0.5\times 10^{-6}\,\mathrm{m/s^{2}}\)）约为月球的一半（虽然太阳质量大，但距离远得多，潮汐 \(\propto M/d^{3}\)）。朔望大潮（日月叠加）和上下弦小潮（部分抵消）就是由此而来。 答：月球潮汐 \(\sim 1.1\,\mathrm{\mu m/s^{2}}\)（约 \(g\) 的千万分之一）。

🎛️ 互动仿真 Interactive Simulation

弯曲时空 — 橡胶膜嵌入图 Curved Spacetime: Rubber-Sheet Embedding. 一个可调质量的中心天体在二维空间切片上「压」出一个凹陷（红色网格），凹陷深度按 \(\Phi(r)=-GM/r\) 缩放可视化。一颗粒子（绿点）在该中心引力下运动，轨迹（紫色拖尾）即为弱场测地线（牛顿极限）。滑块：中心质量 \(M\)、初始切向速度 \(v_0\)、发射角 \(\alpha\)（0°=纯切向，90°=纯径向）、方位角与仰角（手动旋转）。观察要点：(1) 增大 \(M\)，凹陷变深、轨道弯曲加剧；(2) 当 \(v_0=\sqrt{GM/r_0}\)（圆周速度）时轨道闭合为圆；(3) \(v_0=\sqrt{2GM/r_0}\)（逃逸速度）时轨道开放为抛物线；(4) 能量 \(E=v^{2}/2-GM/r\) 与角动量 \(L=r\times v\) 实时显示，判定束缚/非束缚。注意：嵌入图只是直觉辅助，真正的运动通过度规与测地线方程计算；橡胶膜类比本身有循环论证（见误区 5）。

📝 本章小结 Summary & Key Formulas

\[ \boxed{\;\eta\;\equiv\;2\,\frac{|a_1-a_2|}{|a_1+a_2|}\;\lesssim\;10^{-15},\qquad m_g=m_i\;} \]
\[ \boxed{\;c^{2}\,\mathrm{d}\tau^{2}\;=\;g_{\mu\nu}\,\mathrm{d}x^{\mu}\,\mathrm{d}x^{\nu},\qquad \eta_{\mu\nu}=\mathrm{diag}(+1,-1,-1,-1)\;} \]
\[ \boxed{\;\mathrm{d}s^{2}=\Bigl(1-\dfrac{r_s}{r}\Bigr)c^{2}\mathrm{d}t^{2}-\Bigl(1-\dfrac{r_s}{r}\Bigr)^{-1}\!\mathrm{d}r^{2}-r^{2}\mathrm{d}\Omega^{2},\quad r_s=\dfrac{2GM}{c^{2}}\;} \]
\[ \boxed{\;\dv[2]{x^{\mu}}{\tau}+\Gamma^{\mu}_{\;\alpha\beta}\,\dv{x^{\alpha}}{\tau}\dv{x^{\beta}}{\tau}=0,\qquad \Gamma^{\mu}_{\;\alpha\beta}=\tfrac{1}{2}g^{\mu\nu}\bigl(\partial_{\alpha}g_{\beta\nu}+\partial_{\beta}g_{\alpha\nu}-\partial_{\nu}g_{\alpha\beta}\bigr)\;} \]
\[ \boxed{\;\dv[2]{\vec{x}}{t}=-\vec{\nabla}\Phi\;\;(\text{Newtonian limit, }g_{00}=1+2\Phi/c^{2})\;} \]
\[ \boxed{\;\Delta a^{i}=-c^{2}R^{i}_{\;0j0}\,\Delta x^{j};\qquad \Delta a_r=+\dfrac{2GM}{r^{3}}\,\Delta r,\quad \Delta a_\perp=-\dfrac{GM}{r^{3}}\,\Delta r\;} \]
\[ \boxed{\;\frac{\Delta\nu}{\nu}=-\frac{gH}{c^{2}}\quad(\text{引力红移, Pound–Rebka})\;} \]