24小时热门版块排行榜    

查看: 469  |  回复: 3
当前主题已经存档。

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

[交流] 【求助】gromacs中“non-default goups in your mdp file”错误

大家好,请教一下:

我用gromacs运行,mdrun时md.mdp文件如下:
; Berendsen temperature coupling is on in three groups
Tcoupl                  =  v-rescale
tau_t                    =  0.1                   0.1           0.1                 0.1
tc_grps                    =  Protein_A           SIM          SOL         NA+
ref_t                     =  310                  310         310         310

top文件如下:
[ molecules ]
; Compound        #mols
Protein_A             1
SIM                            1
SOL                      46639
NA+                            2

但是在运行时却出现如下错误,请问是哪里有问题了?

Fatal error:
Group Protein_A not found in indexfile.
Maybe you have non-default goups in your mdp file, while not using the '-n' option of grompp.
In that case use the '-n' option.

[ Last edited by lei0736 on 2009-11-13 at 21:46 ]
回复此楼

» 猜你喜欢

个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

armeria

铁虫 (初入文坛)

★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+5,VIP+0):谢谢 新人加油 5-27 09:21
top文件中[ molecules ]中的Compound,比如Protein_A ,并不是index文件中的Group。

Gromacs默认的index文件是每次运行grompp等程序的时候根据数据库里的aminoacids.dat等文件自动生成的,主要是根据残基的代码区分蛋白,非蛋白,溶剂,离子等,总共有十个左右的组。这个默认的index文件只是计算需要,并不保存。如果您不需要特定的组,可以直接用Gromacs默认的组,比如Protein和Non Protein。

如果您想区分特定的组,比如这里的Protein_A等,需要用make_ndx命令生成您所需要的index文件,这样的index文件不仅包含了Gromacs默认的组,还包含了您自己设定的组,这样grompp就可以调用这个生成的index文件去匹配mdp文件中的组名了。

[ Last edited by armeria on 2009-5-27 at 02:17 ]
2楼2009-05-27 02:11:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gwdavid

木虫 (著名写手)

小木虫散兵坑坑长

非常感谢您!
个人博客:[url]http://blog.sina.com.cn/gaiwei[/url]
3楼2009-05-27 08:29:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
gwdavid(金币+4,VIP+0):非常感谢! 6-8 18:25
从你md.mdp文件来看,你所有的group的温度都是一样的,耦合的频率也是一样的,那么就没有必要吧他们都写出来,你这样写就可以搞定:
Tcoupl                  =  v-rescale
tau_t                    =  0.1           
tc_grps                    =  System
ref_t                     =  310   

如果你按照你的方式写,出现的问题可能的原因是你的ITP文件的[ atoms ]中残基组(resid)这一项写的不是Protein_A
4楼2009-06-08 16:58:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gwdavid 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自科面上基金字体 +6 iwuli 2026-03-12 7/350 2026-03-16 21:18 by sculhf
[考研] 环境工程调剂 +6 大可digkids 2026-03-16 6/300 2026-03-16 17:16 by barlinike
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +3 Losir 2026-03-12 3/150 2026-03-14 12:11 by 热情沙漠
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 考研材料与化工,求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 329求调剂 +3 miaodesi 2026-03-12 4/200 2026-03-13 20:53 by 18595523086
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考研] 08食品或轻工求调剂,本科发表3篇sci一区top论文,一志愿南师大食品科学与工程 +3 我是一个兵, 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见