如何修复GROMACS运行时出现的‘No such moleculetype SOL’报错

报错根因 #

1. No such moleculetype SOL报错 #

• 你当前的拓扑文件仅定义了蛋白质的分子类型，没有引入水模型对应的SOL分子类型定义文件（itp格式），GROMACS无法识别[molecules]块中列出的SOL分子。

2. Moleculetype 'protein' has no atoms报错 #

• 你手动添加的SOL的[moleculetype]块位置错误，插入到了protein的[moleculetype]声明和它对应的[atoms]块之间，打断了蛋白质的分子定义结构，导致程序判定protein分子类型下无原子条目。

修正操作 #

• 删掉你手动添加的SOL相关[moleculetype]定义，恢复蛋白质分子定义的结构完整性。
• 在拓扑文件的[atomtypes]块之后、蛋白质[moleculetype]声明之前，添加对应水模型的itp引用，常用的tip3p水对应语句为#include "tip3p.itp"，路径需和你引用力场itp的路径保持一致。

  若后续添加离子时出现NA/CL等分子类型不存在的报错，同步添加#include "ions.itp"即可。

• 确认拓扑文件整体顺序符合GROMACS规范：

  力场itp引入 → 水/离子等小分子itp引入 → 蛋白质[moleculetype]定义 → 蛋白质[atoms]/[bonds]等属性块 → [system]块 → [molecules]块
  

• 重新生成tpr文件后再运行genion命令即可正常执行。

内容的提问来源于stack exchange，提问作者rosanna