作为 S³ 上测地流的 Kepler 问题

Jul 28, 2026 · updated Jul 29, 2026

翻译声明：本文由 LLM 从原文翻译而来，可能存在翻译不准确之处。建议阅读原文以获得最准确的内容。

Kepler 问题通常从它的守恒量入手。能量和角动量将运动约化到一个平面上，随后便可以在极坐标中算出轨道。这是求解该问题的一种有效方式，但它并没有穷尽其中隐藏的几何。

对于负能量，角动量与 Runge–Lenz 矢量揭示出一种 $SO (4)$ 对称性。我们不会在这里推导它们的完整代数关系，而只是将这一对称性视作一条线索： $SO (4)$ 是 $ℝ^{4}$ 的旋转群，也是标准三维球面 $𝕊^{3}$ 的保向等距群。这自然引出了一个问题：负能量 Kepler 问题能否被理解为 $𝕊^{3}$ 上的自由运动？

沿着Moser 正则化的思路，我们将直接构造这一对应。这个构造把正则化的 Kepler 流与 $𝕊^{3}$ 上的测地流等同起来，并让我们从大圆中重新得到 Kepler 椭圆。

三维球面的球极投影

将 $ℝ^{4}$ 中的一点写作 $(𝜉_{0}, 𝜉_{⟂}) \in ℝ \times ℝ^{3}$ ，并令

𝕊^{3} = {(𝜉_{0}, 𝜉_{⟂}) \in ℝ \times ℝ^{3} | 𝜉_{0}^{2} + {‖ 𝜉_{⟂} ‖}^{2} = 1} .

我们将北极记为 $𝑁 = (1, 0)$ 。从 $𝑁$ 出发的球极投影给出一个坐标卡

𝜎_{𝑁} : 𝕊^{3} - {𝑁} \to ℝ^{3}, 𝑥 = 𝜎_{𝑁} (𝜉_{0}, 𝜉_{⟂}) = \frac{𝜉_{⟂}}{1 - 𝜉_{0}} .

其逆映射为

𝜉_{0} = \frac{{‖ 𝑥 ‖}^{2} - 1}{1 + {‖ 𝑥 ‖}^{2}}, 𝜉_{⟂} = \frac{2 𝑥}{1 + {‖ 𝑥 ‖}^{2}} .

因此， $𝕊^{3}$ 是 $ℝ^{3}$ 的单点紧化：被略去的北极出现在 $‖ 𝑥 ‖ \to \infty$ 处。

标准球面度规拉回为 Euclidean 度规的一个共形倍数，

𝑔_{rd} = \frac{4}{{(1 + {‖ 𝑥 ‖}^{2})}^{2}} \sum_{𝑖 = 1}^{3} {(d 𝑥_{𝑖})}^{2},

而它的余度规为

𝑔_{rd}^{- 1} = \frac{{(1 + {‖ 𝑥 ‖}^{2})}^{2}}{4} \sum_{𝑖 = 1}^{3} 𝜕_{𝑥_{𝑖}}^{2} .

对于一个协向量 $𝑦 \in 𝑇_{𝑥}^{*} ℝ^{3}$ ，标准球面上自由粒子的 Hamilton 函数因而为

𝐹 (𝑥, 𝑦) ≔ \frac{1}{2} 𝑔_{rd}^{- 1} (𝑦, 𝑦) = \frac{{(1 + {‖ 𝑥 ‖}^{2})}^{2}}{8} {‖ 𝑦 ‖}^{2} .

单位速测地流位于能级集 $𝐹 = \frac{1}{2}$ 上。

仅有底空间上的投影还不足以给出动力系统之间的映射：Kepler 问题存在于相空间中。我们需要考虑球极投影的余切提升，其局部坐标为 $(𝑥, 𝑦)$ 。这个构造中出人意料的一点在于， $𝑥$ 对应 Kepler 动量，而 Kepler 位置则成为与它共轭的协向量 $𝑦$ 。

负能量 Kepler 能量面

约去质心运动后，三维 Kepler 问题的相空间为

𝒫︀ = 𝑇^{*} (ℝ^{3} - {0})

其中位置 $𝑞 \in ℝ^{3} - {0}$ 、动量 $𝑝 \in ℝ^{3}$ ，Hamilton 函数为

𝐻 (𝑞, 𝑝) = \frac{1}{2} {‖ 𝑝 ‖}^{2} - \frac{1}{‖ 𝑞 ‖} .

碰撞点集 $𝑞 = 0$ 被排除在外，而当一条轨道趋近碰撞时，Hamilton 向量场会变得奇异。

对能量进行缩放后，只需研究归一化能量面 $Σ_{- 1 / 2}$ 。它的能量方程为

- \frac{1}{2} = \frac{1}{2} {‖ 𝑝 ‖}^{2} - \frac{1}{‖ 𝑞 ‖},

或者等价地，

\frac{{(1 + {‖ 𝑝 ‖}^{2})}^{2}}{4} {‖ 𝑞 ‖}^{2} = 1 .

现在引入正则变换

Φ (𝑞, 𝑝) = (- 𝑝, 𝑞) ≕ (𝑥, 𝑦) .

归一化能量条件变为

\frac{{(1 + {‖ 𝑥 ‖}^{2})}^{2}}{4} {‖ 𝑦 ‖}^{2} = 1 .

将这个方程与球面 Hamilton 函数 $𝐹$ 比较，我们立即得到

Φ (Σ_{- 1 / 2}) = {𝐹 = \frac{1}{2}} = 𝑆^{*} (𝕊^{3} - {𝑁}) .

因此，在球极投影的余切坐标中，归一化的 Kepler 能量面恰好就是标准球面的单位速约束。

流的对应

能量超曲面相同，还不足以等同两侧带参数的流。令

𝐾 ≔ 𝐻 \circ Φ^{- 1} = \frac{1}{2} {‖ 𝑥 ‖}^{2} - \frac{1}{‖ 𝑦 ‖} .

由于 $Φ$ 是正则变换，Kepler 向量场的推前就是 $𝐾$ 的 Hamilton 向量场：

Φ_{*} 𝑋_{𝐻} = \sum_{𝑖 = 1}^{3} [\frac{𝑦_{𝑖}}{{‖ 𝑦 ‖}^{3}} 𝜕_{𝑥_{𝑖}} - 𝑥_{𝑖} 𝜕_{𝑦_{𝑖}}] .

对于球面上的自由粒子 Hamilton 函数，利用

d 𝐹 = - 𝜄_{𝑋_{𝐹}} 𝜔, 𝜔 = \sum_{𝑖 = 1}^{3} d 𝑦_{𝑖} \land d 𝑥_{𝑖},

我们得到

𝑋_{𝐹} = \sum_{𝑖 = 1}^{3} [\frac{{(1 + {‖ 𝑥 ‖}^{2})}^{2}}{4} 𝑦_{𝑖} 𝜕_{𝑥_{𝑖}} - \frac{1 + {‖ 𝑥 ‖}^{2}}{2} {‖ 𝑦 ‖}^{2} 𝑥_{𝑖} 𝜕_{𝑦_{𝑖}}] .

在共同的能量超曲面

Σ = {𝐹 = \frac{1}{2}} = {𝐾 = - \frac{1}{2}},

上，约束 $(1 + {‖ 𝑥 ‖}^{2}) ‖ 𝑦 ‖ = 2$ 将这个向量场化为

𝑋_{𝐹} |_{Σ} = \sum_{𝑖 = 1}^{3} [\frac{𝑦_{𝑖}}{{‖ 𝑦 ‖}^{2}} 𝜕_{𝑥_{𝑖}} - ‖ 𝑦 ‖ 𝑥_{𝑖} 𝜕_{𝑦_{𝑖}}] = ‖ 𝑦 ‖ Φ_{*} 𝑋_{𝐻} .

因此，

Φ_{*} 𝑋_{𝐻} = \frac{1}{‖ 𝑦 ‖} 𝑋_{𝐹} on Σ .

若以 $𝑡$ 表示 Kepler 时间、 $𝑠$ 表示球面测地流的弧长参数，那么

d 𝑡 = ‖ 𝑦 ‖ d 𝑠 = ‖ 𝑞 ‖ d 𝑠, d 𝑠 = \frac{d 𝑡}{‖ 𝑞 ‖} .

这就是 Sundman 时间重参数化。它将归一化的 Kepler 流转化为去点球面上的单位速测地流。

碰撞正则化

至此，这个构造给出

Σ_{- 1 / 2} ≅ 𝑆^{*} (𝕊^{3} - {𝑁}) .

在固定能量下，趋近碰撞意味着

𝑞 \to 0, ‖ 𝑝 ‖ \to \infty .

在 $𝑥 = - 𝑝$ 与 $𝑦 = 𝑞$ 之下，这变为 $‖ 𝑥 ‖ \to \infty$ ，恰好对应球极投影坐标卡中缺失的北极。

仅仅把北极加入底空间还不够。在碰撞处，剩余的数据是正则化轨道的方向。自然的补全方式是加入单位协向量球面

𝑆_{𝑁}^{*} 𝕊^{3} ≅ 𝕊^{2}

作为 $𝑁$ 上方的纤维。因此，补全后的能量面为

{\hat{Σ}}_{- 1 / 2} ≅ 𝑆^{*} 𝕊^{3} ≅ 𝑇_{1} 𝕊^{3} .

测地流在这个补全空间上是光滑的。穿过 $𝑁$ 的大圆给出正则化的径向碰撞—弹出轨道。在固定负能量下，它们是 $𝐿 = 0$ 的退化椭圆。

从大圆到 Kepler 轨道

现在可以直接从一条大圆读出 Kepler 轨道。令 $𝜓$ 为大圆的单位速角参数，并平移参数原点，使 $𝜓 = 0$ 对应近心点；再在轨道平面内选取正交单位向量 $𝑒_{1}, 𝑒_{2} \in ℝ^{3}$ 。应用余切球极投影并使用 $𝑞 = 𝑦$ ，得到

𝑞 (𝜓) = (cos 𝜓 - 𝑒) 𝑒_{1} + \sqrt{1 - 𝑒^{2}} sin 𝜓 𝑒_{2}, ‖ 𝑞 (𝜓) ‖ = 1 - 𝑒 cos 𝜓 .

参数 $𝑒$ 度量了大圆与北极的接近程度。消去 $𝜓$ ，得到

{(𝑞 \cdot 𝑒_{1} + 𝑒)}^{2} + \frac{{(𝑞 \cdot 𝑒_{2})}^{2}}{1 - 𝑒^{2}} = 1,

因此，大圆的余切投影是一条偏心率为 $𝑒$ 的 Kepler 椭圆。

由于大圆以单位速运动，我们有 $d 𝜓 = d 𝑠$ 。Sundman 关系于是变为

d 𝑡 = ‖ 𝑞 (𝜓) ‖ d 𝜓 = (1 - 𝑒 cos 𝜓) d 𝜓,

从而

𝑡 - 𝜏 = 𝜓 - 𝑒 sin 𝜓 .

因此，大圆的角参数、正则化时间与偏近点角是同一个变量，至多相差一个加法常数。

至此，我们从三维球面上的测地线出发，重新得到了 Kepler 问题的负能量分支，并完成了它的求解。

Back to all articles