跳到主要内容
Sandbox Physics
从远处力矩走到 Kozai–Lidov 轨道循环

AC09 · L2

从四极势推导偏心率循环

Sandbox Physics ·

把伴星的差分引力势沿内外轨道平均，就得到最简单的四极循环。守恒的平均能量和轴向角动量共同决定偏心率能达到多大；向量方程则直接演化椭圆的形状与方向，避免轨道角度在圆轨道处的奇点。

先减去恒星与行星共同承受的加速度

用外轨向量表示伴星相对于中心恒星的位置，内轨向量表示行星的位置。中心恒星本身也在向伴星加速。在以它为原点、坐标轴不旋转的参考系里，要加入间接势，扣除这份共同加速度。

内轨距离远小于外轨距离时，可以展开伴星势。一次项恰好与间接势抵消，常数项不产生力，最先留下的是关于内轨位置的二次项，也就是四极潮汐势。

再沿圆外轨做平均。伴星方向在一个平面内均匀扫过，其投影平方的平均值，等于内轨位置在外轨平面内分量平方的一半。所得势函数绕外轨法向量轴对称。此时只平均了外轨时钟，内轨仍保留瞬时位置。

U=−Gm2∣R−r∣+Gm2R⋅rR3U=-\frac{Gm_2}{|\boldsymbol R-\boldsymbol r|}+\frac{Gm_2\boldsymbol R\cdot\boldsymbol r}{R^3}
第二项产生的加速度扣除中心恒星朝伴星的加速。
UQ=−Gm22R3[3(r⋅R^)2−r2],⟨UQ⟩2=−Gm24a23[r2−3(r⋅n^2)2]U_Q=-\frac{Gm_2}{2R^3}\left[3(\boldsymbol r\cdot\hat{\boldsymbol R})^2-r^2\right],\qquad \langle U_Q\rangle_2=-\frac{Gm_2}{4a_2^3}\left[r^2-3(\boldsymbol r\cdot\hat{\boldsymbol n}_2)^2\right]
前式保留外轨相位，后式把圆外轨的快速运动平均掉。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1 · Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

沿椭圆积分时，要使用正确的时钟

在内轨平面选两个单位向量，一个指向近心点，另一个与它垂直。偏近点角可以用这两个向量参数化位置。不过，等间隔的偏近点角并不对应等长时间，平均测度中必须保留下面的时间权重。

以距离平方为例，把参数化位置代入后，奇数余弦项积分为零，余弦平方的平均为二分之一，最后会留下偏心率平方修正。对沿外轨轴的投影平方做同样计算，便得到第二个恒等式。

把两个积分代回外轨平均势，再除以 L1 定义的势尺度，整理后就得到无量纲哈密顿量。偏心率向量指向近心点，长度等于偏心率，因此还保存了椭圆朝向；仅用一个偏心率数值无法完整描述这种运动。

r=a1[(cos⁡E−e)p^+1−e2sin⁡E q^],dM=(1−ecos⁡E) dE\boldsymbol r=a_1\left[(\cos E-e)\hat{\boldsymbol p}+\sqrt{1-e^2}\sin E\,\hat{\boldsymbol q}\right],\quad dM=(1-e\cos E)\, dE
时间平均要沿平近点角积分，再除以一整圈。
⟨r2⟩1=a12(1+32e2),⟨(r⋅n^2)2⟩1=a122[1−e2−j∥2+5e∥2]\langle r^2\rangle_1=a_1^2\left(1+\frac32e^2\right),\qquad \left\langle(\boldsymbol r\cdot\hat{\boldsymbol n}_2)^2\right\rangle_1=\frac{a_1^2}{2}\left[1-e^2-j_\parallel^2+5e_\parallel^2\right]
平行偏心率表示偏心率向量沿外轨法向量的投影。
ϕQ=⟨UQ⟩1,2Φ0=18[1−6e2−3j∥2+15e∥2]\phi_Q=\frac{\langle U_Q\rangle_{1,2}}{\Phi_0}=\frac18\left[1-6e^2-3j_\parallel^2+15e_\parallel^2\right]
这里采用单位测试粒子质量的势能符号约定，文献若使用负扰动函数，公式符号会相反。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1

用一个哈密顿量得到两条向量方程

角动量向量与偏心率向量描述同一条椭圆，不能任意独立取值：它们互相垂直，长度平方之和为一。轨道泊松括号给出下面的叉积演化规则；不能把这些受约束向量直接当成普通笛卡尔位置和动量。

对多项式哈密顿量求向量梯度，再代入叉积规则，并以世俗时间归一化，就得到后面两条显式方程。Lab 积分的正是这组方程。

把方程分别与原向量点乘，可以检验约束：向量与包含自己的叉积正交。同样的计算还表明轴向角动量不变，而不显含时间的哈密顿量保持其数值。数值程序保留原始积分向量，用这些关系检查误差，并不强行投影来掩盖漂移。

djdτ=−j×∇jϕ−e×∇eϕ,dedτ=−j×∇eϕ−e×∇jϕ\frac{d\boldsymbol j}{d\tau}=-\boldsymbol j\times\nabla_j\phi-\boldsymbol e\times\nabla_e\phi,\qquad \frac{d\boldsymbol e}{d\tau}=-\boldsymbol j\times\nabla_e\phi-\boldsymbol e\times\nabla_j\phi
轨道向量括号与势函数的符号约定共同固定正负号。
j⋅e=0,j2+e2=1,j⋅n^2=constant,ϕQ=constant\boldsymbol j\cdot\boldsymbol e=0,\quad j^2+e^2=1,\quad \boldsymbol j\cdot\hat{\boldsymbol n}_2=\mathrm{constant},\quad \phi_Q=\mathrm{constant}
这些量分别检验向量几何、轴对称性和平均能量。
j′=34[j∥j×n^2−5e∥e×n^2]\boldsymbol j^{\prime}=\frac34\left[j_\parallel\boldsymbol j\times\hat{\boldsymbol n}_2-5e_\parallel\boldsymbol e\times\hat{\boldsymbol n}_2\right]
撇号表示对无量纲世俗时间求导。
e′=34[2j×e−5e∥j×n^2+j∥e×n^2]\boldsymbol e^{\prime}=\frac34\left[2\boldsymbol j\times\boldsymbol e-5e_\parallel\boldsymbol j\times\hat{\boldsymbol n}_2+j_\parallel\boldsymbol e\times\hat{\boldsymbol n}_2\right]
偏心率向量方向的变化同时描述近心点转动。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1 · Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

找出临界倾角与最大偏心率

考虑初始轨道无限接近圆的极限。轴向角动量趋近于初始互倾角的余弦。在增长支上，偏心率达到最大时，近心点方向与交线相差九十度，此时偏心率向量的轴向分量等于偏心率乘以互倾角正弦。

用守恒的轴向角动量消去倾角，再令峰值时和初始时的哈密顿量相等，就得到下面的因式分解。顺行轨道要超过约三十九点二三度，才会出现非零增长支。

初始倾角为六十度时，最大偏心率约为零点七六三七六三，近心距约为半长轴的零点二三六二三七倍。这是初始近圆极限的解析例子。默认初始偏心率零点零五拥有略微不同的平均能量，因此在六十五度时达到约零点八三八九零九，而零种子极限给出零点八三八零四七。

严格为圆的向量状态在这组方程里会一直保持圆。所谓不稳定，是任意小的扰动能够增长，而增长所需时间也依赖种子大小。程序若偷偷加入偏心率，就改变了初值问题；本页明确显示并保存所用的有限初始偏心率。

0=emax⁡2[9−15cos⁡2i01−emax⁡2],emax⁡=1−53cos⁡2i00=e_{\max}^2\left[9-\frac{15\cos^2 i_0}{1-e_{\max}^2}\right],\qquad e_{\max}=\sqrt{1-\frac53\cos^2i_0}
此式采用近圆初值、测试粒子及圆外轨四极近似。
icrit=arccos⁡35=39.2315∘i_{\mathrm{crit}}=\arccos\sqrt{\frac35}=39.2315^\circ
在这组近似下，超过临界互倾角才会打开增长支。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1

发生共振的是哪个慢角度

近心点幅角描述偏心率向量在内轨平面里相对于升交点的方向。升交点是行星穿过外轨平面、朝其法向量正侧运动的位置。采用这个约定，偏心率向量的轴向分量等于偏心率、互倾角正弦和近心点幅角正弦的乘积。

代回哈密顿量，就得到偏心率与幅角的表达式。固定能量的等值线规定慢变轨迹：有些幅角会循环经过全部方向，有些则绕九十或二百七十度往返摆动。后一种叫作振荡，角度围绕某个偏好值变化，而不是不断转过整圈，这就是通常所说的 Kozai–Lidov 共振行为。

仅凭偏心率大幅循环，还不能区分两者。这里的有限种子与初始近心点方向给出的默认平均轨迹属于幅角循环：十二个世俗时间内，展开后的幅角累计前进约四百零八度。导出的偏心率和角动量向量保留了重建幅角的信息。轨道接近圆时，角度容易失去定义，向量描述仍然可用。

ϕQ=18[1−6e2−3j∥2+15e2(1−j∥21−e2)sin⁡2ω]\phi_Q=\frac18\left[1-6e^2-3j_\parallel^2+15e^2\left(1-\frac{j_\parallel^2}{1-e^2}\right)\sin^2\omega\right]
近心点幅角在内轨平面从升交点起算；守恒轴向投影用于消去互倾角。

对应阅读：Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1 · Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

加入保守的相对论近心点进动

加入一个与角动量长度成反比的负势，就能产生最低阶轨道平均相对论进动。它不施加角动量力矩，而是在轨道平面内转动偏心率向量。高偏心率时，这种转动会加快，破坏此双平均模型里牛顿力矩的持续同相作用。

把默认轨道的无量纲相对论强度设为零点八，平均最大偏心率降至约零点六二九二七八。测试还独立求解了修正后的能量根，结果与向量积分一致。保存的直接轨迹仍采用牛顿引力；琥珀色曲线明确作为牛顿参照，不能读成相对论直接积分。

进动公式要求弱引力场和低速运动，近心点必须远离引力半径。计算会在声明的近星边界前，或引力半径达到近心距百分之一前停止。后一项是保守选定的适用边界，不代表物理突变。耗散潮汐和引力辐射均未加入。放宽双平均后，快变化与进动可以发生更复杂的相互作用，最后的原文导读会讨论这个范围。

ϕ=ϕQ−ϵGRj,e′∣GR=ϵGRj3j×e,ω˙GR=ϵGRtKj2\phi=\phi_Q-\frac{\epsilon_{\mathrm{GR}}}{j},\quad \left.\boldsymbol e^{\prime}\right|_{\mathrm{GR}}=\frac{\epsilon_{\mathrm{GR}}}{j^3}\boldsymbol j\times\boldsymbol e,\quad \dot\omega_{\mathrm{GR}}=\frac{\epsilon_{\mathrm{GR}}}{t_Kj^2}
系数无量纲，它与物理光速的关系由下式固定。
ϵGR=3GMa1c2Mm2(a2a1)3,GMc2q≪1\epsilon_{\mathrm{GR}}=\frac{3GM}{a_1c^2}\frac{M}{m_2}\left(\frac{a_2}{a_1}\right)^3,\qquad \frac{GM}{c^2q}\ll1
质量和长度固定时，独立改变教学用进动系数，等于改变所对应的相对论紧致度。

对应阅读：Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2 · Klein & Hamilton, Short-range forces beyond double averaging, arXiv:2606.19456v1

保存的数值例子与计算方法

合成动力学单位令引力常数、两颗恒星质量和初始内轨半长轴均为一。平均原始向量采用 RK4，最大世俗步长为零点零零一；直接参照为固定 REBOUND 5.2.2 IAS15 状态，保存间隔为零点零零五个世俗时间。

有限初始偏心率与初始倾角e0=0.0500,i0=65.0∘e_0=0.0500,\quad i_0=65.0^\circ
默认平均最大偏心率emax⁡=0.838909e_{\max}=0.838909
初始近圆解析极限的最大值emax⁡(e0→0)=0.838047e_{\max}(e_0\to0)=0.838047
保守平均相对论项下的最大值ϵGR=0.800,emax⁡=0.629278\epsilon_{\mathrm{GR}}=0.800,\quad e_{\max}=0.629278
代表状态的平均哈密顿量ϕQ=0.056315120\phi_Q=0.056315120

代表状态取世俗时刻一点八。有限初始偏心率与近圆解析极限分别列明；可下载分析能复算峰值，并提供实际直接积分状态。

亲手改变轨道交换

  1. 比较偏心率交换与相对论进动预设，观察平均峰值下降，同时留意直接参照仍使用声明的牛顿模型。
  2. 将初始倾角提高到近极轨预设，观察记录在排除边界停止。一次未完成的高偏心率跃升没有完整循环周期。

练习与解释性答案

  1. 为什么均匀采样偏近点角，会得到错误的平均潮汐势？

    展开解释

    它把快速经过的近心点附近赋予过多权重，却低估了行星在远心点附近停留的时间。平近点角测度补回时间权重，偏心率平方项的系数也随之改变。

  2. 把三十度初始倾角代入近圆最大偏心率公式，根号下为负意味着什么？

    展开解释

    说明非零增长支在此处不存在，并不是物理偏心率变成了虚数。这组近似下，轨道处于稳定的低偏心率分支。

  3. 某次积分的偏心率始终有界，但两个轨道向量逐渐失去正交性，可以相信它的最大偏心率吗？

    展开解释

    还不能。这两个向量已经无法对应同一条物理开普勒椭圆。应同时检查两条几何约束、轴向投影和哈密顿量，再缩小步长或修正积分器，恢复它们的精度。只截断一个标量，两个向量仍然不能构成正确的椭圆。

来源与计算范围

  1. Tremaine & Yavetz, Why do Earth satellites stay up?, arXiv:1309.5244v1

    固定为2013年9月20日版本，DOI10.1119/1.4874853。势函数18、向量方程23、条件24和第IV节给出四极模型及近圆轨道临界角极限。

  2. Liu, Muñoz & Lai, Short-range forces in triple systems, arXiv:1409.6717v2

    固定为2015年1月7日版本，DOI10.1093/mnras/stu2396。阅读第2–4节、四极势7、八极势8、GR项32–33/42和图3–4、7。浏览器采用圆外轨、测试粒子四极模型，导读继续阅读更完整的模型。

  3. Klein & Hamilton, Short-range forces beyond double averaging, arXiv:2606.19456v1

    固定阅读2026年6月17日预印本v1，另有八月v2。图1区分平均和物理选择，图2–3展示周期变化，附录B与直接N体例子比较。本页导读原文，浏览器没有复现其中单平均GR及潮汐计算。

反馈这个页面有错误、看不明白，或者缺了什么？告诉我们
反馈类型

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