24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 475  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程302分求调剂 +10 zyx上岸! 2026-04-04 10/500 2026-04-06 22:49 by chenzhimin
[考研] 308求调剂 +13 倘若起风了呢 2026-04-05 13/650 2026-04-06 14:20 by 蒋皓禹
[考研] 315求调剂 +14 欣喜777 2026-04-04 15/750 2026-04-06 06:25 by houyaoxu
[考研] 一志愿同济大学323分(080500)求调剂 +8 yikeniu 2026-04-01 8/400 2026-04-05 18:15 by cql1109
[考研] 286求调剂 +3 草木不言 2026-04-04 3/150 2026-04-04 22:40 by lbsjt
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 278求调剂 +14 范婷娜 2026-04-04 15/750 2026-04-04 22:15 by lqwchd
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 11408,284分,二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[考研] 321求调剂 +17 y-yh 2026-04-01 20/1000 2026-04-03 12:57 by y-yh
[考研] 316求调剂 +14 舟自梗 2026-04-01 18/900 2026-04-03 10:28 by linyelide
[考研] 262求调剂 +6 励志一定发文章 2026-04-02 7/350 2026-04-03 09:54 by linyelide
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 江苏苏北高校诚邀调剂同学 +3 zzll406 2026-03-31 3/150 2026-03-31 16:54 by 及时行乐fan
信息提示
请填处理意见