计算化学公社

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

[GROMACS] 能量最小化时被迫停止,不能完成能量最小化,应该如何解决?

[复制链接 Copy URL]

3

帖子

0

威望

81

eV
积分
84

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 Rororozoroni 于 2023-12-18 10:10 编辑

我在使用官方所给minim.mdp文件对我建立的模型进行能量最小化(设置要跑5000step),但是跑了50多个step就报了以下的提示:
Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 1000 (which may not be possible for your system).
It stopped because the algorithm tried to make a new step whose size was too
small, or there was no change in the energy since last step. Either way, we
regard the minimization as converged to within the available machine
precision, given your starting configuration and EM parameters.

Double precision normally gives you higher accuracy, but this is often not
needed for preparing to run molecular dynamics.
You might need to increase your constraint accuracy, or turn
off constraints altogether (set constraints = none in mdp file)

writing lowest energy coordinates.
我注意到可以在mdp文件设置constraints=none,然后我将这一句写入minim.mdp文件,运行后仍然报了以上的提示,后续进行nvt就不跑了,跑到step0就说One or more water molecules can not be settled.Check for bad contacts and/or reduce the timestep if appropriate.Wrote pdb files with previous and current coordinates,我用vmd观察我的模型没有发现原子堆叠,相互之间是分开的没有重合或交叉,目前来看可能就是能量最小化没有运行成功,没有让体系能量最小化,所以跑不了nvt,请问遇到这种情况应该如何解决?

minim.mdp

992 Bytes, 下载次数 Times of downloads: 26

用于能量最小化的配置文件

system_newbox.gro

288.68 KB, 下载次数 Times of downloads: 3

我建立的模型文件

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

2#
发表于 Post on 2023-12-18 16:37:39 | 只看该作者 Only view this author
你要模拟真空下磷脂双层膜吸附小分子?结构文件里没有水
一般都是在水环境中模拟的

mdp文件没明显问题

此类报错原因该说的下文都说了,没有什么可补充的
http://sobereva.com/soft/Sobtop#FAQ8
北京科音自然科学研究中心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

3

帖子

0

威望

81

eV
积分
84

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-12-18 17:12:35 | 只看该作者 Only view this author
sobereva 发表于 2023-12-18 16:37
你要模拟真空下磷脂双层膜吸附小分子?结构文件里没有水
一般都是在水环境中模拟的

老师您好,我最后的文件中是有水的,但是那个.gro文件太大了放不上来,就放了一个未溶解水的.gro,这个大小刚好可以放上来,我会认真按照您给的链接进行排查纠错,感谢您的指点!

5

帖子

0

威望

23

eV
积分
28

Level 2 能力者

4#
发表于 Post on 2024-4-6 10:30:10 | 只看该作者 Only view this author
Rororozoroni 发表于 2023-12-18 17:12
老师您好,我最后的文件中是有水的,但是那个.gro文件太大了放不上来,就放了一个未溶解水的.gro,这个大 ...

你好,这个问题解决了,可以请教一下吗

5

帖子

0

威望

23

eV
积分
28

Level 2 能力者

5#
发表于 Post on 2024-4-6 10:35:07 | 只看该作者 Only view this author
你好,这个问题结局了吗,可以请教一下吗

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 17:11 , Processed in 0.168894 second(s), 23 queries , Gzip On.

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