计算化学公社

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

[GROMACS] 求助:分子动力学在进行能量极小化步骤时,出现报错:Listed nonbonded...

[复制链接 Copy URL]

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

跳转到指定楼层 Go to specific reply
楼主
各位老师好,我在分子动力学中,进行能量极小化步骤时,出现报错,报错如下:

Listed nonbonded interaction between particles 14786 and 14895 at distance 3.338 which is larger than the table limit 2.200 nm.

==============================================================================

This is likely either a 1,4 interaction, or a listed interaction inside a smaller molecule you are decoupling during a free energy calculation. Since interactions at distances beyond the table cannot be computed, they are skipped until they are inside the table limit again. You will onlyseethis message once, even if it occurs for several interactions IMPORTANT:This should not happen in a stable simulation, so there is D  probablysomething wrong with your system. Only change the table-extension distance in the mdp file if you are really sure that is the reason.
===============================================================================

我有在公社中看到和我出现一样问题的帖子,知道是我的拓扑文件出现了问题,但是我反复检查了几遍,依然没法发现问题所在。所以想来求助一下各位老师。
em.mdp文件我用的是师兄用过的原文件,拓扑文件也在附件里面。
小分子的mol2文件用guassianview加氢后保存,保存后用acpype在线网站(https://www.bio2byte.be/acpype/status/xqxlwVdknA)生成的小分子的itp和gro文件,在acpype网站中,Charge Method选择的是user,Multiplicity是1,Net Charge是0,Atom Type为AMBER。
em.mdp文件我用的是师兄用过的原文件,拓扑文件也在附件里面。

希望各位老师能够帮我解决这个报错的问题。



微信图片_20230313230901.jpg (919 KB, 下载次数 Times of downloads: 9)

微信图片_20230313230901.jpg

MOL.mol2

10.41 KB, 下载次数 Times of downloads: 5

小分子mol2

em.mdp

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

topol.top

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

拓扑文件

MOL_GMX.itp

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

小分子itp文件

MOL_GMX.gro

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

小分子gro文件

1562

帖子

0

威望

5007

eV
积分
6569

Level 6 (一方通行)

2#
发表于 Post on 2023-3-14 09:49:30 | 只看该作者 Only view this author
本帖最后由 牧生 于 2023-3-14 10:53 编辑

大概率是结构不合适,这两个氢原子靠有点近,优化结构的过程就有可能造成它们脱掉,你尝试一下把两个氢原子的朝向不同的方向。别的地方可能还有别的不合理的地方。你可以自己先优化一下结构,看看有没有不合理的地方,然后再用sobtop去搞top文件。


同时,VMD打开gro文件时,有提示ERROR) Suspicious pbc side length (a=0 b=0 c=0). Have you forgotten to set the pbc parameters?
ERROR) Suspicious pbc angle (alpha=0 beta=0 gamma=0).





又菜又爱玩

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2023-3-14 20:15:47 | 只看该作者 Only view this author
牧生 发表于 2023-3-14 09:49
大概率是结构不合适,这两个氢原子靠有点近,优化结构的过程就有可能造成它们脱掉,你尝试一下把两个氢原子 ...

谢谢老师!!已经在修改了

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2023-3-18 11:09:55 | 只看该作者 Only view this author
牧生 发表于 2023-3-14 09:49
大概率是结构不合适,这两个氢原子靠有点近,优化结构的过程就有可能造成它们脱掉,你尝试一下把两个氢原子 ...

老师好,我更改了小分子的结构以后,还是会出现我提问的这个结果,然后我单独跑了蛋白,发现蛋白质也是显示这样的结果,我就推测是蛋白质的初始结构出了问题,但是我不知道蛋白质的初始结构如何更改,想请教一下老师

1562

帖子

0

威望

5007

eV
积分
6569

Level 6 (一方通行)

5#
发表于 Post on 2023-3-18 11:41:32 | 只看该作者 Only view this author
本帖最后由 牧生 于 2023-3-18 11:43 编辑

试了下,我把这个小分子放进了一个大的空盒子里面,虽然没有算电荷,但是这个小分子在真空确实是可以跑起来的。

box.gro

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

MOL.itp

209.69 KB, 下载次数 Times of downloads: 4

MOL.pdb

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

MOL.top

333 Bytes, 下载次数 Times of downloads: 4

又菜又爱玩

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2023-3-18 11:55:31 | 只看该作者 Only view this author
牧生 发表于 2023-3-18 11:41
试了下,我把这个小分子放进了一个大的空盒子里面,虽然没有算电荷,但是这个小分子在真空确实是可以跑起来 ...

那是不是说明是蛋白质的问题,我该如何修改这个问题呢?我看老师之前有用到双精度,我去尝试了,发现也报错...

1562

帖子

0

威望

5007

eV
积分
6569

Level 6 (一方通行)

7#
发表于 Post on 2023-3-18 11:56:24 | 只看该作者 Only view this author
wyxlfw 发表于 2023-3-18 11:55
那是不是说明是蛋白质的问题,我该如何修改这个问题呢?我看老师之前有用到双精度,我去尝试了,发现也报 ...

我用的单精度的,跑起来没有明显问题。。没有必要双精度。你试试我传上来的文件
又菜又爱玩

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

8#
 楼主 Author| 发表于 Post on 2023-3-18 12:36:09 | 只看该作者 Only view this author
牧生 发表于 2023-3-18 11:56
我用的单精度的,跑起来没有明显问题。。没有必要双精度。你试试我传上来的文件

好的!!谢谢老师!!

24

帖子

0

威望

91

eV
积分
115

Level 2 能力者

9#
 楼主 Author| 发表于 Post on 2023-3-18 16:19:16 | 只看该作者 Only view this author
wyxlfw 发表于 2023-3-18 12:36
好的!!谢谢老师!!

老师好,我试了一下老师的文件,发现还是出现这个报错,我怀疑是蛋白质的初始结构有问题,我该如何查看这两个原子之间的关系呢?

43

帖子

0

威望

207

eV
积分
250

Level 3 能力者

10#
发表于 Post on 2023-11-30 11:18:49 | 只看该作者 Only view this author
牧生 发表于 2023-3-14 09:49
大概率是结构不合适,这两个氢原子靠有点近,优化结构的过程就有可能造成它们脱掉,你尝试一下把两个氢原子 ...

你好,我想请教下为什么这俩H离得太近会在优化过程中造成他们脱掉呢,分子力学不是不会涉及键的断裂和生成吗老哥?

本版积分规则 Credits rule

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

GMT+8, 2026-2-26 07:30 , Processed in 0.170604 second(s), 23 queries , Gzip On.

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