MathLabs

数学物理

黑洞的几何

爱因斯坦场方程的史瓦西解与克尔解:视界、光子球与最内稳定圆轨道、光线偏折、奇点定理,以及引力波,并与事件视界望远镜和LIGO相联系。

直观坠入时空中的一口井

把时空想象成不是一个空荡荡的舞台,而是一张会伸缩的橡皮膜:一个质量放在上面会把膜压出一个凹陷,而附近的一切——行星、光,甚至时间本身——都必须穿过那个凹陷运动。把足够多的质量塞进足够小的区域,凹陷就变成一口没有底的井:掉进去的光再也爬不出来。那个无法返回的区域就是黑洞,它的边界叫事件视界。本页将精确算出这口井究竟有多深,什么东西会掉进去、什么东西只是被它弯曲绕行,以及我们如今如何能看到这口井投下的影子,并听到两口这样的井相撞时发出的涟漪。

一个漏斗形曲面,远离中心处向外展开并趋于平坦,靠近中心处在一个狭窄的喉部附近急剧向下弯曲,表示空间在黑洞事件视界附近如何被拉伸。
弗拉姆抛物面 z(r)=2rs(r−rs)z(r) = 2\sqrt{r_s(r-r_s)}:将质量周围的空间几何嵌入到普通三维空间中的一种方式,展示了径向距离在 r=rsr=r_s 附近如何被拉伸。逐渐变陡的漏斗只是一种图示,并非真实的四维几何——但 r=rsr=r_s 处的喉部是真实的:那正是事件视界。

中学这口井有多大?史瓦西半径

定义: 逃逸速度与史瓦西半径

质量为 MM、半径为 rr 的球体的牛顿逃逸速度为 vesc=2GM/rv_{\text{esc}}=\sqrt{2GM/r}:飞离到无穷远所需的最小速度。令 vesc=cv_{\text{esc}}=c(光速)并解出 rr,得到一个特殊半径——史瓦西半径 rs=2GMc2r_s = \dfrac{2GM}{c^2}——在此半径处连光也无法逃脱。这个牛顿式的捷径给出了正确的数值,尽管光无法逃脱的真正原因(时空曲率,而不是光必须对抗的某种"引力")完全是相对论性的;下面给出严谨的推导。

rs=2GMc2r_s = \dfrac{2GM}{c^2}

这里 GG 是牛顿引力常数,MM 是质量,cc 是光速。对太阳而言,rs≈2.95 kmr_s\approx 2.95\text{ km}——太阳需要被压缩到一个小镇大小才能变成黑洞;对地球而言,rs≈8.9 mmr_s\approx 8.9\text{ mm},比一颗葡萄还小。

大学史瓦西度规

定义: 史瓦西度规

1916年,卡尔·史瓦西找到了爱因斯坦场方程唯一的球对称真空解。在坐标 (t,r,θ,ϕ)(t,r,\theta,\phi) 中,时空间隔为 ds2=−(1−rsr)c2 dt2+(1−rsr)−1dr2+r2 dΩ2ds^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Ω2=dθ2+sin⁡2θ dϕ2d\Omega^2 = d\theta^2+\sin^2\theta\,d\phi^2 是单位球面上的度规,rs=2GM/c2r_s=2GM/c^2 同上。它描述任意球对称、不自转质量——恒星、行星或黑洞——外部的时空,并在 r→∞r\to\infty 时退化为平直的闵可夫斯基时空。

ds2=−(1−rsr)c2 dt2+(1−rsr)−1dr2+r2 dΩ2ds^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

有两个半径看起来是奇异的:gttg_{tt} 和 grrg_{rr} 都发散的 r=0r=0,以及 gtt→0g_{tt}\to0 且 grr→∞g_{rr}\to\infty 的 r=rsr=r_s。这其中只有一个是真正的物理奇点。克雷奇曼标量 K=RαβγδRαβγδ=48 G2M2c4r6K = R_{\alpha\beta\gamma\delta}R^{\alpha\beta\gamma\delta} = \dfrac{48\,G^2M^2}{c^4 r^6} 是由完整黎曼张量构造出的、与坐标无关的曲率度量;它在 r=0r=0 处发散(真实、不可避免的无穷曲率),但在 r=rsr=r_s 处却完全有限。因此 r=rsr=r_s 处的发散只是坐标奇点——史瓦西坐标系的产物,而非时空本身的缺陷——正如地球仪经线在北极的"奇异性"只是经纬度坐标的产物,而非球面本身的缺陷。改用埃丁顿-芬克尔斯坦坐标(将 tt 替换为内向零坐标 v=t+r∗v=t+r_*,r∗=r+rsln⁡∣r/rs−1∣r_*=r+r_s\ln|r/r_s-1|)或完全扩展的克鲁斯卡尔-塞凯赖什坐标,可使所有度规分量在 r=rsr=r_s 处都保持正则,揭示出它是一个完全光滑的零曲面:事件视界。

史瓦西黑洞周围的关键半径
量公式以 rsr_s 为单位
事件视界rs=2GM/c2r_s=2GM/c^21 rs1\,r_s
光子球(不稳定圆形光轨道)rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2}1.5 rs1.5\,r_s
最内稳定圆轨道(ISCO)rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2}3 rs3\,r_s
影子半径(EHT成像所得)rshadow=332rs≈2.6 rsr_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s≈2.6 rs\approx2.6\,r_s

大学黑洞周围的轨道:光子球与ISCO

无论是光线还是有质量粒子,在史瓦西几何中运动时都沿其测地线守恒两个量——一个类似能量的常数和一个类似角动量的常数——这与开普勒问题完全一样,但有效势中带有一个牛顿理论没有对应物的相对论性扭曲,正是它产生了两个具有里程碑意义的半径。

对于史瓦西黑洞,光子的圆轨道仅存在于 rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2} 处,且总是不稳定的;有质量粒子的圆轨道只有在 r≥rISCOr \ge r_{\text{ISCO}} 时才稳定,在 rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} 处边缘稳定,在 3rs/2<r<3rs3r_s/2<r<3r_s 区间内不稳定。

为什么成立?

把径向运动想象成一个球在一维势阱 Veff(r)V_{\text{eff}}(r) 中滚动:圆轨道位于势阱平坦之处(dVeff/dr=0dV_{\text{eff}}/dr=0),且只有在势阱向上弯曲(d2Veff/dr2>0d^2V_{\text{eff}}/dr^2>0)而非向下弯曲之处才稳定。由于相对论性势中多出了一个牛顿引力所没有的 1/r31/r^3 项,势阱在靠近黑洞处会出现一个极大值——在该半径以内根本不存在任何稳定圆轨道,无论粒子转动多快;物质只会直接坠落进去。

证明

第一步(守恒量与径向方程)。 沿限制在赤道面 θ=π/2\theta=\pi/2 上的任意测地线(类时 κ=−c2\kappa=-c^2,或零测地线 κ=0\kappa=0),度规不依赖于 tt 和 ϕ\phi 给出两个守恒量 E~=(1−rs/r)c2t˙\tilde E = (1-r_s/r)c^2\dot t 和 L~=r2ϕ˙\tilde L = r^2\dot\phi(点 =d/dτ=d/d\tau)。将其代入归一化条件 gμνx˙μx˙ν=κg_{\mu\nu}\dot x^\mu\dot x^\nu=\kappa 并化简,对类时运动得到 (drdτ)2=E~2−(1−rsr)(c2+L~2r2)≡E~2−Veff2(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)

第二步(圆轨道条件)。 展开得到 Veff2(r)=c2−c2rsr+L~2r2−L~2rsr3V^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} 圆轨道要求 rr 恒定,即任意时刻都有 dr/dτ=0dr/d\tau=0,这要求 E~2=Veff2(r)\tilde E^2=V^2_{\text{eff}}(r) 且(为使 rr 不发生漂移)dVeff2/dr=0dV^2_{\text{eff}}/dr=0:ddrVeff2=c2rsr2−2L~2r3+3L~2rsr4=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

第三步(解出圆轨道的角动量)。 两边乘以 r4r^4 并解出 L~2\tilde L^2,得到维持半径 rr 处圆轨道所需的角动量:L~2(r)=c2rsr22r−3rs\tilde L^2(r) = \frac{c^2 r_s r^2}{2r-3r_s} 这已经包含了光子球的信息:当 r→(3/2)rsr\to (3/2)r_s 从上方趋近时,分母 2r−3rs→0+2r-3r_s\to0^+,L~2→+∞\tilde L^2\to+\infty——没有任何有限角动量能在如此靠近处维持圆轨道,这正是下面达到的无质量(光子)极限。

第四步(边缘稳定性——ISCO)。 当 rr 从无穷远减小时,L~2(r)\tilde L^2(r) 先减小,达到极小值,然后在 r=3rs/2r=3r_s/2 处发散;圆轨道稳定当且仅当进一步缩小 rr 需要增大 L~2\tilde L^2(即稳定性由 dL~2/drd\tilde L^2/dr 的符号决定)。对第三步的结果求导,dL~2dr=2c2rs r(r−3rs)(2r−3rs)2\frac{d\tilde L^2}{dr} = \frac{2c^2 r_s\,r(r-3r_s)}{(2r-3r_s)^2} 它在 r=0r=0(排除)、r=3rsr=3r_s 处为零,在 r>3rs/2r>3r_s/2 范围内再无其他零点。这个唯一的内部根就是边缘稳定半径:rISCO=3rs=6GMc2r_{\text{ISCO}} = 3r_s = \dfrac{6GM}{c^2} 当 r>3rsr>3r_s 时轨道稳定(dL~2/dr>0d\tilde L^2/dr>0,与开普勒问题一样,rr 越大需要的角动量越大);当 3rs/2<r<3rs3r_s/2<r<3r_s 时轨道不稳定。

第五步(由零测地线得到光子球)。 对零测地线,令 κ=0\kappa=0 作同样代换,得到 (drdλ)2=E2−L2r2(1−rsr)≡E2−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),其中 Vγ2(r)=L2(1r2−rsr3)V^2_\gamma(r) = L^2\left(\frac{1}{r^2}-\frac{r_s}{r^3}\right)。令括号内表达式的导数为零,ddr(1r2−rsr3)=−2r3+3rsr4=0  ⟹  r=3rs2\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} 直接给出圆形光子轨道的唯一半径,rph=32rs=3GMc2r_{\text{ph}} = \dfrac{3}{2}r_s = \dfrac{3GM}{c^2},与第三步得到的极限一致;它是不稳定的,因为括号 1/r2−rs/r31/r^2-r_s/r^3 在该处取局部极大值而非极小值——任何向内的扰动都会使光子螺旋坠入视界,任何向外的扰动都会使它逃逸到无穷远。

以碰撞参数 b≫rsb\gg r_s(最近点远在视界之外)掠过质量 MM 的光线,在最低阶上被偏折角度 Δϕ=4GMc2b\Delta\phi = \dfrac{4GM}{c^2 b}——恰好是朴素粒子式计算所预言的牛顿值 2GM/(c2b)2GM/(c^2b) 的两倍。

为什么成立?

牛顿引力若被当作作用在高速粒子上的力来处理,会预言一定的偏折,但它只考虑了时间的曲率(质量附近时钟变慢)。广义相对论额外加入了来自空间曲率的等量贡献(质量附近的尺子沿径向收缩),两种效应相加,使偏折角加倍——正是这个数值因子让1919年日食观测队得以区分爱因斯坦理论与牛顿理论。

证明

第一步(光的轨道方程)。 对限制在赤道面上的零测地线,用与光子球相同的守恒量 E=(1−rs/r)c2t˙E=(1-r_s/r)c^2\dot t 和 L=r2ϕ˙L=r^2\dot\phi,代入 gμνx˙μx˙ν=0g_{\mu\nu}\dot x^\mu\dot x^\nu=0,并将对 τ\tau 的导数转换为 u≡1/ru\equiv1/r 对 ϕ\phi 的导数(利用 r˙=−L du/dϕ\dot r = -L\,du/d\phi),再求一次导后得到标准的光线弯曲轨道方程 d2udϕ2+u=3rs2u2,u≡1r\frac{d^2u}{d\phi^2}+u = \frac{3r_s}{2}u^2, \qquad u\equiv\frac1r

第二步(零阶:直线)。 去掉右边项(令 rs=0r_s=0,平直时空)得到 u0′′+u0=0u_0''+u_0=0,其在 ϕ=π/2\phi=\pi/2 处以最近距离 bb 通过的解是直线 u0(ϕ)=sin⁡ϕbu_0(\phi) = \frac{\sin\phi}{b}——这是极坐标形式下未偏折的路径。

第三步(一阶微扰)。 写出 u=u0+rsu1u=u_0+r_s u_1 并代入第一步的方程,只保留关于 rsr_s 的线性项;u0′′+u0=0u_0''+u_0=0 部分抵消,剩下 u1u_1 满足的受迫线性振子方程:u1′′+u1=32b2sin⁡2ϕu_1''+u_1 = \frac{3}{2b^2}\sin^2\phi 一个与最近点两侧对称弯曲相匹配的特解是 u1(ϕ)=12b2(1+cos⁡2ϕ)u_1(\phi) = \frac{1}{2b^2}\left(1+\cos^2\phi\right),因此完整的一阶轨迹为 u(ϕ)=sin⁡ϕb+rs2b2(1+cos⁡2ϕ)u(\phi) = \frac{\sin\phi}{b} + \frac{r_s}{2b^2}\left(1+\cos^2\phi\right)

第四步(提取总偏折角)。 平直空间中 u→0u\to0 恰好发生在 ϕ=0\phi=0 和 ϕ=π\phi=\pi 处;有了 rsr_s 修正后,u→0u\to0 发生在 ϕ=−δ1\phi=-\delta_1 和 ϕ=π+δ2\phi=\pi+\delta_2 处,δ1,δ2\delta_1,\delta_2 为小量。将 u(−δ1)=0u(-\delta_1)=0 对 δ1\delta_1 和 rsr_s 展开到一阶,得到 −δ1/b+rs/b2=0-\delta_1/b+r_s/b^2=0,故 δ1=rs/b\delta_1=r_s/b;对另一端做相同计算得 δ2=rs/b\delta_2=r_s/b。出射渐近线偏离与入射渐近线严格反向的总角度为 Δϕ=δ1+δ2=2rsb=4GMc2b\Delta\phi = \delta_1+\delta_2 = \frac{2r_s}{b} = \frac{4GM}{c^2 b},代回 rs=2GM/c2r_s=2GM/c^2 即得所要证明的公式。

进阶奇点不可避:彭罗斯定理

史瓦西黑洞内部的 r=0r=0 奇点,可能看起来只是求解爱因斯坦方程时所假设的严格球对称性带来的特殊现象——或许对称性较低、更贴近现实的坍缩能够避开它。1965年,罗杰·彭罗斯证明了相反的结论:一旦形成囚陷面(一个闭合曲面,其上向外和向内两束光线都在会聚,而不像远离任何黑洞处那样只有向内的一束会聚),那么对于任何满足合理能量条件的物质,无论对称性如何,奇点都不可避免。这是第一个现代意义上的奇点定理,开启了如今主导数学相对论的整体、几何方法。

若时空包含一个非紧的柯西曲面,满足零能量条件,并包含一个闭合囚陷面,则时空在未来方向上是零测地不完备的:至少存在一条光线无法向未来无限延伸——这就是彭罗斯意义下的奇点。

为什么成立?

在广义相对论中,引力总是使光线会聚而非发散(这正是能量条件所刻画的内容);因此在某个闭合曲面上已经两侧都在会聚的光线束会持续会聚,而几何上会聚的光线必然相交——但在一个性质良好、整体双曲的时空内部,两条光线相交(共轭点)若不为该光线的存在制造一个边界,就无法自洽。该定理把这种会聚的直觉转化为一个严格的不等式(雷查德胡里方程)和一个整体拓扑障碍。

证明

第一步(用符号表述囚陷面)。 设 SS 是一个闭合的类空二维曲面,θ\theta 表示从 SS 正交出发的零测地线束的膨胀率——即光束横截面积增长的相对速率。远离任何质量处,外向光束膨胀(θout>0\theta_{\text{out}}>0),内向光束收缩(θin<0\theta_{\text{in}}<0)。若在整个 SS 上反而有 θout<0,θin<0\theta_{\text{out}}<0, \qquad \theta_{\text{in}}<0,则称 SS 是囚陷的:即使是向外的光也在被拖向内部。

第二步(雷查德胡里方程)。 沿切向量为 kak^a 的零测地同余族,对膨胀率的定义 θ=∇aka\theta=\nabla_a k^a 求导,并利用黎曼张量的定义交换导数次序,得到精确的运动学恒等式 dθdλ=−θ22−σabσab−Rabkakb\frac{d\theta}{d\lambda} = -\frac{\theta^2}{2} - \sigma_{ab}\sigma^{ab} - R_{ab}k^ak^b,其中 σab\sigma_{ab} 是该同余族的剪切张量。这个方程纯粹是几何性的——在假设任何物理之前,它由曲率的定义得出,说明任何一族光线必然如何弯曲。

第三步(能量条件迫使进一步会聚)。 物理上合理的物质满足零能量条件 Rabkakb≥0R_{ab}k^ak^b \ge 0,且剪切项 σabσab≥0\sigma_{ab}\sigma^{ab}\ge0 恒成立(它是平方和);从第二步恒等式右边去掉这两个非负项,得到不等式 dθdλ≤−θ22\frac{d\theta}{d\lambda} \le -\frac{\theta^2}{2}

第四步(里卡蒂不等式迫使在有限时间内发散)。 在囚陷面上,λ=0\lambda=0 时 θ0<0\theta_0<0。只要 θ\theta 保持为负,将不等式除以 θ2>0\theta^2>0 并整理,得到 ddλ(1θ)≥12  ⟹  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} 由于右边随 λ\lambda 增大而无界增长,而 1/θ1/\theta(为负)必须满足此下界,1/θ1/\theta 被迫趋向 0−0^-,即在某个有限仿射参数 λ∗≤2/∣θ0∣\lambda_*\le 2/|\theta_0| 处 θ→−∞\theta\to-\infty:与 SS 正交的零测地同余族必须在有限仿射时间内出现共轭点(散焦面,相邻光线在此相交)。

第五步(从局部会聚到整体不完备)。 因果结构理论的一个标准结果(此处引用但不重新推导)表明,整体双曲时空中的零测地线在经过共轭点之后不能保持因果边界性(不能停留在某集合未来的边界上)。结合非紧的柯西曲面,这就排除了测地线回环或以通常方式终止的可能,从而迫使至少一条与 SS 正交的零测地线未来不完备:它无法延伸到任意大的仿射参数。根据定义,这种不完备性正是彭罗斯所说的奇点——该定理表明,只要囚陷面形成,无需任何对称性假设,奇点就必然出现。

大学实际应用与典型例题

黑洞几何学在2010年代不再纯属理论。事件视界望远镜(EHT)将一个行星尺度的射电天线阵列作为单一干涉仪工作,拍摄黑洞投射在周围发光气体上的阴影——一个接近光子球相关半径 rshadow≈2.6 rsr_{\text{shadow}}\approx2.6\,r_s 的暗盘,而非视界本身。而LIGO–Virgo–KAGRA探测到两个黑洞螺旋靠近并合并时发出的啁啾引力波,其质量信息编码在波的频率扫描和极其微小的应变振幅 h∼10−21h\sim10^{-21} 之中。

例题: 计算人马座A*的阴影大小

人马座A*是我们星系中心的超大质量黑洞,质量 M≈4.0×106 M⊙M\approx4.0\times10^6\,M_\odot,距离 D≈8 kpcD\approx8\text{ kpc}。请估算其阴影的角直径,并与事件视界望远镜2022年测得的约 51.8 μas51.8\,\mu\text{as} 相比较。

解答

第一步(史瓦西半径)。 对一个太阳质量,rs(M⊙)=2GM⊙/c2≈2.95 kmr_s(M_\odot)=2GM_\odot/c^2\approx2.95\text{ km}。按质量线性缩放,rs≈4.0×106×2.95 km≈1.18×107 km=1.18×1010 mr_s\approx4.0\times10^6\times2.95\text{ km}\approx1.18\times10^7\text{ km}=1.18\times10^{10}\text{ m}。

第二步(阴影直径)。 阴影半径为 rshadow=332rs≈2.6 rsr_{\text{shadow}} = \dfrac{3\sqrt3}{2}r_s \approx 2.6\,r_s,因此直径为 2rshadow=2×2.6 rs≈5.2 rs≈6.14×1010 m2r_{\text{shadow}}=2\times2.6\,r_s\approx5.2\,r_s\approx6.14\times10^{10}\text{ m}。

第三步(换算为角度)。 小角近似下,角大小 == 物理大小 // 距离。取 D=8 kpc=8×3.086×1019 m=2.47×1020 mD=8\text{ kpc}=8\times3.086\times10^{19}\text{ m}=2.47\times10^{20}\text{ m},角度为 6.14×1010/2.47×1020≈2.49×10−10 rad6.14\times10^{10}/2.47\times10^{20}\approx2.49\times10^{-10}\text{ rad}。

第四步(把弧度换算为微角秒)。 1弧度等于 2.063×1011 μas2.063\times10^{11}\,\mu\text{as},因此阴影直径约为 2.49×10−10×2.063×1011 μas≈51.3 μas2.49\times10^{-10}\times2.063\times10^{11}\,\mu\text{as}\approx51.3\,\mu\text{as}——与EHT测得的 51.8 μas51.8\,\mu\text{as} 相差在 1%1\% 以内,有力地证实了人马座A*的行为与史瓦西几何(更准确地说是自旋接近零的克尔几何)的预言完全一致。

例题: GW150914的啁啾质量与应变数量级

LIGO首次探测到的GW150914,来自两个质量约为 m1≈36 M⊙m_1\approx36\,M_\odot 和 m2≈29 M⊙m_2\approx29\,M_\odot 的黑洞在距离 D≈410 MpcD\approx410\text{ Mpc} 处并合。计算控制波频率扫描的啁啾质量 M\mathcal M,并估算地球处应变 hh 的数量级。

解答

第一步(啁啾质量)。 M=(m1m2)3/5(m1+m2)1/5\mathcal M = \dfrac{(m_1 m_2)^{3/5}}{(m_1+m_2)^{1/5}} 由 m1m2≈1044 M⊙2m_1m_2\approx1044\,M_\odot^2 及 m1+m2≈65 M⊙m_1+m_2\approx65\,M_\odot:(m1m2)3/5≈10440.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},(m1+m2)1/5≈650.2≈2.3(m_1+m_2)^{1/5}\approx65^{0.2}\approx2.3,得 M≈100/2.3 M⊙≈28 M⊙\mathcal M\approx100/2.3\,M_\odot\approx28\,M_\odot——非常接近LIGO报告的约 28 M⊙28\,M_\odot。

第二步(源的长度尺度)。 总质量对应的史瓦西半径为 rs,tot=2G(m1+m2)/c2≈65×2.95 km≈1.9×105 mr_{s,\text{tot}}=2G(m_1+m_2)/c^2\approx65\times2.95\text{ km}\approx1.9\times10^5\text{ m}:这是并合临近时系统的尺寸。

第三步(数量级应变估计)。 临近并合时轨道速度接近光速 cc 的相当一部分;取 v/c∼0.5v/c\sim0.5,用标度关系 h∼rs,totD(vc)2h \sim \dfrac{r_{s,\text{tot}}}{D}\left(\dfrac{v}{c}\right)^2,并取 D=410 Mpc≈1.26×1025 mD=410\text{ Mpc}\approx1.26\times10^{25}\text{ m}:h∼(1.9×105/1.26×1025)×0.25≈3.8×10−21h\sim(1.9\times10^5/1.26\times10^{25})\times0.25\approx3.8\times10^{-21}。

第四步(比较)。 这一数量级估计与LIGO实际测得的应变 h∼10−21h\sim10^{-21} 相差不过几倍——这个微小的数字并不意味着源处效应微弱(临近并合时那里的时空弯曲程度几乎达到了物理所允许的极限),而只是引力辐射在数百兆秒差距的距离上按 1/D1/D 被稀释的结果。

进阶旋转黑洞:克尔几何与能层

真实的天体物理黑洞诞生于自转恒星或并合过程,带有角动量 JJ,可用自旋长度 a=J/(Mc)a=J/(Mc) 参数化。罗伊·克尔于1963年找到了唯一的稳态轴对称真空解。自转把单一的史瓦西曲面 r=rsr=r_s 分裂成两个不同的曲面:位于 r+=GM/c2+(GM/c2)2−a2r_+ = GM/c^2 + \sqrt{(GM/c^2)^2-a^2} 的外事件视界,以及其外侧更大、呈扁球状的静止极限面 rergo(θ)=GMc2+(GMc2)2−a2cos⁡2θr_{\text{ergo}}(\theta) = \dfrac{GM}{c^2} + \sqrt{\left(\dfrac{GM}{c^2}\right)^2 - a^2\cos^2\theta}。该曲面在两极(θ=0,π\theta=0,\pi)与视界相切,在赤道(θ=π/2\theta=\pi/2)向外鼓出至 r=rsr=r_s,它所围成的区域就是能层——一个位于视界之外、但参考系拖曳极端强烈的区域,在其中任何观测者若不超光速运动,都无法相对于远处恒星保持固定的 (r,θ,ϕ)(r,\theta,\phi):时空本身正沿着黑洞自旋方向被拖着旋转,能层内的一切都必须随之同向转动。由于能层位于视界之外,粒子进入后仍可逃逸到无穷远;彭罗斯(1969)证明,让粒子在能层内分裂,可使逃逸出的碎片带走比入射粒子更多的能量——从而开采黑洞的旋转能量。

进阶时空的涟漪:线性化引力波

在远离源、时空近乎平直之处,度规可写成闵可夫斯基度规的小微扰 gμν=ημν+hμνg_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}(∣hμν∣≪1|h_{\mu\nu}|\ll1)。在对迹反转微扰 hˉμν=hμν−12ημνh\bar h_{\mu\nu}=h_{\mu\nu}-\tfrac12\eta_{\mu\nu}h 取调和(洛伦兹)规范 ∂μhˉμν=0\partial^\mu\bar h_{\mu\nu}=0 时,爱因斯坦场方程线性化为标准波动方程 □hˉμν=−16πGc4Tμν\Box \bar h_{\mu\nu} = -\dfrac{16\pi G}{c^4}T_{\mu\nu}(其中 □=−c−2∂t2+∇2\Box=-c^{-2}\partial_t^2+\nabla^2):度规的微扰以光速 cc 传播,由能量-动量张量 TμνT_{\mu\nu} 作为源,正如电磁波由电流作为源一样。由于质量与动量守恒,孤立系统的单极矩和偶极矩不能辐射;主导辐射来自质量四极矩 QijQ_{ij} 的时间二阶导数——在功率中则是时间三阶导数——由此给出爱因斯坦关于辐射功率的四极公式:P=G5c5⟨Q...ijQ...ij⟩P = \dfrac{G}{5c^5}\left\langle \dddot Q_{ij}\dddot Q^{ij}\right\rangle 前置因子 G/c5G/c^5 极其微小(∼3×10−53 W−1\sim3\times10^{-53}\text{ W}^{-1}),这就是为什么只有像并合黑洞或中子星那样具有相对论性、太阳质量或更大尺度的系统,才能辐射出足以被探测到的强度。

研究开放问题与前沿方向

若黑洞质量 MM 增大到原来的三倍,其史瓦西半径 rsr_s、光子球半径 rphr_{\text{ph}} 以及最内稳定圆轨道半径 rISCOr_{\text{ISCO}} 如何变化?

为什么我们知道史瓦西度规中的 r=rsr=r_s 是坐标奇点,而非物理上的曲率奇点?

事件视界望远镜在M87和人马座A周围拍到的暗“阴影”,其大小是由什么决定的?

尽管双黑洞并合在瞬间释放的功率超过可观测宇宙中所有恒星的总和,为什么LIGO在地球上测得的引力波应变却如此微小(h∼10−21h\sim10^{-21})?

参考文献

  1. Charles W. Misner, Kip S. Thorne, John Archibald Wheeler (1973). Gravitation
  2. Roger Penrose (1965). Gravitational Collapse and Space-Time Singularities · DOI:10.1103/PhysRevLett.14.57
  3. B. P. Abbott et al. (LIGO Scientific Collaboration and Virgo Collaboration) (2016). Observation of Gravitational Waves from a Binary Black Hole Merger · DOI:10.1103/PhysRevLett.116.061102
  4. Event Horizon Telescope Collaboration (2022). First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way · DOI:10.3847/2041-8213/ac6674
  5. Sergiu Klainerman, Jérémie Szeftel (2021). Kerr stability for small angular momentum · arXiv:2104.11857 [预印本,未经同行评审]
  6. Geoffrey Penington (2019). Entanglement Wedge Reconstruction and the Information Paradox · arXiv:1905.08255 [预印本,未经同行评审]