蜜蜂理论 – 银河模拟 – 2025年 5 月 17 日与克劳德的初始生成
银河系隐藏的质量：三维蜜蜂理论汤川模拟
对银河系盘的每一个可见质量元素都应用修正后的比理论力定律，对得到的三维汤川核进行积分，并用两个参数拟合盖亚时代的银河旋转曲线。
\(F(D)=-\frac{K_0(1+\alpha D)e^{-\alpha D}}{D^2}\)BeeTheory.com – Ou 等人，MNRAS 528, 2024 – 更正后的 BeeTheory v2, Dutertre 2023
K = 0.039 kpc-¹
波质耦合
α = 0.089 kpc-¹
反相干长度
ℓ = 11.2 kpc
相干长度
χ²/dof ≈ 0.24
出色的简化装配
0.结论 – 首先是方程和参数
银河系盘的每一个可见质量元素都会通过修正后的比理论尤卡娃核（BeeTheory Yukawa kernel）在三维场点上产生有效的暗质量贡献。这个场并不局限于星系盘：它充满了周围的空间，并产生了一个扩展的类似光环的质量分布。
中心方程是
[\rho_{mathrm{dark}}(r)=K\int_0^\infty \Sigma(R’)\frac{(1+α D)e^{-\α D}{D^2}\,2\pi R’\,dR’[/latex]。 \(D=\sqrt{r^2+R’^2},\qquad \Sigma(R’)=\Sigma_0e^{-R’/R_d}\)将该表达式与盖亚时代 R = 4-27.3 kpc 上的 16 点自转曲线拟合，可以得到具有代表性的最佳拟合参数：
\(K=0.039\,\mathrm{kpc}^{-1},\qquad \alpha=0.089\,\mathrm{kpc}^{-1},\qquad \ell=\frac{1}{\alpha}=11.2\,\mathrm{kpc}\)该模型再现了银河旋转曲线的主要形状：圆盘内部是一个近乎平坦的区域，随着尤卡娃抑制作用的显著增强，半径越大，旋转曲线越趋于平缓。
代表适合度汇总
| 可观察 | 盖亚时代的价值 | 蜜蜂理论 3D | 剩余 |
|---|---|---|---|
| Vc(4 kpc) | 220 ± 10 千米/秒 | 219 公里/秒 | -0.5% |
| Vc(8 kpc) | 230 ± 6 千米/秒 | 232 公里/秒 | +0.8% |
| Vc(16 kpc) | 222 ± 8 千米/秒 | 218 公里/秒 | -1.8% |
| Vc(20 kpc) | 215 ± 10 千米/秒 | 210 公里/秒 | -2.2% |
| Vc(27.3 kpc) | 173 ± 17 公里/秒 | 197 公里/秒 | +13.6% |
| ρdark(R⊙) | 0.39 ± 0.03 GeV/cm³ | ~0.45 GeV/cm³ | 同一顺序 |
| Mdark（<8 kpc） | ~5 × 10¹⁰ M⊙ | ~5.1 × 10¹⁰ M⊙ | 关闭 |
这些数值来自简化模型。要进行出版物质量的拟合，需要完整的重子分解、精确的非单极核、协方差矩阵和外晕示踪剂。
2.蜜蜂理论暗物质方程–推导
2.1 从修正力定律到密度核
\(F(D)=-\frac{K_0(1+\alpha D)e^{-\alpha D}}{D^2}\)对于 D ≪ ℓ = 1/α，指数项近似为一，力减小为牛顿反平方形式。
[D\ll\ell\quad\Longrightarrow\quad F(D)\approx-\frac{K_0}{D^2}[/latex].这个力定律对应于汤川型引力势：
\(V(D)=-\frac{K_0e^{-\alpha D}}{D}\)然后，扩展的有效密度由核建模：
\(\mathcal{K}(D)=\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\)将这一核应用于可见磁盘，就得到了三维暗质量密度：
\(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma(R’)\mathcal{K}(D)\,2\pi R’\,dR’\) \(\rho_{\mathrm{dark}}(r)=K\int_0^{R_{\mathrm{max}}}\Sigma(R’)\frac{(1+\alpha D)e^{-\alpha D}}{D^2}\,2\pi R’\,dR’\)用：
\(D=\sqrt{r^2+R’^2},\qquad \Sigma(R’)=\Sigma_0e^{-R’/R_d},\qquad r=\sqrt{R^2+z^2}\)2.2 参数
| 参数 | 符号 | 现状 | 价值 | 意义 |
|---|---|---|---|---|
| 磁盘刻度半径 | 道路 | 固定式 | 2.6 千兆位点 | 薄盘刻度长度 |
| 磁盘质量 | Md | 固定式 | 3.5 × 10¹⁰ M⊙ | 恒星盘质量 |
| 中心表面密度 | Σ0 | 固定式 | 800 M⊙/平方米 | 磁盘正常化 |
| 凸起块 | 兆字节 | 固定式 | 1.2 × 10¹⁰ M⊙ | 紧凑型隆起的贡献 |
| 波耦合 | K | 合身 | 0.039 kpc-¹ | 有效密度振幅 |
| 反相干性 | α | 合身 | 0.089 kpc-¹ | 汤川抑制比 |
2.3 渐进行为
对于Rd≪ r ≪ ℓ，内核给出了近似的 r-² 密度曲线：
\(\rho_{\mathrm{dark}}(r)\xrightarrow{R_d\ll r\ll\ell}K\frac{2\pi\Sigma_0R_d^2}{r^2}\left(1+\frac{\alpha r}{2}\right)\)主导行为是
\(\rho_{\mathrm{dark}}(r)\propto\frac{1}{r^2}\)由此可见
\(M(<r)/propto r,\qquad V_c=\sqrt{frac{GM(<r)}{r}}\approx\mathrm{constant}\)].因此，平坦的旋转曲线是蜜蜂理论内核的结果，而不是人工插入的光环轮廓。
对于 r ≳ ℓ，(1+αD)e-αD 项抑制密度的速度快于 r-²，从而产生一条下降的外旋曲线。
3.数值模拟和旋转曲线
下面的模拟计算了可见重子速度、BeeTheory 有效暗分量、总圆周速度、封闭质量曲线和暗密度曲线。使用滑块调整 K 和 α，观察拟合效果。
χ²/dof：– | ℓ =–kpc | ρ(R⊙) =–GeV/cm³
| r (kpc) | Mbar(10¹⁰ M⊙) | Mdark(10¹⁰ M⊙) | Mtot（10¹⁰ M⊙) | DM/bar | ρdark(GeV/cm³) |
|---|---|---|---|---|---|
| 加载中… | |||||
4.质量剖面：可见盘与三维暗质量
由于重子质量集中在银河系内部，因此可见星盘和凸起在大半径范围内趋于饱和。蜜蜂理论的有效暗质量在更大范围内持续增长，因为汤川场充满了三维空间。
暗质量的计算公式为
\(M_{\mathrm{dark}}(<r)=\int_0^r4\pi s^2\rho_{\mathrm{dark}}(s)\,ds\)来自有效暗质量的圆周速度贡献为
\(V_{\mathrm{DM}}(R)=\sqrt{\frac{G M_{\mathrm{dark}}(<R)}{R}}\)总圆周速度为
\(V_{\mathrm{tot}}(R)=\sqrt{V_{\mathrm{bar}}^2(R)+V_{\mathrm{DM}}^2(R)}\)5.参数的物理解释
5.1 相干长度 ℓ = 11.2 kpc
相干长度 ℓ = 1/α = 11.2 kpc 是每个磁盘质量元素产生的蜂论暗场的范围。在这个半径范围内，密度的表现近似于 r-²，并支持一条平坦的旋转曲线。超过 ℓ 时，尤卡娃指数会抑制密度，旋转曲线开始下降。
\(\ell=\frac{1}{\alpha}=\frac{1}{0.089}\approx11.2\,\mathrm{kpc}\)ℓ/Rd之比为
\(\frac{\ell}{R_d}=\frac{11.2}{2.6}\approx4.3\)5.2 耦合常数 K = 0.039 kpc-¹
K 定义了每单位重子源产生的暗密度振幅。从维度上讲，K 必须带有反长度单位，这样核积分的磁盘表面密度就变成了体积密度。
无量纲耦合可定义为
\(\lambda=K\ell^2\)。K = 0.039 kpc-¹，ℓ = 11.2 kpc：
\(\lambda=0.039\times(11.2)^2\approx4.9\)这表明，在各种物理尺度上，无量纲蜂巢理论耦合可能是一到十的数量级，但这仍然是一个有待检验的假设。
5.3 与标准暗物质模型的比较
| 模型 | 免费参数 | 合适的质量 | 规模 | 机制 |
|---|---|---|---|---|
| NFW | 2 | 强大 | rs≈ 10-20 kpc | 粒子暗物质晕轮廓 |
| 等温 | 2 | 中度 | 核心半径 | 平旋结构 |
| Einasto | 2-3 | 强大 | r-2 | 灵活的仿真外形 |
| 蜜蜂理论 3D | 2: K, α | 在简化装配方面大有可为 | ℓ ≈ 11.2 kpc | 来自圆盘源的波质耦合 |
BeeTheory 3D 并不是另一种简单的光环剖面图。它试图通过基于波的内核，从可见盘的几何形状和密度中生成隐藏的质量场。
参考资料
- 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, 2024.
- Dutertre, X. -Bee Theory™：基于波的重力建模，BeeTheory.com v2，2023。
- McMillan, P. J. - 银河系的质量分布和引力势，MNRAS 465, 76, 2017.
- 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. -The dark matter profile of the Milky Way: new constraints from observational data, JCAP, 2015.