求助，模拟二氧化硅吸附甲烷气体报错，单独二氧化硅可以跑，整体跑不了

hyq12

求助
用GROMACS模拟二氧化硅吸附甲烷气体，跑完能量最小化之后，跑动力学模拟报错如下，periodic-molecules = yes开启了，单独二氧化硅试过了可以跑，但是整体没法跑，请问这是哪里出错了呢？
甲烷用的amber力场，二氧化硅使用sobtop例子里面的力场，水分子用tip3p力场
Fatal error:
There are inconsistent shifts over periodic boundaries in a molecule type
consisting of 288 atoms. The longest distance involved in such interactions is
1.869 nm which is above half the box length. Either you have excessively large
distances between atoms in bonded interactions or your system is exploding

牧生

原则上，sobtop自带的SiO2是指晶体内部的各个原子参数。而你这个模拟，应该将其中一个面用羟基饱和，然后再在这个被饱和的表面放置气体分子。被饱和的羟基表面，需要加入相应的O和H的类型以及参数。

http://bbs.keinsci.com/thread-31431-1-1.html

sobereva

这种事，加的分子数由少到多反复测试
一次加一堆自然不方便排查原因