蜜蜂理论 – 科学文章 – 2025

银河系的隐藏质量：
基于波的推导和数值拟合

从 “蜜蜂理论”（BeeTheory）的假设出发，即每个质量元素都会发射出一个以$e^{-D/\ell}$衰减的引力波场，我们分析推导出了三维暗质量分布，并将其与盖亚2024的旋转曲线拟合，找到了模型的两个基本参数。

BeeTheory.com – Technoplane S.A.S. – Based on Ou et al.

$\rho_0 = 1.14\;\text{GeV/cm}^3$

中心暗密度

$r_s = 9.6\;\text{kpc}$

波相干尺度

$\chi^2/\text{dof} = 0.44$

拟合优度

$0.41\;\text{GeV/cm}^3$

预测的 $\rho_\text{dark}(R_\odot)$

$\sim 5\times 10^{11}\,M_\odot$

200 kpc以内的暗质量总量

0.结论–结果第一

基于 “蜜蜂理论 “的波模型–其中每个可见质量元素$dV$都会产生一个引力场，该引力场在三维空间中以$e^{-D/ell}$的指数形式衰减–预言了一个暗质量密度曲线，该曲线在银河系圆盘上进行积分后，趋近于NFW形式。

仅使用两个自由参数拟合银河系的盖亚2024旋转曲线，该模型达到了$chi^2/\mathrm{dof} = 0.44$。

最佳拟合参数是：中心暗密度 $\rho_0 = 1.14\,\mathrm{GeV/cm}^3$和相干尺度半径 $r_s = 9.6\,\mathrm{kpc}$。这些直接映射到 “蜜蜂理论 “的两个参数：波耦合常数 $\lambda$ 和相干长度 $\ell = r_s\sqrt{2}$ 。\大约13.6,mathrm{kpc}$。

该模型预测的本地暗物质密度为：$\rrh_\text{dark}(R_\odot = 8\,\mathrm{kpc}) = 0.41\,\mathrm{GeV/cm}^3$ – 在测量值$0.39 \pm 0.03\,\mathrm{GeV/cm}^3$的5%以内。200 kpc范围内的总暗物质质量为$\sim 7.1 \times 10^{11}\,M_\odot$，与最近的卫星运动学测量结果一致。

中心暗密度

$\rho_0 = 1.14\;\frac\{text{GeV}}\{text{cm}^3}$

相当于3.0times10^7,M_odot,text{kpc}^{-3}$。在 $r=0$ 时的波场振幅。

波相干尺度

$r_s = 9.6\;text\{kpc}$

波场从内部机制过渡到外部机制的尺度。

拟合优度

$chi^2/\text{dof} = 0.44$

拟合效果极佳。16 个数据点中有 15 个在 $1\sigma$范围内。

可观测	盖亚2024测量	蜜蜂理论预测	残差
$V_c(R_odot = 8\\text{kpc})$	$230 /pm 6\;\text{km/s}$	$231\\text{km/s}$	$+0.4\%$
$V_c(20\,\text{kpc})$	$215 /pm 10\;\text{km/s}$	$208\\text{km/s}$	$-3.3\%$
$V_c(27.3\,\text{kpc})$	$173pm 17\\text{km/s}$	$199\\text{km/s}$	$+15\%$, $1.5\sigma$
$rho_\text{dark}(R_\odot)$	$0.39 \pm 0.03\text{GeV/cm}^3$	$0.41\;\text{GeV/cm}^3$	$+5\%$
$M_\text{dark}(<8\,\text{kpc})$	$\sim 5\times10^{10}\,M_\odot$	$5.1\times10^{10}\,M_\odot$	$+2\%$
$M_\text{dark}(<200\,\text{kpc})$	$\sim(5\text{–}9)\times10^{11}\,M_\odot$	$7.1\times10^{11}\,M_\odot$	在范围内

1.蜜蜂理论假设：质量辐射波

经典万有引力和相对论万有引力描述了万有引力的作用方式，但没有描述万有引力存在的原因。蜜蜂理论提出了一种机制：每个质量元素 $dV$ 都是量子波场的源头，它在三维空间中向外传播，并随着与源头的欧几里得距离 $D$ 呈指数衰减。

这个波场携带着有效的引力能量–确切地说，它是 “隐藏的质量”。

蜜蜂理论波质量公设 $$d\rho_\text{wave}(\mathbf{r}) = \frac{\lambda}{\ell}\;\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV$$

这里，$\lambda$是波-质耦合常数，$\ell$是相干长度，$\rrho_\text{vis}$是可见重子质量密度，$D = \mathbf{r}-\mathbf{r}’|$是源到场点的欧几里得距离。

任意点的总暗物质密度 $\mathbf{r}$ 是星系中每个可见质量元素的波场的叠加：

暗物质总密度-叠加积分 $$\rho_\text{dark}(\mathbf{r}) = \frac\{lambda}{\ell}\（星系）\rho_\text{vis}(\mathbf{r}’)\;\exp\!\left(-\frac{|\mathbf{r}-\mathbf{r}’|}{\ell}\right) dV’$$

这是可见质量分布与指数核的三维卷积。

根据假设，暗质量不是球形对称的。它反映了源的几何形状。
暗质量充满了整个三维空间，而不仅仅是银河系平面。
一旦知道重子分布，两个参数$(\lambda,\ell)$就能完全决定暗质量分布。

2.可见源：指数盘

银河系的恒星盘可以很好地用指数表面密度来描述：

星盘表面密度 $$\Sigma(R) = \Sigma_0\,e^{-R/R_d}, \qquad \Sigma_0 = 800\,M_\odot\,\text{pc}^{-2}, \quad R_d = 2.6\,\text{kpc}$$

与半径相比，圆盘的厚度可以忽略不计，因此其体积密度可以用圆盘表面密度和垂直三角函数来表示。

来自薄圆盘的暗密度–精确双积分 $$\rho_\text{dark}(R,z) = \frac{\lambda}{\ell}\int_0^\infty\!\int_0^{2\pi} \Sigma(R’)\,\exp\!\left(-\frac{\sqrt{R^2+R’^2-2RR’\cos\phi+z^2}}{\ell}\right) R’\,d\phi\,dR’$$

场点 $(R,z)$ 位于圆盘平面的圆柱半径 $R$ 及其上方高度 $z$。设定 $r = \sqrt{R^2+z^2}$，我们使用单极近似法进行方位角积分：

单极核-方位角平均值 $$K_\phi(r,R’) \equiv \int_0^{2\pi} e^{-D/\ell}\,d\phi \;\approx\; \frac{2\pi\ell}{r}\,\sinh\!\left(\frac{r}{\ell}\right)\exp\!\left(-\frac{r+R’}{\ell}\right)$$

这就把双积分降到了一维：

一维主积分 $$\rho_\text{dark}(r) = \frac{\lambda\Sigma_0}{\ell}\int_0^\infty R’\,e^{-R’/R_d}\cdot\frac{2\pi\ell}{r}\,\sinh\!\!\left(\frac{r}{\ell}\right)e^{-(r+R’)/\ell}\,dR’$$

2.1 分析结果 – NFW 的出现

分析$R’$积分可以得到

闭式蜜蜂理论暗密度 $$\rho_\text{dark}(r) = \frac{2\pi\lambda\Sigma_0\,R_d^2}{r}\cdot \frac{\ell^2}{(R_d+\ell)^2}\cdot \sinh\!\!left(\frac{r}\{ell}\right) e^{-r/\ell}$$

在内部系统中

内部机制 $$sinh(r/\ell)\,e^{-r/\ell}#frac{r}{\ell}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto r^{-1}$$

在外层系统中

外部制度 $$sinh(r/\ell)\,e^{-r/\ell}\approx \tfrac{1}{2}\quad\Longrightarrow\quad \rho_\text{dark}(r)\propto \frac{e^{-r/\ell}}{r}$$

这两种状态之间的过渡发生在 $r（sim/ell$）处。这是可以将蜂论波谱与 NFW 尺度行为进行比较的区域。

关键理论结果

类似于NFW的暗物质剖面是根据应用于指数盘源的蜜蜂理论波-质量公设分析得出的。在这种解释中，NFW 参数并不是任意的光环参数；它们与蜜蜂理论波参数和圆盘几何相关联。

2.2 蜂论-NFW 词典

蜂论到 NFW 参数的映射 $$r_s = \ell, \qquad \rho_0^text{NFW} = \frac{2\pi\lambda\Sigma_0 R_d^2}{r_s}\cdot\frac{1}{(1+R_d/r_s)^2}$$

蜜蜂理论解释的拟合参数 \ell = r_s = 9.6\,\text{kpc}, \qquad \lambda = \frac\{rho_0 (R_d+r_s)^2}{2\pi\Sigma_0 R_d^2}$$

蜂论数值参数 $$\boxed{ell = 9.6\,\text{kpc}, \qquad \lambda = \frac{3.0\times10^7 \times (12.2)^2}{2\pi \times 8\times10^8 \times 6.76} = 0.132}$$

蜜蜂理论参数	物理意义	最佳拟合值	来自
$\ell$	相干长度。等于 NFW 尺度半径 $r_s$.	$9.6\,\text{kpc}$	$V_c(R)$ 下降的形状
$\lambda$	无量纲波质耦合	$0.132$	绝对速度尺度
$\rho_0$	在 $r=0$ 时的峰值暗质量密度。	$1.14\,\text{GeV/cm}^3$	由 $\lambda$ 和 $\ell$ 计算得出
$r_s$	密度斜率之间的过渡半径.	$9.6\text{kpc}$	与 $\ell$ 相同

3.从质量缺失到旋转曲线

3.1 质量缺失问题

牛顿动力学要求

封闭总质量 $$M_\text{tot}(<R) = M_\text{bar}(<R) + M_\text{dark}(<R), \qquad M_\text{dark}(<R) = \frac{V_c^2 R}{G} – M_\text{bar}(<R)$$

半径 $R$ 内的重子质量有两个部分：指数盘和紧凑的隆起。

重子封闭质量 $$M_text{disk}(<R) = 2\pi\Sigma_0 R_d^2\!\left[1 – \left(1+frac{R}{R_d}\right)e^{-R/R_d}\right], \qquad M_text{bulge} = 1.2\times10^{10}\,M_\odot$$

3.2 圆周速度分解

圆周速度分解 $$V_c^2(R) = V_\text{disk}^2(R) + V_\text{bulge}^2(R) + V_\text{dark}^2(R)$$

磁盘贡献使用的是修正贝塞尔函数的弗里曼公式：

弗里曼圆盘速度 $$V_text{disk}^2(R) = \frac{2\,G\,M_d}{R_d}\,y^2\!\left[I_0(y)\,K_0(y) – I_1(y)\,K_1(y)\right], \quad y = \frac{R}{2R_d}$$

隆起使用的是 Hernquist 曲线：

亨奎斯特凸起贡献 $$V_\text{bulge}^2(R)=\frac{G\,M_b\,R}{(R+a)^2},\qquad a=0.6\,\text{kpc}$$

包含在 $R$ 中的 NFW 暗质量有一个解析形式：

NFW 包裹暗质量 $$M_\text{dark,NFW}(<R) = 4\pi\,\rho_0\,r_s^3\!\left[\ln\!\left(1+\frac{R}{r_s}\right) -\frac{R/r_s}{1+R/r_s}\right]$$

为什么仅重子预测的速度是错误的

当 $M_d = 3.5times10^{10}\,M_odot$和 $M_b = 1.2times10^{10}\,M_odot$时，重子模型预测的速度约为162text{km/s}$，接近$8text{kpc}$，低于观测到的速度$sim230text{km/s}$。

4.数值模拟和参数拟合

4.1 输入数据 – Gaia 2024

来自 Ou 等人（2024 年）的 16 个数据点，跨度为 $R=4$-$27.3\,\text{kpc}$：

const OBS_R = [4, 5, 6, 7, 8, 9, 10, 11, 12, 14, 16, 18, 20, 22, 24, 27.3]；
const OBS_V = [220,228,232,231,230,229,228,227,226,224,222,219,215,208,200,173]；
const OBS_ERR = [10,8,7,7,6,6,6,7,7,8,9,10,11,13,17]；

4.2 算法

计算 $V_text{bar}(R)$

使用弗里曼圆盘 + 亨奎斯特凸起。贝塞尔函数使用多项式近似计算。

评估 NFW 暗速度

使用闭合形式的 NFW 封闭质量，计算 $V_\text{dark}(R)$。

计算总速度和 $\chi^2$

$V_\text{tot}(R)=\sqrt{V_\text{bar}^2+V_\text{dark}^2}$.

最小化 $\chi^2(\rho_0,r_s)$

在 $\rho_0$ 和 $r_s$ 上使用两道网格。

关键单元说明

牛顿常数的正确值是 kpc-km-s-$M_\odot$ 单位：

$$G = 4.302\times10^{-6}\,\text{kpc}\,\text{km}^2\,\text{s}^{-2}\,M_\odot^{-1}$$

使用$4.302\times10^{-3}$是一个常见的单位误差，它给出的速度太大了。

4.3 交互式旋转曲线

银河旋转曲线 – Gaia 2024 与 BeeTheory模型对比

仅重子蜂巢理论 $V_\text{total}$ 仅暗物质盖亚2024数据

参数资源管理器 – 调整 $\rho_0$ 和 $r_s$

$\rho_0$ 30 $10^6\,M_\odot\,\text{kpc}^{-3}$

$r_s$ 9.6 kpc

$chi^2/\text{dof}$： – | $\rho_\text{dark}(8\,\text{kpc})$: – GeV/cm³

4.4 结果–三维质量曲线

半径为 $r$ 的圆球内的暗物质质量在 $r_s$ 内陡然上升，在 $r_s$ 外则呈对数增长。

质量分布：可见盘 vs 全部 vs 暗物质

可见盘 + 隆起暗物质总质量

$r$	$M_\text{bar}(<r)$	$M_\text{dark}(<r)$	$M_\text{tot}(<r)$	DM/bar 比率	$V_c$
5 kpc	$3.2\times10^{10}\,M_\odot$	$2.6\times10^{10}\,M_\odot$	$5.7\times10^{10}\,M_\odot$	0.81	229 km/s
8 kpc	$4.0\times10^{10}\,M_\odot$	$5.1\times10^{10}\,M_\odot$	$9.0\times10^{10}\,M_\odot$	1.28	231 km/s
15 kpc	$4.5\times10^{10}\,M_\odot$	$1.1\times10^{11}\,M_\odot$	$1.56\times10^{11}\,M_\odot$	2.44	216 km/s
30 kpc	$4.6\times10^{10}\,M_\odot$	$2.2\times10^{11}\,M_\odot$	$2.66\times10^{11}\,M_\odot$	4.78	196 km/s
100 kpc	$4.6\times10^{10}\,M_\odot$	$5.1\times10^{11}\,M_\odot$	$5.54\times10^{11}\,M_\odot$	11.1	154 km/s
200 kpc	$4.6\times10^{10}\,M_\odot$	$7.1\times10^{11}\,M_\odot$	$7.56\times10^{11}\,M_\odot$	15.4	128 km/s

5.两个参数的物理解释

5.1 相干长度 $\ell = r_s = 9.6\,\text{kpc}$

$\ell$ 是每个质量元素发射的引力波场保持相位的范围。在这个半径内，波的干涉是建设性的，暗密度会缓慢下降。在这个半径之外，破坏性干扰会导致密度下降得更快。

$\ell = 9.6\,\text{kpc}\approx 3.7R_d$ 这个值有一个自然的解释：相干长度由圆盘尺度半径乘以一个数量级为一的因子设定。

5.2 耦合常数 $\lambda = 0.132$

$\lambda$ 决定了每相干长度下每单位可见质量产生多少波质量。

从 $\lambda$ 得出的局部暗质量与可见光质量之比 $$\frac{\rho_\text{dark}(R_\odot)}{\rho_\text{vis}(R_\odot)} \approx \lambda\cdot\frac{\pi\ell}{R_\odot}\cdot\frac{R_d^2}{(R_d+\ell)^2/\ell}\大约 4.2$$

200 kpc范围内的全局暗质量与重子质量之比大约为$M_\text{dark}/M_\text{bar}\approx15$，这与巨大的隐质量成分是一致的。

蜜蜂理论预测：光环形状

由于暗质量是通过三维卷积从盘源中产生的，因此光环并不是完美的球形。精确的非单极计算预测银河系暗晕的小轴与主轴比约为$q=c/a/approx0.82$。

6.总结与展望

蜜蜂理论从一个单一的物理假设出发–即每个可见质量元素都会产生一个引力波场，在三维空间中以$e^{-D/\ell}$的形式衰减–提供了一个基于波的暗物质类密度曲线推导。

该模型拟合了盖亚2024年的银河旋转曲线，在两个自由参数的作用下达到了$chi^2/\text{dof}=0.44$：

最佳拟合蜜蜂理论参数–银河系 $$\ell = r_s = 9.6\,\text{kpc},\qquad \lambda = 0.132$$ $$\Longrightarrow\quad \rho_\text{dark}(R_\odot)=0.41\,\text{GeV/cm}^3,\qquad M_\text{dark}(<200\,\text{kpc})=7.1\times10^{11}\,M_\odot$$

除了旋转曲线之外，该模型还做出了三个可检验的预测：

光环形状：暗质量是圆盘扁平的，轴比为$q（约0.82$）。
参数普遍性：同样的$(\lambda,\ell)$关系应该适用于具有已知圆盘参数的外部星系。
相干缩放：$\ell\approx3.7R_d$表明圆盘大小与暗晕尺度半径之间存在缩放关系。

参考文献

Ou, X., Eilers, A.-C., Necib, L., Frebel, A. –The dark matter profile of the Milky Way inferred from its circular velocity curve, MNRAS 528, 693-710 (2024)
Navarro, J. F., Frenk, C. S., White, S. D. M. –A Universal Density Profile from Hierarchical Clustering, ApJ 490, 493 (1997)
Freeman, K. C. –On the disks of spiral and S0 galaxies, ApJ 160, 811 (1970)
Pato, M., Iocco, F., Bertone, G. –银河系暗物质分布的动态约束，JCAP 12, 001 (2015)
McMillan, P. J. – 银河系的质量分布和引力势能，《英国国家科学院院刊》（MNRAS）465, 76 (2017)
Abramowitz, M., Stegun, I. A. –Handbook of Mathematical Functions, Dover (1972)

银河系的隐藏质量：基于波的推导和数值拟合