求助用LAMMPS计算MOFs对水分子的吸附能的问题

筱婷397

我想快速计算出一批MOFs对水分子的吸附数据，为了便于后续比较，我用LAMMPS计算了MOFs吸附一个水分子后体系的pe（UFF4MOF力场和TIP4P_2005_long水模型）、纯MOF的pe（UFF4MOF力场），基于吸附能计算公式，现在我想得到单个水分子的能量。我多次尝试在一个很大的周期性盒子里放置一个水分子，计算出的pe都是约等于0。
我现在的问题是：
1. 我想计算的吸附能，是用pe相减吗？
2. 是因为经典水模型在分子内用fix处理键/角，所以水分子的pe≈0吗？那我可以直接用 Eads = E_pe(MOF+H2O) − E_pe(MOF) 计算吸附能数据吗？
3. 还有什么更好的方法，可以快速得到大量MOFs的吸附数据吗？

以上问题请教各位老师，谢谢！