计算化学公社

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

[GROMACS] 求助:在能量最小化时提示错误信息One or more water molecules can not be settled

[复制链接 Copy URL]

55

帖子

0

威望

771

eV
积分
826

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 芮啦 于 2021-4-7 08:56 编辑

各位老师,我定义了一个尺寸为4nm的立方体盒子,在盒子的中心放置适配体-铅离子,在CHARMM27力场下用Tip3p水解。进行到能量最小化这一步一直报错如下:
step 15: 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

step 21: 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

step 23: 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

......
step 492: 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

我只用金属离子进行分子动力学模拟,只用适配体进行分子动力学模拟,缩小盒子的尺寸等操作后依然出现同样的信息,提示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。我检查了我的pdb文件和坐标文件没有问题,麻烦各位老师能够指明一下解决的方向。

tp.pdb

41.25 KB, 下载次数 Times of downloads: 7

tp.gro

23.63 KB, 下载次数 Times of downloads: 10

topol.top

1.18 KB, 下载次数 Times of downloads: 18

t_ions.7z

399.41 KB, 下载次数 Times of downloads: 8

topol_DNA_chain_A.itp

148.75 KB, 下载次数 Times of downloads: 6

topol_Ion_chain_A2.itp

956 Bytes, 下载次数 Times of downloads: 7

em.mdp

718 Bytes, 下载次数 Times of downloads: 80

23

帖子

0

威望

225

eV
积分
248

Level 3 能力者

31#
发表于 Post on 2025-4-10 21:44:38 | 只看该作者 Only view this author
请问解决了吗

2

帖子

0

威望

27

eV
积分
29

Level 2 能力者

30#
发表于 Post on 2024-4-12 00:44:29 | 只看该作者 Only view this author
本帖最后由 juqing 于 2024-4-12 00:46 编辑
芮啦 发表于 2021-4-8 09:41
谢谢sob老师,我的初始结构和拓扑文件应该没有问题,在建模的过程中系统也没有警告和注释,这个问题困扰 ...

请问解决了吗,如果在linux下遇到一样的问题应该怎么做呢?

1560

帖子

0

威望

4995

eV
积分
6555

Level 6 (一方通行)

29#
发表于 Post on 2023-12-19 07:29:44 | 只看该作者 Only view this author
xiaohaha123 发表于 2023-12-18 22:26
请问能量最小化成功的标准是啥

一般情况下,如果在.top中,各个分子出现顺序正常,力场合适,正常结束em后,能跑起来MD,基本上就正确。
又菜又爱玩

2

帖子

0

威望

31

eV
积分
33

Level 2 能力者

28#
发表于 Post on 2023-12-18 22:26:49 | 只看该作者 Only view this author
牧生 发表于 2023-7-20 21:31
我现在的意见,只要em以后,跑md报错,大大大大概率是分子拓扑不正确

请问能量最小化成功的标准是啥

1560

帖子

0

威望

4995

eV
积分
6555

Level 6 (一方通行)

27#
发表于 Post on 2023-7-20 21:31:34 | 只看该作者 Only view this author
chuxuedexiaobai 发表于 2023-7-20 21:07
您好,我现在的模拟也是能量最小化一直降不下来,但没有报错,然后反复跑em也不行只要跑NVT还没开始就崩 ...

我现在的意见,只要em以后,跑md报错,大大大大概率是分子拓扑不正确
又菜又爱玩

55

帖子

0

威望

195

eV
积分
250

Level 3 能力者

26#
发表于 Post on 2023-7-20 21:07:05 | 只看该作者 Only view this author
牧生 发表于 2021-6-19 14:45
是的。如果反复多次还不能最小化,就考虑用双精度,
如果还不行,就考虑是自己做的模型是否合理。

您好,我现在的模拟也是能量最小化一直降不下来,但没有报错,然后反复跑em也不行只要跑NVT还没开始就崩溃了,怎么检查我的模型有没有问题呢,我是双相体系,两相单独分开能量最小化就可以,但合起来就不行了
,可以讲讲您解决的步骤吗,非常感谢您。

28

帖子

0

威望

457

eV
积分
485

Level 3 能力者

25#
发表于 Post on 2022-7-27 20:35:06 | 只看该作者 Only view this author
chengb917 发表于 2022-6-1 16:00
请问问题解决了吗?我现在碰到同样的问题。

解决了,就是结构和mdp的问题,需要认真检查。

15

帖子

0

威望

161

eV
积分
176

Level 3 能力者

24#
发表于 Post on 2022-6-1 16:00:57 | 只看该作者 Only view this author
木林森 发表于 2021-6-21 22:49
我现在就是在linux下,同样的问题一直报错

请问问题解决了吗?我现在碰到同样的问题。

15

帖子

0

威望

217

eV
积分
232

Level 3 能力者

23#
发表于 Post on 2021-11-16 12:26:25 | 只看该作者 Only view this author
木林森 发表于 2021-10-13 09:30
我的是初始结构出了小问题

您好,我也出现了同样的问题,但是报Fmax的原子都是加进去的水,挺奇怪的

28

帖子

0

威望

457

eV
积分
485

Level 3 能力者

22#
发表于 Post on 2021-10-13 09:30:51 | 只看该作者 Only view this author
蔡亚磊 发表于 2021-9-28 19:23
补充一下,我是跑蛋白的时候遇到的这个问题,前面最小化,nvt.npt都没问题。大概跑了两百多万步就出现上 ...

我的是初始结构出了小问题

4

帖子

0

威望

95

eV
积分
99

Level 2 能力者

21#
发表于 Post on 2021-9-28 19:23:06 | 只看该作者 Only view this author
木林森 发表于 2021-9-9 10:39
回复的晚了。
我之后就是检查结构,检查mdp文件解决的。
多数情况可能是结构的问题,不麻烦的话重建也 ...

补充一下,我是跑蛋白的时候遇到的这个问题,前面最小化,nvt.npt都没问题。大概跑了两百多万步就出现上述问题(大概是整体时间的一半),而且我换了一个蛋白也一样。

4

帖子

0

威望

95

eV
积分
99

Level 2 能力者

20#
发表于 Post on 2021-9-28 19:11:15 | 只看该作者 Only view this author
木林森 发表于 2021-9-9 10:39
回复的晚了。
我之后就是检查结构,检查mdp文件解决的。
多数情况可能是结构的问题,不麻烦的话重建也 ...

你好,我也遇到同样的问题。请问一下你是那一步的结构出问题了。

17

帖子

0

威望

71

eV
积分
88

Level 2 能力者

19#
发表于 Post on 2021-9-14 17:13:32 | 只看该作者 Only view this author
木林森 发表于 2021-9-9 10:39
回复的晚了。
我之后就是检查结构,检查mdp文件解决的。
多数情况可能是结构的问题,不麻烦的话重建也 ...

好的,感谢感谢

28

帖子

0

威望

457

eV
积分
485

Level 3 能力者

18#
发表于 Post on 2021-9-9 10:39:00 | 只看该作者 Only view this author
sdf 发表于 2021-8-11 20:45
您好,请问您的这个问题解决了吗?我是在linux下跑,也是一直出错,想请教一下您。Steepest Descents con ...

回复的晚了。
我之后就是检查结构,检查mdp文件解决的。
多数情况可能是结构的问题,不麻烦的话重建也可以。

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

GMT+8, 2026-2-20 01:10 , Processed in 0.213259 second(s), 30 queries , Gzip On.

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