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

牧生

非好汉 发表于 2023-3-20 12:44
大佬，还是觉得得明白哪里错了，不然下次还会犯这种错，还请大佬指教。

pdb的盒子信息，你试试把小数点加上去，比如CRYST1   30.000  28.000   75.000  90.00  90.00  90.00 P 1           1

非好汉

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

不好意思，找不到昨天的回复贴了，昨天您让我尝试把小数点加上CRYST1   28.000   24.000   69.000  90.00  90.00  90.00 P 1           1，但是还是出现了很多pdb文件，然后在28，24，23的长方体内放600个水分子是不是有点多，根据阿伏伽德罗常数和水的密度来算，这个长方体能放515左右个水分子

牧生

把盒子稍微改大一点就好了啊，反正NPT的时候，会自发收缩到合理的大小。用下面这个参数就可以啊，不要试图自己手动去改一个“恰到好处”的大小。

CRYST1   30.000  28.000   75.000  90.00  90.00  90.00 P 1           1

非好汉

牧生 发表于 2023-3-21 11:15
把盒子稍微改大一点就好了啊，反正NPT的时候，会自发收缩到合理的大小。用下面这个参数就可以啊，不要试图 ...

好的，大佬

liudinghao

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

老师，对冰扩出真空区，gmx solvate加水之后，把一些水分子加入到我的冰结构当中了，这该怎么处理呢？

sobereva

liudinghao 发表于 2025-4-25 16:51
老师，对冰扩出真空区，gmx solvate加水之后，把一些水分子加入到我的冰结构当中了，这该怎么处理呢？

利用VMD去除
参考北京科音分子动力学与GROMACS培训班（http://www.keinsci.com/KGMX）里模拟冰的例子用到的做法