相互作用能计算

本项目计算了一个分子在固体表面上的相互作用能曲线，用于确定最优吸附距离、吸附强度与脱附趋势。这张图展示的是典型的 Lennard-Jones 型相互作用：近距排斥、中距势阱、远距 van der Waals 尾。

这张吸附曲线，藏着三个信息

曲线最低点就是分子在表面上的”最佳站位”。坑深 0.38 eV 是吸附强度——低于 0.5 eV 一般算弱物理吸附，室温下能可逆脱附；坑底对应的距离 3.48 Å 比化学键长得多，说明主要靠范德华力粘着。曲线尾巴按 r⁻⁶ 衰减，正是色散作用的指纹。看一张图，吸附强不强、靠多近、靠什么力，全有了。

计算逻辑

相互作用能定义为：

E_int(r) = E_{molecule+surface}(r) − E_{molecule} − E_{surface}

本项目沿垂直于表面的方向，以 0.1 Å 为步长从 2.5 Å 扫描到 9 Å，计算每个距离下的体系总能量。扫描覆盖了排斥区、势阱区与长程尾区，保证曲线形状完整。

关键参数

• 泛函/方法：DFT-D3 或半经验色散矫正方法，必须描述长程 van der Waals 作用。
• 表面模型：足够厚的 slab，底层固定，表层允许弛豫。
• 分子取向：保持分子初始构象，仅改变质心到表面的垂直距离。
• 扫描步长：1 Å，势阱附近可加密到 0.05 Å。

vdW 方法怎么选

相互作用能曲线形状几乎完全由色散项决定，泛函选错全盘皆输：

本项目用 PBE-D3：在可接受成本内给出了 r₀ = 3.48 Å、ε = 0.38 eV 的合理势阱。若关心吸附构型随覆盖度的变化，应升级到 vdW-DF 系列。

结果判读

• ε = 0.38 eV：物理吸附强度，室温下部分可逆脱附。
• r_0 = 3.48 Å：吸附距离，可用于估算吸附层厚度。
• r⁻⁶ 尾：远距离行为符合色散相互作用，验证了方法对弱相互作用的描述能力。

可借鉴经验

• 不描述色散的泛函（如纯 PBE）会把吸附能算成接近 0 甚至排斥，必须加 D3/D4 或直接用 vdW-DF。
• 扫描时分子不能自由旋转，否则势阱位置会被平均化；应固定取向做单变量扫描。
• 吸附能 < 0.2 eV 通常认为室温下不稳定，> 1 eV 则倾向化学吸附。

吸附能结果自检

得到 ε = 0.38 eV、r₀ = 3.48 Å 后，本项目用四步确认它不是数值假象：

1. 色散开关对比：关掉 D3 重算，若吸附能跌到接近 0 或变正，说明确实是 vdW 主导，结果可信。
2. 扫描范围检查：曲线必须覆盖到 E_int 回升到 0 附近的长程区，否则势阱深度被截断低估。
3. slab 厚度测试：把 slab 从 4 层加到 6 层，吸附能变化应 < 0.03 eV，否则表面未充分弛豫。
4. 与实验对照：物理吸附能在1–0.5 eV 区间是常识范围，超出要警惕单位或模型错误。

局限

单点扫描只给出垂直方向的势能曲线，真实的吸附构型还需要在表面上做二维平移与旋转搜索。分子相互作用能相关方法可参考本站的分子对接计算栏目。更多服务见科研学术网首页。