使用RASPA和lammps计算MOF对气体吸附量计算结果错误

偶有

目前我在研究一种柔性MOF对苯和环己烷的吸附，目前我已经计算了MOF的DDEC6，用于MOF的力场中。MOF使用DREIDING（C、H、O、N）和UFF（Zr)力场，苯和环己烷使用TraPPE力场，MOF、苯、环己烷都是使用刚性分子去模拟的。我的MOF材料在吸附苯的过程中会逐渐从闭口相变为开口相。我用Zeo+分析得到的三个相的几何苯吸附量上限如下Phase1 0.49 mmol/g，Phase2 4.43 mmol/g，Phase3 10.14 mmol/g，但我现在GCMC计算的phase1在123pa时苯的吸附量为6.56mmol/g，在12154.916pa下phase3苯吸附量为18.42mmol/g，苯的饱和蒸汽压为12700pa，出现模拟的吸附量远大于实验值和理论值是什么原因（对比数值如下图），有大佬遇到过这种问题吗
之前使用lammps用MOF-FF力场在GCMC/MD模拟时phase3吸附量也在15.6mmol/g左右，但是化学势对应的压力远大于饱和蒸汽压，最后没有使用lammps去模拟

yuzc

1. 力场不够精准，尤其是带多价金属的。
2. 有相变，但是你没用osmotic ensemble。
3. GCMC的原理，没法让分子自由扩散进孔道，所以如果PLD非常小，GCD非常大，程序会按照GCD尺寸为上限填满吸附相。导致高估。前人对这种方法有一个妙招，就是在不该进入气体的孔道里面塞一个大刚球。