星系本征加速度

——从牛顿点质量场到星系整体本征运动的动力学假说

摘要：本文提出星系整体本征加速度的候选动力学公式，将星系重子总质量M与空间加速度直接联系起来，而不预先引入暗物质晕及其密度分布。由星系旋转动力学方程GRoDE导出有限核心的径向结构，再提出独立的场强猜想A=β√M，得到α = −β√M·r/(r²+c²)。其中r表示由星系中心指向测点的位置矢量，r=|r|，c为核心尺度。该公式产生从中心上升、向外趋平的旋转曲线，并给出V_flat⁴=β²M。重子Tully–Fisher关系支持外围质量平方根标度，但有限核心结构、β的普适性和盘外三维适用性仍须独立检验。

关键词：GRoDE；本征加速度；本征力学；星系旋转曲线；重子Tully–Fisher关系；暗物质

一、牛顿点质量加速度与星系本征加速度

牛顿点质量场由源质量M和距离r直接决定。为避免任何字体或公式编辑器造成矢量丢失，本文统一使用下标“_vec”表示矢量：

α_N,vec = −GM·r_vec/r³    （1）

这里G为万有引力常数，r_vec为位置矢量。相应圆周运动满足V²=GM/r，外围速度随半径下降。

本文提出另一种星系本征加速度：

α_G,vec = −β√M·r_vec/(r²+c²)    （2）

式（2）中的M专指恒星和气体等重子物质的总质量，而不包含为解释运动而反推的暗晕质量。β是待检验的质量—场强耦合系数(本征常数)，c是特征长度。两式均直接给出加速度大小及其指向中心的方向，但质量标度与径向衰减不同。

必须强调：牛顿理论也可以通过分布式质量产生非平方反比的局域加速度。因此式（2）的意义不是宣称牛顿数学绝不容纳平坦旋转，而是提出一种无需事先指定暗晕质量分布的物理假说。

二、GRoDE提供径向结构

星系旋转动力学方程GRoDE写成：

r·dV/dr = V(1 − V²/A)    （3）

当A>0时，其非零解族为：

V²(r) = A·r²/(r²+c²)    （4）

对圆周运动，向心加速度大小为α(r)=V²(r)/r，于是：

α(r) = A·r/(r²+c²)    （5）

若把这个径向关系延拓为中心矢量场：

α_vec(r) = −A·r_vec/(r²+c²)    （6）

因此，式（2）的分母r²+c²来自GRoDE方程的解析解，而非在拟合完成后随意添加。但从盘面圆周速度推断三维球对称场，是额外的物理假设，不属于观测数据已经直接证实的内容。

三、质量平方根假说

GRoDE本身规定V(r)如何随尺度变化，却不规定场强A如何由星系总质量确定。本文独立提出：

A = β√M    （7）

这是本文的基本猜想，不强求先从牛顿力学推导，也不把Tully–Fisher关系当作定义。它需要由观测判断真伪。将式（7）代入式（6），得到：

α_vec(r) = −β√M·r_vec/(r²+c²)    （8）

这里最核心的命题是：星系整体本征场的强度参数与重子总质量的平方根成正比。本征常数β若跨星系近似不变，则它可能具有普适物理意义。

四、中心、外围和本征运动

当r远小于c时，加速度及圆周速度近似为：

α(r) ≈ (β√M/c²)·r；  V(r) ≈ [√(β√M)/c]·r    （9）

中心区域因而表现出近似线性的上升旋转曲线。当r远大于c时：

α(r) ≈ β√M/r；  V²(r) → β√M    （10）

由加速度对半径求导，还可得到：

dα/dr = β√M·(c²−r²)/(r²+c²)²    （11）

r = c 时，α_max = β√M/(2c)    （12）

因此c不仅是旋转曲线的特征尺度，也是候选本征加速度达到最大值的半径。

在本征力学表述中，可将无额外扰动时的运动写为a_vec=α_vec。于是式（8）可以用于研究非圆轨道和轨道稳定性，而不仅是圆周速度拟合。其中心势可写为：

Φ(r) = (β√M/2)·ln(r²+c²) + 常数    （13）

这个对数势在数学形式上已有先例；本征力学能否提供独立物理内容，取决于质量参数的决定方式以及能否预测额外观测现象。

五、无需预设暗物质晕

常见暗晕模型把总旋转速度分为恒星、气体与暗晕分量，再拟合暗晕密度轮廓。本文直接把式（8）作为总本征加速度的候选定律，圆周运动随之给出：

V²(r) = β√M·r²/(r²+c²)    （14）

计算式（14）无需在模型开始时假设暗晕密度与暗晕总质量。这一点并不等于已经排除暗物质的存在；透镜、星系团、宇宙学与三维垂直动力学仍是不可绕过的独立检验。

六、Tully–Fisher关系与β

由式（14）的远区极限，立刻有：

V_flat⁴ = β²M    （15）

这就是重子Tully–Fisher关系的标度形式。已有相关观测显示，质量与平坦旋转速度四次方之间存在紧密关系。在速度以km/s、质量以太阳质量M⊙计时，参考关系M/M⊙≈47[V_flat/(km/s)]⁴，给出β≈1/√47≈0.146（数值单位为(km/s)²·M⊙^(−1/2)）。

此前的探索性GRoDE计算中，DDO 154、NGC 6503、NGC 7793等星系分别报告β≈0.1435、0.1452、0.1511；但NGC 3521与NGC 2841报告值分别为0.2123、0.2579，显示明显偏离。另有部分数据佐证β的近宇宙常数性。

β约为0.15的线索令人关注，但不能据此宣布本征常数β已经成为新的宇宙常数。特别需要区分：由于β=V_flat²/√M，其近似集中性本来就与BTFR密切相关；真正独立的检验，应是提前固定β，再预测多颗星系的完整旋转曲线与非圆运动。

七、原创性的边界与真正需要回答的问题

已有MOND研究提出过远区加速度与√M/r成正比的形式，经典对数势也使用过相同的径向数学结构。因此，本文不把任何单一数学因子宣称为首次出现。本文试图建立的特色，是把GRoDE的径向动力学、A=β√M的场强猜想及本征运动概念联合起来，作为星系总加速度的可检验候选规律。

未来需要回答：① 是否存在跨星系普适的β？② c能否由独立可测的结构量决定？③ 完整旋转曲线、非圆轨道及盘外三维运动是否服从同一场？④ 与现有暗晕和修正动力学模型相比，是否能给出新的判别性预测？

八、结论

星系本征加速度的核心猜想，可凝练为：

α_vec(r) = −β√M·r_vec/(r²+c²)    （16）

该式把星系重子质量与空间加速度直接联系起来，GRoDE赋予其有限核心径向结构，外围自然产生V_flat⁴=β²M。它提供了一条无需预设暗物质晕的建模路线，但仍须经受完整旋转曲线、三维动力学及其他独立观测的严格检验。它的价值最终取决于自然界是否遵循这一简洁的质量—加速度规律。

参考文献

Milgrom, M. (1983). A modification of the Newtonian dynamics as a possible alternative to the hidden mass hypothesis. Astrophysical Journal, 270, 365.

McGaugh, S. S. (2012). The Baryonic Tully–Fisher Relation of Gas-Rich Galaxies as a Test of ΛCDM and MOND. Astronomical Journal, 143, 40.

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

支持单位：