求助，想用packmol构建冰水混合物，如何添加四点水模型？

非好汉 · 发表于 Post on 2023-3-30 11:25:44

各位大佬，想请问一下，我在用packmol构建冰水混合物时，如何添加四点水模型。我先用genice构建了tip4p冰的pdb，然后不知道如何向盒子中添加四点水模型，如果先用packmol构建不带水的冰盒子再用gromacs solvate就会将整个盒子填满，没有真空区。

Anbling · 发表于 Post on 2023-11-29 15:14:55

牧生发表于 2023-11-27 12:44
哪里奇怪啊，四点水模型不就是这样的？？虚拟点和氧原子有部分重合，看起来就是这样的啊

谢谢您的帮助！已在win下安装成功genice。

Anbling · 发表于 Post on 2023-11-29 15:12:39

非好汉发表于 2023-11-27 10:33
1.这个图片，看起来确实有点怪，但是我用vesta看起来没有问题，不知道会不会是vmd的自动识别成键，这个问 ...

谢谢您的详细解答。我在win下安装成功了，或许是下面大佬说的安装了conda的原因~

牧生 · 发表于 Post on 2023-11-27 12:44:22

本帖最后由牧生于 2023-11-27 13:15 编辑

非好汉发表于 2023-11-27 10:33
1.这个图片，看起来确实有点怪，但是我用vesta看起来没有问题，不知道会不会是vmd的自动识别成键，这个问 ...

哪里奇怪啊，四点水模型不就是这样的？？虚拟点和氧原子有部分重合，看起来就是这样的啊

genice可以在win下安装，不过要用conda。不仅麻烦，还有可能报错。。

在原生linux下，或者WSL，pip install genice即可安装好

非好汉 · 发表于 Post on 2023-11-27 10:33:49

Anbling 发表于 2023-11-24 16:33
楼主您好。
第一种方式的pdb通过VMD得到下图。（1）我有点疑惑这个O原子为何像裂开了一样，是这个模型的水 ...

1.这个图片，看起来确实有点怪，但是我用vesta看起来没有问题，不知道会不会是vmd的自动识别成键，这个问题我不能解释清楚，你需要问下大佬们
2.genice需要在linux系统里安装（windows里我没有安装成功），网页搜索genice（genice2）是可以搜出来的

Anbling · 发表于 Post on 2023-11-24 16:33:57

楼主您好。
第一种方式的pdb通过VMD得到下图。（1）我有点疑惑这个O原子为何像裂开了一样，是这个模型的水分子就是这样吗？

（2）还有个问题是，TIP4P冰的pdb您说是通过genice得到的，请问有什么教程指路吗，通过“genice使用”关键词在论坛或百度都没有搜到相关信息。
（3）最后，感谢楼主提醒tip4p/ice.itp文件！

非好汉 · 发表于 Post on 2023-11-24 13:08:22

Anbling 发表于 2023-11-24 12:27
楼主您好，对于您回复我的两种方式，我都进行了尝试，想让您再帮忙看一下。
（1）用tip4p.gro先删除多余分 ...

我觉得第一种方式的pdb就可以了，第二种方式你把cryst1这行删了就可以了，在网上可以看一下pdb文件的格式。把单个水分子在可视化软件中打开看看，比如vesta。其次就是水分子的拓扑文件引用对就行，在这里应该是tip4p/ice.itp。在我看来，pdb文件就是提供原子坐标，itp文件提供键长键角等。

Anbling · 发表于 Post on 2023-11-24 12:27:39

楼主您好，对于您回复我的两种方式，我都进行了尝试，想让您再帮忙看一下。
（1）用tip4p.gro先删除多余分子单个四点水的gro，然后通过openbabel转化成pdb文件，最后删除盒子信息。

（2）用gmx solvate -cs tip4p.gro -maxsol 1 -box 1 1 1 -o box.gro填充单个分子再处理，抠出单个分子。

由于第一次回复，有些着急，前面重复的问题还请您谅解，希望您在有空时可以再帮我解答一下疑惑。感谢！

Anbling · 发表于 Post on 2023-11-24 11:14:57

通过gmx solvate -cs tip4p.gro -maxsol 1 -box 1 1 1 -o box.gro我得到了下图内容，（1）之后是要删除1、2、7行得到单个分子吗？（2）删除后也只是.gro文件，是也要用openbabel转换为pdb文件吗？还请大佬再帮下忙。

Anbling · 发表于 Post on 2023-11-24 11:14:24

非好汉发表于 2023-11-23 18:27
我是用tip4p.gro先删除多余分子单个四点水的gro，然后通过openbabel转化成pdb文件，最后删除盒子信息。也 ...

感谢您的回复！您给了两种方案，我先尝试了一下第二种。因为回复我不会插入图片所以又在下方发了个帖子。

非好汉 · 发表于 Post on 2023-11-23 18:27:43

Anbling 发表于 2023-11-23 12:57
楼主，看到你上面的问题，与我现在遇到的问题很相似，想请教一下您我的问题。我想构建0.9%nacl溶液的冰水混 ...

我是用tip4p.gro先删除多余分子单个四点水的gro，然后通过openbabel转化成pdb文件，最后删除盒子信息。也可以先转成pdb再删；或者用gmx solvate -cs tip4p.gro -maxsol 1 -box 1 1 1 -o box.gro填充单个分子再处理，抠出单个分子。

Anbling · 发表于 Post on 2023-11-23 12:57:40

楼主，看到你上面的问题，与我现在遇到的问题很相似，想请教一下您我的问题。我想构建0.9%nacl溶液的冰水混合物。目前我的想法是建立一个4.4nm的立方体，立方体的中间层为冰，上下两层为水（水中包含钠离子或氯离子）。因为牵扯到冰，想用TIP4P/ICE模型。看到其它的帖子要用packmol构建自己的模型。目前已有na.pdb，cl.pdb。但是TIP4P/ICE.pdb我一直没找到，看您是通过genice构建了tip4p冰的pdb，通过删除tip4p.gro得到tip4p水的pdb（不知道我理解的对不对）。如果是要水的pdb是只剩图红框中内容并另存为.pdb格式吗？（packmol的输入文件）

非好汉 · 发表于 Post on 2023-4-3 17:34:10

谢谢各位大佬的指导，最终使用gromacs文件夹下的tip4p.gro（删的只剩一个SOL（即四个原子））成功了，然后使用packmol构建的模型

shadowcrystal · 发表于 Post on 2023-4-2 15:15:49

本帖最后由 shadowcrystal 于 2023-4-2 15:16 编辑

非好汉发表于 2023-3-30 18:42
大佬，我主要想用packmol构建一个冰水混合溶液（水中有溶质），但是我想加入四点水模型的冰和水，冰已经 ...

我是这样做的：gmx solvate -cp ice.gro -cs tip4p -radius 0.2 -p topol.top -o ice_water.gro
如果需要四点水的top文件，我放附件里了

Lacrimosa · 发表于 Post on 2023-4-1 11:05:00

非好汉发表于 2023-3-31 21:34
大佬，试过了，没有成功，tip4p.gro是个有很多分子的文件，然后去尝试只留里面的1SOL，然后用openbabel转 ...

gmx editconf不就可以吗？？？建议你先把tip4p.gro转成pdb再删除其他原子

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

[GROMACS] 求助，想用packmol构建冰水混合物，如何添加四点水模型？