计算化学公社

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

[GROMACS] 能量最小化进行不下去

[复制链接 Copy URL]

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

本帖最后由 吧唧爱吃糖 于 2022-6-23 16:06 编辑

求助各位老师,在进行能量最小化的时候未达到设定要求就停止了。
想问一下各位老师,无论怎么跑都不可能到我设定的条件能量最小化算不算已经达到最小了那?

zzzzzzee.png (294.63 KB, 下载次数 Times of downloads: 6)

这是日志文件提醒

这是日志文件提醒

zzzzzxxxxx.png (354.51 KB, 下载次数 Times of downloads: 5)

这个是pdb文件的构型

这个是pdb文件的构型

min.mdp

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

这是mdp文件

1.gro

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

这是gro文件

419

帖子

1

威望

2551

eV
积分
2990

Level 5 (御坂)

2#
发表于 Post on 2022-6-23 16:12:50 | 只看该作者 Only view this author
如果之前用的单精度浮点数版本,尝试使用双精度浮点数版本(gmx -> gmx_d,需自行编译)

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

3#
发表于 Post on 2022-6-23 16:22:48 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-6-23 16:24 编辑

看了一下,你的gro和top中分子出现的顺序是不对应的,所以你选择not water的时候,最下面一层晶体的水就没有了,而上面的乱七八糟的地方,就不能识别为水。

这个gro看起来可能像是一个单胞扩展为超胞得到的,所以解决办法是重复多次单胞中的分子的顺序和数量(首先自己要确定好单胞中的分子顺序)。
比如2X2X2的甲烷水合物超胞,在top中就写成

SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
SOL         46
CH4        8
理解到这个意思以后,就可以自己发挥一下
又菜又爱玩

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2022-6-23 18:28:22 | 只看该作者 Only view this author
牧生 发表于 2022-6-23 16:22
看了一下,你的gro和top中分子出现的顺序是不对应的,所以你选择not water的时候,最下面一层晶体的水就没 ...

啥?师兄我没放top文件啊,师兄是不是回错人了?

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2022-6-23 18:29:03 | 只看该作者 Only view this author
snljty2 发表于 2022-6-23 16:12
如果之前用的单精度浮点数版本,尝试使用双精度浮点数版本(gmx -> gmx_d,需自行编译)

您好,使用的是双精度,我使用的是gmx_mpi_d

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

6#
发表于 Post on 2022-6-23 20:10:57 | 只看该作者 Only view this author
吧唧爱吃糖 发表于 2022-6-23 18:28
啥?师兄我没放top文件啊,师兄是不是回错人了?

  用vmd打开1.gro,你选择not water的时候,最下面一层晶体的水就没有了,而上面的乱七八糟的地方,就不能识别为水。




这个gro看起来可能像是一个胞体扩展为超胞得到的,所以解决办法是重复多次单胞中的分子的顺序和数量(首先自己要确定好单胞中的分子顺序)。
又菜又爱玩

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2022-6-23 20:19:27 | 只看该作者 Only view this author
牧生 发表于 2022-6-23 20:10
用vmd打开1.gro,你选择not water的时候,最下面一层晶体的水就没有了,而上面的乱七八糟的地方,就不 ...

师兄我是先建立了下层界面的模型单胞然后扩胞成3*2*2的然后和上面的相应尺寸的游离分子拼成的整个pdb然后转成了gro。师兄的解决办法我没看懂,是说先用下层水合物模型单胞和相应尺寸的游离分子拼成单胞然后再把整个整体扩成3*2*2的吗?

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

8#
发表于 Post on 2022-6-23 21:43:18 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-6-24 07:16 编辑
吧唧爱吃糖 发表于 2022-6-23 20:19
师兄我是先建立了下层界面的模型单胞然后扩胞成3*2*2的然后和上面的相应尺寸的游离分子拼成的整个pdb然后 ...

大概明白你的步骤了。
你是让3*2*2的超胞在下层,上层无规律排列的是你自己加上去的,想模拟在水合物晶种条件下,看水合物的生长。
用packmol把晶种和上层放进一个盒子,top的顺序就应该为如下四行(具体顺序应该看你packmol里面加入的顺序)

晶种里面水分子或者客体分子      
晶种里面客体分子或者水分子
(或许要需要多重复几次晶体单胞中分子的顺序和个数)
加入的水分子或者客体分子
加入的客体分子或者水分子

如果还是不明白,就看看这个压缩包里面,topol.top和input.inp里面的对应关系,这里面topol.top的写的看似很杂乱,其实就只是在简单的重复单胞中的水分子和CO2分子出现的顺序和个数。多看几遍,你应该就能明白了
https://doi.org/10.5281/zenodo.6558623

又菜又爱玩

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2022-6-24 09:21:58 | 只看该作者 Only view this author
牧生 发表于 2022-6-23 21:43
大概明白你的步骤了。
你是让3*2*2的超胞在下层,上层无规律排列的是你自己加上去的,想模拟在水合物晶 ...

师兄,我看明白了。我用packmol试一下,但是我检查了我的top文件中分子顺序和gro文件中是一致的,因为我用editconf 命令将pdb转化为gro的时候有几个分子丢失了我还特意改过top文件。为什么用gromacs搭建的模型会出现能量最小化提前截止的情况那?

下面是我建模的命令

gmx_mpi_d insert-molecules -ci H-c.pdb -nmol 480 -box 6.1887 5.897 1.9788 -o 480H.pdb
gmx_mpi_d insert-molecules -ci Methylcyclohexane-c.pdb -nmol 8 -f 480H.pdb -o 480H-8M.pdb
gmx_mpi_d insert-molecules -ci H2O-c.pdb -nmol 1060 -f 480H-8M.pdb -o 480H-8M-1060H2O.pdb
gmx_mpi_d editconf -f 480H-8M-1060H2O.pdb -box 6.1887 8.279 1.9788 -center 0 5.3305 0 -o 480H-8M-1060H2O-c.pdb
gmx_mpi_d editconf -f M-H2O-322.pdb -box 6.1887 8.279 1.9788 -center 0 1.191 0 -o M-H2O-322-c.pdb
gmx_mpi_d solvate -cp M-H2O-322-c.pdb -cs 480H-8M-1060H2O-c.pdb -o SH-M-480H-8M-1060H2O-zkc.pdb
这个建模之后显示插入的分子数量是对得上的,但是转成gro之后1060和480都少了,我就改了top文件然后grompp是没问题的,但是mdrun就提示能量最小化提前结束了

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

10#
发表于 Post on 2022-6-24 09:32:33 | 只看该作者 Only view this author
本帖最后由 牧生 于 2022-6-24 09:34 编辑
吧唧爱吃糖 发表于 2022-6-24 09:21
师兄,我看明白了。我用packmol试一下,但是我检查了我的top文件中分子顺序和gro文件中是一致的,因为我 ...
  1. 用editconf 命令将pdb转化为gro的时候有几个分子丢失
复制代码


发生这种情况,已经就有错误了。

你的操作中,不需要反复用insert-molecules,也不要用solvate去添加水,否则会偶尔有几个水强行加到冰晶里面,导致分子间距离太近而爆炸。。只需要一步packmol以后就可以了,而且绝对不会错


好好看那个压缩包里面的步骤
又菜又爱玩

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2022-6-24 10:20:52 | 只看该作者 Only view this author
牧生 发表于 2022-6-24 09:32
发生这种情况,已经就有错误了。

你的操作中,不需要反复用insert-molecules,也不要用solvate去 ...

好的师兄,我去尝试一下packmol谢谢师兄的指点

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2022-6-24 12:46:33 | 只看该作者 Only view this author
本帖最后由 吧唧爱吃糖 于 2022-6-24 13:23 编辑
牧生 发表于 2022-6-24 09:32
发生这种情况,已经就有错误了。

你的操作中,不需要反复用insert-molecules,也不要用solvate去 ...

师兄中午好,我使用Packmol尝试建立一个大小为61.887 70.0 19.788的盒子将我的水合物相固定在61.887 23.82 19.788范围内,但是我发现生成的pdb文件中盒子是只有1.2的。我将我的文件附在下面,希望师兄指点。

new_1.inp

175 Bytes, 下载次数 Times of downloads: 2

M-H2O-322.pdb

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

mix.pdb

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

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

13#
发表于 Post on 2022-6-24 13:18:51 | 只看该作者 Only view this author
吧唧爱吃糖 发表于 2022-6-24 12:46
师兄中午好,我使用Packmol尝试建立一个大小为61.887 70.0 19.788的盒子将我的水合物相固定在61.887 23.8 ...

1、我不太明白你把水合物相固定在某个范围是为了什么。
2、我让你看的那个压缩包里面的input.inp,已经非常清晰明了。但我感觉你并没有掌握到packmol的精髓。你现在需要做的是,自己在一张纸上,画一个立方体,把其中一个点定位0,0,0的原点,然后找出顶点的坐标,将这一段放晶体,再把立方体延伸,再找各个顶点的坐标,把延伸这一段放水和其余分子。如果还不明白,那就在纸上把input.inp里面的坐标画成立方体看看。
3、61.887 23.82 19.788这三个坐标没必要那么精确,可以写成62 24 20
4、肉眼来看,这个像是I型水合物,那么水分子形成的笼子的尺寸,足够放入一个甲苯分子吗??
又菜又爱玩

54

帖子

0

威望

131

eV
积分
185

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2022-6-24 14:35:59 | 只看该作者 Only view this author
本帖最后由 吧唧爱吃糖 于 2022-6-24 14:51 编辑
牧生 发表于 2022-6-24 13:18
1、我不太明白你把水合物相固定在某个范围是为了什么。
2、我让你看的那个压缩包里面的input.inp,已经 ...

我仔细看了师兄推荐的文件,然后修改了我的inp文件但是不知道为什么进行不下去。提示错误如下,麻烦师兄指点,十分感谢师兄。


师兄我在inp文件中在水合物相加fixed之后成功了,但是水合物相和游离分子相之间有真空层应该怎么解决那?

1111111.png (88.05 KB, 下载次数 Times of downloads: 5)

1111111.png

3AY_N(1J41658@AA5RY}UXO.png (512.11 KB, 下载次数 Times of downloads: 3)

3AY_N(1J41658@AA5RY}UXO.png

H.pdb

458 Bytes, 下载次数 Times of downloads: 1

H2O.pdb

556 Bytes, 下载次数 Times of downloads: 1

Methylcyclohexane.pdb

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

new_1.inp

395 Bytes, 下载次数 Times of downloads: 3

Me-SH-zj-322.pdb

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

1480

帖子

0

威望

4550

eV
积分
6030

Level 6 (一方通行)

15#
发表于 Post on 2022-6-24 14:59:59 | 只看该作者 Only view this author
吧唧爱吃糖 发表于 2022-6-24 14:35
我仔细看了师兄推荐的文件,然后修改了我的inp文件但是不知道为什么进行不下去。提示错误如下,麻烦师兄 ...

基本合理了,自己再琢磨一下就行了。。如果还不行,尝试把盒子再略加大点。真空层没关系,NPT的时候会自发调节到合理的位置。
又菜又爱玩

本版积分规则 Credits rule

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

GMT+8, 2025-8-17 01:46 , Processed in 0.215966 second(s), 23 queries , Gzip On.

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