求助，模拟冰水时，如何获得top和itp文件，如何避免冰与水之间的真空区

非好汉 · 发表于 Post on 2023-3-16 11:53:05

1.刚开始，从晶体库获得冰晶文件并进行扩胞，用gauss view构建了单个水分子，然后利用packmol进行混合，本来想用sobtop得到top和itp文件，但是看了一些帖子，发现不能整个体系去生成，然后想去尝试用sobtop去得到扩胞后的冰晶的top和itp文件，但是因为涉及用力场还是用hessian，所以在群里问了大家，sob老师说直接用tip4p/ice水模型，我不是很懂他的意思，是想把tip4p/ice的top和itp当成冰的top和itp吗？2.用packmol构建盒子的时候，会出现冰晶胞放不进预设的放冰晶胞的地方（外面再包一层水），调大放冰晶胞的盒子，冰晶胞与水分子之间会有比较大的真空区。
想问一下老师，像这样模拟冰水时，如何获得top和itp文件，又如何避免冰与水之间的真空区？

Lacrimosa · 发表于 Post on 2023-3-16 13:02:49

冰不需要单独构建力场，纯水在低于凝固点下做模拟即可获得冰。因此选用合适的水模型，设置合理的温度压力就足够了。在选择水模型之前要了解该水模型的相图，确定你所需冰结构存在于该模型的相图中，并且找到相应的温度压力。
以ice 1为例，水模型选tip4p/ice就很合适。tip4p水模型的ice 1的晶体结构可以用genice直接生成，而且是gro格式可以直接用。top文件里设置好水分子的数量，在凝固点下先跑NPT平衡一段时间，如果结构能维持住就说明可以了。若要获得固液混合体系，可以给冰中的一部分施加position restraint维持住这部分的结构，然后在凝固点以上模拟，将其余部分融化，由此得到固液混合体系。

牧生 · 发表于 Post on 2023-3-16 18:21:18

你先看懂这个例子
https://zhuanlan.zhihu.com/p/75785937

然后水的参数就抄作业
http://bbs.keinsci.com/thread-32973-1-1.html
真空区不用管，NPT自发会合理接触的

sobereva · 发表于 Post on 2023-3-17 03:53:57

模拟冰和模拟水盒子仅有两个差别：
(1)需要用能够描述冰比较好的水模型的itp文件代替常用的诸如spce.itp
(2)初始结构用冰的结构而非液态水的结构

对于模拟冰-水界面，对冰扩出真空区，gmx solvate加水（不建议packmol，加水明显不如gmx solvate致密）

非好汉 · 发表于 Post on 2023-3-17 11:27:57

Lacrimosa 发表于 2023-3-16 13:02
冰不需要单独构建力场，纯水在低于凝固点下做模拟即可获得冰。因此选用合适的水模型，设置合理的温度压力就 ...

好的，谢谢。但是这个genice，我在Windows上用pip和conda装过好多次了，也用了镜像，一直装不上去，换了别人的电脑也不行，不知道是不是只有linux可以装。

非好汉 · 发表于 Post on 2023-3-17 11:28:32

牧生发表于 2023-3-16 18:21
你先看懂这个例子
https://zhuanlan.zhihu.com/p/75785937

好的，谢谢，已经在尝试了。

非好汉 · 发表于 Post on 2023-3-17 11:30:24

sobereva 发表于 2023-3-17 03:53
模拟冰和模拟水盒子仅有两个差别：
(1)需要用能够描述冰比较好的水模型的itp文件代替常用的诸如spce.itp
...

好的，谢谢老师，我再试试。

牧生 · 发表于 Post on 2023-3-17 11:31:33

本帖最后由牧生于 2023-3-17 13:18 编辑

非好汉发表于 2023-3-17 11:27
好的，谢谢。但是这个genice，我在Windows上用pip和conda装过好多次了，也用了镜像，一直装不上去，换了 ...

在centos下，安装miniconda，然后用pip只需要一步命令即可安装上genice。

非好汉 · 发表于 Post on 2023-3-20 11:09:01

牧生发表于 2023-3-17 11:31
在centos下，安装miniconda，然后用pip只需要一步命令即可安装上genice。

大佬，这个centos下的genice2已经安装好了，结果我又卡在了那篇知乎文章上了，在这一步gmx grompp -f em.mdp -c icebox.pdb -p icebox.top -o em.tpr -maxwarn 2总会错，还请看一下我的文件

冰水.rar (77.81 KB, 下载次数 Times of downloads: 46)

牧生 · 发表于 Post on 2023-3-20 11:23:23

非好汉发表于 2023-3-20 11:09
大佬，这个centos下的genice2已经安装好了，结果我又卡在了那篇知乎文章上了，在这一步gmx grompp -f em. ...

[ molecules ]
; Molecule nmols
SOL 1464

非好汉 · 发表于 Post on 2023-3-20 11:31:39

牧生发表于 2023-3-20 11:23
[ molecules ]
; Molecule nmols
SOL 1464

谢谢大佬，刚刚尝试了一下，生成em.tpr了

非好汉 · 发表于 Post on 2023-3-20 11:37:54

牧生发表于 2023-3-20 11:23
[ molecules ]
; Molecule nmols
SOL 1464

但是进行gmx mdrun -v -deffnm em这一步时，出现了好多pdb文件，

，请问这是怎么回事？

牧生 · 发表于 Post on 2023-3-20 11:40:57

主要问题在于 icebox.pdb中的盒子信息不太对。。
不过我试了一下，问题不大，得到的em.gro打开是没有问题的。

非好汉 · 发表于 Post on 2023-3-20 11:57:41

牧生发表于 2023-3-20 11:40
主要问题在于 icebox.pdb中的盒子信息不太对。。
不过我试了一下，问题不大，得到的em.gro打开是没有问题 ...

好的

非好汉 · 发表于 Post on 2023-3-20 12:44:12

牧生发表于 2023-3-20 11:40
主要问题在于 icebox.pdb中的盒子信息不太对。。
不过我试了一下，问题不大，得到的em.gro打开是没有问题 ...

大佬，还是觉得得明白哪里错了，不然下次还会犯这种错，还请大佬指教。

icebox.pdb (351.96 KB, 下载次数 Times of downloads: 21)

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GROMACS] 求助，模拟冰水时，如何获得top和itp文件，如何避免冰与水之间的真空区

浏览过的版块