计算化学公社

标题: 使用RASPA和lammps计算MOF对气体吸附量计算结果错误 [打印本页]

作者
Author: 偶有 时间: 3 day ago
标题: 使用RASPA和lammps计算MOF对气体吸附量计算结果错误
本帖最后由偶有于 2026-4-10 11:16 编辑

目前我在研究一种柔性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去模拟

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)