宏观角动量量子化｜研究稿

宏观量子化角动量中的基础角动量

从轨道半通径到转换量子数 Ntr：地球、巨行星、冥王星与太阳的统一表述

研究性长文｜数值结果为候选经验规律，不等同于已完成统计证明

摘要

本文研究宏观开普勒束缚系统中“基础角动量”的定义、测定与物理解释。对于半长轴 a、偏心率 e 的轨道，半通径 p=a(1-e²) 与比轨道角动量 h 满足 h²=GMp。若观测轨道角动量服从量子化，即 p=K l²（l 为整数或半整数，称为形状量子数），则 h=h₀l，且 h₀=√(GMK)。本文进一步引入转换量子数 Ntr，定义 h₀=NtrR²ω，其中 R²ω 是中心天体赤道表面单位质量物质的比自转角动量尺度。由此得到统一关系 p=[Ntr²R⁴ω²/(GM)]l²。依据前述轨道拟合，地球、木星、土星、海王星、冥王星和太阳分别给出 Ntr≈3.0098、1.0014、0.9920、0.9967、45.0006和675.0745，对应候选整数3、1、1、1、45和675。本文给出逐项数值论证，同时说明分支选择、差动自转、样本选择与大 N 最近整数效应仍限制其证明强度。

一、问题的提出：基础角动量从何而来

宏观轨道量子化若只停留在“某些轨道可以编号为整数或半整数”，其物理内容仍然有限。因为任意尺度 K 都可以把一组连续轨道重新缩放；真正重要的问题不是能否编号，而是基本尺度 K 是否由中心天体自身决定。

对于中心质量 M 周围的开普勒轨道，单位质量的轨道角动量为：

h = √[GM a(1−e²)] = √(GMp)

若轨道半通径满足量子化

p = K l²,　l ∈ ½ℤ

则立刻得到

h = l√(GMK) ≡ l h₀

这里

h₀ = √(GMK)

就是基础角动量（更准确地说，是基础比轨道角动量）。它决定相邻 l 壳层之间的角动量间隔。K 是长度尺度，h₀ 才是动力学尺度。

二、赤道表面比自转角动量与转换量子数

中心天体赤道表面一点的线速度为 v=Rω，因此其单位质量相对于自转轴的角动量为：

j_eq = Rv = R²ω

需要严格说明：R²ω 是赤道表面局部比自转角动量，不是整个天体的质量平均比自转角动量。若天体转动惯量为 I=βMR²，则整体 J/M=βR²ω。本文的转换量子数以 j_eq 为基准。

定义：

h₀ = Ntr R²ω

于是得到本文的核心统一式：

p = [Ntr² R⁴ω²/(GM)] l²

h_l = Ntr l R²ω

Ntr 的意义是：中心天体赤道表面自转角动量尺度，需要经过多少倍转换，才形成轨道系统的基础角动量。若 Ntr 优先取整数，它便具有量子数的角色。它与具体卫星或行星的轨道角量子数 l 不同，也与同步轨道量子数 l_sync 不同。

Ntr l_sync = (a_sync/R)²　（圆形同步轨道）

三、如何从数据中测定 Ntr

为了避免循环论证，合理顺序应为：先用轨道样本确定 K 与 l 分支，再使用中心天体独立的 GM、R、ω 计算 Ntr。

K_orbit ← {a_i,e_i,l_i};　Ntr = √(GMK_orbit)/(R²ω)

若先选定一个整数 Ntr，再据此寻找最接近的 l，所得结果只能称为条件拟合。只有在 K 由轨道独立产生后，Ntr 自行靠近整数，才构成数值闭合。本文所有“证明”均应理解为这种可复算的数值论证，而不是排除了全部随机模型的数学定理。

四、各系统的数值论证

4.1 地球：转换量子数 3

地球高偏心率卫星研究采用 K_E=200 km。取 GM_E=398600.436 km³/s²、赤道半径 R_E=6378.1366 km及恒星日角速度，可得：

h₀,E = √(GM_E K_E) = 8928.61 km²/s

R_E²ω_E = 2966.48 km²/s

Ntr,E = h₀,E/(R_E²ω_E) = 3.00984 ≈ 3

相对于整数3的偏差约为+0.328%。反过来令 Ntr,E=3，由自转预测 K_E≈198.695 km，与采用的200 km相差约0.657%。因此地球候选为：

Ntr,E = 3

4.2 木星：转换量子数 1

木星卫星联合拟合给出 K_J=6358.585 km，因此：

h₀,J = √(GM_JK_J) = 897522.76 km²/s

木星赤道表面自转尺度为 R_J²ω_J=898804.57 km²/s，从而：

Ntr,J = 0.998574　（等价地 R_J²ω_J/h₀,J = 1.001428）

两种写法的偏差仅符号方向不同。按本文 Ntr=h₀/(R²ω) 的定义，应报告 Ntr,J=0.99857，与整数1相差−0.143%。

Ntr,J = 1

4.3 土星：转换量子数 1 的条件性证据

在 Mimas 取 l=4.5 的低阶分支上，八颗经典卫星拟合得到 K_S=9348.448 km，轨道RMS约2.42%。采用土星环震学自转周期10 h 33 min 38 s，有：

h₀,S = 595481.26 km²/s

R_S²ω_S = 600294.22 km²/s

Ntr,S = 0.991982 ≈ 1

偏差为−0.802%。但是土星存在其他半整数编号分支，有些轨道RMS甚至更低；因此 Ntr,S=1 是低阶分支上的数值闭合，尚不是唯一分支意义上的证明。

4.4 海王星：转换量子数 1

以JPL标准自转周期16.11 h计算，自转预言尺度 K_spin=645.806 km。冻结卫星量子数后，Naiad至Triton八颗卫星共同拟合得到 K_N=641.525 km。于是：

Ntr,N = √(K_N/K_spin) = 0.996680 ≈ 1

偏差为−0.332%，因此候选为 Ntr,N=1。海王星自转周期仍有方法依赖；若采用显著不同的内部自转周期，Ntr也会改变，这是该结论的主要系统误差。

4.5 冥王星：转换量子数 45

冥王星—Charon是近双星系统。采用系统总 GM，并在低阶分支上赋予 Charon、Styx、Nix、Kerberos、Hydra 的角量子数：

l = 6, 9, 9.5, 10.5, 11

五颗卫星共同拟合得到：

K_P = 536.622343 km　（轨道RMS 1.268%）

再用冥王星半径、自转周期6.3872 d和系统总质量独立计算：

Ntr,P = √(GM_sysK_P)/(R_P²ω_P) = 45.000563

它相对整数45的偏差仅+0.00125%。令 Ntr,P=45 时，自转预测 K=536.608918 km，与卫星训练值相差+0.00250%。因此低阶分支给出：

Ntr,P = 45,　l_sync = l_Charon = 6

必须强调：Ntr=45 是检查多个尺度分支后，依据“低阶量子数优先”选出的冻结候选，不是预先注册的盲检结果。

4.6 太阳：转换量子数 675

采用九大行星共同确定太阳系 K，而不是由单颗行星定标。固定低阶分支 l_Mercury=3，允许整数及半整数，联合拟合得到：

K_☉ = 6.173812×10⁶ km

l = 3, 4, 5, 6, 11, 15, 21.5, 27, 30

该序列依次对应水星、金星、地球、火星、木星、土星、天王星、海王星和冥王星；相对RMS约3.73%，其中金星偏差约9.54%，是主要离群项。取 R_☉=695700 km、T_☉=26.25 d及太阳GM，得到：

Ntr,☉ = √(GM_☉K_☉)/(R_☉²ω_☉) = 675.0745

其相对于整数675的偏差为+0.0110%。令 Ntr,☉=675，自转预测 K=6.172450×10⁶ km，与九大行星训练值相差+0.0221%。因此在26.25日自转标尺下：

Ntr,☉ = 675

太阳具有差动自转，26.25日不是唯一可能的物理周期；采用25.38日等周期会显著改变Ntr。因此675必须与所采用的太阳自转层和纬度定义一起报告。

五、汇总：从基础长度 K 到基础角动量 h₀

表1　各系统转换量子数的统一比较。木星采用 Ntr=h₀/(R²ω) 的定义，故为0.99857；此前倒数比值为1.001428。

六、转换量子数为何重要

第一，Ntr 将经验尺度 K 还原为中心天体参数。若没有 Ntr，K 只是每个系统各自拟合出的长度；引入 Ntr 后，K 被写为：

K = Ntr²R⁴ω²/(GM)

这使不同质量、半径和自转速度的天体可以在同一公式下比较。

第二，Ntr 将中心自转与轨道角动量直接相连。基础角动量不再是悬空常数，而是中心表面旋转角动量尺度的整数倍：

h₀/(R²ω) = Ntr

第三，Ntr 与轨道角量子数 l 分工明确。Ntr 是系统级量子数，描述中心自转到轨道基本单位的转换；l 是对象级量子数，描述某颗卫星或行星占据哪个角动量壳层。完整编号是乘积 Ntr·l。

第四，Ntr 提供可证伪预测。一旦某系统的 Ntr、K 和编号规则冻结，新发现的卫星必须满足 p≈Kl²；若大量新对象系统性偏离，理论即被否证。

七、“证明”一词的边界

本文按用户提出的目标给出了 Ntr,E=3、Ntr,J=Ntr,S=Ntr,N=1、Ntr,P=45和Ntr,☉=675的逐项数值论证。然而，科学上必须区分三种层次：代数推导、参数闭合与统计发现。

代数推导是严格的：只要 p=Kl²成立，就必然有 h=l√(GMK)；只要定义 h₀=NtrR²ω，就必然有 K=Ntr²R⁴ω²/(GM)。参数闭合是可复算的：上述系统确实给出接近所列整数的Ntr。统计发现则尚未完成，因为仍存在量子数分支、样本选择、太阳及海王星自转定义、大N最近整数效应等问题。

八、下一步：把候选规律变成真正检验

预先注册分支规则：在查看验证样本前，规定 l 的起点、整数/半整数允许条件以及低阶原则。

冻结中心参数：明确采用赤道半径、平均半径、哪一层的内部自转周期，并传播测量误差。

进行随机化检验：以相同数量、相同 l 范围的连续轨道样本建立最近格点残差分布，并校正扫描多个分支的“别处效应”。

使用新增对象盲检：冥王星的新卫星、外太阳系新TNO以及未来更精确轨道应在不修改K和Ntr的条件下接受检验。

寻找动力学机制：必须给出自转如何产生或选择 h₀=NtrR²ω 的方程，而不能只停留在数值对应。

结语

基础角动量 h₀ 是宏观角动量量子化模型的核心。它把轨道几何 p=a(1−e²)、开普勒动力学 h²=GMp 与中心天体自转 R²ω 连接起来。转换量子数 Ntr 的引入，使这一结构从“每个系统都有一个拟合K”提升为“中心自转通过整数倍率生成轨道角动量单位”的可检验命题。

当前候选图景为：地球3；木星、土星、海王星1；冥王星45；太阳675。这组整数跨越人造卫星、巨行星卫星、双矮行星系统和整个太阳系，具有明确的统一美感。但它的科学价值最终不取决于整数有多漂亮，而取决于规则是否预先冻结、新数据能否通过盲检、统计显著性是否经得住选择效应，以及是否能建立产生这些整数的动力学机制。

h_l = Ntr l R²ω

如果这一关系经得住上述检验，那么基础角动量将不再是宏观轨道拟合中的辅助参数，而可能成为描述自转束缚系统的一种新的系统级不变量。

数据来源与说明

JPL Solar System Dynamics：Planetary Satellite Mean Elements，https://ssd.jpl.nasa.gov/sats/elem/sep.html

JPL Solar System Dynamics：Planetary Physical Parameters，https://ssd.jpl.nasa.gov/planets/phys_par.html

JPL Solar System Dynamics：Planetary Satellite Physical Parameters，https://ssd.jpl.nasa.gov/sats/phys_par/

JPL Small-Body Database：小天体轨道根数与历元说明，https://ssd.jpl.nasa.gov/tools/sbdb_lookup.html

Mankovich et al.：Cassini Ring Seismology as a Probe of Saturn’s Interior I: Rigid Rotation，ApJ 871 (2019).

IAU 2015 Resolution B3：名义太阳半径 R☉=695700 km。

支持单位：