计算化学公社

标题: 求助,模拟二氧化硅吸附甲烷气体报错,单独二氧化硅可以跑,整体跑不了 [打印本页]

作者
Author:
hyq12    时间: 2024-12-8 17:54
标题: 求助,模拟二氧化硅吸附甲烷气体报错,单独二氧化硅可以跑,整体跑不了
求助
用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



作者
Author:
牧生    时间: 2024-12-8 18:44
本帖最后由 牧生 于 2024-12-8 18:45 编辑

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

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

作者
Author:
sobereva    时间: 2024-12-9 03:46
这种事,加的分子数由少到多反复测试
一次加一堆自然不方便排查原因





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3