lammps计算结合能

项目过程：MOF-5 吸附 CO₂、CH₄、N₂ 的结合能与吸附热

客户研究 MOF-5 在气体分离中的应用，需要定量比较三种小分子在 MOF-5 孔道中的结合强度，并给出等量吸附热 Q_st，用来和实验吸附等温线做对比。

我用 LAMMPS 做分子动力学采样，结合 GCMC 计算吸附负载。力场方面，MOF 骨架用 UFF/DREIDING 混合参数，气体分子用 TraPPE。LJ 截断 8.5 Å，长程静电用 Ewald。每个体系跑 5 ns NVT MD，温度 298 K，然后用 GCMC 统计不同 fugacity 下的吸附量，再通过 Clausius-Clapeyron 关系计算 Q_st。

结合能分解为 van der Waals 和 Coulomb/LJ 两部分，用来判断哪种相互作用主导。

案例分析：CO₂ 结合能最强，静电贡献不可忽视

第一张案例图是三种气体的结合能分解柱状图。CO₂ 的总结合能约 −28 kJ·mol⁻¹，其中 vdW 贡献约 −20 kJ·mol⁻¹，Coulomb/LJ 贡献约 −8 kJ·mol⁻¹；CH₄ 约 −19 kJ·mol⁻¹，主要来自 vdW；N₂ 约 −12 kJ·mol⁻¹，也是 vdW 主导。CO₂ 的强结合来自它的四极矩与骨架金属位点的静电相互作用。

第二张案例图是等量吸附热 Q_st 随吸附负载的变化。CO₂ 在低负载时 Q_st 约 26.5 kJ·mol⁻¹，随负载增加逐渐降到 21.5 kJ·mol⁻¹；CH₄ 从约 19 降到 15.5 kJ·mol⁻¹。图中还画了实验点，模拟曲线和实验趋势吻合较好。

总结教训：力场参数决定气体分离预测的可靠性

第一，MOF 的静电势对 CO₂ 吸附影响很大。如果去掉骨架电荷，CO₂ 结合能会低估约 8 kJ·mol⁻¹，Q_st 也会明显偏低。

第二，GCMC 的加载量对化学势敏感。我对比了 298 K 下几个 fugacity 点，确保 Q_st 是在 Henry 区到中等负载区计算的，避免高压下吸附饱和带来的非线性影响。

第三，LJ 截断 8.5 Å 对 MOF 体系够用。我测试了 10 Å，结合能差异小于 1 kJ·mol⁻¹，但计算量增加约 30%。

项目经验：结合能分解比总能更直观

对客户来说，−28 kJ·mol⁻¹ 不如“CO₂ 主要靠静电、CH₄ 和 N₂ 靠色散”这句话好理解。做气体分离项目时，我会把结合能分解和选择性曲线一起给，客户可以直接判断该优化 MOF 的哪个化学位点来提升选择性。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：lammps计算结合能。