跑npt 时,指令如下
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -o npt.tpr
遇到的问题:
Fatal error:Group water referenced in the .mdp file was not found in the index file.Group names must match either [moleculetype] names or custom index group names, in which case you must supply an index file to the '-n' option of grompp.
致命错误:中引用的水组。在索引文件中找不到mdp文件。组名必须与[moleculetype]名或自定义索引组名匹配,在这种情况下,您必须为grompp的'-n '选项提供一个索引文件。
跑的模拟没有往盒子里面加水,而是加葡萄糖和柠檬酸
求详细解答
参考GPT和自己的思路:
根据您的问题描述,错误提示是由于在 .mdp 文件中引用了一个名为 "water" 的分组,但在索引文件中找不到这个分组导致的。由于您的模拟系统中没有水分子,因此需要更新 .mdp 文件中的分组信息。
首先,您需要用命令 gmx make_ndx 来创建一个新的索引文件,其中包含您的葡萄糖和柠檬酸分子。在命令行中输入以下命令以创建一个新的索引文件:
gmx make_ndx -f nvt.gro -o index.ndx
在输入以上命令后,将会打开交互式编辑器。您需要使用以下命令来创建一个名为 "System" 的新组,并将葡萄糖和柠檬酸添加到组中:
a 1-XXXX ;(XXXX为葡萄糖在文件中的原子编号)
a 1-XXXX ;(XXXX为柠檬酸在文件中的原子编号)
name 2 System
然后您需要在 npt.mdp 文件中更新水分组的名称。找到以下行:
tc-grps = Protein Non-Protein SOL
并将其更新为以下内容:
tc-grps = Protein Non-Protein System
然后再次执行 gmx grompp 命令,将新的索引文件和更新后的 .mdp 输入文件提供给该命令:
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -n index.ndx -o npt.tpr
这样应该能够成功运行您的 NPT 模拟,而无需使用水分子。
参考GPT和自己的思路:
这个问题的原因是你在 MDP 文件中引用了名为 "water" 的组,但是该组在索引文件中不存在。因此,GROMACS 报错了。
解决方法是在索引文件中添加名为 "water" 的组。你可以使用 make_ndx 命令手动添加此组,或者在生成索引文件时,在 [ molecules ] 部分中添加水分子信息,例如:
[ molecules ]
; Compound #mols
Protein 1
GLC 100
ACN 100
其中 "GLC" 和 "ACN" 分别表示葡萄糖和柠檬酸,你可以根据你的具体模拟系统进行修改。
补充说明:在使用 make_ndx 命令手动生成索引文件时,可以使用 "a SOL" 命令来将所有的水分子选中为一个组。
该回答引用GPTᴼᴾᴱᴺᴬᴵ
这个错误提示是因为在你的 npt.mdp 文件中引用了一个叫做 water 的组,但是在你的拓扑文件 topol.top 中并没有定义这个组,所以 grompp 执行时会出错。
解决这个问题有两种方法:
1.将 npt.mdp 文件中的 water 组改成拓扑文件中存在的组名。
2.在 topol.top 文件中定义一个叫做 water 的组,可以使用 make_ndx 工具生成自定义的索引文件 index.ndx,再将其用于 grompp 命令中的 -n 选项,例如:
gmx make_ndx -f topol.top -o index.ndx
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -n index.ndx -o npt.tpr
在执行 make_ndx 命令时,会出现一个交互式命令行界面,可以在其中选择需要定义的组,例如选择所有的水分子,命令行输入:
> a OW
> q
然后在 topol.top 文件中添加以下行:
[ water ]
OW
然后在 grompp 命令中使用 -n 选项引用该索引文件即可。
gmx grompp -f npt.mdp -c nvt.gro -r nvt.gro -t nvt.cpt -p topol.top -n index.ndx -o npt.tpr
把mdp文件中的 tc-grps =water Non-water 改为 Protein Non-Protein,下面几行出现的唯一的 water 也改为 Protein