带结晶水的GROMACS模拟后结晶水断键

模小助

整理一个 ITP/TOP 相关问题，重点看参数和连接方式是否合理。

近一段时间在研究带结晶水的蛋白配体相互作用的动力学模拟问题，结晶水涉及到与配体 蛋白的氢键问题故保留。蛋白结晶水一起用pdb2gmx生成拓扑文件，配体用acpype生成拓扑文件，同时生成配体限制势并将Multiwfn拟合的RESP电荷放了进去

主要想确认：这种拓扑/力场处理是否合理，有没有更稳妥的设置方式？

admin

拓扑问题建议优先检查 include 顺序、[ moleculetype ] 名称、[ molecules ] 数量以及 atomtype 是否重复/冲突。很多问题不是某一行参数值错了，而是两套参数被同时 include。

可以让 gmx grompp -pp processed.top 生成展开后的拓扑，再从 processed.top 里查目标分子和 atomtype，通常能很快定位实际被读取的内容。