← 返回 资料库 › 数学物理 › 宇宙学与天体物理学 数学物理
黑洞的几何 爱因斯坦场方程的史瓦西解与克尔解:视界、光子球与最内稳定圆轨道、光线偏折、奇点定理,以及引力波,并与事件视界望远镜和LIGO相联系。
直观 坠入时空中的一口井 把时空想象成不是一个空荡荡的舞台,而是一张会伸缩的橡皮膜:一个质量放在上面会把膜压出一个凹陷,而附近的一切——行星、光,甚至时间本身——都必须穿过那个凹陷运动。把足够多的质量塞进足够小的区域,凹陷就变成一口没有底的井:掉进去的光再也爬不出来。那个无法返回的区域就是黑洞 ,它的边界叫事件视界 。本页将精确算出这口井究竟有多深,什么东西会掉进去、什么东西只是被它弯曲绕行,以及我们如今如何能看到这口井投下的影子,并听到两口这样的井相撞时发出的涟漪。
弗拉姆抛物面 z ( r ) = 2 r s ( r − r s ) z(r) = 2\sqrt{r_s(r-r_s)} z ( r ) = 2 r s ( r − r s ) :将质量周围的空间几何嵌入到普通三维空间中的一种方式,展示了径向距离在 r = r s r=r_s r = r s 附近如何被拉伸。逐渐变陡的漏斗只是一种图示,并非真实的四维几何——但 r = r s r=r_s r = r s 处的喉部是真实的:那正是事件视界。 中学 这口井有多大?史瓦西半径 定义: 逃逸速度与史瓦西半径
质量为 M M M 、半径为 r r r 的球体的牛顿逃逸速度为 v esc = 2 G M / r v_{\text{esc}}=\sqrt{2GM/r} v esc = 2 GM / r :飞离到无穷远所需的最小速度。令 v esc = c v_{\text{esc}}=c v esc = c (光速)并解出 r r r ,得到一个特殊半径——史瓦西半径 r s = 2 G M c 2 r_s = \dfrac{2GM}{c^2} r s = c 2 2 GM ——在此半径处连光也无法逃脱。这个牛顿式的捷径给出了正确的数值,尽管光无法逃脱的真正原因(时空曲率,而不是光必须对抗的某种"引力")完全是相对论性的;下面给出严谨的推导。
r s = 2 G M c 2 r_s = \dfrac{2GM}{c^2} r s = c 2 2 GM 这里 G G G 是牛顿引力常数,M M M 是质量,c c c 是光速。对太阳而言,r s ≈ 2.95 km r_s\approx 2.95\text{ km} r s ≈ 2.95 km ——太阳需要被压缩到一个小镇大小才能变成黑洞;对地球而言,r s ≈ 8.9 mm r_s\approx 8.9\text{ mm} r s ≈ 8.9 mm ,比一颗葡萄还小。
大学 史瓦西度规 定义: 史瓦西度规
1916年,卡尔·史瓦西找到了爱因斯坦场方程唯一的球对称真空解。在坐标 ( t , r , θ , ϕ ) (t,r,\theta,\phi) ( t , r , θ , ϕ ) 中,时空间隔为 d s 2 = − ( 1 − r s r ) c 2 d t 2 + ( 1 − r s r ) − 1 d r 2 + r 2 d Ω 2 ds^2 = -\left(1-\dfrac{r_s}{r}\right)c^2\,dt^2 + \left(1-\dfrac{r_s}{r}\right)^{-1}dr^2 + r^2\,d\Omega^2 d s 2 = − ( 1 − r r s ) c 2 d t 2 + ( 1 − r r s ) − 1 d r 2 + r 2 d Ω 2 ,其中 d Ω 2 = d θ 2 + sin 2 θ d ϕ 2 d\Omega^2 = d\theta^2+\sin^2\theta\,d\phi^2 d Ω 2 = d θ 2 + sin 2 θ d ϕ 2 是单位球面上的度规,r s = 2 G M / c 2 r_s=2GM/c^2 r s = 2 GM / c 2 同上。它描述任意球对称、不自转质量——恒星、行星或黑洞——外部的时空,并在 r → ∞ r\to\infty r → ∞ 时退化为平直的闵可夫斯基时空。
d s 2 = − ( 1 − r s r ) c 2 d t 2 + ( 1 − r s r ) − 1 d r 2 + r 2 d Ω 2 ds^2 = -\left(1-\dfrac{r_s}{r}\right)c^2\,dt^2 + \left(1-\dfrac{r_s}{r}\right)^{-1}dr^2 + r^2\,d\Omega^2 d s 2 = − ( 1 − r r s ) c 2 d t 2 + ( 1 − r r s ) − 1 d r 2 + r 2 d Ω 2 有两个半径看起来是奇异的:g t t g_{tt} g tt 和 g r r g_{rr} g r r 都发散的 r = 0 r=0 r = 0 ,以及 g t t → 0 g_{tt}\to0 g tt → 0 且 g r r → ∞ g_{rr}\to\infty g r r → ∞ 的 r = r s r=r_s r = r s 。这其中只有一个是真正的物理奇点。克雷奇曼标量 K = R α β γ δ R α β γ δ = 48 G 2 M 2 c 4 r 6 K = R_{\alpha\beta\gamma\delta}R^{\alpha\beta\gamma\delta} = \dfrac{48\,G^2M^2}{c^4 r^6} K = R α β γ δ R α β γ δ = c 4 r 6 48 G 2 M 2 是由完整黎曼张量构造出的、与坐标无关的曲率度量;它在 r = 0 r=0 r = 0 处发散(真实、不可避免的无穷曲率),但在 r = r s r=r_s r = r s 处却完全有限。因此 r = r s r=r_s r = r s 处的发散只是坐标奇点 ——史瓦西坐标系的产物,而非时空本身的缺陷——正如地球仪经线在北极的"奇异性"只是经纬度坐标的产物,而非球面本身的缺陷。改用埃丁顿-芬克尔斯坦坐标 (将 t t t 替换为内向零坐标 v = t + r ∗ v=t+r_* v = t + r ∗ ,r ∗ = r + r s ln ∣ r / r s − 1 ∣ r_*=r+r_s\ln|r/r_s-1| r ∗ = r + r s ln ∣ r / r s − 1∣ )或完全扩展的克鲁斯卡尔-塞凯赖什坐标 ,可使所有度规分量在 r = r s r=r_s r = r s 处都保持正则,揭示出它是一个完全光滑的零曲面:事件视界。
史瓦西黑洞周围的关键半径 量 公式 以 r s r_s r s 为单位 事件视界 r s = 2 G M / c 2 r_s=2GM/c^2 r s = 2 GM / c 2 1 r s 1\,r_s 1 r s 光子球(不稳定圆形光轨道) r ph = 3 2 r s = 3 G M c 2 r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} r ph = 2 3 r s = c 2 3 GM 1.5 r s 1.5\,r_s 1.5 r s 最内稳定圆轨道(ISCO) r ISCO = 3 r s = 6 G M c 2 r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} r ISCO = 3 r s = c 2 6 GM 3 r s 3\,r_s 3 r s 影子半径(EHT成像所得) r shadow = 3 3 2 r s ≈ 2.6 r s r_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s r shadow = 2 3 3 r s ≈ 2.6 r s ≈ 2.6 r s \approx2.6\,r_s ≈ 2.6 r s
大学 黑洞周围的轨道:光子球与ISCO 无论是光线还是有质量粒子,在史瓦西几何中运动时都沿其测地线守恒两个量——一个类似能量的常数和一个类似角动量的常数——这与开普勒问题完全一样,但有效势中带有一个牛顿理论没有对应物的相对论性扭曲,正是它产生了两个具有里程碑意义的半径。
对于史瓦西黑洞,光子的圆轨道仅存在于 r ph = 3 2 r s = 3 G M c 2 r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} r ph = 2 3 r s = c 2 3 GM 处,且总是不稳定的;有质量粒子的圆轨道只有在 r ≥ r ISCO r \ge r_{\text{ISCO}} r ≥ r ISCO 时才稳定,在 r ISCO = 3 r s = 6 G M c 2 r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} r ISCO = 3 r s = c 2 6 GM 处边缘稳定,在 3 r s / 2 < r < 3 r s 3r_s/2<r<3r_s 3 r s /2 < r < 3 r s 区间内不稳定。
为什么成立? 把径向运动想象成一个球在一维势阱 V eff ( r ) V_{\text{eff}}(r) V eff ( r ) 中滚动:圆轨道位于势阱平坦之处(d V eff / d r = 0 dV_{\text{eff}}/dr=0 d V eff / d r = 0 ),且只有在势阱向上弯曲(d 2 V eff / d r 2 > 0 d^2V_{\text{eff}}/dr^2>0 d 2 V eff / d r 2 > 0 )而非向下弯曲之处才稳定。由于相对论性势中多出了一个牛顿引力所没有的 1 / r 3 1/r^3 1/ r 3 项,势阱在靠近黑洞处会出现一个极大值——在该半径以内根本不存在任何稳定圆轨道,无论粒子转动多快;物质只会直接坠落进去。
证明 第一步(守恒量与径向方程)。 沿限制在赤道面 θ = π / 2 \theta=\pi/2 θ = π /2 上的任意测地线(类时 κ = − c 2 \kappa=-c^2 κ = − c 2 ,或零测地线 κ = 0 \kappa=0 κ = 0 ),度规不依赖于 t t t 和 ϕ \phi ϕ 给出两个守恒量 E ~ = ( 1 − r s / r ) c 2 t ˙ \tilde E = (1-r_s/r)c^2\dot t E ~ = ( 1 − r s / r ) c 2 t ˙ 和 L ~ = r 2 ϕ ˙ \tilde L = r^2\dot\phi L ~ = r 2 ϕ ˙ (点 = d / d τ =d/d\tau = d / d τ )。将其代入归一化条件 g μ ν x ˙ μ x ˙ ν = κ g_{\mu\nu}\dot x^\mu\dot x^\nu=\kappa g μν x ˙ μ x ˙ ν = κ 并化简,对类时运动得到 ( d r d τ ) 2 = E ~ 2 − ( 1 − r s r ) ( c 2 + L ~ 2 r 2 ) ≡ E ~ 2 − V eff 2 ( r ) \left(\frac{dr}{d\tau}\right)^2 = \tilde E^2 - \left(1-\frac{r_s}{r}\right)\left(c^2+\frac{\tilde L^2}{r^2}\right) \equiv \tilde E^2 - V^2_{\text{eff}}(r) ( d τ d r ) 2 = E ~ 2 − ( 1 − r r s ) ( c 2 + r 2 L ~ 2 ) ≡ E ~ 2 − V eff 2 ( r )
第二步(圆轨道条件)。 展开得到 V eff 2 ( r ) = c 2 − c 2 r s r + L ~ 2 r 2 − L ~ 2 r s r 3 V^2_{\text{eff}}(r) = c^2 - \frac{c^2 r_s}{r} + \frac{\tilde L^2}{r^2} - \frac{\tilde L^2 r_s}{r^3} V eff 2 ( r ) = c 2 − r c 2 r s + r 2 L ~ 2 − r 3 L ~ 2 r s 圆轨道要求 r r r 恒定,即任意时刻都有 d r / d τ = 0 dr/d\tau=0 d r / d τ = 0 ,这要求 E ~ 2 = V eff 2 ( r ) \tilde E^2=V^2_{\text{eff}}(r) E ~ 2 = V eff 2 ( r ) 且(为使 r r r 不发生漂移)d V eff 2 / d r = 0 dV^2_{\text{eff}}/dr=0 d V eff 2 / d r = 0 :d d r V eff 2 = c 2 r s r 2 − 2 L ~ 2 r 3 + 3 L ~ 2 r s r 4 = 0 \frac{d}{dr}V^2_{\text{eff}} = \frac{c^2 r_s}{r^2} - \frac{2\tilde L^2}{r^3} + \frac{3\tilde L^2 r_s}{r^4} = 0 d r d V eff 2 = r 2 c 2 r s − r 3 2 L ~ 2 + r 4 3 L ~ 2 r s = 0
第三步(解出圆轨道的角动量)。 两边乘以 r 4 r^4 r 4 并解出 L ~ 2 \tilde L^2 L ~ 2 ,得到维持半径 r r r 处圆轨道所需的角动量:L ~ 2 ( r ) = c 2 r s r 2 2 r − 3 r s \tilde L^2(r) = \frac{c^2 r_s r^2}{2r-3r_s} L ~ 2 ( r ) = 2 r − 3 r s c 2 r s r 2 这已经包含了光子球的信息:当 r → ( 3 / 2 ) r s r\to (3/2)r_s r → ( 3/2 ) r s 从上方趋近时,分母 2 r − 3 r s → 0 + 2r-3r_s\to0^+ 2 r − 3 r s → 0 + ,L ~ 2 → + ∞ \tilde L^2\to+\infty L ~ 2 → + ∞ ——没有任何有限 角动量能在如此靠近处维持圆轨道,这正是下面达到的无质量(光子)极限。
第四步(边缘稳定性——ISCO)。 当 r r r 从无穷远减小时,L ~ 2 ( r ) \tilde L^2(r) L ~ 2 ( r ) 先减小,达到极小值,然后在 r = 3 r s / 2 r=3r_s/2 r = 3 r s /2 处发散;圆轨道稳定当且仅当进一步缩小 r r r 需要增大 L ~ 2 \tilde L^2 L ~ 2 (即稳定性由 d L ~ 2 / d r d\tilde L^2/dr d L ~ 2 / d r 的符号决定)。对第三步的结果求导,d L ~ 2 d r = 2 c 2 r s r ( r − 3 r s ) ( 2 r − 3 r s ) 2 \frac{d\tilde L^2}{dr} = \frac{2c^2 r_s\,r(r-3r_s)}{(2r-3r_s)^2} d r d L ~ 2 = ( 2 r − 3 r s ) 2 2 c 2 r s r ( r − 3 r s ) 它在 r = 0 r=0 r = 0 (排除)、r = 3 r s r=3r_s r = 3 r s 处为零,在 r > 3 r s / 2 r>3r_s/2 r > 3 r s /2 范围内再无其他零点。这个唯一的内部根就是边缘稳定半径:r ISCO = 3 r s = 6 G M c 2 r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} r ISCO = 3 r s = c 2 6 GM 当 r > 3 r s r>3r_s r > 3 r s 时轨道稳定(d L ~ 2 / d r > 0 d\tilde L^2/dr>0 d L ~ 2 / d r > 0 ,与开普勒问题一样,r r r 越大需要的角动量越大);当 3 r s / 2 < r < 3 r s 3r_s/2<r<3r_s 3 r s /2 < r < 3 r s 时轨道不稳定。
第五步(由零测地线得到光子球)。 对零测地线,令 κ = 0 \kappa=0 κ = 0 作同样代换,得到 ( d r d λ ) 2 = E 2 − L 2 r 2 ( 1 − r s r ) ≡ E 2 − V γ 2 ( r ) \left(\frac{dr}{d\lambda}\right)^2 = E^2 - \frac{L^2}{r^2}\left(1-\frac{r_s}{r}\right) \equiv E^2 - V^2_{\gamma}(r) ( d λ d r ) 2 = E 2 − r 2 L 2 ( 1 − r r s ) ≡ E 2 − V γ 2 ( r ) ,其中 V γ 2 ( r ) = L 2 ( 1 r 2 − r s r 3 ) V^2_\gamma(r) = L^2\left(\frac{1}{r^2}-\frac{r_s}{r^3}\right) V γ 2 ( r ) = L 2 ( r 2 1 − r 3 r s ) 。令括号内表达式的导数为零,d d r ( 1 r 2 − r s r 3 ) = − 2 r 3 + 3 r s r 4 = 0 ⟹ r = 3 r s 2 \frac{d}{dr}\left(\frac{1}{r^2}-\frac{r_s}{r^3}\right) = -\frac{2}{r^3}+\frac{3r_s}{r^4}=0 \;\Longrightarrow\; r=\frac{3r_s}{2} d r d ( r 2 1 − r 3 r s ) = − r 3 2 + r 4 3 r s = 0 ⟹ r = 2 3 r s 直接给出圆形光子轨道的唯一半径,r ph = 3 2 r s = 3 G M c 2 r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} r ph = 2 3 r s = c 2 3 GM ,与第三步得到的极限一致;它是不稳定的,因为括号 1 / r 2 − r s / r 3 1/r^2-r_s/r^3 1/ r 2 − r s / r 3 在该处取局部极大值 而非极小值——任何向内的扰动都会使光子螺旋坠入视界,任何向外的扰动都会使它逃逸到无穷远。
以碰撞参数 b ≫ r s b\gg r_s b ≫ r s (最近点远在视界之外)掠过质量 M M M 的光线,在最低阶上被偏折角度 Δ ϕ = 4 G M c 2 b \Delta\phi = \dfrac{4GM}{c^2 b} Δ ϕ = c 2 b 4 GM ——恰好是朴素粒子式计算所预言的牛顿值 2 G M / ( c 2 b ) 2GM/(c^2b) 2 GM / ( c 2 b ) 的两倍。
为什么成立? 牛顿引力若被当作作用在高速粒子上的力来处理,会预言一定的偏折,但它只考虑了时间 的曲率(质量附近时钟变慢)。广义相对论额外加入了来自空间 曲率的等量贡献(质量附近的尺子沿径向收缩),两种效应相加,使偏折角加倍——正是这个数值因子让1919年日食观测队得以区分爱因斯坦理论与牛顿理论。
证明 第一步(光的轨道方程)。 对限制在赤道面上的零测地线,用与光子球相同的守恒量 E = ( 1 − r s / r ) c 2 t ˙ E=(1-r_s/r)c^2\dot t E = ( 1 − r s / r ) c 2 t ˙ 和 L = r 2 ϕ ˙ L=r^2\dot\phi L = r 2 ϕ ˙ ,代入 g μ ν x ˙ μ x ˙ ν = 0 g_{\mu\nu}\dot x^\mu\dot x^\nu=0 g μν x ˙ μ x ˙ ν = 0 ,并将对 τ \tau τ 的导数转换为 u ≡ 1 / r u\equiv1/r u ≡ 1/ r 对 ϕ \phi ϕ 的导数(利用 r ˙ = − L d u / d ϕ \dot r = -L\,du/d\phi r ˙ = − L d u / d ϕ ),再求一次导后得到标准的光线弯曲轨道方程 d 2 u d ϕ 2 + u = 3 r s 2 u 2 , u ≡ 1 r \frac{d^2u}{d\phi^2}+u = \frac{3r_s}{2}u^2, \qquad u\equiv\frac1r d ϕ 2 d 2 u + u = 2 3 r s u 2 , u ≡ r 1
第二步(零阶:直线)。 去掉右边项(令 r s = 0 r_s=0 r s = 0 ,平直时空)得到 u 0 ′ ′ + u 0 = 0 u_0''+u_0=0 u 0 ′′ + u 0 = 0 ,其在 ϕ = π / 2 \phi=\pi/2 ϕ = π /2 处以最近距离 b b b 通过的解是直线 u 0 ( ϕ ) = sin ϕ b u_0(\phi) = \frac{\sin\phi}{b} u 0 ( ϕ ) = b sin ϕ ——这是极坐标形式下未偏折的路径。
第三步(一阶微扰)。 写出 u = u 0 + r s u 1 u=u_0+r_s u_1 u = u 0 + r s u 1 并代入第一步的方程,只保留关于 r s r_s r s 的线性项;u 0 ′ ′ + u 0 = 0 u_0''+u_0=0 u 0 ′′ + u 0 = 0 部分抵消,剩下 u 1 u_1 u 1 满足的受迫线性振子方程:u 1 ′ ′ + u 1 = 3 2 b 2 sin 2 ϕ u_1''+u_1 = \frac{3}{2b^2}\sin^2\phi u 1 ′′ + u 1 = 2 b 2 3 sin 2 ϕ 一个与最近点两侧对称弯曲相匹配的特解是 u 1 ( ϕ ) = 1 2 b 2 ( 1 + cos 2 ϕ ) u_1(\phi) = \frac{1}{2b^2}\left(1+\cos^2\phi\right) u 1 ( ϕ ) = 2 b 2 1 ( 1 + cos 2 ϕ ) ,因此完整的一阶轨迹为 u ( ϕ ) = sin ϕ b + r s 2 b 2 ( 1 + cos 2 ϕ ) u(\phi) = \frac{\sin\phi}{b} + \frac{r_s}{2b^2}\left(1+\cos^2\phi\right) u ( ϕ ) = b sin ϕ + 2 b 2 r s ( 1 + cos 2 ϕ )
第四步(提取总偏折角)。 平直空间中 u → 0 u\to0 u → 0 恰好发生在 ϕ = 0 \phi=0 ϕ = 0 和 ϕ = π \phi=\pi ϕ = π 处;有了 r s r_s r s 修正后,u → 0 u\to0 u → 0 发生在 ϕ = − δ 1 \phi=-\delta_1 ϕ = − δ 1 和 ϕ = π + δ 2 \phi=\pi+\delta_2 ϕ = π + δ 2 处,δ 1 , δ 2 \delta_1,\delta_2 δ 1 , δ 2 为小量。将 u ( − δ 1 ) = 0 u(-\delta_1)=0 u ( − δ 1 ) = 0 对 δ 1 \delta_1 δ 1 和 r s r_s r s 展开到一阶,得到 − δ 1 / b + r s / b 2 = 0 -\delta_1/b+r_s/b^2=0 − δ 1 / b + r s / b 2 = 0 ,故 δ 1 = r s / b \delta_1=r_s/b δ 1 = r s / b ;对另一端做相同计算得 δ 2 = r s / b \delta_2=r_s/b δ 2 = r s / b 。出射渐近线偏离与入射渐近线严格反向的总角度为 Δ ϕ = δ 1 + δ 2 = 2 r s b = 4 G M c 2 b \Delta\phi = \delta_1+\delta_2 = \frac{2r_s}{b} = \frac{4GM}{c^2 b} Δ ϕ = δ 1 + δ 2 = b 2 r s = c 2 b 4 GM ,代回 r s = 2 G M / c 2 r_s=2GM/c^2 r s = 2 GM / c 2 即得所要证明的公式。
进阶 奇点不可避:彭罗斯定理 史瓦西黑洞内部的 r = 0 r=0 r = 0 奇点,可能看起来只是求解爱因斯坦方程时所假设的严格球对称性带来的特殊现象——或许对称性较低、更贴近现实的坍缩能够避开它。1965年,罗杰·彭罗斯证明了相反的结论:一旦形成囚陷面 (一个闭合曲面,其上向外和向内两束光线都 在会聚,而不像远离任何黑洞处那样只有向内的一束会聚),那么对于任何满足合理能量条件的物质,无论对称性如何,奇点都不可避免。这是第一个现代意义上的奇点定理,开启了如今主导数学相对论的整体、几何方法。
若时空包含一个非紧的柯西曲面,满足零能量条件,并包含一个闭合囚陷面,则时空在未来方向上是零测地不完备的:至少存在一条光线无法向未来无限延伸——这就是彭罗斯意义下的奇点。
为什么成立? 在广义相对论中,引力总是使光线会聚而非发散(这正是能量条件所刻画的内容);因此在某个闭合曲面上已经两侧都在会聚的光线束会持续会聚,而几何上会聚的光线必然相交——但在一个性质良好、整体双曲的时空内部,两条光线相交(共轭点)若不为该光线的存在制造一个边界,就无法自洽。该定理把这种会聚的直觉转化为一个严格的不等式(雷查德胡里方程)和一个整体拓扑障碍。
证明 第一步(用符号表述囚陷面)。 设 S S S 是一个闭合的类空二维曲面,θ \theta θ 表示从 S S S 正交出发的零测地线束的膨胀率 ——即光束横截面积增长的相对速率。远离任何质量处,外向光束膨胀(θ out > 0 \theta_{\text{out}}>0 θ out > 0 ),内向光束收缩(θ in < 0 \theta_{\text{in}}<0 θ in < 0 )。若在整个 S S S 上反而有 θ out < 0 , θ in < 0 \theta_{\text{out}}<0, \qquad \theta_{\text{in}}<0 θ out < 0 , θ in < 0 ,则称 S S S 是囚陷的 :即使是向外的光也在被拖向内部。
第二步(雷查德胡里方程)。 沿切向量为 k a k^a k a 的零测地同余族,对膨胀率的定义 θ = ∇ a k a \theta=\nabla_a k^a θ = ∇ a k a 求导,并利用黎曼张量的定义交换导数次序,得到精确的运动学恒等式 d θ d λ = − θ 2 2 − σ a b σ a b − R a b k a k b \frac{d\theta}{d\lambda} = -\frac{\theta^2}{2} - \sigma_{ab}\sigma^{ab} - R_{ab}k^ak^b d λ d θ = − 2 θ 2 − σ ab σ ab − R ab k a k b ,其中 σ a b \sigma_{ab} σ ab 是该同余族的剪切张量。这个方程纯粹是几何性的——在假设任何物理之前,它由曲率的定义得出,说明任何一族光线必然如何弯曲。
第三步(能量条件迫使进一步会聚)。 物理上合理的物质满足零能量条件 R a b k a k b ≥ 0 R_{ab}k^ak^b \ge 0 R ab k a k b ≥ 0 ,且剪切项 σ a b σ a b ≥ 0 \sigma_{ab}\sigma^{ab}\ge0 σ ab σ ab ≥ 0 恒成立(它是平方和);从第二步恒等式右边去掉这两个非负项,得到不等式 d θ d λ ≤ − θ 2 2 \frac{d\theta}{d\lambda} \le -\frac{\theta^2}{2} d λ d θ ≤ − 2 θ 2
第四步(里卡蒂不等式迫使在有限时间内发散)。 在囚陷面上,λ = 0 \lambda=0 λ = 0 时 θ 0 < 0 \theta_0<0 θ 0 < 0 。只要 θ \theta θ 保持为负,将不等式除以 θ 2 > 0 \theta^2>0 θ 2 > 0 并整理,得到 d d λ ( 1 θ ) ≥ 1 2 ⟹ 1 θ ( λ ) ≥ 1 θ 0 + λ 2 \frac{d}{d\lambda}\left(\frac{1}{\theta}\right) \ge \frac12 \;\Longrightarrow\; \frac{1}{\theta(\lambda)} \ge \frac{1}{\theta_0}+\frac{\lambda}{2} d λ d ( θ 1 ) ≥ 2 1 ⟹ θ ( λ ) 1 ≥ θ 0 1 + 2 λ 由于右边随 λ \lambda λ 增大而无界增长,而 1 / θ 1/\theta 1/ θ (为负)必须满足此下界,1 / θ 1/\theta 1/ θ 被迫趋向 0 − 0^- 0 − ,即在某个有限仿射参数 λ ∗ ≤ 2 / ∣ θ 0 ∣ \lambda_*\le 2/|\theta_0| λ ∗ ≤ 2/∣ θ 0 ∣ 处 θ → − ∞ \theta\to-\infty θ → − ∞ :与 S S S 正交的零测地同余族必须在有限仿射时间内出现共轭点(散焦面,相邻光线在此相交)。
第五步(从局部会聚到整体不完备)。 因果结构理论的一个标准结果(此处引用但不重新推导)表明,整体双曲时空中的零测地线在经过共轭点之后不能保持因果边界性(不能停留在某集合未来的边界上)。结合非紧的柯西曲面,这就排除了测地线回环或以通常方式终止的可能,从而迫使至少一条与 S S S 正交的零测地线未来不完备 :它无法延伸到任意大的仿射参数。根据定义,这种不完备性正是彭罗斯所说的奇点——该定理表明,只要囚陷面形成,无需任何对称性假设,奇点就必然出现。
大学 实际应用与典型例题 黑洞几何学在2010年代不再纯属理论。事件视界望远镜 (EHT)将一个行星尺度的射电天线阵列作为单一干涉仪工作,拍摄黑洞投射在周围发光气体上的阴影 ——一个接近光子球相关半径 r shadow ≈ 2.6 r s r_{\text{shadow}}\approx2.6\,r_s r shadow ≈ 2.6 r s 的暗盘,而非视界本身。而LIGO–Virgo–KAGRA 探测到两个黑洞螺旋靠近并合并时发出的啁啾引力波,其质量信息编码在波的频率扫描和极其微小的应变振幅 h ∼ 10 − 21 h\sim10^{-21} h ∼ 1 0 − 21 之中。
例题: 计算人马座A*的阴影大小
人马座A*是我们星系中心的超大质量黑洞,质量 M ≈ 4.0 × 10 6 M ⊙ M\approx4.0\times10^6\,M_\odot M ≈ 4.0 × 1 0 6 M ⊙ ,距离 D ≈ 8 kpc D\approx8\text{ kpc} D ≈ 8 kpc 。请估算其阴影的角直径,并与事件视界望远镜2022年测得的约 51.8 μ as 51.8\,\mu\text{as} 51.8 μ as 相比较。
解答 第一步(史瓦西半径)。 对一个太阳质量,r s ( M ⊙ ) = 2 G M ⊙ / c 2 ≈ 2.95 km r_s(M_\odot)=2GM_\odot/c^2\approx2.95\text{ km} r s ( M ⊙ ) = 2 G M ⊙ / c 2 ≈ 2.95 km 。按质量线性缩放,r s ≈ 4.0 × 10 6 × 2.95 km ≈ 1.18 × 10 7 km = 1.18 × 10 10 m r_s\approx4.0\times10^6\times2.95\text{ km}\approx1.18\times10^7\text{ km}=1.18\times10^{10}\text{ m} r s ≈ 4.0 × 1 0 6 × 2.95 km ≈ 1.18 × 1 0 7 km = 1.18 × 1 0 10 m 。
第二步(阴影直径)。 阴影半径为 r shadow = 3 3 2 r s ≈ 2.6 r s r_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s r shadow = 2 3 3 r s ≈ 2.6 r s ,因此直径为 2 r shadow = 2 × 2.6 r s ≈ 5.2 r s ≈ 6.14 × 10 10 m 2r_{\text{shadow}}=2\times2.6\,r_s\approx5.2\,r_s\approx6.14\times10^{10}\text{ m} 2 r shadow = 2 × 2.6 r s ≈ 5.2 r s ≈ 6.14 × 1 0 10 m 。
第三步(换算为角度)。 小角近似下,角大小 = = = 物理大小 / / / 距离。取 D = 8 kpc = 8 × 3.086 × 10 19 m = 2.47 × 10 20 m D=8\text{ kpc}=8\times3.086\times10^{19}\text{ m}=2.47\times10^{20}\text{ m} D = 8 kpc = 8 × 3.086 × 1 0 19 m = 2.47 × 1 0 20 m ,角度为 6.14 × 10 10 / 2.47 × 10 20 ≈ 2.49 × 10 − 10 rad 6.14\times10^{10}/2.47\times10^{20}\approx2.49\times10^{-10}\text{ rad} 6.14 × 1 0 10 /2.47 × 1 0 20 ≈ 2.49 × 1 0 − 10 rad 。
第四步(把弧度换算为微角秒)。 1弧度等于 2.063 × 10 11 μ as 2.063\times10^{11}\,\mu\text{as} 2.063 × 1 0 11 μ as ,因此阴影直径约为 2.49 × 10 − 10 × 2.063 × 10 11 μ as ≈ 51.3 μ as 2.49\times10^{-10}\times2.063\times10^{11}\,\mu\text{as}\approx51.3\,\mu\text{as} 2.49 × 1 0 − 10 × 2.063 × 1 0 11 μ as ≈ 51.3 μ as ——与EHT测得的 51.8 μ as 51.8\,\mu\text{as} 51.8 μ as 相差在 1 % 1\% 1% 以内,有力地证实了人马座A*的行为与史瓦西几何(更准确地说是自旋接近零的克尔几何)的预言完全一致。
例题: GW150914的啁啾质量与应变数量级
LIGO首次探测到的GW150914,来自两个质量约为 m 1 ≈ 36 M ⊙ m_1\approx36\,M_\odot m 1 ≈ 36 M ⊙ 和 m 2 ≈ 29 M ⊙ m_2\approx29\,M_\odot m 2 ≈ 29 M ⊙ 的黑洞在距离 D ≈ 410 Mpc D\approx410\text{ Mpc} D ≈ 410 Mpc 处并合。计算控制波频率扫描的啁啾质量 M \mathcal M M ,并估算地球处应变 h h h 的数量级。
解答 第一步(啁啾质量)。 M = ( m 1 m 2 ) 3 / 5 ( m 1 + m 2 ) 1 / 5 \mathcal M = \dfrac{(m_1 m_2)^{3/5}}{(m_1+m_2)^{1/5}} M = ( m 1 + m 2 ) 1/5 ( m 1 m 2 ) 3/5 由 m 1 m 2 ≈ 1044 M ⊙ 2 m_1m_2\approx1044\,M_\odot^2 m 1 m 2 ≈ 1044 M ⊙ 2 及 m 1 + m 2 ≈ 65 M ⊙ m_1+m_2\approx65\,M_\odot m 1 + m 2 ≈ 65 M ⊙ :( m 1 m 2 ) 3 / 5 ≈ 1044 0.6 M ⊙ 1.2 ≈ 100 M ⊙ 1.2 (m_1m_2)^{3/5}\approx1044^{0.6}\,M_\odot^{1.2}\approx100\,M_\odot^{1.2} ( m 1 m 2 ) 3/5 ≈ 104 4 0.6 M ⊙ 1.2 ≈ 100 M ⊙ 1.2 ,( m 1 + m 2 ) 1 / 5 ≈ 65 0.2 ≈ 2.3 (m_1+m_2)^{1/5}\approx65^{0.2}\approx2.3 ( m 1 + m 2 ) 1/5 ≈ 6 5 0.2 ≈ 2.3 ,得 M ≈ 100 / 2.3 M ⊙ ≈ 28 M ⊙ \mathcal M\approx100/2.3\,M_\odot\approx28\,M_\odot M ≈ 100/2.3 M ⊙ ≈ 28 M ⊙ ——非常接近LIGO报告的约 28 M ⊙ 28\,M_\odot 28 M ⊙ 。
第二步(源的长度尺度)。 总质量对应的史瓦西半径为 r s , tot = 2 G ( m 1 + m 2 ) / c 2 ≈ 65 × 2.95 km ≈ 1.9 × 10 5 m r_{s,\text{tot}}=2G(m_1+m_2)/c^2\approx65\times2.95\text{ km}\approx1.9\times10^5\text{ m} r s , tot = 2 G ( m 1 + m 2 ) / c 2 ≈ 65 × 2.95 km ≈ 1.9 × 1 0 5 m :这是并合临近时系统的尺寸。
第三步(数量级应变估计)。 临近并合时轨道速度接近光速 c c c 的相当一部分;取 v / c ∼ 0.5 v/c\sim0.5 v / c ∼ 0.5 ,用标度关系 h ∼ r s , tot D ( v c ) 2 h \sim \dfrac{r_{s,\text{tot}}}{D}\left(\dfrac{v}{c}\right)^2 h ∼ D r s , tot ( c v ) 2 ,并取 D = 410 Mpc ≈ 1.26 × 10 25 m D=410\text{ Mpc}\approx1.26\times10^{25}\text{ m} D = 410 Mpc ≈ 1.26 × 1 0 25 m :h ∼ ( 1.9 × 10 5 / 1.26 × 10 25 ) × 0.25 ≈ 3.8 × 10 − 21 h\sim(1.9\times10^5/1.26\times10^{25})\times0.25\approx3.8\times10^{-21} h ∼ ( 1.9 × 1 0 5 /1.26 × 1 0 25 ) × 0.25 ≈ 3.8 × 1 0 − 21 。
第四步(比较)。 这一数量级估计与LIGO实际测得的应变 h ∼ 10 − 21 h\sim10^{-21} h ∼ 1 0 − 21 相差不过几倍——这个微小的数字并不意味着源处效应微弱(临近并合时那里的时空弯曲程度几乎达到了物理所允许的极限),而只是引力辐射在数百兆秒差距的距离上按 1 / D 1/D 1/ D 被稀释的结果。
常见错误. 事件视界并不是 一个物理表面,r = r s r=r_s r = r s 也不是 奇点。对一个自由下落穿过 r = r s r=r_s r = r s 的观测者而言,局部上不会发生任何特殊现象——没有墙壁,没有力的骤增,没有辐射爆发(这正是“坐标奇点”的含义)。真正会把一切碾碎的奇点位于 r = 0 r=0 r = 0 ,一旦穿过视界,它完全位于观测者的未来 ,而不在视界本身。把这两者混为一谈——把视界当作奇点本身,或当作物质可以“撞上”的实体表面——是关于黑洞最常见的误解。 进阶 旋转黑洞:克尔几何与能层 真实的天体物理黑洞诞生于自转恒星或并合过程,带有角动量 J J J ,可用自旋长度 a = J / ( M c ) a=J/(Mc) a = J / ( M c ) 参数化。罗伊·克尔于1963年找到了唯一的稳态轴对称真空解。自转把单一的史瓦西曲面 r = r s r=r_s r = r s 分裂成两个不同的曲面:位于 r + = G M / c 2 + ( G M / c 2 ) 2 − a 2 r_+ = GM/c^2 + \sqrt{(GM/c^2)^2-a^2} r + = GM / c 2 + ( GM / c 2 ) 2 − a 2 的外事件视界 ,以及其外侧更大、呈扁球状的静止极限面 r ergo ( θ ) = G M c 2 + ( G M c 2 ) 2 − a 2 cos 2 θ r_{\text{ergo}}(\theta) = \dfrac{GM}{c^2} + \sqrt{\left(\dfrac{GM}{c^2}\right)^2 - a^2\cos^2\theta} r ergo ( θ ) = c 2 GM + ( c 2 GM ) 2 − a 2 cos 2 θ 。该曲面在两极(θ = 0 , π \theta=0,\pi θ = 0 , π )与视界相切,在赤道(θ = π / 2 \theta=\pi/2 θ = π /2 )向外鼓出至 r = r s r=r_s r = r s ,它所围成的区域就是能层 ——一个位于视界之外 、但参考系拖曳 极端强烈的区域,在其中任何观测者若不超光速运动,都无法相对于远处恒星保持固定的 ( r , θ , ϕ ) (r,\theta,\phi) ( r , θ , ϕ ) :时空本身正沿着黑洞自旋方向被拖着旋转,能层内的一切都必须随之同向转动。由于能层位于视界之外,粒子进入后仍可逃逸到无穷远;彭罗斯(1969)证明,让粒子在能层内分裂,可使逃逸出的碎片带走比入射粒子更多的能量——从而开采黑洞的旋转能量。
进阶 时空的涟漪:线性化引力波 在远离源、时空近乎平直之处,度规可写成闵可夫斯基度规的小微扰 g μ ν = η μ ν + h μ ν g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu} g μν = η μν + h μν (∣ h μ ν ∣ ≪ 1 |h_{\mu\nu}|\ll1 ∣ h μν ∣ ≪ 1 )。在对迹反转微扰 h ˉ μ ν = h μ ν − 1 2 η μ ν h \bar h_{\mu\nu}=h_{\mu\nu}-\tfrac12\eta_{\mu\nu}h h ˉ μν = h μν − 2 1 η μν h 取调和(洛伦兹)规范 ∂ μ h ˉ μ ν = 0 \partial^\mu\bar h_{\mu\nu}=0 ∂ μ h ˉ μν = 0 时,爱因斯坦场方程线性化为标准波动方程 □ h ˉ μ ν = − 16 π G c 4 T μ ν \Box \bar h_{\mu\nu} = -\dfrac{16\pi G}{c^4}T_{\mu\nu} □ h ˉ μν = − c 4 16 π G T μν (其中 □ = − c − 2 ∂ t 2 + ∇ 2 \Box=-c^{-2}\partial_t^2+\nabla^2 □ = − c − 2 ∂ t 2 + ∇ 2 ):度规的微扰以光速 c c c 传播,由能量-动量张量 T μ ν T_{\mu\nu} T μν 作为源,正如电磁波由电流作为源一样。由于质量与动量守恒,孤立系统的单极矩和偶极矩不能辐射;主导辐射来自质量四极矩 Q i j Q_{ij} Q ij 的时间二阶导数——在功率中则是时间三阶导数——由此给出爱因斯坦关于辐射功率的四极公式 :P = G 5 c 5 ⟨ Q ... i j Q ... i j ⟩ P = \dfrac{G}{5c^5}\left\langle \dddot Q_{ij}\dddot Q^{ij}\right\rangle P = 5 c 5 G ⟨ Q ... ij Q ... ij ⟩ 前置因子 G / c 5 G/c^5 G / c 5 极其微小(∼ 3 × 10 − 53 W − 1 \sim3\times10^{-53}\text{ W}^{-1} ∼ 3 × 1 0 − 53 W − 1 ),这就是为什么只有像并合黑洞或中子星那样具有相对论性、太阳质量或更大尺度的系统,才能辐射出足以被探测到的强度。
历史注记
卡尔·史瓦西于1915年12月——爱因斯坦发表场方程仅仅一个月后——在第一次世界大战俄国前线服役期间找到了他的度规;几个月后他因自身免疫疾病去世。此后数十年间,曲面 r = r s r=r_s r = r s 被广泛称为“史瓦西奇点”并被怀疑不具物理实在性,直到大卫·芬克尔斯坦(1958)以及马丁·克鲁斯卡尔和乔治·塞凯赖什(1960)构造出跨越该处仍正则的坐标系,证明它是一层单向膜——事件视界——而非几何的崩溃。罗伊·克尔于1963年找到旋转解,罗杰·彭罗斯于1965年证明了第一个奇点定理(该工作获2020年诺贝尔物理学奖),而2015年9月LIGO探测器记录到了GW150914——来自双黑洞并合的第一个直接引力波信号。
大卫·希尔伯特
研究 开放问题与前沿方向 研究前沿 截至 2026 年
截至2026年的三个活跃前沿:(1) 克尔黑洞的非线性稳定性。 从数学上证明受微小扰动的克尔黑洞会弛豫回另一个克尔度规——而非产生裸奇点或辐射掉其结构——是数学相对论中最重大的偏微分方程问题之一。米哈利斯·达费尔莫斯、古斯塔夫·霍尔泽格尔、伊戈尔·罗德尼扬斯基和马丁·泰勒证明了不自转史瓦西族的完全非线性稳定性(2021年预印本),而塞尔日·克莱内曼、热雷米·谢夫特尔和埃莱娜·焦尔吉完成了慢速旋转 克尔黑洞(∣ a ∣ / M ≪ 1 |a|/M\ll1 ∣ a ∣/ M ≪ 1 )非线性稳定性的里程碑式证明(2021–2022年预印本,2023–2024年发表);涵盖天体物理中常见的快速自转黑洞的完整亚极端区间 ∣ a ∣ < G M / c 2 |a|<GM/c^2 ∣ a ∣ < GM / c 2 仍然开放。(2) 黑洞信息悖论与岛公式(island formula)。 自2019年起(杰弗里·彭宁顿;艾哈迈德·阿尔姆海里、内塔·恩格尔哈特、唐纳德·马洛夫与亨利·麦克斯菲尔德),利用“量子极值面”与岛 ——黑洞内部在佩奇时间之后其量子态被编码进远处霍金辐射的区域——所做的半经典引力计算,在玩具模型和高维模型中重现了蒸发黑洞的幺正佩奇曲线;将这些复本虫洞(replica wormhole)计算推广到无需精细调节假设的真实四维渐近平直黑洞,仍是活跃的研究前沿。(3) 低频引力波。 在LIGO–Virgo–KAGRA于 10 10 10 –10 3 Hz 10^3\text{ Hz} 1 0 3 Hz 频段观测恒星级质量并合的同时,脉冲星计时阵列(NANOGrav、EPTA、PPTA、CPTA)于2023年报告了来自超大质量双黑洞的纳赫兹随机引力波背景的首批证据,而空间引力波探测器LISA (2024年获ESA正式立项,计划于2030年代发射)将打开毫赫兹频段,高精度测绘极端质量比旋近系统的克尔几何。
若黑洞质量 M M M 增大到原来的三倍,其史瓦西半径 r s r_s r s 、光子球半径 r ph r_{\text{ph}} r ph 以及最内稳定圆轨道半径 r ISCO r_{\text{ISCO}} r ISCO 如何变化?
三者都变为原来的三倍,因为每一个都与 M M M 成正比。 r s r_s r s 变为三倍,而 r ph r_{\text{ph}} r ph 和 r ISCO r_{\text{ISCO}} r ISCO 按 3 2 = 9 3^2=9 3 2 = 9 增长。由于逃逸速度带有平方根,r s r_s r s 按 3 \sqrt{3} 3 增长。 r s r_s r s 变为三倍,而 r ph r_{\text{ph}} r ph 和 r ISCO r_{\text{ISCO}} r ISCO 是与 M M M 无关的普适常数。为什么我们知道史瓦西度规中的 r = r s r=r_s r = r s 是坐标奇点,而非物理上的曲率奇点?
因为克雷奇曼标量 K = 48 G 2 M 2 / ( c 4 r 6 ) K=48G^2M^2/(c^4r^6) K = 48 G 2 M 2 / ( c 4 r 6 ) 等曲率不变量在 r = r s r=r_s r = r s 处保持有限,并且正则坐标系(如埃丁顿-芬克尔斯坦坐标或克鲁斯卡尔-塞凯赖什坐标)可以光滑地跨越该处。 因为 K K K 在 r = 0 r=0 r = 0 以外处处恒为零。 因为光线会在 r = r s r=r_s r = r s 处反弹,而不是穿过它。 因为彭罗斯定理禁止在非零半径处出现奇点。 事件视界望远镜在M87和人马座A 周围拍到的暗“阴影”,其大小是由什么决定的?
掠过光子球(r ph = 1.5 r s r_{\text{ph}}=1.5\,r_s r ph = 1.5 r s )的光子的临界碰撞参数,经引力透镜效应放大到 r shadow = 3 3 2 r s ≈ 2.6 r s r_{\text{shadow}}=\tfrac{3\sqrt3}{2}r_s\approx2.6\,r_s r shadow = 2 3 3 r s ≈ 2.6 r s 。 半径为 r s r_s r s 的事件视界本身的几何截面。 位于 r ISCO = 3 r s r_{\text{ISCO}}=3r_s r ISCO = 3 r s 处的吸积盘物理内缘,不含透镜修正。 黑洞的康普顿波长。 尽管双黑洞并合在瞬间释放的功率超过可观测宇宙中所有恒星的总和,为什么LIGO在地球上测得的引力波应变却如此微小(h ∼ 10 − 21 h\sim10^{-21} h ∼ 1 0 − 21 )?
因为应变按 h ∼ ( r s / D ) ( v / c ) 2 h\sim (r_s/D)(v/c)^2 h ∼ ( r s / D ) ( v / c ) 2 标度:时空非常“坚硬”(巨大的功率在近区也只转化为有限的曲率),而且波随后在数百兆秒差距的距离上按 1 / D 1/D 1/ D 稀释。 因为引力波在传向地球的途中被星际尘埃吸收了。 因为引力波随距离按指数 h ∼ e − D / r s h\sim e^{-D/r_s} h ∼ e − D / r s 衰减。 因为单极辐射和偶极辐射几乎完全抵消了四极波。