计算化学公社

 找回密码 Forget password
 注册 Register

求助用Gromacs进行蛋白-配体复合物的动力学模拟在加离子后用grompp时报错

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

正文摘要:

能量最小化时,gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr指令报错。 能量最小化前的添加抗衡离子步骤中,我用6个CL替换掉了6个溶剂分子,使体系达到电荷平衡,并得到更新后的top文件; ...

回复 Reply

Frozen-Penguin 发表于 Post on 2022-10-18 09:55:19
LvYang 发表于 2022-10-18 00:28
感谢老师!问题已解决

ions.itp文件中的117行氯原子一开始的写法是CLA,但是gromacs中输入指令识别的 ...

这种情况一般应该改top文件而不是itp文件,不过对氯离子来说影响应该不大
LvYang 发表于 Post on 2022-10-18 00:28:18
Frozen-Penguin 发表于 2022-10-17 23:27
在力场文件中的ions.itp里面找定义CL的部分
如果没有就重新下载力场文件,如果有应该不会报错

感谢老师!问题已解决

ions.itp文件中的117行氯原子一开始的写法是CLA,但是gromacs中输入指令识别的原子格式是CL,删除A,报错解除。
Frozen-Penguin 发表于 Post on 2022-10-17 23:27:58
LvYang 发表于 2022-10-17 20:19
力场文件charmm36.ff和top文件是处在同一文件目录下的,但是报错的原因我无法理解,平衡离子的过程中我确 ...

在力场文件中的ions.itp里面找定义CL的部分
如果没有就重新下载力场文件,如果有应该不会报错

CL.PNG (10.38 KB, 下载次数 Times of downloads: 7)

CL.PNG
LvYang 发表于 Post on 2022-10-17 20:19:32
Frozen-Penguin 发表于 2022-10-16 13:54
你这样写的话,力场文件charmm36.ff需要放在top文件同目录下,力场文件的目录中需要包含所有需要的itp文 ...

力场文件charmm36.ff和top文件是处在同一文件目录下的,但是报错的原因我无法理解,平衡离子的过程中我确实添加了六个CL,为什么在能量最小化时指令却报错说top文件没有CL分子类型呢?
还是说我需要自己手动创建一个6CL的itp文件放在力场文件charmm36.ff下,如果是,这种文件该怎样创建呢?
Frozen-Penguin 发表于 Post on 2022-10-16 13:54:50
LvYang 发表于 2022-10-16 13:16
感谢您的解惑,我在top文件中找到了引入ions参数的指令

但是依旧不能解决该问题,引入ions.itp文件后 ...

你这样写的话,力场文件charmm36.ff需要放在top文件同目录下,力场文件的目录中需要包含所有需要的itp文件
LvYang 发表于 Post on 2022-10-16 13:16:07
Frozen-Penguin 发表于 2022-10-16 12:29
去力场文件目录下找ions.itp在哪
或者直接用pdb2gmx产生一个top文件,参考里面是怎么引入ions参数的

感谢您的解惑,我在top文件中找到了引入ions参数的指令

但是依旧不能解决该问题,引入ions.itp文件后,进行gmx grompp -f em.mdp -c solv_ions.gro -p topol.top -o em.tpr指令,依旧会有上述同样的报错。
Frozen-Penguin 发表于 Post on 2022-10-16 12:29:54
LvYang 发表于 2022-10-16 12:25
好的,sob老师,我下次发帖一定注意标题格式!

去力场文件目录下找ions.itp在哪
或者直接用pdb2gmx产生一个top文件,参考里面是怎么引入ions参数的
LvYang 发表于 Post on 2022-10-16 12:25:46
sobereva 发表于 2022-10-16 00:00
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此 ...

好的,sob老师,我下次发帖一定注意标题格式!
sobereva 发表于 Post on 2022-10-16 00:00:51
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示,求助帖必须在帖子标题明确体现出此帖内容是求助或提问,并清楚、准确反映出帖子具体内容,避免有任何歧义和含糊性,仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “利用Gromacs进行蛋白-配体复合物的分子动力学模拟” 改了,以后务必注意,下次将删帖+扣分处理

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

GMT+8, 2026-2-24 17:59 , Processed in 0.178349 second(s), 25 queries , Gzip On.

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