跳到主要内容
Sandbox Physics
旋转黑洞与玻色云

守恒能流、物理单位与第一次增长估算

Sandbox Physics ·

克尔黑洞周围的波动方程守恒入射、反射、透射振幅的某个组合,单凭这一点就能证明:波被放大,当且仅当 0<ω<mΩH0<\omega<m\Omega_H。把黑洞质量放回去,以 1/M1/M 计的增长率就成了秒,耦合就成了以电子伏特计的粒子质量。一条简短的类氢公式给出云增长快慢的第一次估算,再拿严格解一对,就知道这个估算能信到什么程度。

打开黑洞超辐射Lab

本章沿用波、自转与克尔黑洞的尺度的结果:单位 G=c=M=1G=c=M=1,视界半径 r+r_+,视界转速 ΩH\Omega_H。

径向方程

克尔背景上质量为 μ\mu(以 1/M1/M 计)的标量场,可以分离成时间因子 e−iωte^{-i\omega t}、角向因子 eimφS(θ)e^{im\varphi}S(\theta) 和径向函数 R(r)R(r)。径向函数满足

ddr(ΔdRdr)+[K2Δ−λ−μ2r2]R=0,\frac{\mathrm d}{\mathrm dr}\left(\Delta\frac{\mathrm dR}{\mathrm dr}\right)+\left[\frac{K^2}{\Delta}-\lambda-\mu^2r^2\right]R=0,

其中 Δ=r2−2r+a2\Delta=r^2-2r+a^2,K=(r2+a2)ω−amK=(r^2+a^2)\omega-am,λ=Λ+a2ω2−2amω\lambda=\Lambda+a^2\omega^2-2am\omega。Λ\Lambda 来自角向方程,黑洞不转时等于 l(l+1)l(l+1),有自转时略有偏移。Dolan,《克尔时空中有质量克莱因—戈尔登场的不稳定性》第II节

Δ\Delta 在内外两个视界 r±=1±1−a2r_\pm=1\pm\sqrt{1-a^2} 上为零。靠近 r+r_+ 时 K2/ΔK^2/\Delta 一项占主导,而 K(r+)=2r+(ω−mΩH)K(r_+)=2r_+(\omega-m\Omega_H)。放大条件里的所有量,都装在这一个边界值里。

一个不会变的朗斯基量

取无质量的波 μ=0\mu=0,频率为实数。这时径向方程的系数全是实数,复共轭 R∗R^\ast 也满足同一个方程。用 R∗R^\ast 乘 RR 的方程,减去 RR 乘共轭方程,势能项正好抵消:

ddr[Δ(R∗R′−RR∗′)]=0.\frac{\mathrm d}{\mathrm dr}\left[\Delta\left(R^\ast R^\prime-RR^{\ast\prime}\right)\right]=0 .

方括号里是朗斯基量,在每个半径上都取同一个值,所以紧贴视界处的值必须等于远处的值。

换成 ψ=r2+a2 R\psi=\sqrt{r^2+a^2}\,R,并用乌龟坐标 r∗r_\ast 表示,计算更干净。r∗r_\ast 在视界处趋于 −∞-\infty,在远处趋于 +∞+\infty。物理上,视界里不能有东西出来,所以视界附近的波只往里掉:ψ→T e−ikHr∗\psi\to T\,e^{-ik_Hr_\ast},kH=ω−mΩHk_H=\omega-m\Omega_H。远处是入射波和反射波之和:ψ→I e−iωr∗+R eiωr∗\psi\to I\,e^{-i\omega r_\ast}+R\,e^{i\omega r_\ast},这里的 RR 指反射振幅。让两端的守恒量相等,得到

ω(∣I∣2−∣R∣2)=(ω−mΩH) ∣T∣2.\omega\left(|I|^2-|R|^2\right)=(\omega-m\Omega_H)\,|T|^2 .

Brito、Cardoso、Pani,《超辐射》第4.4节式(4.53)

左边是入射波与反射波之间少掉的能流。ω>0\omega>0 时,反射波比入射波强,即 ∣R∣>∣I∣|R|>|I|,当且仅当 ω<mΩH\omega<m\Omega_H。这个证明不用任何近似,也不需要小参数。增益用

Z=∣R∣2∣I∣2−1Z=\frac{|R|^2}{|I|^2}-1

来衡量。

例题1。 我们的程序取 a=0.99a=0.99、Mω=0.3M\omega=0.3,从视界向外积分径向函数。在 1.14M1.14M 到 267M267M 之间取六个半径,朗斯基量始终是 0.3053595840.305359584,相对起伏只有 3.6×10−103.6\times10^{-10}。自转0.99这一列的全部100个频率上,能流平衡两边的差都小于 4.2×10−114.2\times10^{-11}。任何放大程序都该先过这一关。完整的复算过程见复算一份公开的超辐射放大表。

从几何单位换回秒和电子伏特

换回普通单位:无量纲时间乘 GM/c3GM/c^3,长度乘 GM/c2GM/c^2,能量乘 ℏc3/(GM)\hbar c^3/(GM)。

GMc3=4.93 μs MM⊙,GMc2=1.48 km MM⊙,ℏc3GM=1.34×10−10 eV M⊙M.\frac{GM}{c^3}=4.93\,\mu\mathrm s\,\frac{M}{\mathrm M_\odot},\qquad \frac{GM}{c^2}=1.48\,\mathrm{km}\,\frac{M}{\mathrm M_\odot},\qquad \frac{\hbar c^3}{GM}=1.34\times10^{-10}\,\mathrm{eV}\,\frac{\mathrm M_\odot}{M}.

最后一个因子就是 ℏ\hbar 除以时间单位,6.58×10−16 eV s6.58\times10^{-16}\,\mathrm{eV\,s} 除以 4.93 μs4.93\,\mu\mathrm s。于是云的耦合 α=Mμ\alpha=M\mu 定下了粒子的静能:

mbc2=α×1.34×10−10 eV M⊙M,m_bc^2=\alpha\times1.34\times10^{-10}\,\mathrm{eV}\,\frac{\mathrm M_\odot}{M},

增长率 MωIM\omega_I,即 MωM\omega 的虚部,定下了场振幅增大到 ee 倍所需的时间:

τ=1ωI=GM/c3MωI.\tau=\frac{1}{\omega_I}=\frac{GM/c^3}{M\omega_I}.

云里的粒子数正比于振幅的平方,按 e2ωIte^{2\omega_It} 增长。

表 2.1 同一团云,两个黑洞:a=0.99a=0.99,α=0.421\alpha=0.421,MωI=1.50×10−7M\omega_I=1.50\times10^{-7},取自下一章的连分数解。
黑洞质量玻色子静能振幅e倍增长时间
10 M⊙10\,\mathrm M_\odot5.6×10−12 eV5.6\times10^{-12}\,\mathrm{eV}328 s328\,\mathrm s
4×106 M⊙4\times10^6\,\mathrm M_\odot1.4×10−17 eV1.4\times10^{-17}\,\mathrm{eV}4.2 yr4.2\,\mathrm{yr}

两行用的是同一个无量纲解。考察更重的黑洞,就是在考察更轻的粒子。

观测者能把这些数和什么联系起来

上面这些量都不能直接观测。黑洞自转要从吸积双星的X射线谱或并合引力波的波形里推断,各有各的模型和误差。如果云里的玻色子成对湮灭,会发出几乎单色的引力波,频率约为玻色子频率的两倍,fGW≈2mbc2/hf_{\mathrm{GW}}\approx2m_bc^2/h。mbc2=2.67×10−12 eVm_bc^2=2.67\times10^{-12}\,\mathrm{eV} 时约为 1.29 kHz1.29\,\mathrm{kHz},落在地面探测器的频段里。Baryakhtar、Galanis、Lasenby、Simon,《自相互作用标量场的黑洞超辐射》图13

要把增长率和这样的测量连起来,光靠本章的计算不够。下面的增长率只描述线性阶段的单个模式,它说明小云长得多快;云最后长到多大,取决于自转下降的历史。

第一次估算:类氢增长率

α\alpha 很小时,云离黑洞很远,可以把类氢波函数和视界附近的解拼接起来求增长率。l=m=1l=m=1 基态的结果是

MωI≈148(a−2αr+)α9.M\omega_I\approx\frac{1}{48}\left(a-2\alpha r_+\right)\alpha^9 .

Brito、Cardoso、Pani,《超辐射》第5.6.2节式(5.27)

括号其实就是超辐射条件。取 ω≈μ=α\omega\approx\mu=\alpha、ΩH=a/(2r+)\Omega_H=a/(2r_+),条件 ω<mΩH\omega<m\Omega_H 就写成 2αr+<a2\alpha r_+<a。陡峭的 α9\alpha^9 因子说明,又大又稀的云只有极少一部分能碰到视界。

例题2。 取 a=0.99a=0.99,r+=1.141r_+=1.141,α=0.1\alpha=0.1。括号等于 0.99−0.228=0.7620.99-0.228=0.762,于是 MωI≈0.762×10−9/48=1.59×10−11M\omega_I\approx0.762\times10^{-9}/48=1.59\times10^{-11}。严格的连分数解是 7.39×10−127.39\times10^{-12},小了2.1倍。对10倍太阳质量,估算给出的e倍增长时间是36天,严格解是77天。

表 2.2 类氢估算除以严格增长率。
α\alphaa=0.5a=0.5a=0.7a=0.7a=0.9a=0.9a=0.99a=0.99
0.050.051.151.151.291.291.441.441.531.53
0.100.101.081.081.411.411.861.862.152.15
0.200.200.790.792.232.233.543.54
0.300.301.071.074.674.67

空格处的耦合下,严格解里的云已经不再增长。在整张表里,估算的频率实部误差都小于0.4%,增长率却只准到几倍以内。截止点附近,两条曲线降向零时会碰巧相交,a=0.9a=0.9、α=0.3\alpha=0.3 时的1.07这种看似吻合就是这样来的。

本章每个数背后的假设

  1. 场是试验场:云足够轻,不改变克尔几何。
  2. 场是自由的,不计轴子那样的自相互作用。
  3. 只追踪一个模式,l=m=1l=m=1、径向无节点。
  4. 黑洞是孤立的,没有吸积,也没有伴星。
  5. 增长是线性的,自己永远停不下来。只有云拿走足够多的自转、让 ΩH\Omega_H 降到云的频率,增长才会停,而这段历史这里没有计算。

只作量级参考:α=0.42\alpha=0.42 时,装下10倍太阳质量黑洞百分之一质量的云,约含 2×10762\times10^{76} 个玻色子。从一个玻色子长到这么多,粒子数要经过 ln⁡(2×1076)≈176\ln(2\times10^{76})\approx176 个e倍,相当于88个振幅e倍增长时间,大约8小时。之后怎样,由饱和、自转下降和云自身的辐射决定。

Lab任务

  1. 在Lab里打开高级选项,把黑洞质量从10倍改成100倍太阳质量。先预言e倍增长时间和玻色子静能怎样变,再看读数。
  2. 选“恰在门槛上”预设,看放大图在 ω=mΩH\omega=m\Omega_H 附近的走向,确认曲线正好在那里穿过零,与能流恒等式的要求一致。
  3. 把自转设为0.99,比较小耦合和大耦合下的实线与虚线增长曲线,找出虚线约为实线两倍的耦合。

习题与解析

  1. 证明 K(r+)=2r+(ω−mΩH)K(r_+)=2r_+(\omega-m\Omega_H)。
  2. 为什么 ω\omega 为复数时,朗斯基论证不再成立?
  3. M=4×106 M⊙M=4\times10^6\,\mathrm M_\odot、α=0.2\alpha=0.2 时,玻色子静能是多少?
  4. a=0.99a=0.99、α=0.2\alpha=0.2 的云 MωI=1.61×10−9M\omega_I=1.61\times10^{-9}。对10倍太阳质量,它的e倍增长时间是多少?
  5. 一团云的 ω=0.38/M\omega=0.38/M。自转降到多少以下,它就不再增长?
解析
  1. 视界上 Δ=0\Delta=0,所以 r+2+a2=2r+r_+^2+a^2=2r_+。于是 K(r+)=2r+ω−am=2r+(ω−m a/(2r+))=2r+(ω−mΩH)K(r_+)=2r_+\omega-am=2r_+\left(\omega-m\,a/(2r_+)\right)=2r_+(\omega-m\Omega_H)。
  2. ω\omega 为复数时 K2K^2 是复数,R∗R^\ast 不再满足同一个方程,势能项也就抵消不了。增长或衰减的模式不是稳态,它的能流随时间变化。
  3. mbc2=0.2×1.34×10−10 eV/(4×106)=6.7×10−18 eVm_bc^2=0.2\times1.34\times10^{-10}\,\mathrm{eV}/(4\times10^6)=6.7\times10^{-18}\,\mathrm{eV}。
  4. τ=4.93×10−5 s/1.61×10−9=3.06×104 s\tau=4.93\times10^{-5}\,\mathrm s/1.61\times10^{-9}=3.06\times10^4\,\mathrm s,约8.5小时。
  5. m=1m=1 时,增长在 ΩH=a/(2r+)\Omega_H=a/(2r_+) 降到0.38时停止。代入 r+=1+1−a2r_+=1+\sqrt{1-a^2} 解出 a=4ω/(1+4ω2)=1.52/1.578=0.963a=4\omega/(1+4\omega^2)=1.52/1.578=0.963,这正是Lab在这个频率下显示的停止自旋。

类氢公式是向导,不是答案。下一章计算严格的增长率,并说明怎样确认它算对了:可复现地求出束缚云的增长率。

在Lab里试一试

打开高级选项改变黑洞质量,核对e倍增长时间和玻色子静能是否按你预言的比例变化。 打开黑洞超辐射Lab

来源

  1. Dolan,《克尔时空中有质量克莱因—戈尔登场的不稳定性》第II节
  2. Brito、Cardoso、Pani,《超辐射》第4.4节式(4.53)
  3. Baryakhtar、Galanis、Lasenby、Simon,《自相互作用标量场的黑洞超辐射》图13

继续阅读

这一页直接连着哪些内容

旋转黑洞与玻色云 · 阅读路径

从被放大的波到长大的云

  1. 文章一个旋转黑洞,为什么可能把周围的波越放越大?
  2. 文章波、自转与克尔黑洞的尺度
  3. 当前页面守恒能流、物理单位与第一次增长估算
  4. 文章可复现地求出束缚云的增长率
  5. 文章复算一份公开的超辐射放大表
  6. 文章自相互作用怎样改变玻色云的自旋限制与引力波信号
这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策