GRODE在仙女星系的惊艳表现

——一次从星系盘走向百千秒差距外晕的跨尺度检验记录

研究纪实稿

记录版本：2026年10月

仙女星系 M31 是我们检验 GRODE（Galaxy Rotation Dynamical Equation，星转方程）过程中一个具有特殊意义的对象。

真正令人印象深刻的，并不是 GRODE 能够在某一段旋转曲线上得到一个漂亮拟合。任何具有两个或三个自由参数的模型，都可能在有限区间内取得不错的拟合结果。

M31 给我们信心的地方在于另一件事：

内区确定 GRODE 参数 (A, c)  →  冻结参数  →  向几十至上百 kpc 的外晕推进

而 GRODE 在这个过程中没有迅速崩溃。这正是此次 M31 检验最值得记录的地方。

一、GRODE 的基本形式

GRODE的基本形式：

rdV/dr=V(1-V²/A)

经过积分，可采用的暗晕旋转速度形式为

V_h²(r) = A r² / (r² + c²)

V_h(r) = √A · r / √(r² + c²)

其中 A 规定暗晕远区的速度平方尺度，c 规定从内区向平坦区过渡的特征半径。

当 r ≪ c 时，

V_h ≈ (√A / c) r

表现为近似线性上升；当 r ≫ c 时，

V_h → √A

因此远区自然得到

M_h(<r) ≈ (A/G) r

ρ_h(r) ∝ r⁻²

也就是说，GRODE 的真正考验不只是能不能拟合星系盘，而是：由内区得到的 A、c 能不能继续支撑真正的外晕动力学。M31 恰好为这个问题提供了罕见的实验场。

二、为什么选择仙女星系

M31 是距离我们最近的大型外星系之一。它与银河系相比有一个巨大优势：我们站在它的外面。

因此可以同时观察星系盘旋转、HI 气体、恒星晕、球状星团、巨型恒星流以及卫星星系。尤其重要的是，M31 的外晕球状星团已经被观测到 R_proj 约 140 kpc，部分对象的三维距离甚至接近约 200 kpc。

已有外晕球状星团研究显示，30 kpc 以外的球状星团仍具有系统性的运动学结构，其速度弥散随半径增加而下降；相关动力学分析给出的 M31 在 200 kpc 内质量约为 1.2–1.6×10¹² M☉，同时存在模型与系统误差。

这使 M31 成为一个难得的实验：同一个动力学模型可以从星系盘一直推到百 kpc 外晕。

三、第一次真正让人吃惊的参数

在这一轮 M31 检验中，我们采用约 6 < r < 20 kpc 作为主要训练区。在重子贡献扣除之后，得到的一组候选 GRODE 参数为

√A = 149.2 km s⁻¹

c = 1.99 kpc

因此

A = (149.2)² ≈ 2.23×10⁴ (km s⁻¹)²

这里真正重要的，不是这两个数字本身，而是从这一刻开始，我们采取了一条严格规则：A、c 冻结。

不再因为外晕数据不好看而重新调整参数；不在 50 kpc 调一次，也不在 100 kpc 调一次，更不能到了 140 kpc 又换一组参数。接下来只问一个问题：这一组在约 6–20 kpc 得到的参数，还能走多远？

四、从 20 kpc 一步跨向百 kpc 外晕

训练区的上限约为 20 kpc，而真正的外晕球状星团资料可以一直走到 100–140 kpc。也就是说，我们不是在做 20→25 kpc 这样的小外推，而是在尝试

20 kpc  →  100–140 kpc

尺度扩大了约 5–7 倍。如果从训练区最内侧约 6 kpc 算起，整个检验实际上覆盖了二十倍以上的径向尺度。若把盘内运动学、HI、外晕 GC 和卫星体系放在一起看，它已经接近我们一直追求的：同一组动力学参数跨越几十倍空间尺度。

五、真正令人惊艳的是球状星团速度弥散

外晕球状星团不是圆盘粒子，因此不能简单写成 V_GC = V_h。必须把 GRODE 给出的引力势放入 Jeans 方程，再投影得到视线速度弥散。

d(νσ_r²)/dr + 2βνσ_r²/r = -ν GM(<r)/r²

M31 外晕 GC 的观测事实非常明确：速度弥散随着半径增加显著下降，已有研究发现这种下降可近似描述为约 -0.5 的幂律。我们当时用冻结的 GRODE 参数作了一次粗 Jeans 外推，得到下列代表性结果。

这里必须特别指出：30 kpc 附近并不好。这是一个真实的压力点，而不是应该隐藏的数据。但是随着半径真正进入外晕，50→100→125→140 kpc，理论与观测反而越来越靠近。尤其约 100–140 kpc 范围内，77↔79、73↔70、64↔68 km/s 这样的数量级一致性，是这次 M31 检验中真正令人吃惊的部分。

说明：上述数值属于当时的 Jeans 粗预测和读图级比较，还不是对逐个 GC 做完整似然分析所得的正式结果。

六、为什么 100–140 kpc 的表现比 30 kpc 更重要

30 kpc 仍然处在一个非常复杂的过渡区域：恒星盘仍有影响；核球和盘质量模型的误差仍会传播；恒星流和非平衡结构很多；GC 的空间分布和各向异性也更复杂。

而到了 100–140 kpc，我们已经真正进入 M31 的远外晕。此时能够继续改变理论曲线的“人为自由度”反而越来越少，参数早已冻结。所以在这个区域如果理论仍然没有跑飞，意义远比再把内盘拟合得漂亮重要。

GRODE 可能真的抓住了暗晕的一种整体动力学结构。

七、质量剖面也没有立即崩溃

冻结参数以后，GRODE 给出的远区质量具有近线性增长：

M_h(<r) ≈ (A/G) r

我们当时得到的一组粗略总质量序列约为：

M(<50 kpc) ≈ 3.5×10¹¹ M☉

M(<100 kpc) ≈ 6.1×10¹¹ M☉

M(<125 kpc) ≈ 7.4×10¹¹ M☉

M(<200 kpc) ≈ 1.1×10¹² M☉

这组数字尤其值得与独立观测比较。巨型南部恒星流动力学曾给出 M(<125 kpc) ≈ 7.5×10¹¹ M☉，其不确定范围约为 +2.5/−1.3×10¹¹ M☉。较新的 M31 综合运动学研究利用上万天体，将旋转曲线构造到约 125 kpc：约 25 kpc 内维持在约 220 km/s，随后向外逐渐下降到约 170 km/s，并得到约 1.14×10¹² M☉ 的维里质量尺度。

因此，我们当时冻结 GRODE 后得到的 M(125)≈7.4×10¹¹ M☉ 以及 M(200)≈1.1×10¹² M☉，并没有随着尺度扩大而出现数量级灾难。

八、这不是普通意义上的“拟合成功”

如果只是在 6–20 kpc 范围里得到一个很好看的旋转曲线，我不会把 M31 称作 GRODE 的“惊艳表现”。真正令我印象深刻的是：用小尺度信息确定参数，然后把同一套动力学结构送入完全不同的示踪体系。

从星系盘走向 HI，再走向恒星晕和球状星团，最后走向卫星和百 kpc 外晕。这些对象的动力学性质完全不同：盘是旋转支撑，GC 和恒星晕主要是压力支撑，卫星又是稀疏的离散示踪体。

同一组 (A, c) 能不能跨示踪体继续工作？

M31 第一次让我们看到：这种可能性至少没有被数据轻易杀死。

九、最漂亮的地方：它原本很容易失败

GRODE 的远区极限非常刚性：

V_h² → A

M_h(<r) → (A/G) r

也就是说，一旦内区把 A 定下来，远区质量斜率几乎已经被确定。它不像某些经验密度模型那样，可以到了外区再加入新的尺度、浓度参数或截断半径不断调整。

因此 M31 实验本来极容易失败。如果冻结参数以后 M(100 kpc) 已经偏离两三倍，或者 σ_GC(100–140 kpc) 完全跑向错误方向，那么实验基本就可以结束。可是并没有发生这种事情。相反，真正进入 100–140 kpc 以后，结果仍处在正确数量级，而且相当接近观测。

十、仙女星系给了 GRODE 什么？

它没有证明 GRODE。但它完成了更加重要的一步：把 GRODE 从一个旋转曲线公式，推到了一个值得进行跨尺度动力学检验的候选规律。

一个旋转曲线公式可以在十几个 kpc 范围内漂亮；一个暗晕动力学规律则必须经得住 10、50、100、200 kpc 不断变化的尺度和示踪体。M31 是我们第一次强烈感受到：GRODE 似乎有可能走这么远。

十一、同时必须留下三个警告

第一，A、c 仍然是候选冻结参数。这里记录的 √A=149.2 km/s、c=1.99 kpc 来自我们当时的工作解，还需要用统一的 M31 重子模型、统一的数据集重新进行正式联合平差。

第二，外晕 GC 比较属于 Jeans 粗预测。尤其 30、50、100、125、140 kpc 的速度弥散比较，还没有完成逐星似然分析，也没有充分边缘化 β(r) 即速度各向异性的不确定性。因此应称为“强烈的初步跨尺度信号”，而不能称作最终统计证明。

第三，M31 外晕并非严格平衡。M31 外晕拥有大量恒星流和吸积遗迹，许多远外晕 GC 与恒星子结构相关，有些甚至属于动力学冷结构。因此任何 Jeans 分析都必须对非平衡成分保持警惕。

恰恰因为这些困难存在，GRODE 还能在如此长的尺度上维持正确的数量级，才更加值得继续追查。

十二、这一战真正留下来的东西

如果以后忘记了所有具体数字，只需要记住这一幅图景：

6–20 kpc：确定一组候选 GRODE 参数，然后封盘。

30 kpc：出现压力。没有重新拟合。继续。

50 kpc：差距缩小。继续。

100 kpc：理论与外晕 GC 已经非常接近。继续。

125 kpc：仍然没有倒。继续。

140 kpc：仍在正确数量级，而且相当接近观测。

这就是这次实验最令人难忘的地方：我们没有追着数据调整 GRODE，而是让被冻结的 GRODE 自己向外走。它竟然走出了星系盘，走过几十 kpc，一直走进仙女星系真正的百 kpc 外晕。

十三、结语：仙女星系给我的信心

科学中最值得珍惜的时刻，不一定是得到一个完美拟合的时候。有时候，是一个原本很容易死亡的想法，被扔进越来越苛刻的数据中以后，竟然还活着。M31 给 GRODE 的，正是这样一次经历。

它并没有让我们宣布胜利。相反，它让问题变得更加严肃：为什么一组在星系盘尺度确定的参数，能够继续描述百 kpc 外的动力学？

如果这只是巧合，那么更多星系会很快把它击碎。如果不是巧合，那么 GRODE 所表达的就可能不是某条旋转曲线的经验形状，而是暗晕内部某种更深的尺度动力学。

“GRODE 已被证明。”——不是这次实验应得出的结论。

仙女星系第一次让我们认真相信：GRODE 可能真的具有跨尺度的暗晕动力学意义。

这，就是 GRODE 在仙女星系中那次令人难忘的惊艳表现。

参考文献与数据线索

Veljanoski, J. et al. (2014), The outer halo globular cluster system of M31. Monthly Notices of the Royal Astronomical Society, 442, 2929–2950.

Fardal, M. A. et al. (2006/相关研究), M31 巨型南部恒星流动力学及质量约束。

2024 年 M31 综合运动学研究（上万天体样本），用于 125 kpc 级旋转曲线与维里质量比较。

本文中的 GRODE 参数、Jeans 粗预测和跨尺度比较为本研究过程中的工作解，尚需统一数据集上的正式联合平差与逐星似然分析。

支持单位：