楼主，您好，我在做MOF吸附氢气的模拟，遇到了困难，可以帮忙看一下吗

我想复现这篇文章哎（10.1021/acs.jpcc.8b09033），但是对分子模拟一窍不通，大佬可以给点自学的学习资料和学习路径吗？

水的等温线一直都是不好算的那种。在298K，对于水来说算是低温，最好要用足够的cycle，不然你的模拟就完全没有平衡。看文献那个等温线，这种z字型线，说明有可能有capillary condensation（一阶相变的一种，气液相变），这种情况下一定要把模拟跑得足够长。

HIll 发表于 2024-1-31 17:32
找到问题了，力场不对。gaff力场对COF 水蒸气吸附误差非常大。
还想请教大家另一个问题，raspa用coulomb g ...

Coulomb grid本质上就是算q_i*q_j/r_ij。当你的结构是不动的情况下，可以通过提前计算，来省去实时计算库仑力能量的时间。就是使用cubic spline来做interpolation。
但是这种interpolation不一定准，如果你的盒子是triclinic，并且tilt很大，可能就不准。还有种可能就是你提前算的时候选点选的太稀松，没有足够的点来代表这个结构。这个warning就是RASPA在模拟刚开始的时候，做了个interpolation准确度的测试，然后发现不够准。
David Dubbeldam建议使用0.15 Å一下的间隔，不然就不一定够准。

HIll 发表于 2024-2-19 10:51
效果会好很多，但是在非常低压的情况下会有一定的吸附量，之后的曲线和文献基本一致。
感觉还是电荷上存 ...

好的 谢谢楼主
好像DREIDING力场会好很多
电荷的话，Qeq DDEC我用的比较多
以后多讨论~

yopoxia 发表于 2024-2-4 16:28
楼主好，按照新的力场，RASPA模拟的准了么。我试了几次raspa的模拟，用的uff，模拟了一些气体吸附，感觉很 ...

效果会好很多，但是在非常低压的情况下会有一定的吸附量，之后的曲线和文献基本一致。
感觉还是电荷上存在一些误差。
我之后准备做一些气体吸附的测试，欢迎讨论，之后再回复你。

楼主好，按照新的力场，RASPA模拟的准了么。我试了几次raspa的模拟，用的uff，模拟了一些气体吸附，感觉很不准啊

找到问题了，力场不对。gaff力场对COF 水蒸气吸附误差非常大。
还想请教大家另一个问题，raspa用coulomb grid的时候出现Coulomb energy interpolation table is wrong该怎么办呢？