---
title: "小黑洞绕自转黑洞的极端质量比旋近互动实验｜Sandbox Physics"
description: "绕上几十万圈的轨道，能读出大黑洞周围的什么？在大质量黑洞附近，轨道绕一圈回不到原处：最近点每绕一圈都会向前转，黑洞自转决定转多快。"
canonical_url: "https://sandboxphysics.com/zh/labs/emri-spacetime-mapping/"
language: "zh-CN"
page_type: "lab"
alternate_languages:
  en: "https://sandboxphysics.com/labs/emri-spacetime-mapping/"
last_modified: "2026-11-14"
---

# 小黑洞绕自转黑洞的极端质量比旋近

**绕上几十万圈的轨道，能读出大黑洞周围的什么？** 在大质量黑洞附近，轨道绕一圈回不到原处：最近点每绕一圈都会向前转，黑洞自转决定转多快。这些转动速率决定引力波的频率。小天体落入之前要绕几十万圈，累积起来的波相位有望测出黑洞的质量、自转和形状。

拖动小天体绕自转黑洞的轨道，看近心点进动、引力波谐波，以及到最内稳定轨道前还剩多少圈。

> 交互仿真：在浏览器中打开本页，即可调整参数并运行模型。

模型假设

单位取引力常数、光速和大黑洞质量为一，坐标为 Boyer–Lindquist 坐标。轨道是精确的束缚赤道面 Kerr 测地线。自转带正负号：取负值表示轨道逆着黑洞自转方向运行。轨道常数按 Glampedakis 与 Kennefick（2002）计算，径向和方位周期由径向相位积分得到，并与直接哈密顿积分一致。引力波谱用四极近似拼接：把测地线放进平直空间，再套用平直空间四极公式，它不是 Teukolsky 波形。对圆轨道，它的能流在六个引力半径处是 Fujita–Tagoshi Teukolsky 表值的一半到一点二倍，十个引力半径以外相差不到百分之十五。剩余波周数和剩余时间来自另一条圆轨道绝热旋近：从半通径出发，用同一四极能流和精确圆轨道能量，停在最内稳定圆轨道外侧。这里省略视界吸收、高阶多极、偏心率衰减、轨道倾角、共振，以及绝热阶以外的全部自力修正。它是参考曲线，不是 LISA 波形模板。画布把半径和方位角画成平直极坐标。

### 这条轨道的谐波

模拟，四极近似

$2\Omega_\phi+n\Omega_r$ $n\Omega_r$ 整数倍

### 到最内稳定轨道前的频率上升

圆轨道绝热参考

$f_{\rm GW}=\Omega_\phi/\pi$

这两张图怎样算出，可信到什么程度

左图：把轨道的质量四极矩放在平直空间里，分解到频率 $2\Omega_\phi+n\Omega_r$ 和 $n\Omega_r$ 上。灰线是径向频率的整数倍，也就是没有进动时谱线所在的位置。右图：从同一 $p$ 出发的圆轨道，按四极能流损失能量，一直到最内稳定轨道外侧。

四极近似的圆轨道能流，除以 Fujita 与 Tagoshi（2004）表中到达无穷远的 Teukolsky 能流（多极到六阶，不含视界能流）：

| $r/M$ | $a=-0.9$ | $a=0$ | $a=0.9$ |
| --- | --- | --- | --- |
| $6$ | $0.506$ | $0.879$ | $1.19$ |
| $10$ | $0.889$ | $1.04$ | $1.15$ |
| $20$ | $1.03$ | $1.07$ | $1.1$ |
| $50$ | $1.03$ | $1.04$ | $1.05$ |

靠近黑洞时，按它数出的波周数误差可达几十个百分点，远处要小得多。LISA 数据分析需要 Teukolsky 能流、偏心和倾斜轨道，以及一阶自力修正。

## 从轨道读出它绕着的黑洞

**背景** 恒星质量黑洞绕着百万倍太阳质量的黑洞运行，就是极端质量比旋近。它每一圈都几乎沿着大黑洞时空的测地线走，所以轨道把这片时空的形状带进了引力波。

**为什么值得理解** 行星轨道几乎不进动；靠近黑洞时，近心点每圈可以转过一整圈以上，黑洞自转还会改变转动的快慢。

- **核心问题**: 轨道和它发出的波，哪些特征取决于黑洞自转？
- **一句话直觉**: 束缚赤道面测地线有两个频率：$\Omega_r$ 描述一进一出的径向运动，$\Omega_\phi$ 描述绕转。每条引力波谱线都落在两者的组合上，而自转同时改变这两个频率。

### 核心数学模型

#### 圆轨道与最后一条稳定圆轨道

$$
\Omega_\phi=\frac{1}{r^{3/2}+a},\qquad \Omega_r^2=\Omega_\phi^2\left(1-\frac{6}{r}+\frac{8a}{r^{3/2}}-\frac{3a^2}{r^2}\right)
$$

单位取 G、c、M 为一。径向频率降到零的位置就是最内稳定圆轨道；顺着自转的轨道会让这个位置向内移。

#### 两个频率决定全部谐波

$$
\omega_{n}=2\Omega_\phi+n\,\Omega_r,\qquad \Delta\phi_{\rm peri}=2\pi\left(\frac{\Omega_\phi}{\Omega_r}-1\right)
$$

转动的四极矩在这些频率上辐射。没有进动时两个频率相等，谱线正好落在径向频率的整数倍上。

#### 能量损失推动缓慢旋近

$$
\frac{dr}{dt}=-\frac{\mathcal F}{dE/dr},\qquad N_{\rm GW}=\int\frac{\Omega_\phi}{\pi}\,dt\propto\frac{1}{q}
$$

能流与单位质量能量都由轨道决定，所以波周数与质量比成反比。这是绝热阶近似；一阶自力修正还会贡献约一弧度量级的相位，对测量来说仍不可忽略。

### 常见难点

#### 这里的进动并不慢

**容易误解** 近心点进动总是像水星那样的小修正。

**正确理解** 靠近分界轨道时，每圈进动可以超过一整圈，形成急冲回旋轨道。对比“急冲与回旋”和“远处，接近开普勒”两个预设。

#### 谱线不是轨道频率的整数倍

**容易误解** 偏心轨道只在一个轨道频率的整数倍上辐射。

**正确理解** 有进动时，谱线相对灰色整数倍刻度整体偏移，偏移量是两个频率之差的两倍。这个偏移本身就是对时空的测量。

#### 波周数属于圆轨道参考

**容易误解** 剩余波周数跟着屏幕上的偏心轨道走。

**正确理解** 它是从同一半通径出发的圆轨道参考，用四极能流计算，没有计入偏心率衰减、视界吸收和自力修正。

### 动手实验

1. 01

   #### 先看玫瑰线

   保持默认设置，看橙色近心点一圈圈向前移。

   **你会看到：** 每个径向周期结束时，近心点都更靠前；读数给出进动角度。

2. 02

   #### 把自转反过来

   保持默认轨道，把自转滑块从 0.7 拖到零，再继续往下拖。

   **你会看到：** 每圈进动变大，粉色的最内稳定轨道圈向外移。自转一变成负值，这条轨道就不再稳定，画面改为显示最近的稳定轨道。

3. 03

   #### 找到急冲回旋的边界

   把橙色近心点手柄向内拖，直到出现提示。

   **你会看到：** 靠近分界轨道时，小天体会在近心点附近多绕几圈，再冲回远处。

4. 04

   #### 数一数剩余波周数

   把质量比改变十倍。

   **你会看到：** 剩余波周数和剩余时间同样变化十倍，轨道和谱线不变。

物理结论

轨道的两个频率，以及它们在旋近中的漂移，携带了黑洞的质量和自转。探测器只要跟住几十万个周期，就能把这种漂移变成测量。

## 相关实验与文章

极端质量比旋近 · 阅读路径

### 小黑洞为大黑洞测绘

这条阅读路径里用到本实验的文章，后面跟着它们推荐的其他实验。

1. [文章 **绕上几十万圈的轨道，能读出大黑洞的什么**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/what-hundreds-of-thousands-of-orbits-reveal/)
2. [文章 **自转黑洞周围的轨道怎样进动，又在哪里失去稳定**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/geodesics-precession-and-the-last-stable-orbit/)
3. [文章 **能量流失怎样决定最后稳定轨道前还剩多少圈**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/radiation-reaction-and-the-adiabatic-inspiral/)
4. [文章 **偏心轨道的引力波为什么是一排谱线**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/fundamental-frequencies-and-harmonics/)
5. [文章 **用已发表的能流表复算一条 EMRI 频率演化**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/recompute-a-frequency-evolution/)
6. [文章 **怎样读 EMRI 的探测预测、自力、共振与多极矩检验文献**](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/reading-emri-forecasts-self-force-resonances-multipoles/)

运动与动力系统 · 入门路径

## 从一次下落到一条轨道

先测重力，再分解运动，加入接触力与碰撞，最后把同样的定律带到轨道上。

[从这里开始 → **自由落体与重力测量** 采集几次下落，用相机记录的位置拟合重力加速度；再加入空气阻尼，看拟合结果怎样偏离真值。 观察 · 14 分钟](https://sandboxphysics.com/zh/labs/free-fall/)

1. [1 自由落体与重力测量](https://sandboxphysics.com/zh/labs/free-fall/)
2. [2 抛体飞行设计器](https://sandboxphysics.com/zh/labs/projectile-flight/)
3. [3 斜面与摩擦](https://sandboxphysics.com/zh/labs/inclined-plane/)
4. [4 一维碰撞台](https://sandboxphysics.com/zh/labs/collision-bench/)
5. [5 弹簧振子与周期标定](https://sandboxphysics.com/zh/labs/mass-spring/)
6. [6 轨道是怎样形成的](https://sandboxphysics.com/zh/labs/orbit/)

[查看完整路径 →](https://sandboxphysics.com/zh/topics/mechanics/#starter-path)