计算化学公社

标题: 求助：MD平衡相可顺利跑完，但产生相跑一段时间报错 [打印本页]

作者
Author: nku_linz 时间: 2023-11-5 00:27
标题: 求助：MD平衡相可顺利跑完，但产生相跑一段时间报错
跑一个含硼原子的聚合物片段（200原子左右）在水中的MD，能量最小化和升温平衡相可顺利进行，跑产生相跑了5 ns后报错，程序提示如图

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

然后根据社长在sobtop程序主页中的FAQ8排查问题

首先考虑了分子结构和拓扑文件结构是否合理。此处分子结构是先用GFN0-xTB初步优化，再用ORCA的B97-3c基组优化得到。电荷使用RESP电荷。拓扑文件是根据sobtop程序主页例2生成，其中缺失的B原子参数和相关成键参数从文献(Biomacromolecules 2023, 24, 7, 3380–3396)中获取。

把分子放到真空中模拟，使用的.mdp文件为em_nopbc.mdp和md_nopbc.mdp，真空中模拟可顺利跑200ps，观察轨迹未看到异常，初步认为结构和拓扑文件应该合理

接着考虑水相模拟时.mdp文件不合理，升温平衡相可顺利进行，因此尝试增加产生相的控压常数，减小模拟步长，去掉H键约束，仍然报错，观看多条水相中的模拟轨迹，分子均维持正常结构

模拟不是一开始就崩，往往维持一段时间才出问题，自己看不出问题在哪，麻烦各位老师指点

用的程序为gmx2019.6_AVX_CUDA，同时附上B原子的参数

(, 下载次数 Times of downloads: 11) 出处：Biomacromolecules 2023, 24, 7, 3380–3396

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