蜜蜂理论 – 科学推导 – 2025年
两个氢原子的波函数：严格推导和校准
从 “蜜蜂理论”（BeeTheory）的指数波函数假设出发，我们推导出了精确的三维相互作用能量，修正了原来的单极近似值，并用两个参数对已知的 H₂ 分子进行了校准，这两个参数与实验的重现率小于 0.2%。
BeeTheory.com – 基于 BeeTheory v2 (Dutertre, 2023) – 经过扩展和修正
κ = 3.509Eh
波质耦合
αeff= 1.727a0
有效波浪范围
要求= 74.2 pm
与实验对比74.1 pm
De= 4.517 eV
对比实验：4.52 eV
0.结论–结果第一
蜜蜂理论 “基于波的模型用球形波函数表示每个氢原子：
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)当两个原子以 R 的间隔相互作用时，该模型会产生有效的吸引力相互作用能，经过全三维积分后，其精确形式为尤卡瓦型势能：
\(E_{\mathrm{att}}(R)=-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}\)结合以原子为单位的核排斥力，这个双参数模型在校准实验数据后再现了氢₂分子的平衡距离和解离能。
蜜蜂理论》原论文的关键结果得到了证实：波的相互作用产生了吸引力。不过，单极近似在这里得到了修正，因为它失去了 R 依赖性。修正后的模型给出了具有校准系数的汤川形式。
\(E(R)=\underbrace{-\frac{\kappa}{\sqrt{\pi}}e^{-R/\alpha_{\mathrm{eff}}}}_{\text{wave attraction}}+\underbrace{\frac{e^2}{4\pi\varepsilon_0R}}_{\text{nuclear repulsion}}\) \(\kappa=3.509E_h,\qquad \alpha_{\mathrm{eff}}=1.727a_0,\qquad a_0=52.92\,\mathrm{pm},\qquad E_h=27.21\,\mathrm{eV}\)κ = 3.509Eh
相当于 95.5 eV。设置吸引力相互作用的振幅。
αeff= 1.727a0
相当于 91.4 pm。这比玻尔半径大 72.7%。
<0.2% 误差
Req= 74.16 pm 和De= 4.517 eV，与实验相符。
1.波函数：精确的三维形式
1.1 蜜蜂理论的起始假设
每个基本粒子都有一个波函数模型，该波函数从其中心开始在三个空间方向上呈指数衰减。对于处于基态的氢原子来说，这不仅仅是一个假设，而是一个精确的量子力学结果：”蜜蜂理论“的波函数与氢原子的 1s 轨道相吻合。
\(\psi_{1s}(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\exp\left(-\frac{r}{a_0}\right),\qquad r=|\mathbf{r}|\)用 α =1/a0 的紧凑符号表示：
\(\psi(r)=\frac{\alpha^{3/2}}{\sqrt{\pi}}e^{-\alpha r}=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)1.2 标准化 – 精确验证
\(\int_0^\infty|\psi(r)|^2\,4\pi r^2\,dr=\frac{4\alpha^3}{\pi}\cdot\pi\int_0^\infty r^2e^{-2\alpha r}\,dr=\frac{4\alpha^3}{1}\cdot\frac{2}{(2\alpha)^3}=1\)1.3 能量 – 薛定谔方程验证
应用与时间无关的薛定谔方程
\(\hat{H}\psi=E\psi\) \(\hat{H}=-\frac{\hbar^2}{2m_e}\nabla^2+V(r),\qquad V(r)=-\frac{e^2}{4\pi\varepsilon_0r}\)在球面坐标中，exp(-αr) 的精确拉普拉奇值为
\(\nabla^2\left(e^{-\alpha r}\right)=\frac{d^2}{dr^2}\left(e^{-\alpha r}\right)+\frac{2}{r}\frac{d}{dr}\left(e^{-\alpha r}\right)=e^{-\alpha r}\left(\alpha^2-\frac{2\alpha}{r}\right)\)对《蜜蜂理论》论文的更正
最初的近似值∇²f(r) ≈ -3α/RAB摒弃了径向依赖性。精确拉普拉斯有两个项：α²e-αr 和 –2αe-αr/r。修正后的推导保留了这两个项。
以原子为单位，ħ =me= e = 1，a0= 1：
\(\nabla^2\psi=\psi(r)\left(1-\frac{2}{r}\right)\) \(T\psi=-\frac{1}{2}\nabla^2\psi=\psi\left(\frac{1}{r}-\frac{1}{2}\right)\) \(V\psi=-\frac{1}{r}\psi\) \((T+V)\psi=\psi\left(\frac{1}{r}-\frac{1}{2}-\frac{1}{r}\right)=-\frac{1}{2}\psi\) \(E_{1s}=-\frac{1}{2}E_h=-13.6057\,\mathrm{eV}\)2.两个波函数之和 – 精确法
将原子 A 放置在原点，原子 B 放置在 Z 轴上的 R 位置。蜂论叠加的总波函数为
\(\Psi(\mathbf{r})=\psi_A(\mathbf{r})+\psi_B(\mathbf{r})=\frac{1}{\sqrt{\pi a_0^3}}\left[e^{-|\mathbf{r}|/a_0}+e^{-|\mathbf{r}-\mathbf{R}|/a_0}\right]\)2.1 在 B 附近评估 A 的波函数
在原子 B 附近，A 的波的贡献为
\(\psi_A\big|_{\mathrm{near}\ B}=\frac{1}{\sqrt{\pi a_0^3}}e^{-|\mathbf{R}+\mathbf{r}|/a_0}\approx\underbrace{\frac{1}{\sqrt{\pi a_0^3}}e^{-R/a_0}}_{C_A(R)}e^{-r/a_0}\)CA(R)的振幅随着原子间距的增加呈指数衰减。它是从原子 A 传到原子 B 的蜂论信号。
| R | CA(R)/N =e-R/a₀ | 物理意义 |
|---|---|---|
| 0.5a0 | 0.607 | 强重叠、排斥机制 |
| 1.0a0 | 0.368 | 在玻尔半径处 |
| 1.4a0 | 0.247 | 近 H₂ 键长度 |
| 2.0a0 | 0.135 | 仍然重要 |
| 3.0a0 | 0.050 | 弱相互作用机制 |
| 5.0a0 | 0.007 | 互动几乎为零 |
2.2 应用于交叉项的哈密顿方程
在 B 附近，有效的局部波为
\(\Psi_{\mathrm{local}}(r)\approx[C_A(R)+N]e^{-r/a_0}\)将动力学算子应用于 A 部分，可以得出
\(\hat{T}\left[C_A(R)e^{-r}\right]=-\frac{1}{2}C_A(R)\nabla^2(e^{-r})\) \(=C_A(R)e^{-r}\left(\frac{1}{r}-\frac{1}{2}\right)\)动能算子的 1/r 项与库仑势成对，产生有效吸引力。
\(\langle\psi_B|e^{-r}/r|\psi_B\rangle=\frac{4}{9}\) \(\langle\psi_B|e^{-r}|\psi_B\rangle=\frac{8}{27}\) \(E_{\mathrm{BT,kin}}(R)=C_A(R)\left[\frac{4}{9}-\frac{1}{2}\cdot\frac{8}{27}\right]=C_A(R)\frac{8}{27}\)3.从动力学耦合到相互作用潜能
3.1 完整的蜜蜂理论互动
原子 A 和原子 B 之间的蜂论相互作用来自 A 的波场与 B 的电子密度的动能耦合。结合核排斥力，总的相互作用能的形式为
\(E_{\mathrm{BT}}(R)=-\kappa\frac{e^{-R/\alpha_{\mathrm{eff}}}}{\sqrt{\pi}}+\frac{1}{R}\)负项是吸引力，1/R 项是核排斥力。两个参数控制着相互作用：κ 和 αeff。
3.2 与原论文的比较
原始近似值
\(\nabla^2f\approx-\frac{3\alpha}{R_{AB}}\)这就失去了相互作用的 R 依赖性，无法产生平衡距离。
校正后的精确拉普拉斯
\(\nabla^2e^{-r}=e^{-r}\left(1-\frac{2}{r}\right)\)这保留了完全的 r 依赖性，并产生了尤卡瓦相互作用。
3.3 为何是汤川势能而非库仑势能
系数e-R/αeff来自 A 的波在 B 的位置上的振幅。在距离较大时，相互作用呈指数衰减。这使得原子尺度的蜜蜂理论相互作用成为一种有限范围的尤卡娃势能。
\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\)在 H₂ 键长度上，吸引力和排斥力达到平衡。
4.校准：两个条件，两个参数
正好有两个自由参数：κ 和 αeff，以及两个来自 H₂ 分子的实验约束。
| 制约因素 | 物理意义 | 数学条件 | 实验值 |
|---|---|---|---|
| 要求 | 键合长度 | dE/dR = 0 | 74.14 pm = 1.401a0 |
| 德 | 解离能 | E(∞) – E(Req) =De | 4.520 eV = 0.1660Eh |
4.1 分析解决方案
条件 1：
\(\frac{dE}{dR}=0\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}\alpha}=\frac{1}{R_{\mathrm{eq}}^2}\)条件 2：
\(E(\infty)-E(R_{\mathrm{eq}})=D_e\quad\Longrightarrow\quad\frac{\kappa e^{-R_{\mathrm{eq}}/\alpha}}{\sqrt{\pi}}=\frac{1}{R_{\mathrm{eq}}}+D_e\)用条件 2 除以条件 1：
\(\alpha=R_{\mathrm{eq}}+D_eR_{\mathrm{eq}}^2\)Req= 1.4014a0，De= 0.1660Eh：
\(\alpha_{\mathrm{eff}}=1.4014+0.1660(1.4014)^2=1.7274a_0\)那么
\(\kappa=\left(\frac{1}{R_{\mathrm{eq}}}+D_e\right)\sqrt{\pi}e^{R_{\mathrm{eq}}/\alpha_{\mathrm{eff}}}=3.509E_h\) \(\boxed{\kappa=3.509E_h=95.5\,\mathrm{eV},\qquad \alpha_{\mathrm{eff}}=1.727a_0=91.4\,\mathrm{pm}}\)4.2 参数的物理解释
| 参数 | 价值 | 蜂论中的物理意义 |
|---|---|---|
| κ | 3.509Eh | 波质耦合振幅。 |
| αeff | 1.727a0 | 相互作用的有效衰变长度。 |
| αeff/a0 | 1.727 | 蜜蜂理论杂交率 |
5.势能曲线及与实验的比较
建议的图表：比较蜜蜂理论、海特勒-伦敦理论和实验参考数据的 H₂ 势能曲线。
备选文字：H₂势能曲线H₂ 势能曲线，横轴为距离 R（以埃为单位），纵轴为能量（以电子伏特为单位）。蜜蜂理论曲线在 R = 0.74 Å 附近达到最小值-4.52 eV，与实验中的₂键距离和解离能相吻合。
| R (a0) | R（下午） | 电波 | 教育 | EBT | EBT(eV) | 现状 |
|---|---|---|---|---|---|---|
| 0.50 | 26.5 | -1.482 | +2.000 | +0.518 | +14.09 | 歪瓜裂枣 |
| 0.80 | 42.3 | -1.246 | +1.250 | +0.004 | +0.11 | 接近于零 |
| 1.00 | 52.9 | -1.110 | +1.000 | -0.110 | -2.98 | 招人喜欢 |
| 1.20 | 63.5 | -0.988 | +0.833 | -0.155 | -4.22 | 招人喜欢 |
| 1.401 | 74.1 | -0.880 | +0.714 | -0.166 | -4.517 | 最低限度 |
| 1.60 | 84.7 | -0.784 | +0.625 | -0.159 | -4.33 | 浅井 |
| 2.00 | 105.8 | -0.622 | +0.500 | -0.122 | -3.32 | 上升的 |
| 3.00 | 158.8 | -0.349 | +0.333 | -0.015 | -0.42 | 接近于零 |
| 5.00 | 264.6 | -0.110 | +0.200 | +0.090 | +2.46 | 斥力尾 |
蜜蜂理论根据校准结构，Req= 74.2 pm，De= 4.52 eV。
Heitler-London:预测键长较大，解离能较低。
实验：Req= 74.14 pm，De= 4.520 eV。
6.完整方程 – 随时可用
6.1 波函数
\(\psi(r)=\frac{1}{\sqrt{\pi a_0^3}}e^{-r/a_0}\)6.2 精确拉普拉斯
\(\nabla^2\psi(r)=\psi(r)\left(\frac{1}{a_0^2}-\frac{2}{a_0r}\right)\)6.3 互动总能量
\(E(R)=-\frac{\kappa}{\sqrt{\pi}}\exp\left(-\frac{R}{\alpha_{\mathrm{eff}}}\right)+\frac{e^2}{4\pi\varepsilon_0R}\) \(E(R)=-\frac{3.509}{\sqrt{\pi}}e^{-R/1.727}+\frac{1}{R}\) \(E(R)=-\frac{3.509E_h}{\sqrt{\pi}}\exp\left(-\frac{R}{1.727a_0}\right)+\frac{e^2}{4\pi\varepsilon_0R}\)6.4 两个氢原子之间的作用力
\(F(R)=-\frac{dE}{dR}=-\frac{\kappa}{\sqrt{\pi}\alpha_{\mathrm{eff}}}e^{-R/\alpha_{\mathrm{eff}}}+\frac{1}{R^2}\) \(F(R)=-\frac{3.509}{\sqrt{\pi}\times1.727}e^{-R/1.727}+\frac{1}{R^2}\)6.5 参数汇总表
| 符号 | 名称 | 价值 | 如何确定 |
|---|---|---|---|
| a0 | 玻尔半径 | 52.918 下午 | 氢量子力学 |
| Eh | 哈特里 | 27.211 eV | 原子单位定义 |
| α | 波浪衰减常数 | 1/a0 | 氢 1s 轨道 |
| κ | 波质耦合 | 3.509Eh | 按要求和去校准 |
| αeff | 有效衰变长度 | 1.727a0 | 根据 H₂ 进行校准 |
| 要求 | 平衡键长度 | 74.14 下午 | 实验 |
| 德 | 解离能 | 4.520 eV | 实验 |
7.未决问题和下一步推导
从 H₂ 到引力–蜂论的缩放问题
在原子尺度上，BeeTheory以κ = 3.509 Eh和αeff = 1.727 a0重现了H₂的化学性质。在星系尺度上，BeeTheory 使用的相干长度是以千帕秒为单位的。悬而未决的问题是相干长度如何从原子系统扩展到天体物理系统。
下一步推导：氦原子和多电子原子
对于氦，波函数可近似为
\(\psi_{\mathrm{He}}(r)=Ne^{-\alpha_{\mathrm{He}}r}\)下一步自然是针对 He₂ 范德瓦耳斯相互作用测试 BeeTheory。
扩展：非相同原子
对于具有不同衰变常数的原子 A 和 B，一般的 “蜜蜂理论 “相互作用可写成
\(E(R)=-\kappa_{AB}\frac{e^{-R/\alpha_{AB}}}{\sqrt{\pi}}+\frac{Z_AZ_B}{R}\)参考资料
- Dutertre, X. –Bee Theory™：基于波的重力建模，BeeTheory.com v2，2023。
- Heitler, W., London, F. –Wechselwirkung neutraler Atome und homöopolare Bindung nach der Quantenmechanik, Z. Physik 44, 455, 1927.
- Kolos, W., Wolniewicz, L. –氢分子 X¹Σg⁺、b³Σu⁺ 和 C¹Πu 态的势能曲线，J. Chem.43, 2429, 1965.
- Herzberg, G. –The Dissociation Energy of the Hydrogen Moleule, J. Mol.Spectrosc.33, 147, 1970.
- Slater, J. C. –Atomic Shielding Constants, Phys. Rev. 36, 57, 1930.
- Atkins, P. W., Friedman, R. –Molecular Quantum Mechanics, 5th ed., Oxford University Press, 2011.