双摆动冥族小天体的轨道动力学研究进展与科学意义——基于 Malhotra 与 Ito（2024）研究的文献综述

摘要

柯伊伯带是太阳系形成与演化过程中留下的原始天体宝库。作为其中与海王星存在 3:2 平运动共振的动力学族群，冥族小天体（Plutinos）的轨道特征为揭示巨行星的迁移历史提供了关键证据。近日点幅角的额外摆动（双摆动现象）是冥王星轨道的显著特征，也是部分冥族小天体共有的动力学属性。本文基于 Malhotra 与 Ito（2024）发表于《Astronomical Journal》的研究《The doubly librating Plutinos》，系统梳理了冥族小天体的研究背景、双摆动现象的动力学基础、识别方法与数据来源，深入分析了双摆动冥族小天体的轨道参数分布特征及其动力学机制，探讨了该研究结果对巨行星迁移理论与柯伊伯带演化模型的约束作用，并指出了当前研究存在的局限性与未来的研究方向。研究表明，约 16% 的长期稳定冥族小天体表现出持续的近日点幅角 $g$ 摆动，这些天体在时间平均的偏心率-倾角平面上沿一条近似双曲线弧聚集，该结构对应于圆限制六体问题中的庞加莱第三类周期轨道族。双摆动冥族小天体的观测丰度（约 16%）显著高于基于频率比估算的随机分布预期（约 0.5%），这一“丰度异常”可由长期动力学稳定性效应与巨行星迁移过程中的优先捕获效应的联合作用得到合理解释。

关键词：冥族小天体；双摆动现象；von Zeipel-Lidov-Kozai 振荡；庞加莱第三类周期轨道；巨行星迁移；柯伊伯带演化

一、引言

太阳系的形成与演化是天体力学与行星科学的核心研究课题之一。自 1992 年柯伊伯带首个天体 1992 QB1 被发现以来，柯伊伯带作为太阳系边缘的原始天体聚集区，其动力学结构与物质组成成为研究太阳系早期演化的重要窗口。柯伊伯带天体根据其轨道特征可分为多个动力学族群，其中冥族小天体是最为突出的一类——它们与海王星存在 3:2 平运动共振，即每绕太阳公转 3 圈，海王星恰好公转 2 圈。冥族小天体的存在及其高丰度被认为是太阳系早期巨行星发生轨道迁移的最强有力证据之一（Nesvorný, 2018; Malhotra, 2019; Gladman & Volk, 2021）。

冥王星作为首个被发现的冥族小天体，其轨道具有两个显著的动力学特征：一是临界共振角 $\sigma =3\lambda -2{\lambda }^{\prime }-\varpi$ 围绕 180° 摆动，这是 3:2 平运动共振的标志；二是近日点幅角 $g$（定义为升交点黄经 $\Omega$ 与近日点黄经 $\varpi$ 之差，即 $g=\varpi -\Omega$）围绕 $\pm 90^{\circ }$ 摆动，这一现象被称为 von Zeipel-Lidov-Kozai（vZLK）振荡。长期以来，天文学家一直关注其他冥族小天体是否也存在类似的双摆动现象，以及这种现象的动力学起源与演化意义。早期的研究由于样本量较小（通常不足 100 个），对双摆动冥族小天体的丰度估计存在较大差异（从 8% 到 33% 不等），且对其轨道参数分布的系统性分析较为缺乏。这促使 Malhotra 与 Ito（2024）利用迄今最大的冥族小天体样本，系统性地开展此项研究。

2024 年，Malhotra 与 Ito 在《Astronomical Journal》上发表了题为《The doubly librating Plutinos》的研究论文。该研究基于 441 个长期稳定的冥族小天体（目前最大的样本），通过长达 1 亿年的数值轨道积分，系统识别了双摆动冥族小天体，并首次发现它们在时间平均的偏心率-倾角平面上沿近似双曲线弧聚集的动力学结构。该研究将这一结构与庞加莱在 19 世纪末发现的第三类周期轨道族联系起来，并对双摆动冥族小天体的丰度异常及其科学意义进行了深入探讨。本文将以此为核心，全面综述双摆动冥族小天体的研究进展，分析其动力学机制与演化意义，为后续相关研究提供参考。

二、冥族小天体的研究背景与动力学基础

2.1 柯伊伯带的动力学分类与冥族小天体的定义

柯伊伯带位于海王星轨道之外，距离太阳约 30 至 50 天文单位（AU），是太阳系中最大的小天体聚集区之一。根据轨道半长轴与动力学特征，柯伊伯带天体可分为经典柯伊伯带天体、共振柯伊伯带天体、离散盘天体与塞德娜型天体四大类。其中，共振柯伊伯带天体是指与海王星存在平运动共振的天体，它们的轨道周期与海王星的轨道周期成简单整数比。

冥族小天体是共振柯伊伯带天体中数量最多、研究最为深入的一类，其定义为与海王星存在 3:2 平运动共振的柯伊伯带天体。这一共振关系意味着冥族小天体的轨道半长轴约为 39.4 AU，与冥王星的轨道半长轴（约 39.5 AU）相近。冥族小天体的动力学特征由临界共振角 $\sigma$ 的摆动所定义，其表达式为：
$$\sigma =3\lambda -2{\lambda }^{\prime }-\varpi$$ 其中，$\lambda$ 为冥族小天体的平黄经，${\lambda }^{\prime }$ 为海王星的平黄经，$\varpi$ 为冥族小天体的近日点黄经。对于所有稳定的冥族小天体，$\sigma$ 围绕 180° 摆动，摆动周期约为 $10^4$ 年量级。这一共振机制能够保护冥族小天体免受海王星的近距离散射，使其在太阳系演化的漫长历史中保持轨道稳定。

这张图展示了平均运动共振的相空间结构。图中闭合的"8字形"轨道代表了共振角（critical argument）的周期性摆动（libration）——即冥族小天体在2:3共振中的"双摆动"现象。当共振角在平衡位置附近周期性振荡时（而非持续单向旋转），就形成稳定的共振俘获状态。

2.2 巨行星迁移理论与冥族小天体的起源

冥族小天体的存在及其高丰度无法用太阳系的静态形成模型解释，这一观测事实直接推动了巨行星迁移理论的发展。根据尼斯模型（Nice Model），太阳系早期的巨行星轨道与现今不同：木星位于约 5.5 AU 处，土星位于约 8 AU 处，天王星与海王星则位于更靠近太阳的位置。随着原行星盘气体的消散，巨行星之间的引力相互作用导致其轨道发生迁移——木星向内迁移，而土星、天王星与海王星向外迁移。当海王星向外迁移时，其平运动共振带也随之向外移动，将原本位于共振带内的原始柯伊伯带天体捕获进入共振轨道，形成了现今的冥族小天体族群（Nesvorný, 2018）。

巨行星迁移理论得到了大量观测与数值模拟的支持，而冥族小天体的轨道特征则为约束迁移的速度、时间尺度与迁移模式提供了关键参数。例如，冥族小天体的偏心率与倾角分布能够反映海王星迁移的平滑程度：如果海王星迁移速度较慢，共振捕获效率较高，且能够保留天体的原始偏心率与倾角；反之，如果迁移速度较快，只有偏心率较高的天体能够被捕获进入共振轨道。

2.3 双摆动现象与 vZLK 振荡

除了临界共振角 $\sigma$ 的摆动外，部分冥族小天体还表现出近日点幅角 $g$ 的额外摆动，这种同时存在两个角度摆动的现象被称为双摆动现象。近日点幅角 $g$ 定义为升交点黄经 $\Omega$ 与近日点黄经 $\varpi$ 之差，即 $g=\varpi -\Omega$，它描述了近日点相对于轨道升交点的位置。

近日点幅角的摆动本质上是 vZLK 振荡的表现形式。vZLK 振荡是限制性三体问题中的一种经典动力学现象，由 von Zeipel、Lidov 与 Kozai 分别于 20 世纪初与 20 世纪 60 年代独立发现。在圆限制三体问题中，当一个小质量天体绕中心天体（太阳）公转，并受到一个遥远的大质量摄动体（如海王星）的引力作用时，如果小天体的轨道相对于摄动体轨道平面具有足够大的倾角（通常大于约 39.2°），那么其轨道半长轴 $a$ 将保持恒定，而一个称为 vZLK 常数的量 $C=(1-e^2){\cos }^{2}i$ 近似守恒。在此情况下，小天体的偏心率与倾角将发生周期性的相互转换，同时近日点幅角 $g$ 将围绕 $0^{\circ }$ 或 $180^{\circ }$（对于顺行轨道）或 $90^{\circ }$ 与 $270^{\circ }$（对于逆行轨道）摆动。

然而，对于处于平运动共振中的冥族小天体，经典的 vZLK 理论不再完全适用。这是因为平运动共振的存在打破了轨道半长轴的守恒性，且除了海王星外，木星、土星与天王星的引力摄动也会对冥族小天体的轨道产生显著影响。因此，共振条件下的 vZLK 振荡需要更复杂的轨道平均方法进行分析（Ito & Malhotra, 2023; Lei et al., 2022; Saillenfest, 2020）。研究表明，在 3:2 平运动共振中，近日点幅角 $g$ 仍然可以围绕 $0^{\circ }$、$90^{\circ }$、$180^{\circ }$ 与 $270^{\circ }$ 四个中心摆动，但摆动的稳定性与范围会受到共振条件的调制。

来自IOPscience期刊论文 “LiDO: Exploring the Stable Plutino Parameter Space” 的精准数据图，展示了冥族小天体（Plutino）的φ₃:₂共振角摆动幅度分布。图中横坐标为摆动幅度（度），显示了不同冥族小天体在2:3共振中的稳定摆动范围（通常集中在60°-100°幅度）。

早期对双摆动冥族小天体的研究主要集中在少数几个天体上，样本量较小。例如，Nesvorný 等人（2000）基于 6 个冥族小天体的样本，估计双摆动天体的比例约为 33%；Volk 等人（2016）基于 24 个天体的样本，得到的比例约为 17%；而 Balaji 等人（2023）基于 85 个天体的样本，估计比例约为 8%。这些差异主要源于样本量的大小、积分时间的长短以及双摆动天体识别标准的不同。此外，观测偏差也会对丰度估计产生显著影响，因为高倾角的冥族小天体更容易被观测遗漏，而双摆动天体往往具有较高的倾角。

三、研究方法与数据来源

3.1 数据来源

Malhotra 与 Ito（2024）的研究采用了 Volk 与 Van Laerhoven（2024）发布的最新冥族小天体列表，该列表包含了截至 2024 年初已被确认的所有冥族小天体。为了获取这些天体的精确轨道参数，研究人员从美国国家航空航天局（NASA）的喷气推进实验室（JPL）Horizons 系统下载了所有天体在初始历元的吻切轨道元素，包括半长轴 $a$、偏心率 $e$、轨道倾角 $i$、升交点黄经 $\Omega$、近日点黄经 $\varpi$ 与平黄经 $\lambda$。同时，研究人员还从 JPL Horizons 系统获取了太阳与八大行星的质量以及行星在初始历元的轨道参数，作为数值模拟的输入数据。

3.2 数值轨道积分方法

为了研究冥族小天体的长期轨道动力学，研究人员采用了 SWIFT_RMVS3 正则化混合变量辛积分器进行数值轨道积分。SWIFT_RMVS3 是由 Levison 与 Duncan（1994）开发的专门用于行星轨道传播的积分器，它采用了正则化技术与混合变量辛算法，能够在保证高精度的同时显著提高计算效率，非常适合进行长时间尺度的天体力学数值模拟。

在本次研究中，数值积分的时间跨度为 1 亿年（100 Myr），初始积分步长设置为 30 天。为了减少数据存储量并便于后续分析，研究人员每 500 年记录一次行星与测试粒子的轨道参数。所有冥族小天体均被视为无质量的测试粒子，这一假设在研究小天体的轨道动力学时是合理的，因为冥族小天体的质量远小于巨行星的质量，它们之间的引力相互作用可以忽略不计。

3.3 长期稳定冥族小天体的识别

并非所有被观测到的冥族小天体都能在太阳系演化的漫长历史中保持轨道稳定。一些天体可能会由于混沌扩散或与巨行星的近距离相遇而脱离 3:2 平运动共振，最终被散射出柯伊伯带。因此，在识别双摆动冥族小天体之前，首先需要筛选出长期稳定的冥族小天体。

Malhotra 与 Ito（2024）采用临界共振角 $\sigma$ 的摆动范围作为长期稳定性的判据。具体来说，如果一个冥族小天体在整个 1 亿年的积分时间内，$\sigma$ 模 $360^{\circ }$ 的最大值与最小值之差小于 $358^{\circ }$，则认为该天体是长期稳定的。这一判据能够有效排除那些在积分过程中脱离共振的天体，因为一旦天体脱离共振，$\sigma$ 将不再摆动，而是随时间单调增加，其范围将接近 $360^{\circ }$（通常 > $359^{\circ }$）。设置 $358^{\circ }$ 的阈值可以在排除脱离共振天体的同时，包容共振内的大振幅摆动。

通过这一判据，研究人员从初始的 453 个冥族小天体中筛选出了 441 个长期稳定的天体。初始样本量与稳定样本量之间的差异主要源于本次研究采用了更长的积分时间（100 Myr），而 Volk 与 Van Laerhoven（2024）的研究仅采用了 10 Myr 的积分时间。尽管如此，441 个天体的样本量仍然是迄今为止最大的冥族小天体研究样本，远大于之前研究中不足 100 个天体的样本量。

3.4 双摆动冥族小天体的识别

在筛选出长期稳定的冥族小天体后，研究人员进一步识别出其中表现出持续近日点幅角 $g$ 摆动的双摆动天体。在识别过程中，一个关键的问题是参考平面的选择。JPL Horizons 系统默认采用 J2000 黄道作为参考平面，但在长期动力学时间尺度上，行星的轨道平面会绕太阳系总角动量矢量垂直的平面（即不变平面）进动。由于巨行星是影响冥族小天体长期动力学的主要摄动体，采用不变平面作为参考平面能够更准确地描述冥族小天体的轨道演化。因此，Malhotra 与 Ito（2024）将所有轨道元素转换到太阳系质心的不变平面坐标系下进行计算。

根据经典的（非共振）vZLK 理论以及共振条件下的动力学分析，在 3:2 平运动共振中，近日点幅角 $g$ 可以围绕 $0^{\circ }$、$90^{\circ }$、$180^{\circ }$ 与 $270^{\circ }$ 四个中心摆动。因此，研究人员针对每个长期稳定的冥族小天体，检查其 $g(t)$ 时间序列是否围绕这四个中心中的任何一个持续摆动。识别判据为：在整个 1 亿年的积分时间内，$g-g_c$ 模 $360^{\circ }$ 的最大值与最小值之差小于 $179^{\circ }$，其中 $g_c$ 为摆动中心（$0^{\circ }$、$90^{\circ }$、$180^{\circ }$ 或 $270^{\circ }$），且 $g-g_c$ 模 $360^{\circ }$ 被映射到 $[-180^{\circ },+180^{\circ }]$ 区间内。

通过这一判据，研究人员在 441 个长期稳定的冥族小天体中识别出了 69 个双摆动天体，占比约为 16%。其中，没有发现围绕 $0^{\circ }$ 或 $180^{\circ }$ 摆动的天体；摆动中心为 $90^{\circ }$ 的有 34 个，为 $270^{\circ }$ 的有 35 个，两者近乎均分。这一结果与之前的研究结果大致相符，但由于样本量更大，其统计显著性更高。需要注意的是，这一判据要求天体在整个 1 亿年的积分时间内持续表现出 $g$ 摆动，因此会排除那些表现出间歇性摆动的天体（见下文第 4.3 节）。

四、双摆动冥族小天体的轨道参数分布特征

4.1 整体冥族小天体的轨道参数分布

在分析双摆动冥族小天体的轨道参数分布之前，首先需要了解整体冥族小天体的轨道参数分布特征。Malhotra 与 Ito（2024）绘制了 441 个长期稳定冥族小天体的吻切半长轴、偏心率与倾角的散点图与直方图（见原文图 1）。

结果显示：

半长轴分布：冥族小天体的半长轴集中在 39.4 AU 附近，这与 3:2 平运动共振的理论半长轴一致。半长轴的分布范围较窄，这是因为平运动共振对天体的半长轴有严格的约束，只有半长轴在共振带内的天体才能保持长期稳定。

偏心率分布：冥族小天体的偏心率低于 0.1 的非常罕见，分布峰值出现在约 0.25 处，且在偏心率高于 0.35 时急剧下降。这一分布特征主要由两个因素决定：一是长期动力学稳定性，偏心率较高的冥族小天体在近日点时会更靠近天王星，容易受到天王星的引力摄动而脱离共振；二是共振相空间体积，海王星的 3:2 平运动共振在 $(a,e)$ 平面上的低偏心率区域存在一个狭窄的“颈部”，导致低偏心率区域的相空间体积较小，因此能够容纳的天体数量较少（Lan & Malhotra, 2019; Balaji et al., 2023）。

倾角分布：冥族小天体的倾角主要集中在不变平面附近的几度范围内，但存在一个明显的高倾角尾部，部分天体的倾角可达 $30^{\circ }$ 以上，甚至有少数天体的倾角超过 $50^{\circ }$。这一高倾角尾部一直是柯伊伯带演化理论中的一个难题，因为经典的巨行星迁移模型难以解释如此多的高倾角天体的起源。此外，由于大多数柯伊伯带观测巡天都集中在不变平面附近的天区，高倾角天体的观测偏差较大，实际的高倾角天体数量可能比观测到的更多（Bannister et al., 2016; Bernstein et al., 2004; Petit et al., 2011）。

4.2 双摆动冥族小天体的轨道参数分布差异

通过将双摆动冥族小天体与整体冥族小天体的轨道参数分布进行对比，Malhotra 与 Ito（2024）发现双摆动天体具有显著更高的偏心率与倾角（见原文图 2）。

具体来说，双摆动天体的偏心率平均值约为 0.28，而整体冥族小天体的偏心率平均值约为 0.23；双摆动天体的倾角平均值约为 $15^{\circ }$，而整体冥族小天体的倾角平均值约为 $10^{\circ }$。这一结果与之前的研究结论一致，表明 vZLK 振荡更容易在高偏心率与高倾角的轨道上发生。

这张图直观展示了海王星（内圈蓝轨）与冥王星（外圈红轨）的3:2轨道共振关系。数字标记显示了共振角 $\phi =3{\lambda }_{Neptune}-2{\lambda }_{Pluto}-{\varpi }_{Pluto}$ 的相位关系，正是这个角度在轨道运动中的周期性双摆动（libration，范围约76°-104°）维持了两者的稳定共振，避免了轨道交叉碰撞。

4.3 偏心率-倾角平面上的聚集特征

Malhotra 与 Ito（2024）研究中最重要的发现之一是，双摆动冥族小天体在时间平均的偏心率-倾角平面上沿一条近似双曲线的弧聚集。为了消除轨道短期振荡的影响，研究人员计算了每个天体在整个 1 亿年积分时间内的平均偏心率 $\bar{e}$ 与平均倾角 $\bar{i}$，并绘制了散点图（见原文图 2(b) 和图 3，其中图 3见下图）。结果显示，所有 69 个双摆动天体都紧密地分布在一条双曲线弧附近，而普通冥族小天体则在 $\bar{e}$-$\bar{i}$ 平面上的分布显得弥散且随机。

值得注意的是，部分非双摆动冥族小天体也靠近这条双曲线弧。通过检查这些天体的 $g(t)$ 时间序列，研究人员发现它们表现出间歇性的 $g$ 摆动——有时围绕 $g=+90^{\circ }$ 摆动，有时围绕 $g=-90^{\circ }$ 摆动，有时则表现为循环运动。这表明这些天体处于双摆动与循环运动的边界区域，其轨道动力学具有较强的混沌性，因此未能满足“持续摆动”的严格判据。如果采用更宽松的识别标准，双摆动冥族小天体的比例将会更高。

为了验证这一聚集特征的物理意义，研究人员计算了两个理论模型下的扰动函数极小值轨迹，并与观测数据进行了对比（见原文图 3）。

第一个模型是圆限制三体问题（CR3BP），即由太阳、海王星与一个无质量的冥族小天体组成的系统；第二个模型是圆限制六体问题（CR6BP），即由太阳、四颗巨行星与一个无质量的冥族小天体组成的系统。在两个模型中，行星的轨道都被假设为圆形且共面。研究人员对平均黄经进行了轨道平均，并考虑了 3:2 平运动共振的约束（即 $a/a_{\text{Neptune}}=(3/2{)}^{2/3}$ 且 $\sigma =180^{\circ }$），计算了不同偏心率与倾角下的扰动函数值，并找到了 $g=90^{\circ }$ 处的局部极小值轨迹。

结果显示，CR6BP 模型得到的红色曲线与双摆动冥族小天体的观测分布吻合得非常好，而 CR3BP 模型得到的黑色点划线则位于更高的偏心率区域。这一结果表明，木星、土星与天王星的引力摄动对双摆动冥族小天体的轨道特征有显著影响——它们的平均效应使得扰动函数极小值轨迹向低偏心率方向偏移。

此外，原文还报告了一个例外天体——2015 VR164。该天体的时间平均倾角仅约 $0.22^{\circ }$，偏心率约 $0.206$，显著偏离图 3 中的红色曲线。由于其轨道平面几乎与不变平面重合，作者推测其 $g$ 摆动可能源于参考平面选择相关的运动学效应，而非真正的动力学 vZLK 振荡，详细分析留待未来研究。

五、动力学机制与理论解释

5.1 扰动函数极小值与摆动中心的稳定性

在天体力学中，小天体的轨道运动可以用哈密顿系统来描述，而扰动函数则描述了摄动体对小天体轨道的引力影响。对于处于共振状态的小天体，其哈密顿系统的平衡点对应于扰动函数的极值点，而平衡点的稳定性则决定了小天体是否能够围绕该点摆动。

对于双摆动冥族小天体，其近日点幅角 $g$ 围绕 $\pm 90^{\circ }$ 摆动，这意味着在 $g=\pm 90^{\circ }$ 处存在扰动函数的局部极小值，且这些极小值是稳定的。Malhotra 与 Ito（2024）的计算结果表明，在 3:2 平运动共振条件下，$g=\pm 90^{\circ }$ 处的局部极小值存在于一定的偏心率与倾角范围内，这些极小值的轨迹构成了偏心率-倾角平面上的双曲线弧。只有当天体的平均偏心率与平均倾角位于这条弧附近时，其近日点幅角 $g$ 才能保持稳定的摆动。

这张图来自天文学顶刊 The Astronomical Journal，展示了冥族小天体在* Kozai共振**与2:3共振耦合作用下的相空间轨迹。图中绿色和红色轨道分别代表不同冥族小天体的双摆动轨迹，在($e\cdot cos\omega$, $e\cdot sin\omega$)空间中形成闭合的libration环。*

5.2 庞加莱第三类周期轨道

Malhotra 与 Ito（2024）进一步指出，双摆动冥族小天体在偏心率-倾角平面上的聚集结构对应于圆限制六体问题中的庞加莱第三类周期轨道族。庞加莱在 1892 年出版的《天体力学新方法》一书中，将限制性三体问题中的周期轨道分为三类：第一类周期轨道是指偏心率与倾角都很小的近圆轨道；第二类周期轨道是指偏心率较大但倾角很小的轨道；第三类周期轨道则是指同时具有较大偏心率与较大倾角的轨道。

庞加莱第三类周期轨道的一个重要特征是，它们的近日点幅角保持不变，即 $g=$ 常数。因此，围绕这些周期轨道的小天体将表现出近日点幅角的摆动。在圆限制三体问题中，第三类周期轨道族对应于 $g=\pm 90^{\circ }$ 处的扰动函数极小值轨迹，这与双摆动冥族小天体的摆动中心一致。而在圆限制六体问题中，由于木星、土星与天王星的额外摄动，第三类周期轨道族的位置发生了偏移——向低偏心率方向移动，这与观测结果相符（见原文图 3）。

这一发现具有重要的理论意义，它将双摆动冥族小天体的动力学特征与经典的天体力学理论联系起来，为理解共振条件下的 vZLK 振荡提供了新的视角。同时，它也表明，庞加莱在一个多世纪前发现的周期轨道族在太阳系的实际天体系统中仍然发挥着重要作用。

5.3 双摆动冥族小天体的丰度异常

Malhotra 与 Ito（2024）研究中另一个重要的结果是，双摆动冥族小天体的观测丰度远高于随机分布预期。为了估计随机分布下的双摆动天体比例，研究人员采用了简化的单自由度摆模型。基于遍历假设，一个天体在共振区内被随机观测到时，其处于 $g$ 摆动状态的概率大致等于两种振荡周期之比。在该模型中，共振带内的相空间体积与小振幅摆动的频率成正比，即与摆动周期成反比。

根据数值模拟结果，冥族小天体的 $\sigma$ 摆动周期约为 2 万年（20 kyr），而 $g$ 摆动周期约为 400 万年（4 Myr），这与冥王星本身的摆动周期相近（Malhotra & Williams, 1997）。因此，随机分布下的双摆动天体预期比例约为 $20\text{kyr}/4\text{Myr}=0.005$，即 0.5%。然而，观测到的双摆动天体比例约为 16%，是预期值的 30 倍以上。这一丰度异常表明，冥族小天体的相空间并非随机均匀分布，而是存在明显的局部过密。

为了解释这一丰度异常，Malhotra 与 Ito（2024）提出了两个可能的动力学机制，这两个机制并非相互排斥，而是可能协同作用：

第一个机制：长期动力学稳定性效应。 双摆动冥族小天体的近日点幅角围绕 $\pm 90^{\circ }$ 摆动，这意味着它们的近日点始终位于轨道平面的最高点或最低点，远离不变平面。因此，双摆动天体在近日点时与巨行星的距离更远，受到的引力摄动更小，轨道更稳定。相比之下，普通冥族小天体的近日点幅角是循环变化的，当它们的近日点经过不变平面时，更容易受到巨行星的摄动而脱离共振。在太阳系演化的数十亿年历史中，这种稳定性差异会导致普通冥族小天体的数量逐渐减少，而双摆动冥族小天体的相对比例逐渐增加。这一效应已经在之前的一些研究中得到了证实（Kuchner et al., 2002; Tiscareno & Malhotra, 2009）。

第二个机制：巨行星迁移过程中的优先捕获效应。 根据巨行星迁移理论，冥族小天体是在海王星向外迁移的过程中被捕获进入 3:2 平运动共振的。共振捕获的效率与迁移速度密切相关：迁移速度越慢，捕获效率越高，且越容易捕获进入较弱的共振。$g$ 摆动共振是一种比 $\sigma$ 摆动共振更弱的共振，其摆动周期更长。因此，在巨行星迁移的后期，当迁移速度减慢时，更多的天体能够被同时捕获进入 $\sigma$ 摆动共振与 $g$ 摆动共振，形成双摆动冥族小天体。这一效应也得到了数值模拟研究的支持（Chiang & Jordan, 2002; Nesvorný & Vokrouhlický, 2016; Pike et al., 2017; Balaji et al., 2023）。

Malhotra 与 Ito（2024）指出，这两个机制可能共同作用：优先捕获为双摆动族群提供了初始“种子”，而长期稳定性效应则在随后的数十亿年中对普通族群进行了“筛选”，共同造就了今天观测到的显著丰度异常。未来的数值模拟研究需要同时考虑这两个效应，以更准确地重现双摆动冥族小天体的观测丰度。

六、研究结果的科学意义与争议

6.1 对巨行星迁移理论的约束

双摆动冥族小天体的丰度与轨道分布特征为约束巨行星迁移理论提供了新的观测证据。如前所述，巨行星迁移的速度与时间尺度直接影响双摆动天体的捕获效率：如果海王星迁移速度过快，那么只有极少数天体能够被捕获进入 $g$ 摆动共振；反之，如果迁移速度过慢，那么双摆动天体的比例将会过高。因此，通过将观测到的双摆动天体比例（16%）与数值模拟结果进行对比，可以精确约束海王星迁移的速度与时间尺度。

此外，双摆动冥族小天体在偏心率-倾角平面上的聚集结构（图 3 中的红色曲线与绿色点的高度吻合）也对巨行星迁移的模式提出了约束。如果巨行星的迁移是平滑且缓慢的，那么被捕获进入共振的天体将保留其原始的偏心率与倾角分布，从而形成清晰的聚集结构；反之，如果巨行星的迁移是剧烈且不连续的，那么天体的偏心率与倾角分布将会被打乱，难以形成明显的聚集特征。因此，观测到的清晰聚集结构支持平滑缓慢的巨行星迁移模式。

6.2 对柯伊伯带演化模型的意义

双摆动冥族小天体的研究也为理解柯伊伯带的整体演化提供了新的线索。例如，高倾角双摆动天体的存在表明，柯伊伯带中存在大量的高倾角原始天体，这对传统的柯伊伯带形成模型提出了挑战。传统模型认为，柯伊伯带天体形成于一个薄的原行星盘中，其初始倾角很小，高倾角天体是通过后期的动力学散射形成的。然而，双摆动天体的高倾角是在共振捕获过程中获得的，这表明部分高倾角天体可能在巨行星迁移之前就已经存在。

需要指出的是，Malhotra 与 Ito（2024）的原文中并未提出存在由“约 40 个普通冥族小天体组成的碰撞族”。原文图 2(b) 中靠近双曲线弧的蓝色点（非双摆动天体）已被作者明确解释为具有间歇性 $g$ 摆动的天体，而非独立的碰撞族。因此，在引用该研究时应避免对此进行错误解读。

6.3 观测偏差与研究局限性

尽管 Malhotra 与 Ito（2024）的研究取得了重要进展，但仍然存在一些局限性需要注意。
首先，观测偏差会对双摆动冥族小天体的丰度估计产生显著影响。如前所述，大多数柯伊伯带观测巡天都集中在不变平面附近的天区，高倾角天体更容易被观测遗漏。而双摆动天体往往具有较高的倾角，因此实际的双摆动天体比例可能比观测到的 16% 更高。Lawler 与 Gladman（2013）的模拟结果表明，考虑观测偏差后，真实的双摆动天体比例可能比观测值高出 20%–40%。

其次，数值模拟的积分时间有限。本次研究采用了 1 亿年的积分时间，这对于研究冥族小天体的中期轨道演化是足够的，但太阳系的年龄约为 46 亿年，更长时间尺度的混沌效应可能会对天体的轨道稳定性产生影响。一些在 1 亿年积分时间内表现为稳定的双摆动天体，可能在数十亿年的时间尺度上脱离共振或转变为循环运动。因此，未来的研究需要进行更长时间尺度的数值模拟，以更准确地评估双摆动天体的长期稳定性。

此外，理论模型采用了简化假设。在计算扰动函数极小值轨迹时，研究人员假设行星的轨道是圆形且共面的，但实际行星的轨道具有一定的偏心率与倾角。非圆非共面的行星轨道会引入额外的摄动项，可能会改变扰动函数极小值的位置与稳定性。因此，未来的理论研究需要考虑更真实的行星轨道模型，以提高理论预测的准确性。

最后，本次研究没有考虑天体的质量与碰撞效应。所有冥族小天体都被视为无质量的测试粒子，它们之间的引力相互作用与碰撞效应都被忽略了。然而，在柯伊伯带的演化过程中，碰撞是一个重要的过程，它会改变天体的轨道参数与物理性质。因此，未来的研究需要结合碰撞演化模型，以更全面地理解冥族小天体的轨道分布特征。

七、总结与未来研究方向

7.1 总结

Malhotra 与 Ito（2024）的研究是双摆动冥族小天体研究领域的一项重要成果。它基于迄今为止最大的冥族小天体样本（441 个长期稳定天体），通过长达 1 亿年的数值轨道积分，系统地研究了双摆动冥族小天体的轨道动力学特征。该研究的主要结论包括：

丰度估计：约 16% 的长期稳定冥族小天体表现出持续的近日点幅角 $g$ 摆动，这些天体具有显著更高的偏心率与倾角。

聚集结构：双摆动冥族小天体在时间平均的偏心率-倾角平面上沿一条近似双曲线的弧聚集。这一结构对应于圆限制六体问题中的庞加莱第三类周期轨道族，且由于木星、土星与天王星的平均摄动效应，该弧线相对于圆限制三体问题中的理论轨迹向低偏心率方向偏移。

丰度异常：双摆动冥族小天体的观测丰度（16%）显著高于基于频率比估算的随机分布预期（0.5%）。这一“丰度异常”可由长期动力学稳定性效应与巨行星迁移过程中的优先捕获效应的联合作用得到合理解释。

这些结果不仅深化了我们对共振条件下 vZLK 振荡的理解，也为约束巨行星迁移理论与柯伊伯带演化模型提供了新的观测证据。

7.2 未来研究方向

尽管本次研究取得了重要进展，但双摆动冥族小天体的研究仍然存在许多问题有待解决。未来的研究可以从以下几个方面展开：

更长时间尺度的数值模拟：进行数十亿年尺度的数值轨道积分，研究双摆动冥族小天体在太阳系整个演化历史中的轨道稳定性与演化过程，更准确地评估长期动力学稳定性效应对丰度的影响。

更精确的理论模型：考虑非圆非共面的行星轨道，建立更真实的圆限制六体问题模型，计算更精确的扰动函数极小值轨迹与周期轨道族，提高理论预测与观测数据的吻合度。

更大样本量的观测数据：随着新一代地面与空间望远镜（如 Vera C. Rubin 天文台）的投入使用，将会发现更多的柯伊伯带天体，特别是高倾角天体。更大的样本量将有助于提高双摆动天体丰度估计的准确性，并验证偏心率-倾角平面上聚集结构的统计显著性。

例外天体的深入研究：对 2015 VR164 这类偏离聚集弧的疑似双摆动天体进行详细的动力学分析，厘清其 $g$ 摆动的真正成因（动力学 vs. 运动学效应）。

相空间体积的精确计算：采用更先进的数值方法（如频率映射分析）计算共振带内的相空间体积，特别是 $g$ 摆动区的相空间体积，更准确地评估双摆动天体的丰度异常程度。

多物理过程的耦合模拟：将巨行星迁移、动力学散射与碰撞演化等物理过程耦合起来，建立更全面的柯伊伯带演化模型，重现双摆动冥族小天体的观测丰度与轨道分布特征。

总之，双摆动冥族小天体的研究是一个充满活力的领域，它将继续为我们揭示太阳系早期的形成与演化历史提供重要的线索。随着观测技术的不断进步与理论模型的不断完善，我们对这一神秘天体族群的理解将会越来越深入。

参考文献

[1] Alexandersen, M., Gladman, B., Kavelaars, J. J., et al. 2016, “A Carefully Characterized and Tracked Trans-Neptunian Survey: The Size Distribution of the Plutinos and the Number of Neptunian Trojans”, AJ, 152, 111
[2] Balaji, S., Zaveri, N., Hayashi, N., et al. 2023, “Can the orbital distribution of Neptune‘s 3:2 mean-motion resonance result from stability sculpting?”, MNRAS, 524, 3039
[3] Bannister, M. T., Kavelaars, J. J., Petit, J.-M., et al. 2016, “The Outer Solar System Origins Survey. II. The Brightness-independent Inclination Distribution of Trans-Neptunian Objects”, AJ, 152, 70
[4] Bernstein, G. M., Trilling, D. E., Allen, R. L., et al. 2004, “The Size Distribution of Trans-Neptunian Bodies”, AJ, 128, 1364
[5] Brouwer, D., & Clemence, G. M. 1961, Methods of Celestial Mechanics (New York: Academic Press)
[6] Chiang, E. I., & Jordan, A. B. 2002, “On the Plutinos and Twotinos of the Kuiper Belt”, AJ, 124, 3430
[7] Gladman, B., & Volk, K. 2021, “Transneptunian Space”, ARA&A, 59, 203
[8] Gomes, R. S. 2003, “The origin of the Kuiper Belt high-inclination population”, Icarus, 161, 404
[9] Ito, T., & Malhotra, R. 2023, “Secular dynamics of resonant Kuiper Belt objects”, in LPI Contributions, Vol. 2851, Asteroids, Comets, Meteors Conference, 2391
[10] Ito, T., & Ohtsuka, K. 2019, “The Lidov-Kozai Oscillation and von Zeipel‘s Theorem”, Monographs on Environment, Earth and Planets, 7, 1
[11] Jefferys, W. H., & Standish, E. M., Jr. 1966, “Further periodic solutions of the three-dimensional restricted problem”, AJ, 71, 982
[12] Jefferys, W. H., & Standish, E. M., Jr. 1972, “Further Periodic Solutions of the Three-Dimensional Restricted Problem. II”, AJ, 77, 394
[13] Kozai, Y. 1962, “Secular perturbations of asteroids with high inclination and eccentricity”, AJ, 67, 591
[14] Kozai, Y. 1985, “Secular perturbations of resonant asteroids”, Celes. Mech., 36, 47
[15] Kuchner, M. J., Brown, M. E., & Holman, M. 2002, “Long-Term Dynamics and the Orbital Inclinations of the Classical Kuiper Belt Objects”, AJ, 124, 1221
[16] Lan, L., & Malhotra, R. 2019, “Neptune‘s resonances in the scattered disk”, Celes. Mech. Dyn. Astron., 131, 39
[17] Lawler, S. M., & Gladman, B. 2013, “Plutino Detection Biases, Including the Kozai Resonance”, AJ, 146, 6
[18] Lei, H., Li, J., Huang, X., & Li, M. 2022, “Dynamics of High-inclination Plutinos”, AJ, 164, 74
[19] Levison, H. F., & Duncan, M. J. 1994, “The Long-Term Dynamical Behavior of Short-Period Comets”, Icarus, 108, 18
[20] Li, J., Zhou, L.-Y., & Sun, Y.-S. 2014, “The orbital distribution of trans-Neptunian objects beyond 50 au”, MNRAS, 437, 215
[21] Malhotra, R. 1993, “The origin of Pluto‘s peculiar orbit”, Nature, 365, 819
[22] Malhotra, R. 1995, “The origin of Pluto‘s orbit: implications for the Solar System beyond Neptune”, AJ, 110, 420
[23] Malhotra, R. 2019, “The migration of the giant planets and the formation of the Kuiper belt”, Geoscience Letters, 6, 12
[24] Malhotra, R., & Ito, T. 2025, “The doubly librating Plutinos”, ApJ, in press (arXiv:2501.01586)
[25] Malhotra, R., & Williams, J. G. 1997, “Pluto‘s Heliocentric Orbit”, in Pluto and Charon, ed. S. A. Stern & D. J. Tholen (Tucson: University of Arizona Press), 127
[26] Morbidelli, A. 2002, Modern Celestial Mechanics: Aspects of Solar System Dynamics (London: Taylor & Francis)
[27] Murray, C. D., & Dermott, S. F. 1999, Solar System Dynamics (Cambridge: Cambridge University Press)
[28] Nesvorný, D., Roig, F., & Ferraz-Mello, S. 2000, “Close Approaches of Trans-Neptunian Objects to Pluto Have Left Observable Signatures on Their Orbital Distribution”, AJ, 119, 953
[29] Nesvorný, D. 2018, “Dynamical Evolution of the Early Solar System”, ARA&A, 56, 137
[30] Nesvorný, D., & Vokrouhlický, D. 2016, “Neptune‘s Orbital Migration Was Grainy, Not Smooth”, ApJ, 825, 94
[31] Petit, J.-M., Kavelaars, J. J., Gladman, B. J., et al. 2011, “The Canada-France Ecliptic Plane Survey—Full Data Release: The Orbital Structure of the Kuiper Belt”, AJ, 142, 131
[32] Pike, R. E., Lawler, S., Brasser, R., et al. 2017, “The Structure of the Distant Kuiper Belt in a Nice Model Scenario”, AJ, 153, 127
[33] Poincaré, H. 1892, Les méthodes nouvelles de la mécanique céleste, Tome I: Solutions périodiques. Non-existence des intégrales uniformes. Solutions asymptotiques (Paris: Gauthier-Villars et fils)
[34] Saillenfest, M., Fouchard, M., Tommei, G., & Valsecchi, G. B. 2016, “Long-term dynamics beyond Neptune: secular models to study the regular motions”, Celes. Mech. Dyn. Astron., 126, 369
[35] Saillenfest, M. 2020, “Long-term orbital dynamics of trans-Neptunian objects”, Celes. Mech. Dyn. Astron., 132, 12
[36] Souami, D., & Souchay, J. 2012, “The solar system‘s invariable plane”, A&A, 543, A133
[37] Tiscareno, M. S., & Malhotra, R. 2009, “Chaotic Diffusion of Resonant Kuiper Belt Objects”, AJ, 138, 827
[38] Tremaine, S. 2023, Dynamics of Planetary Systems (Princeton: Princeton University Press)
[39] Volk, K., & Van Laerhoven, C. 2024, “Dynamical Classifications of Multi-opposition TNOs as of 2023 December”, RNAAS, 8, 36