——从太阳邻域的“握手”，到银河盘多半径上的连续统一检验

关键词：GRoDE；银河旋转曲线；垂直动力学；暗晕；Kz；银河盘尺度长；跨方向检验

说明：本文为阶段性研究稿。文中的银河系 GRoDE 参数、垂直力拟合及盘尺度长结果用于总结当前分析思路；正式投稿前仍需统一数据版本、重子模型与误差传播，并对原始观测表进行可复现核验。

阶段性研究稿 · 2026

GRoDE · 银河系径向—垂直统一检验

一、问题：两套银河动力学，能否来自同一个势？

银河系动力学长期沿着两条互补路线发展。一条沿银河盘方向研究恒星和气体围绕银心的旋转，由旋转曲线约束径向引力；另一条垂直于银河盘，通过恒星在 z 方向的运动约束垂直引力。

V_c²(R) = R ∂Φ/∂R

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

前者是旋转动力学，后者是垂直动力学。它们观察的是同一个银河引力场，只是沿着近乎正交的两个方向。真正值得追问的并不是“能否拟合一条旋转曲线”，而是：同一个动力学模型，能否同时解释银河系径向与垂直方向的引力场？

同一个 Φ(R,z)  →  同时决定 ∂Φ/∂R 与 ∂Φ/∂z

这正是 GRoDE 可以接受跨方向检验的地方。

二、GRoDE 从旋转动力学出发

GRoDE 的基本方程写成

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

其积分解为

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

其中 √A 决定外区渐近平坦速度，c 则给出由内区上升向外区平坦过渡的特征尺度。重要的是，这里并没有预先指定某一种暗物质密度剖面。

由旋转速度继续得到径向加速度

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

对径向加速度积分可得到对数势形式

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

在等效球形延拓下，由泊松方程可得晕密度

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

因此，在这一路径中，对数势不是预先假设，而是由 GRoDE 的旋转动力学推导出来的。

三、两个参数一旦锁定，垂直动力学也被锁定

若采用等效球形延拓

r² = R² + z²

则三维势写成

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

于是垂直加速度立即得到

K_{z,h}(R,z) = − A z/(R² + z² + c²)

这一步是整个“拥抱”检验的关键。因为 A 与 c 已经由径向旋转数据确定，所以 K_{z,h} 不再拥有新的自由参数。

径向动力学确定 A,c  ⇒  垂直动力学成为零新增参数预测

四、太阳邻域：从第一次“握手”开始

采用当前阶段银河系旋转分解中使用的一组代表性 GRoDE 参数

√A ≈ 190.9 km s⁻¹,     c ≈ 5.75 kpc

并取太阳半径 R₀ ≈ 8.15 kpc，GRoDE 给出的太阳处暗晕密度约为

ρ_{h,⊙} ≈ 0.0113 M⊙ pc⁻³

这一数值与现代垂直恒星动力学对太阳邻域暗物质密度的估计处在同一量级，并在我们当前采用的数据组合中显示出良好相容性。由此，旋转动力学与垂直动力学首先在太阳附近相遇。

但单点吻合远远不够。真正重要的问题是：离开太阳以后，它们还能不能继续相遇？

五、从“一个点”升级为“一条径向曲线”

GRoDE 在任意银河半径 R 都可以预言垂直力。固定 z = 1.1 kpc，有

|K_{z,h}(R,1.1)|/(2πG) = A·1.1 / {2πG [R² + 1.1² + c²]}

用同一组 A 与 c，可得到如下暗晕垂直力预测：

于是实验问题变得非常直接：

K_{z,obs}(R,z) − K_{z,bar}(R,z)  ?=  K_{z,h}^{GRoDE}(R,z)

如果这种关系只在太阳处成立，那仍可能只是局部巧合；如果它在多个 R 上连续成立，就意味着径向动力学与垂直动力学开始沿银河盘形成一条连续“拥抱曲线”。

六、反向实验：扣掉 GRoDE 晕，剩下来的究竟是什么？

还可以把问题反过来。我们不先规定银河盘长什么样，而是直接从观测总垂直力中扣掉由径向旋转动力学给出的 GRoDE 暗晕：

K_{z,disk}^{inferred} = K_{z,obs} − K_{z,h}^{GRoDE}

然后问：剩下来的函数是否像一个正常银河盘？如果 GRoDE 暗晕严重错误，那么扣除后很容易留下异常的径向结构；相反，如果剩余项自然接近指数盘，而且尺度长与独立银河盘测量一致，就会形成一种跨方向的内部一致性检验。

七、43 个垂直动力学测量点的阶段性联合检验

在 Bovy–Rix 的多半径垂直动力学数据中，可使用 43 个不同恒星化学族群（MAP）对应的 K_Z(R,1.1 kpc) 约束。我们在当前阶段性分析中固定 GRoDE 的 A 与 c，不让它们根据垂直数据重新调整，只允许剩余盘项采用普通指数形式：

K_{z,disk}/(2πG) = B exp[−(R−8)/R_d]

因此拟合形式为

K_{z,obs} = K_{z,h}^{GRoDE} + K_{z,disk}

对全部 43 个测量点的阶段性联合分析得到盘尺度约

R_d ≈ 2.45 kpc

在排除内银河距离系统误差较明显的 R < 5.5 kpc 区域后，结果收敛为

R_d ≈ 2.10 ± 0.16 kpc

进一步改变径向下限，得到的尺度长仍稳定在约 2.0–2.1 kpc。与此同时，独立银河动力学分析给出的恒星质量盘尺度长约为

R_d ≈ 2.15 ± 0.14 kpc

两者中心值只相差约 0.05 kpc。这个结果值得高度关注：旋转动力学先给出暗晕，垂直动力学在扣除这一暗晕以后，剩余项自行显现出一个与真实银河恒星盘相符的径向尺度。

重要边界：这里的数值属于当前工作阶段的可复现目标结果，正式论文版本必须对原始 43 点表、误差条、选择函数、距离系统误差与重子模型进行统一复核，并给出完整协方差与模型比较。

八、“径向给晕，垂直吐出盘”

整个过程可以压缩成一条非常直观的逻辑链：

银河旋转曲线  ⇒  A,c  ⇒  K_{z,h}^{GRoDE}

K_{z,obs} − K_{z,h}^{GRoDE}  ⇒  R_d ≈ 2.1 kpc

也就是说，径向动力学没有把盘尺度长输入给垂直分析；但垂直数据在扣除径向预测的暗晕后，给出了一个接近真实银河盘的尺度长。

这不是对 GRoDE 的最终证明，却是一种非常具体、非常有区分力的跨方向一致性：

径向动力学给出的暗晕一扣掉，垂直动力学重新显现出了银河盘。

九、为什么多半径“拥抱”比太阳单点“握手”更强？

太阳处的单点比较只要求一个位置上的数值相容，而多半径检验要求同一个 A、c 在 R₁、R₂、R₃……上一整段银河半径同时成立。信息量由“一个数”升级为“一条函数”。

ρ_h^{radial}(R)  ?=  ρ_h^{vertical}(R)

Σ_h^{radial}(R)  ?=  Σ_h^{vertical}(R)

K_{z,h}^{radial}(R,z)  ?=  K_{z,h}^{vertical}(R,z)

如果未来从 R≈6 kpc 一直到 R≈12–14 kpc 都能保持这种连续一致性，那么被检验的就不再是一条旋转曲线，而是一个二维银河引力场。

十、这里并不需要预先引入 q

传统轴对称对数晕通常写作

Φ(R,z) = (v₀²/2) ln(R_c² + R² + z²/q²)

其中 q 是用于描述扁率的几何参数。GRoDE 的研究顺序不同：它首先从旋转动力学导出势，再让垂直动力学检验这个势是否足够。

至少在当前内银河与太阳邻域的检验中，并没有出现“必须依靠额外 q 才能挽救模型”的明确要求。独立的银河暗晕形状研究也允许甚至偏好近球形结构。因此，更自然的顺序可能是：先用动力学决定势，再由三维观测判断是否真的需要额外几何自由度。

先定动力学 → 再检验几何；而不是先塞入几何参数 → 再调动力学。

十一、从“握手”到“跨尺度拥抱”

研究可以分成两个层次。第一层次是在太阳邻域，径向与垂直动力学完成一次局部“握手”；第二层次则要求在多个 R 上持续成立，使两套动力学沿银河盘形成连续一致性。

R = 6, 7, 8, 9, 10, 11, 12 kpc …

当这种一致性沿半径展开时，“拥抱”不再只是修辞，它有严格的数学含义：

∂Φ/∂R 与 ∂Φ/∂z 在不同 R、不同 z 上来自同一个 Φ(R,z)

如果这一点成立，那么旋转动力学与垂直动力学就不再是两套互不相干的经验研究，而只是同一个银河动力学的两个投影。

十二、当前仍需保持的科学边界

第一，银河内盘存在棒结构、旋臂与非轴对称运动，R≲6–7 kpc 的垂直平衡假设需要格外谨慎。

第二，银河外盘存在 warp、bending mode 与 breathing mode，静态 Jeans 分析并非处处适用。

第三，重子质量模型仍是主要系统误差来源之一；径向重子扣除与垂直重子扣除必须尽量使用同一套重子模型。

第四，目前的球形延拓首先应理解为等效球形结构。若未来高精度三维数据明确要求非球形结构，应进一步研究这种偏离是否能够从更完整的动力学中自然产生。

因此，当前最稳健的表述不是“GRoDE 已经证明了银河暗晕结构”，而是：

由银河旋转动力学独立确定的 GRoDE 暗晕，正在表现出与多半径垂直动力学相当有趣的跨方向一致性。

十三、下一步：二维银河动力学地图

下一阶段的目标不应再满足于少数半径，而应构造 K_z,obs(R,z) 的二维地图，并逐点比较

K_{z,h}^{GRoDE}(R,z) = − A z/(R² + z² + c²)

同时检验

ρ_h(R,z),     Σ_h(R,z),     K_{z,h}(R,z)

如此，研究将从一维 V(R) 的旋转曲线，扩展到二维 Φ(R,z) 的统一动力学检验。

结语

银河旋转曲线告诉我们银河在盘面方向怎样产生引力；恒星垂直运动告诉我们银河在竖直方向怎样产生引力。经典力学要求，它们最终必须来自同一个势。

GRoDE 从一个极其简单的径向动力学方程出发，却可以一路走向 V(R)、Φ(R,z)、ρ(R,z)、Σ(R,z) 与 K_z(R,z)。如今，径向旋转所确定的暗晕，从独立垂直动力学数据中扣除以后，剩余项重新呈现出一个尺度长约 2.1 kpc 的银河盘。

这还不是终点，却已经把一个非常具体的问题摆在我们面前：银河系的旋转动力学与垂直动力学，是否正在不同半径、不同尺度上，逐渐显露出同一个动力学结构？

如果这种一致性从太阳附近一直保持到银河外盘，那将不再只是几次漂亮的“握手”，而是旋转动力学与垂直动力学在整个银河尺度上的拥抱。

参考文献与数据线索（工作版）

Bovy, J. & Rix, H.-W. (2013). A Direct Dynamical Measurement of the Milky Way’s Disk Surface Density Profile, Disk Mass, and Dark Matter Profile at 4 kpc < R < 9 kpc. The Astrophysical Journal, 779, 115.

Gaia Collaboration 及后续银河垂直动力学、加速度场与暗晕形状研究：用于太阳邻域暗物质密度、K_z(R,z)、三维引力场和晕轴比的独立约束。

本文 GRoDE 参数与多半径拟合结果为当前对话阶段的工作数值；正式稿应附原始数据表、代码、误差传播与重子模型说明。

阶段性研究稿 · 2026

支持单位：