计算化学公社

标题: mdrun做矿物-溶液-离子体系时遇到的一些问题 [打印本页]

作者
Author: 木林森 时间: 2021-7-19 11:40
标题: mdrun做矿物-溶液-离子体系时遇到的一些问题
最近在做矿物-溶液-离子体系，分别做了能量最小化，npt，nvt，在npt和nvt时各遇到了两个不同问题，以下是我的文件（gro文件略大，没上传）：
(, 下载次数 Times of downloads: 2)

(, 下载次数 Times of downloads: 2)

(, 下载次数 Times of downloads: 6)

(, 下载次数 Times of downloads: 3)

(, 下载次数 Times of downloads: 4)

(, 下载次数 Times of downloads: 1)

(, 下载次数 Times of downloads: 3)

问题一：
在进行npt时，用的是能量最小化后生成的gro文件，提交运行后可以跑起来，但过一段时间就会自动停下并报错：Fatal error:
1 of the 73421 bonded interactions could not be calculated because some atoms involved moved further apart than the multi-body cut-off distance (1.66356 nm) or the two-body cut-off distance (1.66356 nm), see option -rdd, for pairs and tabulated bonds also see option -ddcheck
但如果使用须跑命令，它又可以跑起来，不过过一段时间又会出现同样的报错，当我更换了体系中的离子类型后（比如将Na全部替换为Pb），又可以正常运行，不会报错

问题二：
(, 下载次数 Times of downloads: 8) (, 下载次数 Times of downloads: 26)
我在构建体系时，在矿物局部地方同时将两个-OH连接到一个Al原子上（如图中abc和ade），abc中的羟基是通过将水分子去掉一个氢原子实现的，并在MMT010.itp中增加了相应的键长和键角，但在进行能量最小化时abc中的-OH会从Al上脱离进入真空层（如第二张图），而ade另一个羟基就正常连接。

问题一二出现在不同体系中的，请问大家这个两个问题该怎么样解决？

作者
Author: sobereva 时间: 2021-7-19 22:55
2 可能是bond项没起作用，检查拓扑文件写错没有。谐振势如果生效的话绝对不可能跑那么远

作者
Author: 木林森 时间: 2021-7-20 16:13

sobereva 发表于 2021-7-19 22:55
2 可能是bond项没起作用，检查拓扑文件写错没有。谐振势如果生效的话绝对不可能跑那么远

谢谢sob老师，我仔细检查过itp文件，应该没啥问题，也不知道该怎么解决了

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