24小时热门版块排行榜    

Znn3bq.jpeg
查看: 489  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博自荐 +4 Linxia林夏 2026-04-13 4/200 2026-04-19 19:55 by Equinoxhua
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 求调剂 +6 苦命人。。。 2026-04-18 7/350 2026-04-19 16:27 by 中豫男
[论文投稿] 有没有接收比较快的sci期刊呀,最好在一个月之内的,研三孩子求毕业 20+4 之护着 2026-04-16 6/300 2026-04-19 13:00 by Aaron_zyn
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 320求调剂 +5 深郊akm 2026-04-17 5/250 2026-04-18 19:52 by 王珺璞
[考研] 297,工科调剂? +5 河南农业大学-能 2026-04-14 5/250 2026-04-18 15:17 by Equinoxhua
[考研] 一志愿华中农业071010,320求调剂 +17 困困困困坤坤 2026-04-14 19/950 2026-04-17 20:08 by 关一盏灯cd
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 297工科调剂? +14 河南农业大学-能 2026-04-13 15/750 2026-04-15 13:25 by 黑科技矿业
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
信息提示
请填处理意见