——体密度、面密度与垂直引力的连续交叉检验

摘要

GRoDE 最初是一个关于星系旋转的径向动力学方程。其积分解给出暗晕旋转速度，并在牛顿引力的等效球形解释下导出晕密度。利用银河系径向旋转曲线所确定的同一组 GRoDE 参数，我们在太阳邻域连续比较了三个来自垂直方向的物理量：局部暗物质体密度、|z|<1.1 kpc 内的暗物质面密度，以及 z=1.1 kpc 处的总垂直引力。三组比较分别呈现出 0.0113 对 0.0117±0.0035 M☉ pc⁻³、24.7 对 24±4 M☉ pc⁻²，以及 70.7 对 71±6 M☉ pc⁻²（垂直力以 Kz/2πG 的等效面密度表示）的相容结果。第三次“握手”尤其重要：GRoDE 暗晕垂直力贡献约 14.7 M☉ pc⁻²，与独立重子盘垂直力约 56±6 M☉ pc⁻² 相加得到 70.7 M☉ pc⁻²，几乎与观测总垂直力中心值 71 M☉ pc⁻² 重合。本文同时强调，这些结果仍受重子模型、等效球形延拓及银河盘非平衡系统误差影响，因此应视为值得进一步严格检验的交叉一致性，而不是最终证明。

一、GRoDE：一切从旋转动力学开始

GRoDE 的核心并不是先选择一个暗晕密度模型，而是从径向速度的微分方程出发：

R dV/dR = V (1 − V²/A)

当 V²≪A 时，方程给出 V∝R；当 V²→A 时，dV/dR→0，旋转速度趋向稳定值 √A。因此，内区近线性上升与外区趋平是同一动力学方程的两个极限。积分得到：

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

其中 √A 是暗晕速度尺度，c 是由积分常数形成的空间尺度。到这里，GRoDE 仍主要是在描述“怎样转”。真正的转折发生在把速度结构连接到引力场和质量分布之后。

二、由 GRoDE 诞生的晕密度公式

圆周运动满足 ∂Φ/∂R = V_h²/R，因此由 GRoDE 速度律可得盘面径向势：

Φ(R,0) = (A/2) ln(R² + c²) + C

这里必须注意逻辑顺序：不是先假设对数势再拟合旋转曲线，而是 GRoDE → V_h → g_h → Φ。若将这一盘面势作球形牛顿延拓：

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

再使用泊松方程 ∇²Φ=4πGρ，得到等效球形晕密度：

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

这一公式把 GRoDE 从旋转动力学推进到了质量分布动力学。它还给出中心有限密度 ρ_h(0)=3A/(4πGc²)，并在 r≫c 时自动趋于 ρ∝r⁻²。

三、第一次握手：太阳处局部暗物质体密度

采用现代银河系旋转曲线与较轻的重子质量模型，前述 GRoDE 拟合给出代表参数 √A≈190.9 km s⁻¹、c≈5.75 kpc。在太阳半径 R☉≈8.15 kpc 处，由晕密度公式得到：

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

另一方面，基于 Gaia DR3 大样本 K 型矮星的垂直 Jeans 动力学测量给出：

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

前者来自径向旋转数据，后者来自恒星在银盘上下方向的运动。两条动力学道路在太阳处的局部暗物质密度上相遇：

0.0113  ↔  0.0117 ± 0.0035  M☉ pc⁻³

这是第一次握手。它不是同一条旋转曲线的内部重算，而是径向动力学与垂直动力学之间的交叉一致性。

四、第二次握手：太阳附近暗物质面密度

局部体密度只是 z=0 附近的一个点。进一步可以考察太阳半径处从银盘下方 1.1 kpc 到上方 1.1 kpc 的整根垂直柱。对于等效球形晕，r=√(R²+z²)，因此：

Σ_DM(R☉; |z|<z_max) = 2 ∫₀^{z_max} ρ_h(√(R☉²+z²)) dz

取 z_max=1.1 kpc，并代入与第一次握手完全相同的 GRoDE 参数，可得：

Σ_DM^GRoDE(|z|<1.1 kpc) ≈ 24.7 M☉ pc⁻²

而独立的太阳邻域垂直动力学分析给出的暗物质柱面密度约为：

Σ_DM^obs(|z|<1.1 kpc) ≈ 24 ± 4 M☉ pc⁻²

于是第二次握手为：

24.7  ↔  24 ± 4  M☉ pc⁻²

这一检验不再是一个点，而是对整个 2.2 kpc 高的垂直柱进行积分。需要特别说明：这个积分在等效球形晕假设下是严格的；同时由于 z_max/R☉≈0.13，太阳邻域这段垂直范围内 r 的变化不到约 1%，所以结果对局部几何并不极端敏感。

五、从面密度到垂直力：不能把 Σ 与 Kz/2πG 混为一谈

第二次握手之后，一个更深入的问题出现了：真实柱面密度 Σ_DM 与由垂直引力定义的 Kz/2πG 并不严格相等。一般轴对称泊松方程为：

4πGρ = (1/R) ∂/∂R (R ∂Φ/∂R) + ∂²Φ/∂z²

第一项由旋转曲线控制，第二项是垂直曲率。因此，在把 Kz 换算为面密度时必须保留径向修正。对 GRoDE 等效球形晕，虽然真实暗物质柱面密度为约 24.7 M☉ pc⁻²，但其 z=1.1 kpc 处的垂直力等效量并不是 24.7。

六、第三次握手：垂直引力几乎精确会合

GRoDE 等效球形势写成：

Φ_h(R,z) = (A/2) ln(R² + z² + c²)

直接求垂直导数：

K_z,h(R,z) = − A z/(R² + z² + c²)

取 R=R☉=8.15 kpc、z=1.1 kpc，并继续使用同一组 √A≈190.9 km s⁻¹、c≈5.75 kpc，得到暗晕垂直力：

|K_z,h|/(2πG) ≈ 14.7 M☉ pc⁻²

这是第三次握手的第一部分。它完全来自径向旋转数据确定的 GRoDE 参数，经三维等效势直接求 ∂Φ/∂z。

独立太阳邻域动力学研究给出的重子盘垂直力贡献约为：

|K_z,bar|/(2πG) ≈ 56 ± 6 M☉ pc⁻²

因此，GRoDE 暗晕与独立重子盘相加得到：

14.7 + 56 = 70.7 M☉ pc⁻²

而经典太阳邻域总垂直力观测为：

|K_z^obs(1.1 kpc)|/(2πG) ≈ 71 ± 6 M☉ pc⁻²

于是第三次握手的核心数值就是：

70.7  ↔  71 ± 6  M☉ pc⁻²

这里的中心值几乎重合。这比“量级一致”更强：GRoDE 暗晕垂直力约 14.7，加上独立测得的重子盘垂直力约 56，直接重建出观测总垂直力约 71。其物理结构可写成：

K_z,bar + K_z,h^GRoDE ≈ K_z^obs

第三次握手的重要性在于，比较对象已经不再是密度，而是真正的力场。GRoDE 从“给出质量分布”进一步走到了“给出垂直引力”。

七、三次握手汇总

三次握手并非三个完全统计独立的实验，但它们对应三个不同的物理层次：局部体密度、垂直积分质量和垂直引力。更重要的是，它们都由同一套从径向旋转数据确定的 GRoDE 参数向垂直方向延伸而来。

八、为什么第三次握手最有分量

第一次握手说明 GRoDE 的等效晕在太阳处给出了正确量级的局部暗物质密度；第二次握手说明把这一密度沿 z 积分以后，整根垂直柱的暗物质质量仍与观测相容；第三次握手则进一步要求同一个势的垂直导数也正确。

换言之，三次握手形成了从“密度”到“积分质量”再到“力”的连续检验链：

ρ  →  Σ  →  K_z

其中第三项最难伪装，因为它不只是对密度公式做积分，而是直接对三维势求垂直导数。正因为如此，70.7 与 71±6 的会合，是目前这组三次握手中最值得重视的结果。

九、从等效球形晕走向三维银河动力学

从盘面旋转曲线严格能够直接确定的是 Φ(R,0) 或 ∂Φ/∂R(R,0)。将其延拓为 Φ_h(R,z)=(A/2)ln(R²+z²+c²) 相当于选择了一个球形三维延拓。因此，当前的 ρ_h(r) 最稳妥的称呼仍是“等效球形晕密度”。

然而，如果这种球形延拓只是与真实三维晕完全无关的数学替身，那么离开盘面之后的柱面密度与垂直力应更容易出现明显偏差。现在体密度、柱面密度和垂直力三个层面都没有出现灾难性失败，提示这个等效球形结构可能抓住了太阳邻域暗晕的主要三维引力特征。

十、科学边界与下一步

目前结果仍需保持三个限制。第一，GRoDE 的 A、c 来自银河系重子扣除，因此恒星盘和气体模型是主要系统误差之一。第二，垂直动力学本身受银河盘非平衡、速度椭球倾斜项和重子垂直分布影响。第三，当前等效球形延拓并不能唯一确定真实暗晕的三维形状。

因此下一阶段应在同一套银河系重子模型下，对 A、c、ρ_dm,☉、Σ_DM 和 K_z 做统一误差传播，并进一步检验 K_z(z) 的整条曲线，而不仅仅是 z=1.1 kpc 的单点。若同一个 GRoDE 势能够同时解释 V_c(R)、ρ(R☉,0)、Σ(R☉,z) 和 K_z(R☉,z)，那么 GRoDE 才真正完成从旋转动力学向三维银河系动力学的跃迁。

结语

GRoDE 从一个简洁的径向方程出发，先得到旋转速度，再得到势和晕密度；随后，同一套参数在太阳邻域连续遇到了三个来自垂直方向的观测量。第一次握手是局部暗物质体密度，第二次握手是垂直暗物质面密度，第三次握手则是真正的垂直引力。

第一次：0.0113  ↔  0.0117 ± 0.0035  M☉ pc⁻³

第二次：24.7  ↔  24 ± 4  M☉ pc⁻²

第三次：70.7  ↔  71 ± 6  M☉ pc⁻²

第三次握手尤其值得注意：GRoDE 暗晕垂直力约 14.7 M☉ pc⁻²，加上独立重子盘贡献约 56±6 M☉ pc⁻²，得到 70.7 M☉ pc⁻²，与观测总垂直力 71±6 M☉ pc⁻² 的中心值几乎重合。

如果后续严格误差传播和更多数据仍然保持这种一致性，那么 GRoDE 所连接的就不再只是一条银河旋转曲线，而是径向运动、垂直质量分布和垂直引力共同组成的三维银河动力学。

参考资料（简要）

1. Gaia DR3 K 型矮星垂直 Jeans 动力学研究：太阳附近暗物质体密度 ρ_dm,☉≈0.0117±0.0035 M☉ pc⁻³。

2. 太阳邻域垂直动力学 / orbital-torus imaging 相关研究：暗物质面密度约 24±4 M☉ pc⁻²（|z|≈1.1 kpc）。

3. Holmberg & Flynn 等太阳邻域垂直动力学研究：重子盘垂直力贡献约 56±6 M☉ pc⁻²，总 Kz(1.1 kpc) 约 70–74 M☉ pc⁻²。

4. McKee, Parravano & Hollenbach：太阳邻域重子真实面密度约 47.1±3.4 M☉ pc⁻²；注意真实柱面密度与 Kz/2πG 并不严格相同。

GRoDE · 太阳邻域三次握手

支持单位：