计算化学公社

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

[GROMACS] 执行mdrun后第200000多步就停止了,说one or more water can not be settled

[复制链接 Copy URL]

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

各位大神,在执行gmx mdrun -nt 8 -v -deffnm md_0_10 ,第200000多步就停止了 说one or more water can not be settled ,请问有什么解决方法吗?

43

帖子

0

威望

197

eV
积分
240

Level 3 能力者

2#
发表于 Post on 2019-8-30 10:34:48 | 只看该作者 Only view this author
      one or more water can not be settled 在gromacs模拟的时候是很常见的报错,我个人虽然也才接触没多久,但是之前已经遇到过好多次了。
      这个问题我遇到的时候主要原因还是体系构型可能有问题,比如周期性边界条件的时候,边界处可能是不是存在不合理的地方要调整一下。(我承认这点我说的也有点模糊)这个问题往往和体系能量过高一起出现,至少我是这样的。
      虽然我主要是模拟开始的时候就会出现这条报错而你是运行了一段时间后才报错,说明可能是模拟过程中,构型可能出现异常吧。
      总之,其实这个问题,和你的命令关系不大吧看,主要还是要把系统拿出来看看,不然不好确定问题,以及可能的问题。
      不过不晓得别人是不是可以做出更精准判断。

43

帖子

0

威望

197

eV
积分
240

Level 3 能力者

3#
发表于 Post on 2019-8-30 10:40:22 | 只看该作者 Only view this author
      对了,你那边有吗有装高版本点的gromacs?我理解有些课题组可能对低版本的gromacs有一些源代码的修改,而高版本的gromacs因为变成c++语言了,所以所做的修改就不适宜了,所以可能还是停留在低版本。但是其实高版本对于模拟问题的报错其实精确很多的,可以看看相同的体系相同的跑法放到2018或2019版的gromacs,很可能报错就不一样了,而且可能就知道问题所在了。

43

帖子

0

威望

197

eV
积分
240

Level 3 能力者

4#
发表于 Post on 2019-8-30 10:44:58 | 只看该作者 Only view this author
      中间报错终止的模拟,好像都会产生几个pbd文件的,里面是是报错前不久的系统坐标,可以放到VMD里看看,说不定体系已经变形了都不一定

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2019-8-30 13:39:15 | 只看该作者 Only view this author
cavalier 发表于 2019-8-30 10:44
中间报错终止的模拟,好像都会产生几个pbd文件的,里面是是报错前不久的系统坐标,可以放到VMD里看看 ...

真的很感谢 我是刚开始做这个 现在就是照网上进行练习 我打算重做一遍 看看是不是之前哪个步骤出错了

43

帖子

0

威望

197

eV
积分
240

Level 3 能力者

6#
发表于 Post on 2019-8-30 15:03:23 | 只看该作者 Only view this author
chensy 发表于 2019-8-30 13:39
真的很感谢 我是刚开始做这个 现在就是照网上进行练习 我打算重做一遍 看看是不是之前哪个步骤出错了

      其实我也才做了1年,不过刚开始的时候,最常遇见的就是水分子不能settled,总之现在感觉这个报错很可能原因不是他字面上这个意思,既然你只是刚开始,那方便的话最好是能看看你的系统,这样大家应该好找出问题所在,因为我担心你可能自己复查过程都很难查出问题,一般都是很小的细节出错了。

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

7#
发表于 Post on 2019-8-30 21:56:37 | 只看该作者 Only view this author


北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2019-8-31 08:35:21 | 只看该作者 Only view this author

C:\Users\HuiYue\Desktop\图片1.png
谢谢sob老师,我看看坐标有没有什么问题 再看看减少步长可不可以

图片1.png (13.65 KB, 下载次数 Times of downloads: 61)

错误信息

错误信息

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2019-8-31 08:36:31 | 只看该作者 Only view this author
chensy 发表于 2019-8-31 08:35
谢谢sob老师,我看看坐标有没有什么问题 再看看减少步长可不可以

sob老师,因为我又做了一遍还是出现同样的情况 我再重新检查一下坐标文件 看看行不行

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

10#
 楼主 Author| 发表于 Post on 2019-8-31 08:39:41 | 只看该作者 Only view this author
cavalier 发表于 2019-8-30 15:03
其实我也才做了1年,不过刚开始的时候,最常遇见的就是水分子不能settled,总之现在感觉这个报错很 ...

好的 我又重新做了一遍还是出现这样的结果,我再好好看看 谢谢了

图片1.png (13.65 KB, 下载次数 Times of downloads: 31)

图片1.png

43

帖子

0

威望

197

eV
积分
240

Level 3 能力者

11#
发表于 Post on 2019-8-31 10:47:42 | 只看该作者 Only view this author
他这里有说把当前和之前的坐标写成了一个pdb文件,你可以到VMD里看一下,是不是正常,还是说体系扭曲了或怎么的。不行的话最好把pdb的图给大家看下。

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

12#
发表于 Post on 2019-8-31 20:18:52 | 只看该作者 Only view this author
chensy 发表于 2019-8-31 08:35
谢谢sob老师,我看看坐标有没有什么问题 再看看减少步长可不可以

不要光把注意力集中在初始结构上。如果初始结构有严重问题,通常在MD的早期就已经会暴露出来
仔细看整条模拟的轨迹,如果从一开始直到崩溃为止,轨迹上看不出大毛病,那就更不是初始结构的事。从mdp参数和拓扑文件的合理性上考虑
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

33

帖子

0

威望

121

eV
积分
154

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2019-9-3 18:43:39 | 只看该作者 Only view this author
sobereva 发表于 2019-8-31 20:18
不要光把注意力集中在初始结构上。如果初始结构有严重问题,通常在MD的早期就已经会暴露出来
仔细看整条 ...

谢谢 sob老师,我减少了步长,发现是可以的

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 06:46 , Processed in 0.252572 second(s), 23 queries , Gzip On.

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