计算化学公社

 找回密码 Forget password
 注册 Register

已解决:非标准多肽电荷平衡报错

查看数: 2135 | 评论数: 4 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2022-11-29 15:39

正文摘要:

本帖最后由 wyfgg 于 2022-11-29 17:19 编辑 各位老师好,我目前在进行一个人工合成多肽的MD模拟,多肽中包含两个非标准氨基酸,我参考这个帖子《对有非标准残基的蛋白质做蛋白质-配体分子动力学》和Sob老师的htt ...

回复 Reply

wyfgg 发表于 Post on 2022-11-29 17:17:06
Frozen-Penguin 发表于 2022-11-29 17:01
删掉前几行#include的内容也在往前移,所以是15可能是巧合吧,不过这个和具体错误是什么关系不大。
你的 ...

已解决,非常感谢!
问题根源确实出现在forcefield.itp参数上,按照您的提示修改后,gmx报错
  1. ERROR 1 [file topol.top, line 1102]:                                                                  
  2.   Incorrect number of parameters - found 4, expected 3 or 6 for Proper Dih.  (after the function type).
复制代码
然后我检查了top文件和以前运行的蛋白top文件的差异,发现在[ dihedrals ]中多了一列,可能还是我非标准残基rtp文件的问题,我下面再继续研究一下

再次感谢!
Frozen-Penguin 发表于 Post on 2022-11-29 17:01:13
本帖最后由 Frozen-Penguin 于 2022-11-29 17:04 编辑
wyfgg 发表于 2022-11-29 16:46
老师您好!我随机删去top文件的前端几行,gmx报错显示line 15,因此可能还是这个位置

forcefield.itp ...

删掉前几行#include的内容也在往前移,所以是15可能是巧合吧,不过这个和具体错误是什么关系不大。
你的forcefield.itp里面的参数被注释掉了,好像不是原来的内容,手册里面写的是没有前面的分号的,可能是这个导致的错误。
Frozen-Penguin 发表于 Post on 2022-11-29 16:09:35
报错说的25行应该是考虑#include之后的25行,不是top文件文本的第25行
报错内容是combination rule错误,这个参数在[ defaults ]中,一般出现在力场的forcefield.itp中,是力场本来就有的
在手册的这个页面可以看到相关信息,https://manual.gromacs.org/docum ... y-file-formats.html
这个错误和rtp和hdb应该没有关系,这两个文件里面好像没有combination rule参数

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 02:10 , Processed in 0.181264 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list