计算化学公社

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

[GROMACS] 求助:混合膜的模拟中,能量最小化时提示错误信息Energy minimization has stopped.

[复制链接 Copy URL]

2

帖子

0

威望

91

eV
积分
93

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师,我初学混合膜的模拟,组成混合膜的两种磷脂分子是在GROMOS54A8力场下的, 步骤是跟着培训班教程一步一步走的,但是在能量极小化步骤中出现了以下报错信息:
Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 100 (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.

Polak-Ribiere Conjugate Gradients converged to machine precision in 100 steps,
but did not reach the requested Fmax < 100.
Potential Energy  =  4.2020992e+08
Maximum force     =  4.0882774e+08 on atom 3289
Norm of force     =  3.5014327e+06

不知道该怎么解决,请老师们指教,谢谢。
相关文件附上
memwat.rar (171.65 KB, 下载次数 Times of downloads: 7) mixmem.rar (36.83 KB, 下载次数 Times of downloads: 5)


DPO.itp (10.93 KB, 下载次数 Times of downloads: 3)

POP.itp (10.43 KB, 下载次数 Times of downloads: 3) topol.top (1.42 KB, 下载次数 Times of downloads: 12)

em.mdp (389 Bytes, 下载次数 Times of downloads: 5) em.log (19.91 KB, 下载次数 Times of downloads: 1)


75

帖子

0

威望

675

eV
积分
750

Level 4 (黑子)

2#
发表于 Post on 2023-5-18 00:13:19 | 只看该作者 Only view this author
这么大的力一定是拓扑出问题,建议检查结构,能量最小化如果有别的报错应该会输出错误帧的pdb文件。看一下是从哪个原子开始崩掉的。.top别把分子一个个列出来啊,直接写name写数量就好了。

6万

帖子

99

威望

6万

eV
积分
125177

管理员

公社社长

3#
发表于 Post on 2023-5-18 00:15:51 | 只看该作者 Only view this author
pdb文件里DLP残基是什么东西?都没有定义,和[molecules]里明显不符
必须严格确保pdb文件里原子顺序和[molecules]部分的对应关系
北京科音自然科学研究中心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

2

帖子

0

威望

91

eV
积分
93

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2023-5-19 16:16:56 | 只看该作者 Only view this author
谢谢老师  我修改后能量极小化显示粒子间的距离有问题,再试了几次报错信息又是
Energy minimization has stopped, but the forces have not converged to the
requested precision Fmax < 100 (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.
请问这会不会是我的初始结构有问题啊  我初学混合膜模拟  每一步都是按照教程来的  可能细节不太懂  请老师指教  谢谢

DPO.itp

10.89 KB, 下载次数 Times of downloads: 0

em.log

19.9 KB, 下载次数 Times of downloads: 0

em.mdp

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

POP.itp

10.39 KB, 下载次数 Times of downloads: 0

topol.top

1.42 KB, 下载次数 Times of downloads: 2

memwat.rar

117.58 KB, 下载次数 Times of downloads: 1

本版积分规则 Credits rule

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

GMT+8, 2026-2-25 12:59 , Processed in 0.192523 second(s), 23 queries , Gzip On.

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