计算化学公社

 找回密码 Forget password
 注册 Register
Views: 3076|回复 Reply: 1
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 请问在平衡动力学中,为什么会出现如下错误?

[复制链接 Copy URL]

2

帖子

0

威望

27

eV
积分
29

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我的体系是甲苯-水-表面活性剂体系,表面活性剂溶解在水中,油水都是5*5*8nm,packmol建模。然后通过gromacs的命令行,将z方向添加真空层,扩展成30~40nm。xyz周期性体系。
GAFF+Hessian,用sobtop产生的拓扑文件。
能量最小化之后,做平衡动力学mdrun时候,出现如下报错:
Fatal error:
There are inconsistent shifts over periodic boundaries in a molecule type consisting of 15 atoms. The longest distance involved in such interactions is 9.570 nm which is above half the box length. Either you have excessively large distances between atoms in bonded interactions or your system is exploding.


15个原子明显是甲苯,我无论在建模,还是top文件,都是甲苯放在了第一位。

我唯一找到的报错相关资料,是http://bbs.keinsci.com/thread-22373-1-1.html
sob老师说是gro和itp原子顺序出了问题,但是我反复核查,确定是没问题的。
下面有说真空层不够高,我尝试拓展了真空层,没有用。
附件是我的全部源文件。
万分万分感激!实在不知道怎么做了。

span60.itp

69.75 KB, 下载次数 Times of downloads: 6

表面活性剂拓扑

system.top

1.18 KB, 下载次数 Times of downloads: 7

体系拓扑

toluene.itp

10.71 KB, 下载次数 Times of downloads: 4

甲苯拓扑

em.mdp

384 Bytes, 下载次数 Times of downloads: 7

能量最小化mdp

eq.mdp

669 Bytes, 下载次数 Times of downloads: 5

平衡动力学mdp

2

帖子

0

威望

27

eV
积分
29

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2022-4-29 21:10:14 | 只看该作者 Only view this author
解决了,盒子的横向设小了,导致xy方向上周期出了错

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 07:18 , Processed in 0.204209 second(s), 29 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list