DFT 结合能计算

本项目以双原子分子为模型体系，用 DFT 计算沿键长扫描的势能曲线，拟合 Morse 势得到平衡键长、解离能与振动频率。这张图展示的是一个典型结果：结合能随原子间距变化，最低点对应稳定分子态。下文把建模逻辑与拟合过程整理出来，供类似体系参考。

这张图到底在讲什么

读这张图其实不用复杂公式，抓住一个点就够了：曲线最低的那个坑。坑有多深，就是两个原子抱在一起有多牢——这里坑深 4.52 eV，对应平衡核间距 0.74 Å。坑右边往上爬直到归零，意思是把两个原子拉到无穷远需要的能量，也就是结合能。一眼看过去，这就是典型的双原子分子势能曲线。

计算逻辑

结合能定义为：

E_b(r) = E_{AB}(r) − E_A − E_B

本项目用 DFT 计算三个量：孤立原子 A 的能量、孤立原子 B 的能量，以及在不同核间距 r 下分子 AB 的总能量。扫描 r 从 0.5 Å 到 3.5 Å，得到一条势能曲线。曲线最低点越负，说明分子越稳定。

关键参数

• 泛函：PBE 起步，对弱束缚体系加 D3 色散修正。
• 基组/截断能：小分子用高斯基组 def2-TZVP 或平面波截断 ≥ 400 eV。
• 自旋多重度：孤立原子与分子的自旋态需一致，否则能量差无意义。
• 扫描步长：05–0.1 Å，平衡点附近加密。

单原子能量怎么取（决定结合能准不准）

结合能公式里的 E_A、E_B 不是随便一个单点能。本项目实际操作时遵守三条规则，这也是新手最容易翻车的地方：

• 放在大真空盒里算：孤立原子放进边长 ≥ 12 Å 的立方盒子，避免周期性镜像相互作用。盒子太小，原子会和自己的镜像”成键”，能量被人为拉低，结合能随之虚高。
• 自旋态对齐：O 原子基态是三重态（S=1），H 是双重态（S=1/2）。算 E_A−E_B 时两边自旋基底必须一致，否则能量差没有物理意义。
• BSSE 修正：平面波或大基组下，基组叠加误差会系统性高估结合能。本项目对弱相互作用体系额外跑了 Counterpoise 修正，修正量约 02–0.05 eV，弱束缚体系不能省。

Morse 势拟合

扫描得到的离散点用 Morse 函数拟合：

V(r) = D_e [1 − e^{−a(r − r_e)}]² − D_e

其中 D_e 是阱深（结合能），r_e 是平衡键长，a 控制势阱宽度。本项目通过非线性最小二乘拟合，得到图中参数：r_e = 0.74 Å，D_e = 4.52 eV。由 Morse 参数可进一步估算振动频率 ω_e，与实验光谱对照。

这个案例能借鉴什么

• 结合能计算不是单点能，必须做完整扫描才能判断是否有势阱。
• 扫描范围要足够宽，否则右侧解离极限定不准，结合能会被系统性高估。
• 自旋态错误是最常见的坑：孤立原子如果是开壳层，分子优化后必须重新检查自旋。

这套结果能直接用在哪里

结合能曲线不是孤立成果，本项目把它喂给后续三件事：

• 反应能估算：已知 A+B→AB 的结合能，可直接拼出反应放热/吸热，省掉重复大体系计算。
• 振动频率：由 Morse 参数 a 估算 ω_e ≈ a·√(2·D_e·μ)/(πc)，与实验红外峰对照，验证泛函合理性。
• 势能面粗模：把 Morse 势作为反应坐标初猜，代入 NEB 或过渡态搜索。

输入参数速查（VASP 示例）

扫描时只改 POSCAR 原子坐标，其余参数全程固定，保证曲线各点可比。

局限

Morse 势只描述单键振动，对多原子分子或多重键的拟合能力有限；对强关联体系（如 O₂、过渡金属二聚体）需用多参考方法或 HSE06 修正。更多 DFT 案例可参考本站的第一性原理计算栏目。科研计算服务总览见科研学术网首页。