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

Anbling

楼主，看到你上面的问题，与我现在遇到的问题很相似，想请教一下您我的问题。我想构建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的输入文件）

非好汉

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

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

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

Anbling

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

Anbling

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

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

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

非好汉

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

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

Anbling

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

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

非好汉

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

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

牧生

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

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


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

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

Anbling

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

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

Anbling

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

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