大一EPQ核聚变项目：求证氢原子聚变克服库仑势垒的能量方程

Hey there! 作为大一就着手EPQ核聚变研究项目的同学，先给你点个大大的赞——这个方向真的超有意思！针对你问的氢原子聚变克服库仑势垒的能量方程，我来给你梳理下标准的表述和推导思路，也给你一些验证自己推导的小建议：

1. 经典库仑势垒的核心方程 #

对于氢同位素（比如氘-氘、氘-氚）的聚变，两个带正电的原子核需要克服库仑斥力靠近到核力能起作用的距离，此时对应的经典势垒峰值能量（即理论上需要克服的最小能量，不考虑量子隧穿）的标准方程是：

E_b = (k * Z₁ * Z₂ * e²) / r₀

其中各参数的含义：

• k：库仑常数，取值约为 8.988×10⁹ N·m²/C²
• Z₁, Z₂：两个聚变原子核的电荷数（氢同位素的Z都是1，比如氕、氘、氚）
• e：元电荷，取值约为 1.602×10⁻¹⁹ C
• r₀：两个原子核刚好能触发核力作用的临界距离，通常取核半径量级，计算公式为 r₀ ≈ 1.2×10⁻¹⁵ m × (A₁^(1/3) + A₂^(1/3))，其中A是原子核的质量数（比如氘核A=2，所以氘-氘聚变的r₀≈2.7×10⁻¹⁵ m）

代入数值计算的话，氘-氘聚变的经典势垒能量大概在0.5 MeV左右，但实际核聚变中因为量子隧穿效应，粒子不需要达到这么高的能量就能越过势垒——这也是为什么我们能在实验室（比如托卡马克）里实现可控聚变的关键原因。

2. 推导的核心逻辑 #

这个方程其实是从库仑势能的定义推导来的：
库仑势能是将两个点电荷从无穷远移动到距离r处所需做的功，用积分形式表示为：

E(r) = ∫(∞到r) F dr = ∫(∞到r) (kZ₁Z₂e²)/r'² dr' = (kZ₁Z₂e²)/r

当两个原子核靠近到核力能抵消库仑斥力的临界距离r₀时，这个势能值就是势垒的峰值，也就是经典情况下需要克服的最小能量。

3. 验证你自行推导方程的建议 #

如果你已经有了自己的推导结果，可以从这几个维度对比验证：

• 临界距离r₀的取值：别用原子半径！原子半径比核半径大4-5个数量级，用原子半径会导致结果完全错误，必须取核力作用的临界距离（核半径量级）
• 量纲匹配性：检查你的方程量纲是否为能量（ML²T⁻²），如果量纲不对，推导肯定有问题
• 电荷数与质量数的区分：别把质量数A当成电荷数Z来用，氢同位素的Z都是1，但A不同（氕A=1，氘A=2，氚A=3）
• 量子隧穿的修正：如果你考虑了实际聚变的隧穿效应，势垒的“有效克服能量”会比经典值低，这时候需要引入伽莫夫因子（Gamow factor）来修正，方程会更复杂一些

内容的提问来源于stack exchange，提问作者Tom Bushell