星系旋转动力学方程(GRoDE)

——平坦旋转作为尺度动力系统的稳定固定点

摘要

本文提出将星系旋转曲线重新表述为一个一阶尺度动力学问题。核心候选方程为 r dV/dr = V(1−V²/A)，其中 r 为银心距离，V(r) 为圆周速度，A 为特征速度平方。该方程在无量纲变量 u=V/√A、ξ=ln r 下化为 du/dξ=u(1−u²)，因此 V=√A 是稳定固定点。其积分解给出 V(r)=√A·r/√(r²+c²)，而常见的对数势及其加速度形式 α=Ar/(r²+c²) 则作为该尺度动力学方程的积分结果出现。文章的重点不在重新发现某个经验加速度公式，而在于把平坦旋转重新解释为尺度动力系统的稳定态，并提出可直接检验的局部关系 dlnV/dlnr+V²/A=1。前期多颗星系的试验性检验显示，该关系值得进入预先固定公式、样本与质量标准的大样本确认阶段。

一、问题：能否直接写出旋转速度的尺度动力学方程？

传统星系动力学通常从质量分布或引力势出发，再求加速度与旋转曲线。这里采取相反的出发点：不先假定某个复杂质量模型，而直接询问旋转速度 V(r) 随尺度 r 的演化是否满足一个简单的一阶方程。

候选关系为：

r · dV/dr = V(1 − V²/A)        (1)

式（1）才是本文的核心。它不是某条旋转曲线的事后拟合形式，而是一个关于“速度如何随尺度演化”的微分方程。

二、平坦速度是稳定固定点

将式（1）写成对数尺度变量，可得：

dV/dln r = V(1 − V²/A)        (2)

若 V²<A，则右端为正，速度随半径继续增加；若 V²>A，则右端为负，速度随半径回落；当 V²=A 时，右端为零。因此：

V_f = √A        (3)

是方程的固定点。更重要的是，它是稳定固定点。于是“平坦旋转”不再只是经验上看到的一段平台，而可被解释为尺度动力系统自然趋近的稳定态。

三、无量纲化：一个极简的尺度恢复方程

定义无量纲变量：

u = V/√A,     ξ = ln r

式（1）化为：

du/dξ = u(1 − u²)        (4)

再令 y=u²=V²/A，则：

dy/dξ = 2y(1 − y)        (5)

这是一类 Logistic 型尺度恢复方程。对物理上的正速度分支，y=0 为不稳定点，y=1 为稳定点。于是：

V²/A → 1

不是人为规定，而是方程自身的稳定性结果。

四、积分解：旋转曲线从动力学方程中出现

积分式（5）得到：

ln[y/(1−y)] = 2ln r + C

把积分常数重新写成尺度 c，可得：

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

因此：

V(r) = √A · r/√(r²+c²)        (6)

这里 A 属于动力方程本身，而 c 是积分常数。不同星系可以有不同的 c，却仍服从同一个尺度动力学方程。这一区分比把 A、c 都视为经验拟合参数更重要。

五、经典对数势只是积分结果，而不是本文重心

由圆周运动关系 α=V²/r，可由式（6）推出：

α(r) = Ar/(r²+c²)        (7)

这一加速度形式与经典对数势在赤道面上的结果等价，因此它本身不应被表述为新的发现。本文真正要强调的是逻辑顺序的改变：

尺度动力学方程  →  稳定固定点  →  旋转曲线  →  加速度场

也就是说，重点不再是“找一个能产生平坦旋转的势”，而是“是否存在一个更基础的尺度演化规律，使平坦速度成为稳定态”。

六、可直接证伪的局部恒等式

由式（6）可以消去积分常数 c，得到：

dlnV/dlnr + V²/A = 1        (8)

定义：

Q(r) = dlnV/dlnr + V²/A

如果式（1）成立，则必须有：

Q(r) = 1        (9)

这比单纯看旋转曲线拟合残差更严格，因为它同时涉及速度本身和局部速度梯度。式（9）因此是本文提出的主要观测判据。

七、初步观测检验

我们已经使用公开 SPARC 旋转曲线，对式（8）、式（9）进行了第一轮实质性检验。所分析对象包括 NGC 3198、NGC 6503、DDO 154、NGC 1003、IC 2574、NGC 2403、NGC 3521、NGC 7331、NGC 5055、NGC 2903、NGC 2841 等不同质量和形态的星系。

二参数主骨架在多颗规则或近规则星系中可将整体旋转曲线描述到几个百分点量级的归一化残差。更重要的是，局部量 Q(r) 在若干规则星系中的平均值接近 1。前期计算的代表性结果如下。

在这组初步样本中，⟨Q⟩ 的中位数约为 0.998，而局部散布的中位量级约为 0.17。需要强调：A 是从同一条旋转曲线中估计得到的，因此“平均 Q 接近 1”不能单独视作独立证明；真正有判别力的是局部散布、不同星系之间的重复性，以及预先规定分析方案后的大样本表现。

八、与“本征运动 + 扰动”框架的连接

如果把式（1）的积分解视为平滑本征速度场，那么实际观测可以写成：

V_obs(r) = V_0(r) + δV(r)

对应加速度可分解为：

a_obs(r) = α_0(r) + Δα(r)

其中 α_0=Ar/(r²+c²) 是由尺度动力学方程导出的本征骨架，Δα 则表征局部结构对该骨架的偏离。这样，“骨架”和“扰动”可以在数据中被分别研究，而不必让一个越来越复杂的经验函数把所有局部结构全部吞掉。

九、这不是时间动力学，而是尺度动力学

必须特别说明：式（1）的自变量是 r，而不是时间 t。把 ξ=ln r 作为“演化变量”时，讨论的是尺度流或径向流，不是通常意义上的时间演化。因此，“稳定固定点”“恢复”这些术语在本文中指的是随尺度的稳定结构，而非粒子随时间被某种新力直接驱动。

这一限定很重要，因为它避免把一个关于旋转曲线形状的尺度方程误写成新的时间运动定律。

十、确认阶段：固定公式，不再修改

下一阶段最重要的原则不是继续加三阶、四阶修正，而是冻结核心方程：

r · dV/dr = V(1 − V²/A)

然后预先规定样本选择、质量标准、导数估计方法、A 的估计方法、异常结构分类以及通过/失败判据，对几十乃至上百颗高质量星系进行统一检验。真正需要回答的是以下三个问题：

不同星系在 V/√A 与 r/c 归一化后，是否塌缩到同一条母曲线？

Q(r)=dlnV/dlnr+V²/A 是否在规则星系中系统接近 1？

有棒、warp、非圆流等受扰系统的 Q 散布，是否显著大于规则星系？

只有在预先固定规则的大样本检验中继续成立，式（1）才有资格从“漂亮的尺度动力学候选关系”走向具有普适性的经验定律。

十一、讨论：真正可能的新意在哪里？

经典文献早已知道对数势可以产生 V(r)=V_f r/√(r²+c²) 这一类旋转曲线。因此，本文不把对应的势、速度函数或加速度函数本身作为原创性的核心。真正值得研究的是：能否把同一结构重写成自治尺度动力系统，并把平坦速度理解为稳定固定点；能否用 Q=1 这样的局部关系进行统一观测检验；以及局部偏离是否能独立追踪棒、旋臂和非圆运动等结构。

如果这些问题获得肯定答案，那么贡献将不在于“又一个旋转曲线公式”，而在于给出一种不同的组织方式：把星系旋转从势模型的结果，提升为可直接检验的尺度动力学关系。

另外需声明的是：

传统路径是

Φ→α→ V

而我的路径是

{旋转动力学方程}→ V→α→Φ

如果后续大样本验证继续支持这个动力学方程，那么此文最有分量的一句话可能就是：

{经典对数势或许不是一个任意构造，而是更基本的尺度动力学的积分结果。 }

十二、结论

本文提出把星系旋转曲线写成以下尺度动力学方程：

r · dV/dr = V(1 − V²/A)

它具有三个直接结果：

V=√A 是稳定固定点，因此平坦旋转可被理解为稳定尺度态；

积分后得到 V(r)=√A·r/√(r²+c²)，其中 c 是积分常数；

导出可直接检验的局部恒等式 dlnV/dlnr+V²/A=1。

因此，今后的关键不再是继续装饰公式，而是让大样本旋转曲线数据直接决定它的命运。

参考背景

1. Binney, J. (1981). 关于经典 logarithmic potential（对数势）的早期工作。

2. Binney, J. & Tremaine, S. Galactic Dynamics. Princeton University Press.（经典星系动力学教材，对对数势及平坦旋转曲线有系统讨论。）

3. Lelli, F., McGaugh, S. S., & Schombert, J. M. (2016). SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves. The Astronomical Journal, 152, 157.

注：本文中的“大样本确认”仍属下一阶段工作；文中前期检验属于探索性分析，不应被表述为已经完成的普适性证明。

星系动力学公式

支持单位：