暗晕动力模型及其跨星系验证

——从 GRODE 旋转律到径向—垂向独立检验

柳林涛

摘要

本文提出并整理一个两参数暗晕动力模型（GRODE）。其母方程为 R(dV/dR)=V(1−V²/A)，积分后得到暗晕旋转律 V_h²=A R²/(R²+c²)。在球对称假设下，同一组参数 A 与 c 进一步唯一决定径向加速度、引力势、体密度、包围质量、垂向恢复力、中心频率与中心密度，从而构成一组闭合的暗晕动力公式。模型的核心特征不是以多个独立经验函数分别拟合速度、密度和质量，而是由一个一阶动力方程连续生成完整的三维引力场结构。

为检验该模型是否超越单纯的旋转曲线拟合，本文构造“径向定晕—垂向扣晕—剩余吐盘”的跨方向检验：先利用径向旋转数据确定暗晕参数 A,c，再在独立的垂向动力学数据中扣除 GRODE 暗晕贡献，检验剩余项是否自行恢复出与光度学一致的指数恒星盘。基于公开文献参数与数据重建，银河系给出约 2.1 kpc 的动力学盘尺度，与 Bovy & Rix (2013) 的 2.15±0.14 kpc 一致；IC 5249 的初步两点检验给出约 7.0 kpc，与独立光度尺度 7±1 kpc 一致；NGC 1566 采用独立 H I 径向暗晕约束与恒星垂向速度弥散数据，扣晕后得到约 1.47 kpc 的指数尺度，与约 1.59 kpc 的近红外光度尺度相差约 6%。UGC 463 则表现出支持性迹象，但仍需完整的自洽 Jeans 迭代。UGC 7321、UGC 3091 与 NGC 3982 的试算同时暴露了垂向归一化、数据循环、气体质量与盘断裂等系统问题，说明该方法具有实际的可证伪性。

本文认为，现阶段结果支持将球形两参数 GRODE 作为暗晕研究的首要零假设，并将跨星系、跨方向的独立“吐盘”实验作为比单纯旋转曲线拟合更严格的观测检验。所有跨星系数值结果均应视为基于公开文献参数的阶段性再分析，未来仍需使用原始逐点数据、统一距离尺度与完整误差传播进行系统复核。

关键词：暗晕；旋转曲线；GRODE；垂向动力学；指数盘；暗物质；跨星系验证

1  引言

星系暗晕研究长期面临一个基本困难：旋转曲线能够高精度约束总径向引力，但重子盘与暗晕之间存在显著退化。只要给盘质量、暗晕密度律和核心尺度足够多的自由度，许多不同质量分解都可以产生相似的旋转曲线。因此，一个暗晕模型是否具有物理内容，不能只看它能否“拟合 V(R)”，更应看它能否在不增加参数的情况下预言其他方向、其他观测量的动力学。

本文采用一种相反的建模思想：不从密度分布出发猜测旋转曲线，而从一个速度径向演化方程出发，使速度、加速度、势、密度、质量与垂向恢复力成为同一动力结构的不同投影。模型只含两个基本尺度：A 控制速度平方尺度，c 控制核心长度尺度。若两者由径向数据锁定，则垂向力、中心密度和包围质量不再允许自由调整。

由此自然产生一种比旋转曲线拟合更严格的测试：用径向数据确定暗晕，再用完全独立的垂向动力学数据扣除暗晕，观察剩余是否自行恢复出望远镜直接看到的指数恒星盘。本文将这一程序称为“吐盘检验”。如果暗晕扣错，剩余盘的径向斜率和尺度长度将系统偏移；如果多个不同尺度星系都能在不输入光度盘尺度的情况下恢复各自正确的 R_d，则构成跨星系、跨方向的一致性证据。

2  GRODE 暗晕动力模型

2.1  母方程与旋转律

模型的起点为速度随半径的径向演化方程：

其中 A 具有速度平方的量纲。V=√A 是径向演化的渐近定点：当 V<√A 时旋转速度随半径增长，当 V→√A 时增长率趋于零。因此 √A 对应外区的渐近平坦速度尺度。对母方程积分得到

其中 c 为积分产生的长度尺度。在中心 R≪c 时 V_h≈(√A/c)R，表现为近似刚体旋转；在外区 R≫c 时 V_h→√A。由此，一个表达式同时连接中心线性上升区与外区平坦区。

2.2  径向加速度、势、密度与质量

这些关系全部由同一个旋转律生成。中心区域有 ρ→ρ_0、M(<r)∝r³；外区则有 ρ∝r⁻²、M(<r)∝r，从而自然产生 V→常数。这里的“质量线性律”“密度平方反比律”和“平坦旋转”并非三个相互独立的假设，而是同一远区动力结构的等价表达。

2.3  垂向力与中心特征量

在盘面附近 z→0，球形势给出 K_z≈−[A/(R²+c²)]z。因此，径向旋转曲线一旦锁定 A,c，垂向恢复力的绝对幅度也随之确定。这正是本文跨方向检验的基础。需要强调的是，ν=Ω 在 z→0 对任何球对称势均成立，不能作为 GRODE 独有预言；真正具有模型特异性的，是有限高度下的完整 K_z(R,z) 以及由同一 A,c 决定的径向形状和归一化。

3  跨方向“吐盘”检验

3.1  基本思想

传统质量分解通常同时调节恒星盘、气体盘和暗晕。本文采用更严格的反向程序：首先只用径向旋转数据锁定暗晕；然后冻结暗晕参数，在独立垂向动力学数据中扣除暗晕；最后检验剩余项是否形成指数盘。理想情况下，光度盘尺度 R_d 在此之前不参与暗晕定参。

若在多个半径上得到的剩余恒星面密度满足指数律，则可由动力学独立拟合 R_{d,dyn}，最后与光度学 R_{d,phot} 比较。这一检验的关键不是“盘是否存在”，而是暗晕扣除后剩余盘的径向斜率是否正确。

3.2  Jeans 形式

对于恒星垂向分布，更严格的处理应直接使用 Jeans 方程，而不是把某一高度的 K_z 粗略等效为面密度：

若采用指数垂向密度 ρ_*(z)∝exp(−|z|/h_z)，GRODE 暗晕对盘面速度弥散的贡献可写为

                            σ_z,h²(R,0)=∫₀^∞ e^-z/h_z Az/(R²+z²+c²)dz  

在实际应用中，应从观测 σ_z² 中按平方扣除暗晕与气体的动力学贡献，再反演恒星盘。该处理可避免把暗晕“等效面密度”直接相减造成的双重计数或径向形状失真。

4  跨星系验证

4.1  银河系：基准成功案例

银河系提供了最严格的基准，因为径向旋转约束与垂向恒星动力学都非常丰富。Bovy & Rix (2013) 使用 16,269 颗 SEGUE G 型矮星和 43 个单丰度族群，在约 4.5–9 kpc 范围测得垂向力和总表面密度，并动力学得到银河盘质量尺度 R_d=2.15±0.14 kpc。本文此前的 GRODE 分析采用径向数据锁定暗晕后，在同类多半径垂向力数据中扣除暗晕，得到约 2.1 kpc 的剩余盘尺度。就方法而言，这是“径向定晕—垂向吐盘”的最强案例，因为盘尺度来自独立动力学数据而非光度先验。

4.2  IC 5249：长尺度低表面亮度盘

IC 5249 是一个超薄、低表面亮度的边缘盘星系。van der Kruit 等给出径向尺度 R_d=7±1 kpc、垂向尺度约 0.65±0.05 kpc，并指出 H I 旋转曲线在光学盘的大部分范围已经接近平坦。本文利用公开动力学约束进行两点初步重建：扣除 GRODE 暗晕的垂向贡献后，两个代表半径的剩余盘面密度比值给出 R_d≈7.0 kpc，与独立光度尺度一致。由于该结果目前主要依赖少数代表点，因此证据等级应标为“初步成功”，未来需要逐半径 K_z 或 σ_z 数据进行连续拟合。

4.3  UGC 463：支持性案例与自洽迭代问题

UGC 463 的 DiskMass 数据同时包含旋转、气体、恒星速度弥散和盘光度。原研究发现扣除气体后，动力学恒星面密度与 K 波段光度轮廓具有良好的径向一致性，平均 K 波段质光比约 0.22。与此同时，该星系在 R≳R_d 已明显受暗物质控制，球形暗晕可能使 σ_z 增加约 10%–15%，对应传统盘面密度估计约 20%–30% 的归一化偏差。

这说明 GRODE 在 UGC 463 上更可能首先修正盘质量归一化，而不是显著改变指数尺度。现有结果支持盘形状约 3.4–3.6 kpc 与光度一致，但由于 DiskMass 的径向暗晕分量本身使用了垂向动力学推得的恒星盘，严格的“盲吐盘”存在数据循环。因此本文将 UGC 463 记为支持性案例，而非完全独立的验证。

4.4  NGC 1566：独立数据集的第四战

NGC 1566 是本文最重要的外星系独立测试之一。其径向暗晕来自 WALLABY/ASKAP H I 旋转研究，而恒星垂向速度弥散来自更早的独立光谱观测，两套数据在观测路径上相互独立。WALLABY 给出的质量模型允许构造独立暗晕旋转曲线；本文用 GRODE 直接拟合该径向暗晕曲线，得到近似 √A≈123.5 km s⁻¹、c≈3.55 kpc。随后在约 h_z≈0.6 kpc 的垂向结构假设下，将 GRODE K_z 放入 Jeans 积分并从观测 σ_z² 中扣除暗晕贡献。

在约 1.5–6 kpc 的盘区，扣晕后的 σ_{z,disk}² 呈清晰指数下降，拟合得到 R_{d,GRODE}≈1.47 kpc；将独立 H 波段老恒星盘尺度统一到 WALLABY 所用距离后约为 R_{d,phot}≈1.59 kpc，两者相差约 6%。该结果目前仍是基于公开拟合参数的重建，但因径向与垂向观测独立，具有较高的检验价值。

图1  NGC 1566 第四战“吐盘”示意。圆点为扣除 GRODE 暗晕后的盘动力学代理量，实线为指数拟合，虚线为独立光度盘尺度对应的指数曲线。该图基于公开参数的近似重建。

4.5  阶段性汇总

注：上表中的 GRODE 数值属于基于公开文献参数的阶段性再分析，不应替代对原始逐点数据的正式联合拟合。

5  压力测试、失败点与方法学约束

5.1  UGC 7321：垂向归一化问题

UGC 7321 的 H I flaring 数据显示，不同分析对 H I 厚度、自吸收和速度弥散的处理可显著改变垂向力归算。早期研究与后续三维建模甚至得到明显不同的暗晕形状结论。本文的重新检查表明，若直接把径向暗晕归一化带入垂向数据，可能出现约几十个百分点乃至因子约 2 的幅度张力，而径向形状却仍相对一致。这提示在引入非球形参数 q 之前，应首先审计 σ_z、FWHM/尺度高度定义、H I 自吸收、气体自引力和静力平衡假设。

5.2  UGC 3091：数据循环

UGC 3091 的 DiskMass 暗晕旋转分量是通过总旋转减去重子旋转得到，而重子旋转中的恒星盘又来自同一套垂向速度弥散数据。因此，若再用该暗晕回到同一 σ_z 数据中“吐盘”，会形成 σ_z→Σ_*→V_b→V_DM→K_z,h→σ_z 的循环。该对象提醒我们：严格吐盘检验必须要求“径向定晕数据集 A”与“垂向吐盘数据集 B”观测独立。

5.3  NGC 3982：盘断裂与分子气体不确定性

NGC 3982/UGC 6918 虽具有丰富 DiskMass 数据，但其内盘存在 Type-II 断裂，速度弥散斜率也发生变化，同时分子气体质量极高且不确定。该对象的简单扣晕不能恢复单一指数盘，因此不宜被计为 GRODE 的胜负点。它表明样本选择本身是理论检验的一部分：单指数盘、温和气体修正、可靠 h_z 和独立旋转曲线，是更适合的实验条件。

6  关于暗晕扁率 q：为什么现阶段应以球形为零假设

暗晕形状研究常引入轴比 q，但 q 并非直接观测量，而是经过厚度、速度弥散、自吸收、恒星盘和气体盘等多层归算后的模型参数。本文的跨星系实验显示，至少在若干成功或支持性案例中，并不需要预先引入 q 即可恢复正确的盘尺度。因此，在现阶段更合理的策略是把球形 GRODE 作为首要零假设；只有在统一归算后多个星系仍系统、显著地拒绝球形，才应引入额外几何自由度。

这种做法还有统计学上的优势：A,c 两参数已经同时决定径向与垂向场。保持最小自由度，可以让失败具有明确含义；过早引入 q 则可能把数据归一化问题吸收到几何参数中，从而降低模型的可证伪性。

7  讨论

7.1  两参数模型的真正优势：多个独立预言

GRODE 的价值不在于“两参数也能拟合旋转曲线”，而在于两参数一旦由径向数据确定，就同时锁定 Φ、ρ、M 和 K_z。低自由度使模型承担更多观测责任：任何一个方向的独立数据不符合，都可以对整个结构形成压力；反之，如果径向旋转、垂向动力学与光度盘结构在多个星系上同时吻合，其证据强度显著高于单一曲线拟合。

7.2  从二维旋转到三维场的检验逻辑

对于球对称 GRODE，盘面旋转曲线只是三维势的一个投影。真正的动力学检验应当回答：由盘面径向运动锁定的势，是否能在离开盘面的方向产生正确恢复力？“吐盘实验”正是这一问题的具体实现。它把传统 disk–halo degeneracy 从一个联合拟合问题转化为一个预测问题：暗晕先被独立定参，恒星盘随后作为剩余结构被预测出来。

7.3  尚待完成的工作

现阶段的跨星系结果仍有明显局限。第一，多数外星系检验使用公开论文中的拟合参数、图形读数或代表半径，而不是统一处理原始数据。第二，不同文献采用的距离、倾角、垂向密度函数与气体修正不完全一致。第三，部分对象的暗晕径向分量依赖恒星盘动力学，可能导致数据循环。第四，误差尚未通过蒙特卡洛或贝叶斯层级模型完整传播。

下一阶段应建立统一样本：选择具有独立 H I/Hα 径向旋转数据、独立恒星垂向速度弥散或 K_z 数据、可靠 H I/H₂ 面密度及近红外光度的星系；对每个星系只用径向数据拟合 GRODE A,c，再冻结参数进行 Jeans 垂向反演。最终比较 R_{d,dyn}/R_{d,phot} 的样本分布，而不是依赖单个案例。

8  结论

本文建立了一个由一阶速度演化方程出发的两参数暗晕动力模型。其旋转律 V_h²=A R²/(R²+c²) 可以自然连接中心线性旋转与外区平坦速度，并在球对称条件下唯一生成径向加速度、对数型势、有限核心密度、线性增长的远区包围质量以及完整垂向恢复力。两个参数分别对应动力强度尺度和核心长度尺度，但其后果覆盖多个可独立观测的物理量。

比旋转曲线拟合更重要的是跨方向验证。本文构造并实施了“径向定晕—垂向扣晕—剩余吐盘”检验。银河系给出强基准；NGC 1566 在相互独立的径向 H I 与垂向恒星数据下得到约 1.47 kpc 的动力学盘尺度，与约 1.59 kpc 的光度尺度接近；IC 5249 给出约 7 kpc 的初步一致结果；UGC 463 提供支持性迹象但仍需自洽迭代。与此同时，UGC 7321、UGC 3091 和 NGC 3982 暴露出垂向归一化、数据循环和复杂重子结构等问题，说明 GRODE 并非通过任意自由度自动吸收所有数据。

因此，现阶段最谨慎而有力的结论不是“GRODE 已被证明”，而是：GRODE 已从一个旋转曲线公式进入可进行跨星系、跨方向、跨观测手段压力测试的暗晕动力模型阶段。已有结果支持继续以球形两参数形式作为首要零假设，并通过统一原始数据分析检验其能否在更大样本中重复恢复真实恒星盘结构。

参考文献

[1] Bovy, J. & Rix, H.-W. 2013, “A Direct Dynamical Measurement of the Milky Way’s Disk Surface Density Profile, Disk Scale Length, and Dark Matter Profile at 4 kpc ≲ R ≲ 9 kpc”, The Astrophysical Journal, 779, 115. DOI: 10.1088/0004-637X/779/2/115.

[2] van der Kruit, P. C., Jiménez-Vicente, J., Kregel, M., & Freeman, K. C. 2001/2002, “Kinematics and dynamics of the ‘superthin’ edge-on disk galaxy IC 5249”, Astronomy & Astrophysics; arXiv:astro-ph/0109477.

[3] Westfall, K. B., Bershady, M. A., Verheijen, M. A. W., Andersen, D. R., Martinsson, T. P. K., Swaters, R. A., & Schechtman-Rook, A. 2011, “The DiskMass Survey. IV. The Dark-Matter-Dominated Galaxy UGC 463”, The Astrophysical Journal, 742, 18. DOI: 10.1088/0004-637X/742/1/18.

[4] Bershady, M. A., Verheijen, M. A. W., Swaters, R. A., Andersen, D. R., Westfall, K. B., & Martinsson, T. 2010, “The DiskMass Survey. I. Overview”, The Astrophysical Journal; arXiv:1004.4816.

[5] Elagali, A. et al. 2019, “WALLABY early science – III. An H I study of the spiral galaxy NGC 1566”, Monthly Notices of the Royal Astronomical Society, 487, 2797–2817. DOI: 10.1093/mnras/stz1448.

[6] Bottema, R. 1993, stellar kinematics of spiral galaxy disks including NGC 1566, Astronomy & Astrophysics, 275, 16.

[7] O’Brien, J. C., Freeman, K. C., & van der Kruit, P. C. 2010, series on H I flaring and dark halo shape in edge-on disk galaxies, Astronomy & Astrophysics.

说明

本文为阶段性理论与数据再分析稿。文中 GRODE 跨星系验证数值主要根据公开论文中的参数、拟合曲线及代表点重新计算，尚未对所有对象的原始观测数据进行统一管线重处理。正式投稿前建议完成：① 原始逐点数据获取与统一距离尺度；② 对 A,c 的直接拟合及协方差估计；③ Jeans 方程全高度数值求解；④ 气体、厚度、倾角和质光比的联合误差传播；⑤ 对“吐盘”尺度进行盲化分析，以避免光度尺度先验影响。

支持单位：