---
title: "偏心轨道的引力波为什么是一排谱线"
description: "算出 Kerr 偏心轨道的径向与方位频率，说明四极辐射为何是一排谱线，再用 Peters–Mathews 公式检验总功率。"
canonical_url: "https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/fundamental-frequencies-and-harmonics/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/extreme-mass-ratio-inspirals/fundamental-frequencies-and-harmonics/"
last_modified: "2026-11-12"
---

[极端质量比旋近，从轨道读到探测预测](https://sandboxphysics.com/zh/stories/extreme-mass-ratio-inspirals/)

# 偏心轨道的引力波为什么是一排谱线

Sandbox Physics · 2026-11-12

绕自转黑洞的偏心轨道由两个频率描述：$\Omega_r$ 和 $\Omega_\phi$。由轨道算出的任何量都按这两个频率重复，所以引力波里只能出现它们的整数组合。这一章先算出实验室默认轨道的两个频率，再把它的四极辐射拆成一条条谱线。

## 开始前需要什么

前两章、傅里叶级数，以及引力波功率的四极公式。

## 转折点与第三个根

赤道面 Kerr 测地线满足 $r^4(dr/d\tau)^2=R(r)$，$R=P^2-\Delta\,(r^2+x^2)$，其中 $P=E(r^2+a^2)-aL$，$\Delta=r^2-2r+a^2$，$x=L-aE$。$R$ 是 $r$ 的四次式，四个根为 $0$、$r_3$、$r_p$、$r_a$。比较三次因式的常数项，就得到第三个根。

Glampedakis 和 Kennefick 把 $E$ 和 $x^2$ 写成 $p$、$e$、$a$ 的函数。实验室默认轨道 $a=0.7,\ p=7,\ e=0.4$ 给出 $E=0.94461$，$L=3.10553$，$r_3=1.902$，稳稳落在近心点 $5$ 以内。

$$
r_3\,r_p\,r_a=\frac{2x^2}{1-E^2},\qquad r_p r_a=\frac{p^2}{1-e^2}\quad\Longrightarrow\quad r_3=\frac{2(1-e^2)\,x^2}{p^2\,(1-E^2)}
$$

稳定轨道要求第三个根在近心点以内；两者相遇的地方就是分界轨道。

对应阅读：Glampedakis & Kennefick (2002)

## 一个消掉转折点的参数

对 $r$ 积分，两个转折点处会碰到根号为零。换元 $r=p/(1+e\cos\chi)$ 可以消掉它们。可以算出 $(r_a-r)(r-r_p)=r^2e^2\sin^2\chi/(1-e^2)$，$dr/d\chi=r^2e\sin\chi/p$，因子 $\sin\chi$ 相消，下面的方程对所有 $\chi$ 都光滑。

这几个式子是我们从径向势自己推出来的，没有照抄文献里的积分式，并且与 Boyer–Lindquist 坐标下直接积分哈密顿方程的结果做了比较。默认轨道的两个径向周期相差不到百亿分之二。

$$
r^2\frac{d\chi}{d\tau}=\frac{p}{r}\sqrt{\frac{(1-E^2)\,r\,(r-r_3)}{1-e^2}},\qquad r^2\frac{dt}{d\tau}=a\,x+\frac{(r^2+a^2)\,P}{\Delta},\qquad r^2\frac{d\phi}{d\tau}=x+\frac{a\,P}{\Delta}
$$

单位 $\chi$ 的时间和方位角变化率就是这些式子的比值，$r^2$ 相消。

对应阅读：Glampedakis & Kennefick (2002)

## 两个基本频率

$\chi$ 转一整圈就是一个径向周期。把时间和方位角的变化率从 $0$ 积到 $2\pi$，得到周期 $T_r$ 和期间扫过的方位角 $\Delta\phi$。被积函数光滑且周期，梯形公式指数收敛：32个点就有十五位有效数字。

默认轨道：$T_r=226.95\,M$，$\Delta\phi=9.8259$ 弧度，即 $563.0^\circ$，近心点前移 $203^\circ$。$M=10^6\,M_\odot$ 时，$f_r=0.8946\,\mathrm{mHz}$，$f_\phi=1.3990\,\mathrm{mHz}$。

$$
T_r=\int_0^{2\pi}\frac{dt}{d\chi}\,d\chi,\qquad \Delta\phi=\int_0^{2\pi}\frac{d\phi}{d\chi}\,d\chi,\qquad \Omega_r=\frac{2\pi}{T_r},\qquad \Omega_\phi=\frac{\Delta\phi}{T_r}
$$

两者都是坐标时频率，也就是远处探测器在扣除宇宙学红移之前测到的频率。

对应阅读：Glampedakis & Kennefick (2002)

## 为什么波是一把梳子

把轨道放进平直空间，把 $r$ 和 $\phi$ 当作极坐标。它的质量四极矩有两个独立部分：$A=r^2e^{2i\phi}$ 和 $B=r^2$。简单计算可得，无迹四极矩满足 $\sum_{jk}\dddot Q_{jk}^2=\mu^2\bigl(|\dddot A|^2/2+\dddot B^2/6\bigr)$。

方位角等于匀速转动加上一个每个径向周期重复一次的部分：$\phi=\Omega_\phi t+\psi(t)$。于是 $A=e^{2i\Omega_\phi t}\,G(t)$，$G$ 以 $T_r$ 为周期。把 $G$ 展成傅里叶级数，$A$ 的谱线就落在 $2\Omega_\phi+n\Omega_r$。$B$ 只依赖 $r$，谱线在 $n\Omega_r$。

如果轨道不进动，$\Omega_\phi$ 等于 $\Omega_r$，每条线都会落在 $f_r$ 的整数倍上。实验室用灰色虚线画出这些整数倍；金色谱线偏离它们的距离，就是谱图里看到的进动。

对应阅读：Babak et al. (2017)

## 每条线的功率与两项检验

四极公式 $\mathcal F=\tfrac15\langle\dddot Q_{jk}\dddot Q_{jk}\rangle$ 于是拆成对谱线的求和。$A$ 的每条线贡献 $\tfrac1{10}|A_n|^2\omega^6$。$B$ 是实数，$+n$ 与 $-n$ 两条线功率相同，合起来是 $\tfrac1{15}|B_n|^2\omega^6$。

圆轨道检验：只剩 $A_0=r^2$，频率 $\omega=2\Omega_\phi$，求和正好得到 $\tfrac{32}{5}r^4\Omega_\phi^6$。弱场检验：$p=1000$、$e=0.4$ 时，谱线总和是牛顿椭圆 Peters–Mathews 轨道平均功率的1.0008倍，多出的0.08%来自那个距离上轨道的相对论修正。

默认轨道最强的四条线都来自转动部分，呼吸部分只占约0.01%的功率。偏心率越大，功率越往高 $n$ 走：$e=0.7$ 时只取十四条线最多会漏掉五分之一，所以实验室取四十条。

$$
\mathcal F=\frac{1}{10}\sum_{n}|A_n|^2\,\omega_{2,n}^6+\frac{1}{15}\sum_{n\ge1}|B_n|^2\,\omega_{0,n}^6,\qquad \omega_{2,n}=2\Omega_\phi+n\Omega_r,\quad \omega_{0,n}=n\Omega_r
$$

这是按小天体质量平方归一的能流，轨道保持不变。

| 谱线 $n$ | $M=10^6\,M_\odot$ 时的频率 | 占总功率 |
| --- | --- | --- |
| $1$ | $3.693\,\mathrm{mHz}$ | $23.1\,\%$ |
| $2$ | $4.587\,\mathrm{mHz}$ | $30.5\,\%$ |
| $3$ | $5.482\,\mathrm{mHz}$ | $21.5\,\%$ |
| $4$ | $6.376\,\mathrm{mHz}$ | $10.7\,\%$ |

## 这个模型在 EMRI 波形里的位置

Babak 等人介绍了两类低成本模型。数值近似波形沿 Kerr 测地线运动，再套用平直空间的辐射公式，这一章做的就是这件事。他们做预测用的解析近似波形从开普勒轨道出发，用后牛顿规则加上进动和旋近，并分别在史瓦西 ISCO 和 Kerr ISCO 处截断，用两个版本框住建模误差。

我们的谱线频率对测地线是精确的。谱线功率的准确度取决于四极公式；上一章与 Teukolsky 计算对比，靠近黑洞时最多差约50%。默认轨道的谱线总和是同样 $p$、$e$ 下 Peters–Mathews 值的0.834倍。这只是两种近似之比，说明不了哪一种更接近真实。

对应阅读：Babak et al. (2017) · Fujita & Tagoshi (2004)

## 在实验室里试一试

### [读出梳子和它的偏移](https://sandboxphysics.com/zh/labs/emri-spacetime-mapping/)

1. 默认轨道下，找出谱图里最高的金色谱线，以及离它最近的灰色虚线。
2. 把偏心率往零调，看梳子收拢成一条线。
3. 选“远处，接近开普勒”：进动变小，金色谱线离灰色虚线近得多。

谱线频率来自精确测地线；谱线功率用平直空间四极公式。

## 练习与答案

对 $r=p/(1+e\cos\chi)$，证明 $(r_a-r)(r-r_p)=r^2e^2\sin^2\chi/(1-e^2)$。

记 $u=1+e\cos\chi$。则 $r_a-r=p\,(u-1+e)/((1-e)u)$，$r-r_p=p\,(1+e-u)/((1+e)u)$。分子之积是 $e^2-(u-1)^2=e^2\sin^2\chi$，而 $p^2/u^2=r^2$。

默认谱图里有一条 $5.482\,\mathrm{mHz}$ 的线。它的 $n$ 是多少？离 $f_r$ 最近的整数倍有多远？

$2f_\phi+3f_r=2.798+2.684=5.482$，所以 $n=3$。没有进动时它应在 $5f_r=4.473\,\mathrm{mHz}$。进动把每条转动谱线往上推 $2(f_\phi-f_r)=1.009\,\mathrm{mHz}$，超过一格虚线间距，所以最近的灰色虚线是 $6f_r=5.368\,\mathrm{mHz}$，在这条线下方0.114 mHz。

在这个模型里，圆轨道为什么只辐射一条线？

$r$ 不变时 $G=r^2$ 不随时间变化，$A$ 里只剩 $n=0$；$B$ 是常数，不辐射。唯一的线在 $2\Omega_\phi$。这个模型略去的更高阶多极矩，会在其他 $m$ 的 $m\Omega_\phi$ 处加上更弱的线。

## 一手来源，以及每篇读了哪些部分

1. [Babak et al. (2017), Science with the space-based interferometer LISA. V. Extreme mass-ratio inspirals, Phys. Rev. D 95, 103012](https://arxiv.org/abs/1703.09722v2)

   arXiv 第2版，2018年4月11日。摘要；表I，十二个族群模型及其事件率；第四节 AKS 与 AKK 近似波形；图5，自转0.99黑洞上顺行圆赤道旋近的探测视界；图11至13，内禀参数、距离与天区以及四极矩偏离的 Fisher 矩阵误差分布；结论。

2. [Fujita & Tagoshi (2004), New numerical methods to evaluate homogeneous solutions of the Teukolsky equation, Prog. Theor. Phys. 112, 415](https://arxiv.org/abs/gr-qc/0410018v1)

   arXiv 第1版。表VII，七种自转下圆赤道轨道到 l = 6 的引力波光度；表III，r = 6M 处逐个模式的光度。单位归一化以及“到达无穷远的能流”这一解读是我们的推断，已用四极与后牛顿极限检验。

3. [Glampedakis & Kennefick (2002), Zoom and whirl: eccentric equatorial orbits around spinning black holes and their evolution under gravitational radiation reaction, Phys. Rev. D 66, 044002](https://arxiv.org/abs/gr-qc/0203086v1)

   arXiv 第1版。第二节式1至24：用 p、e 表示的赤道测地线常数，关于 L − aE 的四次方程，以及分界条件；附录A。周期的积分式没有沿用文中印刷形式，我们另行推导并做了数值检验。

## 继续阅读

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

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

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/)