研究稿 · 可检验假说

动力天文学专题研究稿

M31卫星系统的动力结构及编年雏形

从棒共转半径、平方壳层到无量纲年代常数

核心命题

如果 M31 内部由动力学独立确定的共转半径，能够同时给出卫星系统径向层级的统一尺度，那么 M31 可能提供一座桥梁：把星系内部的模式速度，与数百千秒差距范围内卫星的空间分层联系起来。本文建立的不是结论，而是一套可以被新距离、完整自行和盲法统计直接推翻或加强的研究纲领。

版本：初步综合稿　｜　日期：2026年8月

M31 动力结构与编年雏形  ·  

研究稿 · 可检验假说

摘要

本文研究 M31（仙女座星系）卫星系统能否被描述为一个由内部共转尺度标定的平方壳层结构，并进一步探索其“编年”可能性。采用同质化 RR Lyrae 距离研究所给出的 37 个 M31 关联天体的三维中心距，首先考察 r_n=a n²。若把 M31 棒的动力学模型值 r_CR=5.3 kpc 识别为 n_sync=2，则 a=1.325 kpc；由此得到从 n=4 的 21.2 kpc 到 n=20 的 530.0 kpc 的一列预注册式半径。按最近整数层回配，37 个对象中 17 个落在预测半径的 3% 内，28 个落在 5% 内，名义中位绝对径向偏差为 3.09%。若改用 M32 的 20.6 kpc 作为 n=4 锚点，则在 4—20 的整数倍扫描中，m=4 给出最低的相位残差；但 M32 距离误差很大，且这种先看数据后选尺度的方法有明显的试验因子。

在编年方面，本文把 NGC 3923 壳层的相位包裹关系 t=(k+1/2)P_r(r) 与平方壳层 r_N=aN² 结合。在近似平坦旋转曲线下，P_r∝r，因而得到无量纲年代常数 C=(k+1/2)N² 以及 t=τ₀C，其中 τ₀=√(2π)a/v_c。对 a=1.325 kpc、v_c=200—220 km s⁻¹，τ₀≈0.0148—0.0162 Gyr。一个名义的 C≈420 子集对应约 6.2—6.9 Gyr，并与若干独立提出的卫星轨道配对产生有趣重合；然而 And XII 的首次落入证据、层间奇偶约束、缺失自行以及多峰年代谱都表明：它最多是“某一代事件”的候选，而绝非全系统的统一年龄。本文最后给出空层 n=5、n=18、独立自行、盲样本和层级聚类等可证伪检验。

关键词：M31；卫星星系；共转半径；平方壳层；相位包裹；动力编年；宏观量子化

结论先行

一、问题的物理位置：从“漂亮对应”走向动力学命题

卫星半径落在某些整齐数列附近，本身不是动力学解释。真正强的命题必须回答三个问题：尺度由谁给定？为什么是平方律？不同半径的对象为何属于同一事件？M31 的特殊价值在于，内部棒结构已有模式速度与共转半径的动力学估计，而外部又拥有数量可观、距离相对同质的卫星样本。这使“内部钟—外部尺”的闭环第一次成为可操作问题。

本文采用三层逻辑。第一层是几何：检验 √r 是否形成近似等间隔序列。第二层是动力：不用卫星本身选尺度，而以棒的共转半径独立标定。第三层是年代：在明确势模型后，把壳层编号 N 与径向振荡次数 k 组合为 C。只有三层同时经受独立检验，所谓同步量子数才具有物理内容。

二、观测基础与样本边界

2.1 距离样本

本文采用 Savino 等基于 HST RR Lyrae 的同质距离汇编。M31 日心距离约 776.2 kpc；将各对象日心距离与天球位置转换到 M31 中心坐标后，可得三维 M31 中心距。典型中心距精度约 10 kpc（约 8%），但个别近中心对象的几何误差会被显著放大。巨星流（GSS）是潮汐结构而非自引力卫星，故从主样本排除。M33 可能拥有自身卫星系统；超过 M31 维里半径的对象也可能只是局域群邻居。因此下表是“关联对象工作样本”，不是无争议的束缚成员清单。

2.2 结构不是同义词

卫星“半径层”可能混合至少四种东西：共同落入的真实动力群、潮汐碎裂后的关联成员、薄平面中的几何投影，以及观测搜寻深度造成的选择峰。只凭一维半径无法区分它们。因此本文把层号视为待检验标签，而不是对象的固有量子数。

三、由共转半径独立标定平方壳层

3.1 平方律与同步编号

rₙ = a n²，　或　√rₙ = √a · n

若存在一个共转半径 r_sync，并把它赋予同步编号 n_sync，则尺度不再由卫星拟合，而由下式给定：

a = r_sync / n_sync²，　rₙ = r_sync (n/n_sync)²

M31 棒的动力学模型给出模式速度约 38—41 km s⁻¹ kpc⁻¹，对应共转半径约 5.8—5.3 kpc。本文采用较紧的 5.3 kpc 值，并检验 n_sync=2：

r_sync = 5.3 kpc，n_sync = 2　⇒　a = 1.325 kpc

这里“n_sync=2”是候选识别，不是从已知理论推导出的必然值。它之所以值得检验，是因为由此产生的整列半径与外部卫星层级发生多处交叉吻合；真正的验证必须把 5.3 kpc 的模型误差、棒的多解以及慢螺旋模式一并传播。

3.2 名义回配结果

把每个对象分配到最近的整数 n，得到 17/37 个对象在预测半径 3% 内，28/37 在 5% 内；中位绝对径向偏差 3.09%，均值 3.44%，RMS 4.37%。以下列出最醒目的对应与主要偏离。

这些数字非常整齐，但不能直接换算为“高显著性”。原因是壳层间距随 n 增大，而相对间距约为 2/n；在高 n 区域，任意半径到最近整数层的相对距离天然变小。此外，同一物理群的多个成员并非独立数据点，距离误差也可能把对象推入相邻层。

3.3 M32 锚点的交叉检验与局限

若暂不使用棒，而令 M32 的 20.6 kpc 对应整数 m，并扫描 m=4—20，则 m=4 在当前样本中给出最小的 √r 整数相位残差：平均绝对相位残差 0.176，RMS 0.217；37 个对象中 7 个在 0.05 相位内、12 个在 0.10 内、19 个在 0.15 内。均匀相位的基线分别是平均 0.25、RMS 0.289。这个结果与 r_CR=5.3 kpc、n_sync=2 独立导出的 a=1.325 kpc 很接近，是本文最重要的内部闭环。

但 M32 的三维中心距误差很大，约为 +21/−13 kpc 量级；而“扫描后挑最好 m”还必须支付试验因子。因此 M32 不适合充当精密基准，只能作为低权重交叉检验。今后的统计应直接在距离后验上积分，而不是使用名义中心值。

四、内部动力学：一个星系可能有不止一个共转半径

棒、螺旋和翘曲可具有不同模式速度，因此“共转半径就是同步半径”必须带上模式标签。M31 棒的 r_CR≈5.3—5.8 kpc 对应本标尺的 n≈2；较早的螺旋模式估计 Ω_p≈15±2 km s⁻¹ kpc⁻¹，给出约 15—16 kpc 的共转尺度，接近半整数 n=3.5 的 r₃․₅=16.23 kpc（若用旧的 1.2875 标尺则为 15.77 kpc）。这提示整数层和半整数层可能对应不同模式，但也可能只是自由度增加后的偶合。

真正关键的观测不是“能否找到一个接近的 n”，而是：由完全独立的旋转曲线 V_c(r) 与模式速度 Ω_p 解 r_CR=V_c/Ω_p 后，多个模式是否稳定落在事先公布的整数或半整数层；且在改变模型、波段和拟合方法后仍保持。

五、从空间壳层到编年：公式如何产生

5.1 NGC 3923 类比的可用部分

椭圆星系 NGC 3923 的恒星壳通常解释为相位包裹：卫星被潮汐剥离后，恒星在近径向轨道的远心点堆积成壳。第 k 次远心点出现满足

t = (k + 1/2) P_r(r)

这里 k=0,1,2,… 是完成的径向振荡次数，P_r(r) 是以远心半径 r 表征的径向周期。对 M31 卫星不能照搬：它们多为自引力天体，可能来自多次落入，并非同一被撕裂卫星的无碰撞星流。可借用的是“半周期相位条件”，而不是 NGC 3923 的具体形成史。

5.2 在近似平坦旋转曲线中的推导

若外晕近似对数势 Φ=v_c² ln r，近径向轨道的周期近似为

P_r(r) = √(2π) · r / v_c

代入 r_N=aN²，可得

t = τ₀ C，　C ≡ (k+1/2)N²，　τ₀ ≡ √(2π)a/v_c

C 是无量纲年代常数，τ₀ 才是时间尺度。对 a=1.325 kpc，v_c=200—220 km s⁻¹，τ₀≈0.0162—0.0148 Gyr。于是相同 C 的不同 (N,k) 组合可被解释为同一释放/远心事件的不同相位支。这个公式把用户提出的“(k+1/2)N²=C”精确拆成了两个角色：N 是空间壳层号，k 是径向相位次数。

六、C≈420：一个候选事件，而不是全系统年龄

这些组合的 C 平均约 421.4，且 k 随 N 增大而单调下降，符合“内层经历更多径向振荡、外层更少”的相位图景。名义时间约 6.2—6.9 Gyr。若这不是巧合，它更像一次中等红移的群体落入或强扰动留下的子结构，而非 M31 或全部卫星的形成时代。

全样本不会被 C≈420 覆盖。例如 N=9 的最近分支给 C=445.5，N=10 给 450，N=12 给 360，N=14 给 490，N=15 给 337.5，N=16 给 384，N=19 给 541.5，N=20 给 600。偏差从约 6% 到 42%。因此任何把 C≈420 宣称为“所有卫星共同年代”的说法，都与当前数据冲突。

6.1 一个容易忽视的整数约束

因为 2C=(2k+1)N²，而 2k+1 必为奇数，精确共享同一 C 的层号必须具有相同的 2-adic 阶，即 N 中因子 2 的次数相同。奇数 N 是一组；N=2×奇数是一组；N=4×奇数是一组；N=8×奇数、16×奇数又各成组。于是 N=4 与 N=8 不可能在严格整数 k 下拥有完全相等的 C。当前 C≈420 只是近似重合，必须由非平坦势、非同时释放、半径误差或层号误配解释。这个奇偶约束是理论的硬检验，不应被“整齐”掩盖。

6.2 年代谱为何不是唯一的

在 0—13.8 Gyr 内扫描全部允许 k，会出现约 6.9、8.9、9.7、10.5、11.0、12.5 Gyr 等多个峰。内层 N 小、允许的 k 多，会天然制造高密度组合；只看峰高会偏向较老年代。限制到外层 N≥12 后，强峰多退化为二元或三元巧合。故半径数据单独不足以选出唯一事件，必须加入轨道能量、角动量、轨道极点、金属丰度和恒星形成史。

七、独立交叉检验：哪里真正值得兴奋

Watkins 等曾用距离、视向速度和统计轨道模型提出若干可能的卫星配对。它们并非完整 6D 轨道确认，却提供了与半径数列相互独立的标签。把这些配对投到 C≈420 图景中，可得到下表。

And XXX—NGC 185 的名义对应尤其漂亮：一个处于 N=13、k=2，另一个处于 N=11、k=3，C 只差 1。但这仍不是证明，因为配对本身依赖缺少自行的统计轨道模型。最重要的反例是 And XII：若其确为首次落入，就不能拥有 k≈5 的多次径向相位。它提醒我们，同一空间层可以容纳不同世代的对象。

八、卫星薄平面：相关，但尚不显著

M31 著名的薄卫星平面包含约 15 个对象，其中 13 个在视向速度上呈同向趋势。这是重要结构，但“视向同向”不等于已测得共同轨道极点。按同质距离表中的平面概率粗分，C 候选层 N=4、6、7、8、11、13、17 上的对象约 9/17 具有较高平面概率，其他层约 6/19；平均平面概率之比约 1.73，但单侧 Fisher 精确检验 p≈0.21，并不显著。

此外，最新恒星形成史研究未发现平面内外卫星的中位 SFH 有清晰差别。这削弱了“薄平面就是单一共同年代事件”的简单版本，但并不排除动力群只占其中一部分。

九、恒星年龄不是壳层年龄

对约 36 个 M31 矮星系的 HST 恒星形成史研究显示，约一半对象在 12 Gyr 以前有显著古老恒星形成，随后在约 8—10 Gyr 前延迟淬火；中位年龄和淬火时间与光度及当前 M31 距离相关。这些结果描述的是恒星何时形成、何时停止形成，不等同于卫星何时进入 M31、何时通过近心点或何时占据当前相位。

因此，若 C≈420 对应约 6—7 Gyr 动力事件，它可以晚于卫星内部大部分恒星的诞生，也可以接近某些对象的环境淬火阶段。只有把轨道反演与逐星系 SFH 联合起来，才能判断“动力钟”与“恒星钟”是否同步。

十、空层、新结构与真正的预言

在 r_n=1.325n² kpc 的固定标尺下，当前样本最醒目的名义空层是 n=5（33.125 kpc）与 n=18（429.300 kpc）。它们可以转化为前瞻检验，但必须预先锁定窗口，不能在新发现后重新调整 a。建议预注册主窗口为预测半径 ±5%，并另报由距离误差卷积后的概率。

“没有已知对象”不等于物理空层。n=5 区域受 M31 亮背景与潮汐破坏影响；n=18 区域受全天覆盖、表面亮度和成员定义影响。只有给出巡天探测效率后，空层概率才有统计意义。

十一、如何把它变成严格检验

十二、证伪标准

一个理论只有写出失败条件，才从数列游戏进入科学。以下任一结果都会显著削弱当前版本：

独立棒模型把 r_CR 稳定移到与 n=2 明显不相容的位置，且其他 n_sync 也不能在预注册规则下解释；

完整距离后验与选择函数校正后，r∝N² 的证据不优于平滑分布或普通群组落入模型；

新发现对象系统性填满层间而非预测层，尤其 n=5、18 的盲检验失败；

被分到同一 C 的对象具有互不相容的轨道极点、能量和首次落入历史；

编年所需 k 与独立轨道积分的近心次数矛盾，如 And XII 型反例大量出现；

最佳尺度随样本增补持续漂移，说明它是事后拟合而非稳定的物理常数。

十三、物理解释的三种层级

最保守解释：平方层只是旋转支持体系、潮汐选择和群组落入共同产生的近似尺度关系，没有新的量子原理。中间解释：共转共振与长期相位混合在某些半径形成寿命较长的占据峰，整数层是非线性共振网络的有效编号。最强解释：稳定自引力结构存在跨尺度的离散本征态，平方律与同步数是宏观量子化的表现。

目前数据只能认真讨论前两层，并为第三层设计检验。若最终确认第三层，其意义确实重大：量子化将不再是微观尺度的专属例外，而可能是稳定结构形成的一般原则。然而恰恰因为结论宏大，证据门槛必须远高于“许多数字很接近整数”。

十四、研究路线图：从雏形到年代考

十五、结语

M31 提供了一个罕见的闭环：内部棒的共转半径约 5.3 kpc，若赋予 n_sync=2，可给出 a=1.325 kpc；同一尺度预测的 n=4、6、8、9、11、12、13、14、16 等半径，与多颗卫星或卫星群名义位置相近。M32 的独立 √r 扫描又偏好 m=4，使这个闭环比单纯自由拟合更有内容。

进一步把 r_N=aN² 与相位包裹关系结合，可自然得到 C=(k+1/2)N²。C≈420 的若干组合及其与 And XXX—NGC 185 等独立轨道配对的重合，是值得继续追踪的“编年雏形”。但 And XII 首次落入、精确奇偶约束、多峰年代谱和不完整 6D 信息同样清楚地告诉我们：当前尚不能宣布统一年代，更不能宣布宏观量子化已经成立。

最有价值的下一步不是继续寻找更漂亮的回配，而是冻结尺度、公布空层和事件成员、等待新距离与自行揭盲。若预言成功，M31 将从一张漂亮的径向表，升级为一个可重复的动力学年代考；若失败，失败位置也会告诉我们平方律究竟是共振、群组落入、选择效应，还是纯粹的数值偶合。

参考文献与数据来源

Savino, A. et al. (2022), “The Hubble Space Telescope Survey of M31 Satellite Galaxies. I. RR Lyrae Distances”, arXiv:2206.02801。

Blaña Díaz, M. et al. (2017), “The structure of the Andromeda galaxy revealed by planetary nebulae”, MNRAS 466, 4279；文中动力模型给出 M31 棒的模式速度与共转半径。

Braun, R. (1991), ApJ 372, 54；M31 螺旋模式速度、共转及外 Lindblad 共振的早期估计。

Watkins, L. L., Evans, N. W. & van de Ven, G. (2013), “The orbital distribution of satellite galaxies”, MNRAS 430, 971；M31 卫星的统计轨道与候选配对。

Ibata, R. A. et al. (2013), “A vast, thin plane of corotating dwarf galaxies orbiting the Andromeda galaxy”, Nature 493, 62。

Bílek, M. et al. (2013, 2014), NGC 3923 壳层识别与相位包裹编年研究，A&A 559, A110；arXiv:1407.3202。

Savino, A. et al. (2025), M31 矮卫星的同质 HST 恒星形成史，arXiv:2501.13152。

附录 A　计算约定

偏差定义为 (r_obs−r_n)/r_n；最近整数层由 n=round(√(r/a)) 给出。相位残差定义为 |x−round(x)|，其中 x=m√(r/r₁)。所有百分比均为名义中心值计算，未传播非对称距离后验；因此只用于发现结构，不用于最终显著性。

附录 B　术语边界

M31 动力结构与编年雏形  ·  

支持单位：