CoderXL's Blog

Back

黑洞视界、逐渐红移的信号与穿过视界的世界线示意黑洞视界、逐渐红移的信号与穿过视界的世界线示意

随笔 →

如果一个宇航员落入黑洞,远处的我们会看到什么?常见的回答是:“他会越来越慢,最后永远冻结在事件视界上。”

这句话抓住了一部分数学,却容易让人想象出一张永久悬停的照片。更准确的图景是:他在有限的自身时间内穿过视界;我们收到的视界外信号不断红移、延迟、变弱,最终难以探测。穿越事件本身的光不会传到外部。 坐标时间的发散与有限固有时并不矛盾。MIT 8.962 的第 21–22 讲 ↗讨论了这个区别。

先换一只钟,做一次实验#

下面让宇航员从某个半径相对当地静止观察者静止释放,随后关闭发动机,径向自由落体。他向外发出固定自身频率的光,我们在 rO=50rsr_O=50r_s 处悬停接收。这个位置已相当远,但仍保留其引力时间膨胀的修正。

THOUGHT EXPERIMENT · 01

最后一束光,什么时候抵达?

史瓦西 · 径向自由落体

宇航员从静止释放,每隔相同的自身时间向外发光。你可以换一只钟,重新观看这段旅程。

正在计算尺度…

r 是球面的面积半径,r₀ 表示到中心的坐标半径;视界外的坐标间隔为 r₀ − rₛ。调整质量或光速时保持 r₀/rₛ 不变。现实中的局部光速恒为 c₀。

落入黑洞的发射事件与远处收到的光 左边表示宇航员所在的发射位置,右边表示该事件的信号。两侧对照同一发射事件,径向刻度经过压缩。 事件视界 r = rₛ 圆环是边界标记 01 / 发射事件 r = 6 rₛ 径向刻度压缩 · 非照片 接收者 50 rₛ 02 / 收到的信号 可见光 波形与颜色为示意

两侧对应同一发射事件,不代表远程同时发生。右侧变暗只用相对 g² 表示能量与到达率的变化。

准备就绪

发射半径 r / rₛ6.0000
宇航员固有时 τ0 s
接收时延 ΔτO0 s
频率比 g = νO / νe—
接收周期 / 发射周期—
收到的波长(发射 550 nm)—
发射事件的坐标时间 t0 s
运动钟率 dτ / dt—
光束相对功率 (g/g₀)²—

动画使用无量纲时间;1× 用约 24 秒展示一次旅程,读数给出对应的物理时间。

红移:引力与运动共同作用

自由落体 悬停对照
频率比随发射者固有时变化;悬停曲线表示在同一半径悬停的另一名发射者 0 0.25 0.5 0.75 1 g01发射者 τ / τH(1 = 穿越视界)

悬停对照 g = √(f(r)/f(rO)),每个点对应一名固定在该半径的发射者;自由落体多出向内运动的多普勒红移。它们不是同一条世界线上的两个时钟。

可以按这个顺序探索:

  1. 用默认的“远处接收者的钟”播放。位置趋近视界,信号的周期拉长,波长移入红外。
  2. 切换到“宇航员随身的钟”。同一个宇航员会穿过视界,继续抵达经典解的奇点。
  3. 重置后增大黑洞质量。无量纲曲线相同,实际半径和秒数变大。
  4. 改变释放半径,比较红移曲线与穿越时间;再试试假想光速,观察 rs∝c−2r_s\propto c^{-2} 与 rs/c∝c−3r_s/c\propto c^{-3}。

“重置”保留当前参数,“恢复默认参数”恢复整个实验。改变物理参数会从释放事件重新开始。动画用约 24 秒展示一段无量纲旅程,实际时间由读数给出;接收者模式的终点是 r/rs−1≈10−8r/r_s-1\approx10^{-8} 的显示精度,不是到达视界。

一张度规,三种时间#

史瓦西解描述什么#

只考虑球对称、不旋转、不带电的黑洞,外部没有吸积物质。宇航员的质量足够小,不改变背景时空。这时可以用史瓦西度规:

ds2=−f(r)c2dt2+dr2f(r)+r2dΩ2,f(r)=1−rsr,rs=2GMc2.ds^2=-f(r)c^2dt^2+\frac{dr^2}{f(r)}+r^2d\Omega^2, \qquad f(r)=1-\frac{r_s}{r}, \qquad r_s=\frac{2GM}{c^2}.

这里 GG 是引力常数,MM 是黑洞质量,dΩ2=dθ2+sin⁡2θ dϕ2d\Omega^2=d\theta^2+\sin^2\theta\,d\phi^2。rr 是面积半径:固定 rr 的球面面积为 4πr24\pi r^2。它不是从中心出发用尺子量出的距离。比如在外部某个 tt 恒定的切片上,径向固有长度元为 dℓ=dr/f(r)d\ell=dr/\sqrt{f(r)}。

视界位于 r=rsr=r_s。采用现实光速 c0=299 792 458 m/sc_0=299\,792\,458\ \mathrm{m/s},有

rs≈2.953 km MM⊙,rsc0≈9.851 μs MM⊙.r_s\approx2.953\ \mathrm{km}\,\frac{M}{M_\odot}, \qquad \frac{r_s}{c_0}\approx9.851\ \mathrm{\mu s}\,\frac{M}{M_\odot}.

这些是实验把无量纲轨迹换成公里和秒的尺度。光速滑块是在固定 GG、MM 时换一个假想模型中的常数 cc;并不表示真实光线在当地变慢。任何局部惯性系测到的真空光速仍是 c0c_0。度规与面积半径的定义可参见 David Tong,General Relativity,第 6 章 ↗。

哪只钟在走#

时间如何定义用来回答什么
宇航员固有时 τ\tau沿他的世界线,ds2=−c2dτ2ds^2=-c^2d\tau^2他自己经历了多久?
史瓦西坐标时间 tt在无穷远静止钟的极限下具有钟读数的含义在这张坐标地图中,事件标在哪里?
接收者固有时 τO\tau_O固定 rOr_O 的钟满足 dτO=f(rO) dtOd\tau_O=\sqrt{f(r_O)}\,dt_O光真正抵达时,接收者的钟走了多久?

尤其要区分:发射事件的 tt 与光抵达时的 tOt_O。光从宇航员走到接收者,还需要传播时间;而且宇航员正在移动,这段传播时间也在变化。

红移不只有“引力让钟变慢”#

先看悬停的发射者#

对固定在视界外某个半径的钟,dr=dΩ=0dr=d\Omega=0,所以

dτstaticdt=f(r).\frac{d\tau_{\mathrm{static}}}{dt}=\sqrt{f(r)}.

若发射者和接收者都悬停,接收频率与发射频率的比值是

gstatic≡νOνe=f(r)f(rO).g_{\mathrm{static}}\equiv\frac{\nu_O}{\nu_e} =\sqrt{\frac{f(r)}{f(r_O)}}.

越靠近视界,gg 越小,波长越长。这个比较要求发射者保持同一个半径;越靠近视界,保持悬停需要的加速度越大,视界上不存在这样的静止钟。静止观察者的红移关系可参见 Harvey Reall,Black Holes,§2.2 ↗。

自由落体还带着多普勒红移#

令 r0r_0 是静止释放半径,定义无量纲守恒能量与向内的径向速度量

E=1−rsr0,q(r)=E2−f(r)=rsr−rsr0.E=\sqrt{1-\frac{r_s}{r_0}}, \qquad q(r)=\sqrt{E^2-f(r)} =\sqrt{\frac{r_s}{r}-\frac{r_s}{r_0}}.

径向类时测地线满足

drdτ=−cq(r),dtdτ=Ef(r).\frac{dr}{d\tau}=-cq(r), \qquad \frac{dt}{d\tau}=\frac{E}{f(r)}.

下面从这些标准测地线关系推导实验中使用的红移。当地悬停观察者测到宇航员的向内速率为 βc\beta c,其中 β=q/E\beta=q/E。宇航员向外发光,运动方向与光相反,于是还要乘上狭义相对论的退行多普勒因子:

g(r)=f(r)f(rO)1−β1+β=f(r)f(rO) [E+q(r)](r>rs).\boxed{ g(r)=\sqrt{\frac{f(r)}{f(r_O)}} \sqrt{\frac{1-\beta}{1+\beta}} =\frac{f(r)}{\sqrt{f(r_O)}\,[E+q(r)]} } \qquad (r>r_s).

红移量常写成 z=g−1−1z=g^{-1}-1。相邻波峰的接收周期拉长为 1/g1/g 倍,收到的波长为 λO=λe/g\lambda_O=\lambda_e/g。这里的周期关系是瞬时频率或极短间隔的极限;较长时间间隔的两束脉冲需要分别计算到达时刻。

注意另一条关系:

dτdt=f(r)E≠g(r).\frac{d\tau}{dt}=\frac{f(r)}{E}\ne g(r).

左边比较宇航员固有时与发射事件的坐标时间;右边比较我们实际收到的频率。它们不能互换。比如默认 r0=6rsr_0=6r_s 时,在 r=2rsr=2r_s,悬停对照的 g≈0.714g\approx0.714,自由落体的 g≈0.339g\approx0.339,而 dτ/dt≈0.548d\tau/dt\approx0.548。三个数分别回答三个问题。

作为另一个边界条件的例子:若从无穷远静止开始落下,E=1E=1;再让接收者也趋于无穷远,公式化为 g=1−rs/rg=1-\sqrt{r_s/r}。实验采用有限半径静止释放,不能直接套用这个特例。

光越晚发出,越晚抵达#

径向光线满足 ds2=0ds^2=0,在视界外有

drdt=±cf(r).\frac{dr}{dt}=\pm c f(r).

积分时自然出现龟坐标(tortoise coordinate):

r∗(r)=r+rsln⁡∣rrs−1∣,dr∗dr=1f(r).r_*(r)=r+r_s\ln\left|\frac{r}{r_s}-1\right|, \qquad \frac{dr_*}{dr}=\frac1{f(r)}.

从 re>rsr_e>r_s 向外发出的光,到达 rOr_O 的坐标时刻是

tO=te+r∗(rO)−r∗(re)c.t_O=t_e+\frac{r_*(r_O)-r_*(r_e)}{c}.

把释放时发出的第一束光的接收时刻设为接收钟的零点,就得到实验中显示的“接收时延”:

ΔτO=f(rO)[te+r∗(r0)−r∗(re)c].\boxed{ \Delta\tau_O=\sqrt{f(r_O)} \left[t_e+\frac{r_*(r_0)-r_*(r_e)}{c}\right] }.

所以读数为零并不是光瞬间传到了远处,而是减去了第一束光的传播时间。

还可以直接求导检查红移公式。沿着下落轨迹,dre/dτ=−cqdr_e/d\tau=-cq,于是

dtOdτ=E+qf(re),dτdτO=f(re)f(rO) (E+q)=g.\frac{dt_O}{d\tau} =\frac{E+q}{f(r_e)}, \qquad \frac{d\tau}{d\tau_O} =\frac{f(r_e)}{\sqrt{f(r_O)}\,(E+q)}=g.

这给出一个直观理解:红移因子就是发射时钟相对于接收时钟的瞬时速率,前提是按同一束光把发射与接收事件对应起来。

若 ϵ=r/rs−1→0+\epsilon=r/r_s-1\to0^+,在这个自由落体模型下,忽略与 ϵ\epsilon 无关的常数,

te∼−rscln⁡ϵ,ΔτO∼−2f(rO)rscln⁡ϵ,g∼ϵ2Ef(rO).t_e\sim-\frac{r_s}{c}\ln\epsilon, \qquad \Delta\tau_O\sim-2\sqrt{f(r_O)}\frac{r_s}{c}\ln\epsilon, \qquad g\sim\frac{\epsilon}{2E\sqrt{f(r_O)}}.

接收时间比发射事件的坐标时间多了一份对数延迟。视界外越接近穿越的连续信号,越迟才收到,频率也越低。动画反解这个到达时间函数,而不是把位置按屏幕时间匀速推向黑洞。

他自己多久能穿过视界#

用 x0=r0/rsx_0=r_0/r_s,可以把静止释放的径向轨迹写成参数解:

r(η)rs=x0cos⁡2η2,τ(η)=rscx03/22(η+sin⁡η),0≤η≤π.\frac{r(\eta)}{r_s}=x_0\cos^2\frac\eta2, \qquad \tau(\eta)=\frac{r_s}{c}\frac{x_0^{3/2}}{2} (\eta+\sin\eta), \qquad 0\le\eta\le\pi.

在 η=0\eta=0 释放;在

ηH=2arccos⁡1x0,τH=τ(ηH)\eta_H=2\arccos\frac1{\sqrt{x_0}}, \qquad \tau_H=\tau(\eta_H)

处穿过视界。τH\tau_H 明确有限。沿这条理想轨迹,到达 r=0r=0 的固有时也有限:

τsing=rscπ2x03/2.\tau_{\mathrm{sing}}=\frac{r_s}{c}\frac\pi2x_0^{3/2}.

固定 r0/rsr_0/r_s 时,这些时间正比于 M/c3M/c^3。这解释了为什么调整质量会改变实际秒数,却不会改变归一化后的世界线形状。

为了让实时解算能平滑穿过视界,实验使用入射 Eddington–Finkelstein 时间

v=t+r∗c,ds2=−f(r)c2dv2+2c dv dr+r2dΩ2.v=t+\frac{r_*}{c}, \qquad ds^2=-f(r)c^2dv^2+2c\,dv\,dr+r^2d\Omega^2.

这张坐标图在 r=rsr=r_s 没有发散。沿宇航员世界线还有一个适合计算的关系:

dvdτ=E−qf(r)=1E+q.\frac{dv}{d\tau}=\frac{E-q}{f(r)}=\frac1{E+q}.

最右边在视界处仍然有限。组件采用上面的解析 r(η)r(\eta)、τ(η)\tau(\eta),数值积分这个平滑关系,并在视界外重建 tt 与到达时间。坐标构造可参见 Tong 第 6 章的史瓦西黑洞部分 ↗。

三张时空图,读同一条世界线#

CAUSAL ATLAS · 02

换一张地图,视界就不再是终点

宇航员世界线
启用 JavaScript 后显示坐标图

纵轴朝上表示未来。史瓦西时间在视界处发散,换成平滑的坐标后,世界线可以继续。

与上方实验联动:改变释放半径会重画世界线,播放或拖动进度会移动圆点。光线按零测地线绘制;箭头指向未来。图中省略两个角度维度。窄屏可在图内横向滚动。

史瓦西图的横轴是 rr,纵轴是 tt。曲线接近视界时不断向上延伸;光线的坐标速率也趋于零。这体现了该坐标在视界处的局限,局部测量的光速并没有变成零。

穿越视界图使用 t~=v−r/c\tilde t=v-r/c,并将释放事件的坐标时间平移为零。入射光线满足 dt~/dr=−1/cd\tilde t/dr=-1/c;另一族光线满足

drdt~=cr−rsr+rs.\frac{dr}{d\tilde t}=c\frac{r-r_s}{r+r_s}.

视界上这族光线沿着视界;视界内部则指向更小的 rr。所有未来因果方向都无法把内部信息带回外部。

彭罗斯图把无穷远压缩到有限图框,同时保留因果关系。它不是普通距离图:不能从纸上的长度读取真实的路程或时间。这里从龟坐标构造零坐标与 Kruskal 坐标,再做反正切压缩:

u=t−r∗c,U=−e−cu/(2rs),V=ecv/(2rs),u=t-\frac{r_*}{c},\qquad U=-e^{-cu/(2r_s)},\qquad V=e^{cv/(2r_s)}, X=arctan⁡V−arctan⁡Uπ,T=arctan⁡V+arctan⁡Uπ.X=\frac{\arctan V-\arctan U}{\pi}, \qquad T=\frac{\arctan V+\arctan U}{\pi}.

上式的 UU 表达式写的是外部;平滑延拓后,未来黑洞内部 U>0U>0。组件对 U,VU,V 做了一次不改变因果结构的常数缩放,使穿越点便于观察。

在这张图里,径向光线是 45∘45^\circ,类时世界线落在光锥内。未来视界是一条光线边界,宇航员穿过它后走向上方的类空奇点。下方的白洞与左侧的第二外部属于永恒史瓦西解的最大延拓,不是通常恒星坍缩形成的黑洞所必需的区域,也不意味着可以穿越虫洞旅行。延拓与因果图见 Reall,§2.11 与 §5.1 ↗。

“看不见穿越”不等于“永久看得见人”#

图中的发射波长为 550 nm。红移会先把它推向可见光的红端,再推入红外;继续红移后,就不再能用“越来越红的可见照片”来理解。人体图标在红外阶段保留的红色是教学提示。

同一束理想的径向光中,单个光子的能量缩小为 gg 倍,光子到达率再缩小为 gg 倍,所以只考虑这两项时,功率比例为 g2g^2。动画用 (g/g0)2(g/g_0)^2 表示相对暗化。真实成像还要考虑几何扩散、逃逸立体角、引力透镜、探测波段和灵敏度,因此图标亮度不等于望远镜测到的精确通量。

所谓“无限延迟”是连续信号在穿越前的极限。若宇航员每隔一个有限固有时间发一个短脉冲,穿越前只发出有限个脉冲;最后一个视界外脉冲可以在有限时刻收到。视界上和内部发出的脉冲则不会抵达外部。连续红移公式描述的是相邻波峰的极限,而不是要求现实中有无穷多个可检测光子。

视界不是墙,奇点才是这个模型的终点#

视界是因果边界,没有实体表面。换一张坐标图能消掉视界处的坐标发散,却消不掉 r=0r=0 的曲率发散。例如 Kretschmann 标量为

K=RμνρσRμνρσ=48G2M2c4r6=12rs2r6.K=R_{\mu\nu\rho\sigma}R^{\mu\nu\rho\sigma} =\frac{48G^2M^2}{c^4r^6} =\frac{12r_s^2}{r^6}.

r=rsr=r_s 处 KK 有限,r→0r\to0 时发散。局部是否会先被潮汐力撕裂,还取决于黑洞质量;径向长度 LL 上的潮汐加速度差量级为 Δa∼2GML/r3\Delta a\sim2GML/r^3,在视界处约为 c2L/rs2c^2L/r_s^2。固定现实光速时,质量越大的黑洞,视界附近的这项潮汐效应越弱。

这个实验把宇航员当作理想测试粒子,忽略自旋、吸积盘、非径向运动、辐射反作用与量子效应。它可以讲清自由落体与信号传播;动画到达经典奇点时停止,不能据此推断奇点之后的物理。

延伸阅读#

这里的图和推导围绕同一个静止释放模型绘制;阅读其他资料时,先核对释放速度、发射方向与接收者的位置,避免把不同边界条件的公式混用。

Comment seems to stuck. Try to refresh?✨