——从一阶尺度动力学到平坦旋转、对数势与暗晕结构

论文初稿

摘要

星系暗物质研究通常从假定暗晕密度分布或引力势出发，再通过泊松方程和运动方程计算旋转曲线。本文提出一种相反的研究路线：首先建立暗晕旋转速度随径向尺度演化的一阶动力学方程——GRODE（Galactic Rotation Ordinary Differential Equation（星转方程）），再由动力学依次推出旋转速度、径向加速度、引力势、等效累计质量和等效密度。GRODE 写为 r(dV/dr)=V(1−V²/A)，其中 A 是方程唯一的动力学参数。给定边界条件后，方程积分产生尺度常数 c，得到 V²=A r²/(r²+c²)，继而推出 g=A r/(r²+c²) 与 Φ=(A/2)ln(r²+c²)+C。因此，广泛用于星系动力学研究的 cored logarithmic potential 在本框架中不是先验假设，而是动力学方程在给定边界条件后的积分结果。本文还总结了基于公开 SPARC 旋转曲线数据对 NGC 3198、NGC 2841、NGC 5055、NGC 2403、NGC 6503 和 NGC 2903 六个星系进行的统一初步检验，并讨论 Pota 等（2015）在 NGC 1407 上采用 cored logarithmic halo 的工作。六个样本在具有足够径向杠杆臂的情况下均显示 GRODE 函数形式具有良好描述能力；其中若干样本的 1/V_h²–1/r² 线性化关系表现出较高线性度。本文据此提出：对数势可能不是暗晕研究的起点，而是暗晕动力学的结果。

关键词： GRODE；暗晕动力学；星系旋转曲线；对数势；SPARC；NGC 1407

1. 引言

星系旋转曲线是暗物质问题最直接的动力学证据之一。在普通牛顿重子质量分布下，外盘圆周速度一般应随半径下降；然而大量螺旋星系表现出近似平坦的外区旋转曲线。

传统暗晕研究通常先选择密度轮廓或势函数，例如 NFW、Einasto、Burkert、pseudo-isothermal 或 logarithmic halo，然后利用观测数据确定其参数。这种路线在经验上十分成功，但留下一个更上游的问题：为什么暗晕动力学应当选择这种密度或这种势？

本文尝试倒转逻辑方向：先建立动力学，再依次得到速度、加速度、势、累计质量和等效密度。本文把这一一阶径向尺度动力学方程称为 GRODE。

2. GRODE：暗晕的一阶尺度动力学

GRODE 的基本形式为：

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

其中 r 为银心距离，V(r) 为暗晕动力学对应的圆周速度，A>0 为方程唯一的动力学参数，√A 具有速度量纲。方程的一个重要特点是半径只以 r(d/dr)=d/d(ln r) 的形式出现，因此它本质上描述的是径向尺度空间中的动力学流。

令 u=V/√A，x=ln r，则式（1）化为：

du/dx = u (1 − u²)    (2)

无量纲化以后，方程中不再显含任何参数。因此，所有满足 GRODE 的体系共享同一条无量纲尺度流；不同星系之间的差别进入速度尺度和边界条件，而不是进入无量纲方程本身。

3. 稳定点与平坦旋转速度

方程 du/dx=u(1−u²) 存在三个不动点：u=0、u=+1、u=−1。对正常正向旋转只考虑 u>0。当 0<u<1 时，du/dx>0，因此速度随尺度增加；在 u=1 处，du/dx=0，且该状态为稳定点。

V(r) → √A    (3)

因此，平坦旋转速度在 GRODE 中不是人为附加的边界条件，而是动力系统的稳定状态。另一方面，在内区 V²≪A 时，式（1）近似为 r(dV/dr)≈V，从而得到：

V ∝ r    (4)

同一个一阶动力学方程因此同时生成“内区近似线性上升—外区趋于平坦”的连续过渡。

4. GRODE 的解析解与边界尺度

令 y=V²，则式（1）可写为 r(dy/dr)=2y(1−y/A)。分离变量并积分，可得 ln[y/(A−y)]=2 ln r+C₁。把积分常数重新写成一个长度尺度 c，得到：

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

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

这里必须区分方程参数与解的积分常数。A 属于 GRODE 方程本身，而 c 是积分以后、给定边界条件才被确定的尺度。若给定 V(r₀)=V₀，则有：

c² = r₀² (A/V₀² − 1)    (7)

因此，从方程层级看，GRODE 是一个“单动力学参数 + 边界条件”的体系，而不是把 A 与 c 并列作为两个基本动力学参数。

5. 从 GRODE 到径向加速度

圆周运动满足 V²=rg。代入式（5），得到：

g(r) = A · r/(r² + c²)    (8)

r≫c 时：g(r) ≃ A/r    (9)

在内区 r≪c 时，g(r)≈(A/c²)r；在远区 r≫c 时，g(r)≈A/r。于是 GRODE 产生从中心近线性加速度到外区 1/r 加速度的连续过渡。

6. 对数势不是假设，而是 GRODE 的积分结果

由 g=dΦ/dr，将式（8）积分得到：

Φ(r) = (A/2) ln(r² + c²) + C    (10)

这正是经典的 cored logarithmic potential。数学上，从已知对数势当然可以反推出同一速度解并重新写成 GRODE；但理论层级不同。在传统 LOG 模型中，势函数通常被直接选作参数化模型，而在 GRODE 中，对数势是动力学方程加边界条件后的积分结果。

7. GRODE 的等效暗晕质量

若用牛顿球对称质量表示同一径向加速度，g=GM_h(<r)/r²，则：

M_h(<r) = (A/G) · r³/(r² + c²)    (11)

r≫c 时：M_h(<r) ≃ (A/G) r    (12)

因此，GRODE 自动产生外晕近线性累计质量律。这里的 M_h 可首先理解为与 GRODE 动力学场对应的牛顿等效质量；它是否等同于某种具体粒子暗物质的真实累计质量，需要额外物理假设与观测判据。

8. GRODE 的等效密度

由球对称关系 ρ=(1/4πr²)dM/dr，得到：

ρ_h(r) = [A/(4πG)] · (r² + 3c²)/(r² + c²)²    (13)

r≫c 时：ρ_h(r) ≃ A/(4πG r²)    (14)

于是同一条动力学方程同时生成 V→常数、g∝r⁻¹、M(<r)∝r 和 ρ∝r⁻²。它们不是彼此独立的四个假设，而是 GRODE 的四种等价表现。

9. 一个直接的线性化观测判据

由式（5）取倒数：

1/V_h² = 1/A + (c²/A) · 1/r²    (15)

定义 X=1/r²、Y=1/V_h²，则 GRODE 预言 Y=b+mX，其中 b=1/A，m=c²/A。因此：

A = 1/b，    c = √(m/b)    (16)

这提供了一个简单而强的诊断：如果重子贡献被合理扣除，那么同一颗星系不同半径的数据应近似落在同一条 1/V_h²–1/r² 直线上。该线性化关系适合作为辅助检验，但在重子主导的内区，V_h² 是两个较大数之差，误差会被取倒数操作放大，因此主检验仍应在原始速度空间中进行，并严格传播观测与重子质量误差。

10. 重子质量与盘系数 γ_d

实际观测测得的是总圆周速度。对无明显核球的盘星系，可写为：

V_obs² = γ_d V_d² + V_g² + A · r²/(r² + c²)    (17)

其中 V_d 是由盘光度分布计算得到的参考恒星贡献，γ_d 表示盘质量归一化或等效质光比缩放。星系光度分布可以较好观测，但恒星质量仍依赖恒星年龄、金属丰度、初始质量函数与星族历史，因此 γ_d 不宜对所有星系机械固定为同一个常数。

在本文六个星系的阶段性试算中，γ_d、A 与 c 被联合估计，而不强制规定 γ_d=0.5。这样可以避免把某个预设的恒星质量归一化直接传递到暗晕参数 c 的估计中。

11. 六个独立星系的初步统一检验

本文按照同一 GRODE 函数形式，对六个径向覆盖较长的公开旋转曲线样本进行了阶段性统一试算。筛选特别关注 r_max/c>5，因为只有观测半径真正跨越数倍 c，才有机会同时看到 r~c 的转折区与 r≫c 的渐近区。

上述参数与线性度属于当前阶段基于公开旋转曲线进行的统一试算结果，不应被误写成既有文献的独立参数测量。六个样本最重要的共同点，不是 A 或 c 具有相同数值，而是结构和质量明显不同的星系均能表现出 GRODE 所要求的“上升—转折—渐近平坦”暗晕速度结构。

11.1 NGC 3198：经典长平坦曲线

NGC 3198 是暗晕模型检验中的经典对象。阶段性反演得到 √A≈136 km/s、c≈4.7 kpc，且 r_max/c≈9.4，说明数据同时覆盖了转折区和较深的渐近区。其 1/V_h²–1/r² 线性度约为 R²≈0.986，是六个样本中最清晰的线性化案例之一。

11.2 NGC 2841：高质量压力样本

NGC 2841 具有显著核球和复杂重子结构，最大观测半径约 64 kpc。试算得到 c≈3 kpc，因此 r_max/c≈21，是六个样本中最深的渐近区检验之一。对如此复杂的高质量星系，同一 GRODE 函数形式仍能跨越约二十个 c 的径向尺度，值得进一步精确分析。

11.3 NGC 5055：下降型总旋转曲线

NGC 5055 的总旋转曲线在外区缓慢下降，但这并不要求暗晕分量下降。因为 V_obs²=V_bar²+V_h²，只要重子盘贡献下降得更快，总速度可以下降而 V_h 仍单调趋向 √A。该样本说明“总旋转曲线下降”与“GRODE 暗晕平台”并不矛盾。

11.4 NGC 2403：较低质量晚型盘

NGC 2403 表明 GRODE 候选适用范围并不局限于高质量大盘星系。去除中心可能受非圆运动影响的区域后，阶段性得到 √A≈119 km/s、c≈3.2 kpc，并有 r_max/c≈6.5。

11.5 NGC 6503：动力学漂亮但存在质光比压力

NGC 6503 的试算给出 r_max/c≈16，已经深入渐近区，且速度空间拟合与线性化反演得到的 A、c 相互接近。不过其最佳恒星盘质量归一化偏低，因此它应被保留为重要压力样本，而不是简单称作“完美案例”。

11.6 NGC 2903：强棒星系的检验

NGC 2903 中心受到棒结构与非圆运动影响，不宜将最中心区域强行纳入简单轴对称圆周模型。主要使用较规则外盘区时，阶段性得到 γ_d≈0.38、√A≈165 km/s、c≈3.7 kpc，线性化 R²≈0.988。这个结果说明：已知的非圆运动区需要单独处理，而不应直接归因于暗晕动力学失效。

12. Pota 等（2015）的 NGC 1407 工作

2015 年，Vincenzo Pota、Aaron J. Romanowsky、Jean P. Brodie、Jorge Peñarrubia、Duncan A. Forbes 等人在 Monthly Notices of the Royal Astronomical Society 发表了对椭圆星系 NGC 1407 的多示踪体动力学研究。他们联合使用恒星、蓝球状星团和红球状星团，通过 Jeans 动力学建模约束可见物质与暗物质分布。

除了 generalized NFW 等模型外，该工作明确采用了一种 cored logarithmic（LOG）模型，其暗物质累计质量写为：

M_d(r) = (1/G) · v₀² r³/(r₀² + r²)    (18)

Pota 等把 r₀ 定义为 core radius，把 v₀ 定义为 asymptotic circular velocity，并将该模型描述为暗物质分布的一种“alternative parametrization”。他们给出的最大似然参数约为 r₀=44^{+19}_{−17} kpc、v₀=483^{+69}_{−67} km/s，并指出 LOG 与 gNFW 均能较好拟合恒星与两类球状星团的动力学数据。

这项工作非常重要，因为它说明与 GRODE 积分解同形的函数形式，早已在真实椭圆星系的多示踪体动力学中表现出良好经验适用性；但在 Pota 等人的理论语境里，LOG 仍是众多可比较暗晕参数化之一，而不是一条更基本动力学方程的积分结果。

13. Pota LOG 模型与 GRODE 的严格对应

GRODE 导出的累计质量式为：

M_h(r) = (A/G) · r³/(r² + c²)    (19)

与 Pota 等的 LOG 模型逐项比较可得严格对应：

A = v₀²，    c = r₀    (20)

因此，Pota 2015 在 NGC 1407 上成功使用的 LOG halo 与 GRODE 的积分解在数学函数形式上完全一致，但理论路径不同：Pota 等从候选暗晕参数化出发，再做 Jeans 数据拟合；GRODE 则从一阶动力学方程出发，再推出 V(r)、g(r)、Φ(r)、M(r) 与 ρ(r)。

14. GRODE 与传统对数势的理论层级

传统 cored logarithmic potential 通常直接写为 Φ=(v₀²/2)ln(r²+r₀²)，其中 v₀ 与 r₀ 都作为势模型参数进入。GRODE 则在方程层级只含 A；尺度 c 是积分常数，由边界条件确定。

所以从参数结构看，可以把二者概括为：传统 LOG 是“两参数势模型”；GRODE 是“一个动力学参数 + 一个边界条件”。进一步无量纲化之后，GRODE 变成 du/dx=u(1−u²)，无量纲动力学本身甚至是零参数的。

在这一意义下，GRODE 并不是简单把 LOG 势换一种写法，而是在尝试回答 LOG 势本身没有回答的问题：为什么对数势会出现？GRODE 的回答是：因为平坦速度是径向尺度动力学的稳定点，对数势是这一稳定流积分后的结果。

15. 对暗物质概念的重新排序

传统逻辑通常是“暗物质密度→引力场→旋转速度”。GRODE 允许另一种排序：“暗晕动力学→引力场→等效质量与密度”。因此，由泊松关系反演出的 ρ_h(r) 首先应被理解为与 GRODE 动力学场等效的牛顿源密度；仅凭这一反演本身，并不能自动证明该密度必然对应同空间分布的某种微观粒子暗物质。

16. GRODE 的远区问题

传统 LOG halo 的经典疑问是 ρ∝r⁻² 会导致 M(<r)∝r，从而累计等效质量不收敛。若 GRODE 是基础动力学，则应区分“动力学等效质量”与“实际粒子暗物质总质量”是否必然等同。是否需要人为截断 GRODE，不能仅由传统 LOG halo 的粒子暗物质解释决定，而应交给百 kpc 级乃至更远的外晕动力学观测。

17. 结论与下一步检验

GRODE 以一阶径向尺度动力学为起点，解析地产生从内区近线性旋转到外区平坦速度的过渡，并连续推出加速度、对数势、等效累计质量和等效密度。其重要理论区别在于：对数势不再作为先验势函数输入，而是动力学方程加边界条件后的积分结果。六个 SPARC 长径向样本的阶段性统一试算显示，不同结构和质量的星系均能表现出 GRODE 候选所要求的暗晕速度结构；Pota 等（2015）对 NGC 1407 的工作则提供了独立历史证据，说明同形 LOG halo 在真实多示踪体动力学中具有良好经验表现。

目前这些结果仍不足以把 GRODE 宣布为普适定律。下一阶段需要在完整 SPARC 样本上预先定义筛选规则，统一传播距离、倾角、气体质量与质光比不确定性，进行 MCMC 后验分析和与 NFW/Burkert/Einasto/LOG 等模型的信息准则比较；同时应利用 M31、Centaurus A 等百 kpc 级外晕数据检验同一动力学是否能够跨越星系盘与真正外晕。

如果这些跨样本、跨尺度检验持续成立，那么 GRODE 的意义将不再只是一个旋转曲线拟合形式，而可能成为描述暗晕动力学的一种更基础表述。

参考文献

1. Pota, V., Romanowsky, A. J., Brodie, J. P., Peñarrubia, J., Forbes, D. A., Napolitano, N. R., Foster, C., Walker, M. G., Strader, J., & Roediger, J. C. (2015). The SLUGGS survey: multipopulation dynamical modelling of the elliptical galaxy NGC 1407 from stars and globular clusters. Monthly Notices of the Royal Astronomical Society, 450, 3345–3358. DOI: 10.1093/mnras/stv831.

2. 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. DOI: 10.3847/0004-6256/152/6/157.

3. Begeman, K. G. (1989). H I rotation curves of spiral galaxies. I - NGC 3198. Astronomy & Astrophysics, 223, 47–60.

4. Rusli, S. P., et al. (2013). Stellar-dynamical modelling of massive early-type galaxies, including logarithmic-halo modelling of NGC 1407. Monthly Notices of the Royal Astronomical Society.

阶段性说明

本文中的六星系 GRODE 参数、线性度及部分重子归一化参数属于基于公开旋转曲线数据的阶段性统一计算结果。正式投稿前，应使用原始公开数据重新执行完整的加权/MCMC 拟合，并统一公布样本选择、异常区剔除规则、误差传播方法、先验条件和代码，以避免后验筛选与模型比较偏差。

支持单位：