GRoDE 暗晕密度的第二个交叉检验

摘要

前文讨论了太阳位置暗物质体密度的一次“握手”：由银河系径向旋转曲线，经 GRoDE 动力学得到的等效晕密度，与 Gaia 垂直动力学得到的太阳邻域暗物质密度相容。本文进一步考察第二个量——太阳附近 |z|<1.1 kpc 范围内的暗物质面密度。

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

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

若将 GRoDE 等效球形晕密度沿垂直方向积分，可得到太阳附近暗物质柱面密度。采用此前由径向旋转数据确定的 GRoDE 参数，得到约 24.7 M☉ pc⁻²；而现代垂直动力学研究给出的暗物质面密度约为 24±4 M☉ pc⁻²。两条来源不同的动力学路线再次在太阳附近相遇。

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

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

一、第一次握手之后，还能再检验什么？

GRoDE 的母方程是：

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

其积分解为：

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

在牛顿引力的等效球形解释下，由泊松方程得到：

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

第一次握手比较的是太阳所在盘面附近的局部体密度，是一个“点”的检验。第二次握手则把视野扩展到银盘上下各 1.1 kpc，检验整个 2.2 kpc 高垂直柱中的暗物质总量。若等效球形晕只在盘面上偶然有效，那么沿 z 方向积分后更容易暴露偏差。

二、什么叫暗物质面密度？

在太阳的银河半径 R₀ 处，取一个很小的底面积 dA，向银盘上方和下方各延伸 z_max。圆柱内部的暗物质质量为：

dM = dA ∫₋zmax^+zmax ρ(R₀,z) dz

因此定义暗物质面密度：

Σ_DM(R₀;|z|<z_max) = dM/dA = ∫₋zmax^+zmax ρ(R₀,z) dz

这里的“面密度”不是圆柱侧表面的密度，而是把整根垂直柱中的质量沿 z 方向压缩到银盘平面后，每单位底面积所对应的质量。因此单位由体密度 M☉ pc⁻³ 乘一个长度 pc，变成 M☉ pc⁻²。

三、GRoDE 如何给出太阳附近的暗物质面密度？

GRoDE 当前给出的是一个等效球形晕。对于固定银河半径 R₀，在离开银盘高度 z 的位置，与银河中心的三维距离为：

r = √(R₀²+z²)

因此该位置的晕密度写成：

ρ_h(R₀,z) = ρ_h(√(R₀²+z²))

若上下对称，则太阳附近的暗物质面密度为：

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

取 z_max=1.1 kpc，并代入此前由银河系径向旋转数据确定的 GRoDE 参数，可以得到：

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

这里没有重新拟合垂直动力学数据。也就是说，这个数来自“径向旋转曲线 → 重子扣除 → GRoDE 参数 A,c → 晕密度 → 垂直积分”这一整条链。

四、垂直动力学如何测出同一个量？

太阳附近恒星会在银盘上下振荡。它们的高度分布、垂直速度和速度弥散共同约束垂直引力场：

K_z = −∂Φ/∂z

再结合 Jeans 方程和泊松方程，可以估计太阳附近一定高度范围内的总面密度 Σ_tot。扣除恒星与气体等重子面密度后，剩下的就是暗物质面密度。现代 Gaia/APOGEE 垂直动力学分析给出的典型结果约为：

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

这一结果主要来自恒星在 z 方向的运动，而不是银河旋转曲线。

五、第二次握手

于是，两条不同的动力学道路再次相遇。由径向旋转动力学出发，经 GRoDE 得到：

Σ_DM^GRoDE ≈ 24.7 M☉ pc⁻²

由垂直恒星动力学出发得到：

Σ_DM^obs ≈ 24 ± 4 M☉ pc⁻²

二者在误差范围内高度相容。这就是“旋转动力学与垂直动力学的第二次握手”。

六、为什么第二次握手比第一次更强？

第一次握手比较的是一个局部体密度 ρ(R⊙,0)，第二次比较的是 −1.1 kpc<z<+1.1 kpc 整个垂直范围内的积分质量。

如果 GRoDE 的等效球形晕一离开盘面就迅速偏离真实三维晕，那么误差会在积分过程中累积，面密度就更容易偏离观测。现在积分结果仍然与垂直动力学相容，说明至少在太阳邻域和约 kpc 的垂直尺度上，GRoDE 等效球形晕没有表现出明显失效。

但这里仍不能把“相容”夸大为“已经证明真实晕严格球形”。更稳妥的说法是：等效球形表示在太阳邻域不仅能描述盘面局部密度，也能给出合理的垂直柱质量。

七、两次握手并列

现在可以把两次交叉检验放在一起。

八、等效球形晕应如何理解？

从盘面旋转曲线直接得到的首先是三维势在盘面上的径向截面 Φ(R,0)，并不能唯一恢复真实三维势 Φ(R,z)。因此 GRoDE 当前推导出的球形晕，更严格地应称为“等效球形晕”。

第二次握手增加了一个新的信息：如果这种球形延拓只是纯数学替代，而真实晕的垂直结构完全不同，那么沿 z 积分的面密度未必仍能与观测相容。现在结果仍然落入独立测量区间，因此等效球形晕可能不仅在径向引力场意义上有效，在太阳邻域的垂直质量分布上也可能是一个很好的近似。

九、仍需保持的科学边界

第一，GRoDE 的 A,c 来自银河系旋转曲线的重子扣除，因此恒星盘和气体质量模型会影响结果。第二，观测得到的暗物质面密度也依赖重子面密度、动力学平衡、Jeans 模型以及速度椭球倾斜项等系统误差。第三，面密度仍然是积分量，不能单独唯一确定暗晕的完整三维形状。

因此目前最合适的结论是：GRoDE 的等效晕与太阳邻域垂直动力学，在局部体密度和垂直柱面密度两个层面均表现出相容性。

十、从“两次握手”走向三维银河系动力学

下一步可以继续检验完整垂直结构，包括 ρ(R⊙,z) 与 K_z(R⊙,z)。如果同一组 GRoDE 参数能够同时解释 V_c(R)、ρ_dm,⊙、Σ_DM(z) 和 K_z(z)，那么 GRoDE 将从旋转曲线动力学进一步进入三维银河系动力学。

结语

GRoDE 最初研究的是星系怎样旋转。从母方程出发，它给出旋转速度，再给出等效晕密度。第一次握手表明太阳处的局部体密度与垂直动力学相容；第二次握手进一步表明，在上下总计 2.2 kpc 的垂直柱中，GRoDE 给出的暗物质面密度也与独立垂直动力学结果相容。

从一个点，到一整根垂直柱；从局部体密度，到积分面密度。两条不同方向的动力学道路再次来到同一密度尺度。值得继续追问的问题因此变得更加清楚：为什么一个从径向旋转动力学诞生的等效晕，能够连续通过垂直方向的质量检验？

参考说明

文中数值基于本轮对话中采用的银河系 GRoDE 参数与太阳邻域垂直动力学代表性结果。正式论文投稿前应统一使用同一套银河系重子模型、明确误差传播，并对原始观测文献进行逐条引用。

GRoDE 暗晕密度研究

支持单位：