GRoDE研究札记｜晕密度公式

晕密度公式

——从 GRoDE 基本微分方程到暗物质晕密度公式及初步观测检验

柳林涛

关键词： GRoDE；暗物质晕；旋转曲线；晕密度；对数势；泊松方程；SPARC

一、GRoDE 基本方程：晕密度公式的真正源头

本文必须从 GRoDE 本身讲起。晕密度公式并不是从某个预设的暗物质密度剖面出发，也不是先假设一个对数势再由泊松方程反推而来。它的源头是 GRoDE 的基本径向动力学方程：

r · dV_h/dr = V_h (1 − V_h²/A)

这是本文的重中之重。它直接规定暗晕剩余旋转速度 V_h 随半径 r 的演化。A 具有速度平方量纲，√A 是该方程的稳定速度尺度：当 V_h≪√A 时，右端近似为 V_h，因而 dlnV_h/dlnr≈1，得到V_h∝r；当 V_h→√A 时，右端趋于零，速度增长自然停止并进入平坦区。

因此，平坦旋转并不是额外塞入的边界条件，而是 GRoDE 微分方程自身的稳定状态。下面只需对这一方程积分，就能得到通常用于拟合旋转曲线的代数形式。令 y=V_h²，则

r · dy/dr = 2y(1 − y/A)

积分得：  ln[y/(A−y)] = 2 ln r + C

把积分常数写成 c²，可得  y/(A−y)=r²/c²，于是

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

这就是说，熟悉的两参数速度律不是本文的经验起点，而是 GRoDE 基本微分方程的积分解。其中 A=V_{h,∞}² 决定外区速度尺度，c 是从内区近线性上升向外区平坦过渡的核心尺度。晕密度公式正是从这一 GRoDE 解继续向下推导出来的。

二、从 GRoDE 积分解到晕密度：由运动学进入动力学

由 GRoDE 基本方程积分得到 V_h²=A r²/(r²+c²) 后，圆周运动立即给出暗晕对应的径向引力加速度：

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

在球对称径向动力学下，有 g_h=dΦ_h/dr。对半径积分，得到势函数：

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

这一步非常关键：对数势不是本文的先验假设，而是由 GRoDE 速度律必然积分得到的结果。随后再把这一派生势代入球对称泊松方程：

∇²Φ = (1/r²) d/dr [ r² dΦ/dr ] = 4πGρ

即可得到由 GRoDE 诞生的晕密度公式：

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

因此，完整逻辑顺序应写成：GRoDE 基本微分方程 → GRoDE 速度积分解 → 引力加速度 → 势函数 → 泊松方程 → 晕密度。这个顺序不能颠倒。对数势和晕密度都是 GRoDE 的派生结果，而不是 GRoDE 的输入假设。

三、公式的物理结构

1. 有限中心密度　r→0 时，ρ_h(0)=3A/(4πGc²)。因此模型没有 NFW 型中心尖峰，而具有有限密度核。

2. 外区自然趋于 r⁻²　r≫c 时，ρ_h(r)≈A/(4πGr²)。这正对应外区近似平坦的旋转速度。

3. 两个参数具有直接动力学意义　A=V_{h,∞}² 决定速度尺度，c 决定核心尺度；两者共同给出中心密度。

4. 可写成无量纲统一曲线　令 x=r/c、ρ₀=3A/(4πGc²)，则 ρ/ρ₀=(x²+3)/[3(1+x²)²]。一旦按 c 和 ρ₀ 归一化，整个密度形状不再含自由参数。

ρ_h/ρ₀ = (x² + 3) / [3(1 + x²)²],    x=r/c

四、与传统对数势的关系：GRoDE 为什么不需要先引入 q

经典轴对称 logarithmic halo 常写成

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

其中 q 是等势面的轴比：q=1 为球形，q<1 为扁球，q>1 为长球。这个参数在研究三维晕形状、垂直运动或非球对称轨道时具有明确的几何用途。

但本文的路线不同。GRoDE 首先给出单一径向变量 r 的动力学速度律，并由这一速度律直接导出势与密度；在这一推导链中并没有出现必须额外引入的几何扁率参数 q。对当前基于盘面旋转曲线的径向动力学检验，q=1 的球形形式已经足以闭合地得到速度、势与密度，并在目前检验的星系中给出良好描述。

因此，本文不把 q 作为 GRoDE 晕密度公式的基本参数。更准确地说：传统轴对称对数势中的 q 是为了描述三维等势面形状而加入的附加几何自由度；而 GRoDE 的晕密度公式由径向动力学直接产生，在当前检验层次只需要 A、c 两个参数。若未来独立三维观测确实要求非球形结构，应当再由三维动力学去解释这种偏离，而不是在基本公式中预先加入 q。

五、如何检验晕密度公式

检验不应只停留在“旋转曲线拟合得好”。密度包含速度曲线的一阶导数，因此是更严格的形状检验。本文采用的基本步骤如下：

1. 从 SPARC 等高质量旋转曲线数据取得 V_obs、气体贡献和恒星盘/核球贡献。

2. 扣除重子贡献，构造暗晕剩余速度：V_h²=V_obs²−V_gas²−Υ_dV_disk²−Υ_bV_bulge²。

3. 仅用 V_h(r) 拟合 GRoDE 的 A 和 c。

4. 由拟合所得 A、c 计算 GRoDE 预测密度 ρ_GR(r)。

5. 另一方面，不采用任何暗晕参数化密度模型，直接由剩余速度反演：ρ_np(r)=[4πGr²]⁻¹ d[rV_h²]/dr。

6. 比较 ρ_np 与 ρ_GR，并用 Δρ=RMS{log10[ρ_np/ρ_GR]} 量化误差。

ρ_np(r) = [1/(4πGr²)] · d[rV_h²(r)]/dr

Δρ = RMS{ log₁₀[ρ_np(r)/ρ_GR(r)] }

六、四个星系的初步检验

本轮讨论中，我们先后对四个性质差异明显的星系进行了初步检验。下表列出的数值用于展示当前结果的量级；由于密度反演涉及重子质光比、平滑与求导选择，正式论文中仍需统一算法、误差传播与盲化/预注册规则后重新计算。

这四个对象形成了一个很有价值的梯度：从暗物质主导的矮星系，到普通大螺旋，再到强重子、结构复杂的系统，同一个两参数密度律并未出现灾难性失效。尤其在 DDO 154、NGC 3198 和 NGC 5055 中，非参数反演密度与 GRoDE 密度的总体差异处在约 0.02–0.08 dex 的量级。

NGC 6946 的误差较大并不应被掩盖。它提示当前检验的主要瓶颈之一是重子扣除和数值求导：当 V_h² 是 V_obs² 与较大的恒星/气体项之差时，质光比的小变化会被密度求导进一步放大。因此，密度公式的最终准确度必须在统一误差模型下评估。

七、这一结果的真正意义

这个密度函数在数学形式上与经典球形 logarithmic halo 的密度—势对相同，但它在本文中的来源不同：本文没有从 logarithmic potential 出发。这里首先存在的是 GRoDE 速度动力学；对数势只是 GRoDE 的积分结果，晕密度又是这一派生势经泊松方程得到的结果。

因此，本文的第一层研究价值是：GRoDE 给出了经典球形 logarithmic halo 密度律的一条动力学生成路径。第二层研究价值是：既然密度律是由 GRoDE 推导出来的，就可以把它当作一个明确的观测预言，用非参数反演的真实晕密度去检验其准确度。

若后续在更大样本、统一重子质量模型和完整误差传播下仍能保持较小的 Δρ，那么最值得强调的将不是“重新发现了对数晕”，而是：一个从动力学出发的两参数结构，能够同时解释旋转曲线和暗晕密度。

八、下一步：把“漂亮”变成“坚固”

1. 扩大到 20–100 个以上 SPARC 星系，并固定同一套样本选择标准。

2. 对距离、倾角、气体质量、恒星质光比和旋转速度误差进行蒙特卡洛传播。

3. 预先固定平滑/求导算法，避免根据结果调整窗口导致后验选择偏差。

4. 把 GRoDE 晕密度与 Burkert、pISO、NFW、coreNFW 等模型做同等自由度或信息准则下的公平比较。

5. 对暗物质主导低表面亮度星系优先做“干净样本”，再单独分析强重子星系。

6. 寻找独立于盘面旋转曲线的检验，例如垂直动力学、恒星流或透镜，以真正检验三维晕形状。

九、结语

GRoDE 的基本微分方程、它的速度积分解，以及由此诞生的晕密度公式可以写成一条连续链：

r · dV_h/dr = V_h (1 − V_h²/A)

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

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

最根本的一式是 r·dV_h/dr=V_h(1−V_h²/A)。速度律 V_h²=A r²/(r²+c²) 是它的积分解，而晕密度公式则由这一积分解经引力场、势函数和泊松方程继续推出。对当前径向动力学检验，不需要预先引入传统轴对称对数势中的 q 参数；A 与 c 已构成一个最简、可检验的两参数体系。

因此，这项工作的核心不是从对数势“借来”一个晕密度，而是从 GRoDE 的基本微分方程出发，让速度律、引力场、势函数和晕密度依次自然诞生。也正因为如此，GRoDE 从旋转曲线的运动学描述，进一步进入了可以直接讨论暗物质质量分布的动力学层面。下一步的关键，就是用更大样本检验这条由 GRoDE 推导出的晕密度律究竟能保持多高的观测精度。

参考文献

Binney, J. 1981, MNRAS, 196, 455.（经典 logarithmic halo / 非球形对数势模型）

Binney, J. & Tremaine, S. 2008, Galactic Dynamics, 2nd ed., Princeton University Press.

Evans, N. W. & Williams, A. A. 2014, MNRAS, 443, 791–798.（讨论球形 Binney logarithmic potential 及对应密度）

Lelli, F., McGaugh, S. S. & Schombert, J. M. 2016, AJ, 152, 157. SPARC: Mass Models for 175 Disk Galaxies with Spitzer Photometry and Accurate Rotation Curves.

支持单位：