——GRoDE 晕密度公式的一次跨方向动力学检验

摘要

星系旋转曲线研究的是银盘内的径向动力学，而太阳邻域暗物质密度还可以通过恒星相对于银盘的垂直运动独立测量。两类观测在几何方向、数据来源和分析方法上不同。

GRoDE 从一个简洁的旋转动力学方程出发，先给出暗晕旋转速度，再在牛顿引力的等效球形解释下导出晕密度。利用银河系径向旋转数据和现代重子分解，GRoDE 在太阳位置给出的局部晕密度约为 0.0113 M☉·pc⁻³；而独立的 Gaia 垂直动力学测量给出约 0.0117 ± 0.0035 M☉·pc⁻³。两个来源不同的结果在太阳位置相遇，构成一次值得重视的“径向动力学—垂直动力学”交叉检验。

一、从 GRoDE 开始，而不是从暗晕模型开始

这项研究的出发点不是选择 NFW、Burkert、伪等温球或 logarithmic halo 等现成暗晕模型，而是 GRoDE 的动力学母方程：

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

当 V² ≪ A 时，方程近似为 r dV/dr ≈ V，因此 V ∝ r，内区呈近似线性上升；当 V² → A 时，dV/dr → 0，速度趋向稳定值 V∞ = √A。也就是说，平坦旋转速度是方程自身的稳定状态，而不是额外附加的条件。

积分得到：

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

其中 c 是积分产生的长度尺度，决定旋转曲线由上升区向平坦区过渡的位置。至此，GRoDE 首先解决的是“怎样转”的问题。真正的关键，是继续从速度走向质量分布。

二、从旋转速度到晕密度

圆周运动满足：

dΦ/dr = V_h²/r

将 GRoDE 速度律代入：

dΦ/dr = A r/(r² + c²)

积分得到径向势：

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

这里的逻辑顺序非常重要：不是先假设对数势，再得到 GRoDE；而是 GRoDE → 速度 → 引力场 → 势。对数形式是 GRoDE 的派生结果。

如果把这个径向势作球形牛顿延拓，再用泊松方程 ∇²Φ = 4πGρ，就得到：

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

这就是 GRoDE 的等效球形晕密度公式。它立即给出有限中心密度 ρ₀ = 3A/(4πGc²)，并在 r ≫ c 时自然过渡到 ρ_h ∝ r⁻²，对应外区近似平坦的旋转速度。

三、这里的“球形”首先是等效球形

从旋转曲线直接得到的是银盘上的径向引力：

V_c²(R) = R · ∂Φ/∂R |_(z=0)

因此严格说，旋转曲线确定的是三维势 Φ(R,z) 在盘面 z=0 上的径向截面 Φ(R,0)，而不是完整的三维势。很多不同的三维势都可以拥有相同的盘面径向截面。

因此，由球对称延拓得到的 ρ_h(r) 最稳妥的物理解释是：产生同一盘面径向引力场的“等效球形晕密度”。它不必等同于真实暗晕的严格三维形状。

四、径向动力学给出的太阳处晕密度

对银河系而言，可以先从观测总旋转速度中扣除恒星盘、气体和核球的重子贡献：

V_h²(R) = V_obs²(R) − V_bar²(R)

然后只用剩余的暗晕速度拟合 GRoDE：

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

采用现代较轻的银河系重子模型与径向旋转数据进行拟合时，我们此前得到的代表性参数约为 √A ≈ 190.9 km·s⁻¹，c ≈ 5.75 kpc；在太阳半径 R☉ ≈ 8.15 kpc 处，对应 V_h(R☉) ≈ 156 km·s⁻¹。

把这些参数直接代回晕密度公式，不再调参，得到：

ρ_GRoDE(R☉) ≈ 0.0113 M☉·pc⁻³

这个数来自径向旋转数据，而不是来自太阳附近恒星上下运动。

五、另一条道路：垂直动力学

太阳附近的恒星不仅绕银河中心旋转，还在银盘上下振荡。垂直方向的恢复力由

K_z(R,z) = −∂Φ/∂z

控制。通过大量恒星的离盘高度 z、垂直速度 v_z 以及恒星数密度随 z 的变化，可以利用 Jeans 方程和泊松方程反推出垂直引力场与当地总质量密度。扣除恒星、气体等重子成分以后，就得到太阳附近的暗物质密度。

Gaia DR3 的 K 型矮星垂直动力学研究给出的代表性结果为：

ρ_dm,☉ = 0.0117 ± 0.0035 M☉·pc⁻³

这类测量主要来自垂直方向的恒星运动，而不是银河旋转曲线，因此与 GRoDE 的径向预测具有明显不同的数据来源和动力学方向。

六、两条动力学道路在太阳处握手

现在出现了本文最值得关注的会合：

ρ_GRoDE(R☉) ≈ 0.0113 M☉·pc⁻³

ρ_vertical(R☉) = 0.0117 ± 0.0035 M☉·pc⁻³

中心值相差约几个百分点。考虑到垂直动力学本身仍有统计和系统误差，这不能被解释为“已经验证到百分之几精度”；真正重要的是，两种不同方向、不同数据和不同分析方法得到的结果落在同一个局部暗物质密度尺度上。

可以把两条道路概括为：

径向旋转曲线 → 重子扣除 → GRoDE(A,c) → ρ_h(R☉)

垂直恒星运动 → K_z(z) → Jeans + Poisson → ρ_dm,☉

最终，它们在太阳位置相遇。这就是“旋转动力学与垂直动力学的握手”。

七、为什么这比旋转曲线内部检验更重要

如果 A、c 由同一条旋转曲线确定，再从同一条旋转曲线做球对称密度反演，那么这主要是内部闭合检验。太阳处的交叉检验更进一步：GRoDE 主要从 ∂Φ/∂R 出发，而垂直动力学主要约束 ∂Φ/∂z。

一个研究盘面内的径向运动，一个研究盘面上下的运动。若二者都指向同一个暗物质密度尺度，其物理意义明显更强。

八、仍需保持的科学边界

首先，GRoDE 对太阳处密度的预测依赖银河系重子质量模型。若恒星盘质量取得过大，扣除后留下的 V_h 会过小，进而低估局部暗物质密度。因此，重子模型是这一检验的重要系统误差来源。

其次，当前 GRoDE 晕密度首先是等效球形密度。旋转曲线本身不能唯一恢复真实暗晕的三维几何。未来仍需要结合恒星流、卫星轨道、透镜和更完整的垂直动力学，判断等效球形密度与真实三维暗晕之间究竟有多接近。

九、值得继续追问的问题

真正值得追问的已经不只是“GRoDE 能不能拟合旋转曲线”，而是：为什么一个从径向旋转动力学诞生的晕密度，能够预测垂直动力学测得的局部暗物质密度？

如果这种一致性只是偶然，大样本和更多独立观测会很快将其打破；如果它能够持续存在，那么 GRoDE 的意义就会从旋转曲线描述进一步上升为连接旋转运动、引力场、势和暗晕密度的一套可检验动力学框架。

结语

GRoDE 从一个极其简洁的旋转动力学方程出发：

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

它先给出

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

再给出

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

现在，银河系提供了一个特殊检验机会：从银盘内的径向旋转出发，我们得到太阳处的 GRoDE 晕密度；从恒星在银盘上下的运动出发，Gaia 垂直动力学也给出了当地暗物质密度。

两条道路最终来到同一个地方。这个结果还不是终点，也不能代替更完整的误差传播和独立三维验证，但它已经构成了一幅非常清楚的物理图景：在太阳附近，径向旋转动力学与垂直恒星动力学发生了一次握手，而连接它们的，正是晕密度。

参考说明

1. Gaia DR3 K 型矮星垂直动力学研究给出的局部暗物质密度：ρ_dm,☉ ≈ 0.0117 ± 0.0035 M☉·pc⁻³。

2. 现代银河系旋转曲线和重子分解用于构造 V_h(R)，GRoDE 拟合只使用径向动力学信息；太阳处垂直动力学密度不参与 A、c 的拟合。

3. 本文中的 GRoDE 晕密度应首先理解为“等效球形晕密度”，真实三维暗晕形状需要更多独立观测确定。

支持单位：