24小时热门版块排行榜    

查看: 581  |  回复: 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 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:10 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:59 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:32 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 by Z7qfCsdIiD1T
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 01:13 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 01:02 by Z7qfCsdIiD1T
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 22:39 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 5sJFDKp2kwuI 2026-11-06 14/700 2026-11-09 11:52 by GwrQVVR6Oenx
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-08 23:01 by l48hFd8Nwae2
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
信息提示
请填处理意见