24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 586  |  回复: 7
当前主题已经存档。

wangyibo20

金虫 (小有名气)

[交流] 【求助】gromacs出错

【求助】我加完水盒子后预处理grompp出错:

Fatal error:
[ file "/usr/local/gromacs/share/gromacs/top/spce.itp", line 42 ]:
             Atom index (1) in settles out of bounds (1-0)
是怎么回事?
命令为:grompp_d -f em.mdp -c 1_box.pdb -p 1.top -o 1_pp.tpr

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

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★
lei0736(金币+2,VIP+0):谢谢 3-20 09:43
你使用的是什么力场,有添加amber力场吗,可能是你的include的次序有问题,能不能看下你的top文件呢
2楼2009-03-19 23:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyibo20

金虫 (小有名气)

我用的Amber03力场

多谢了,top文件如下:
3楼2009-03-20 08:10:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 3-21 19:23
你的top中有包含spce.itp,这个与amber连用可能出问题,你写个新的水的itp文件,命名为amber-spc.itp把它包含在你的top中,把原来包含那个spce.itp注销掉,新的itp内容如下:
[ moleculetype ]
; molname   nrexcl
SOL      2

[ atoms ]
;   nr   type  resnr residue  atom   cgnr     charge       mass
     1  amber99_54   1  SOL     OW      1      -0.82  15.99940
     2  amber99_55   1  SOL    HW1      1       0.41   1.00800
     3  amber99_55   1  SOL    HW2      1       0.41   1.00800

#ifdef FLEXIBLE
[ bonds ]
; i   j  funct length   force.c.
1  2  1  0.1   345000   0.1     345000
1  3  1  0.1   345000   0.1     345000

[ angles ]
; i   j  k  funct angle force.c.
2  1  3  1  109.47   383   109.47   383
#else
[ settles ]
; OW  funct doh   dhh
1  1  0.1   0.16333

[ exclusions ]
1  2  3
2  1  3
3  1  2
#endif


你调下试试
4楼2009-03-20 10:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyibo20

金虫 (小有名气)

又出现别的错误了。我的配体用MOL表示,可很多时候都提示没有这个残基类型,那我应该用什么表示配体?

另外我重做时用Amber99sb力场,tip3p水模型时,在grompp这一步提示
原子类型OWT3找不到,又是怎么回事?

[ Last edited by wangyibo20 on 2009-3-20 at 19:01 ]
5楼2009-03-20 18:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nkyx

铁杆木虫 (正式写手)

★ ★
lei0736(金币+2,VIP+0):谢谢 3-24 09:33
配体残基类型要和你配体的top文件的一致。
直接用文本编辑器打开Amber99sb力场文件看看里面,水的原子类型;或者直接找到水的各种类型文件看看。
6楼2009-03-23 10:49:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyibo20

金虫 (小有名气)

关键是我没找到哪个文件里有这种水原子呀
7楼2009-03-23 19:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)

★ ★ ★ ★
lei0736(金币+4,VIP+0):谢谢 3-26 08:50
是关键词的问题,残基名,类型名,这类是程序的关键词,如果你的top里的与你调用力场里的不一样,程序就会不识别,你的top文件里有包含多个*.itp文件,这些文件要嘛在你的本任务的目录下,要嘛在gromacs的安装目录下,你可以在安装目录的/share/top/里看到很多力场文件
8楼2009-03-25 00:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyibo20 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学363调剂求助 +4 fanzhang6666 2026-04-06 5/250 2026-04-07 10:12 by dongzh2009
[考研] 292求调剂 +4 lilllllxccc 2026-04-05 5/250 2026-04-07 09:29 by 纺大杨老师
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 081200-11408-367学硕求调剂 +4 1_2_3111 2026-04-06 4/200 2026-04-07 08:13 by jp9609
[考研] 考研调剂 +7 15615482637 2026-04-04 7/350 2026-04-06 22:56 by chenzhimin
[考研] 调剂 +3 mcbbc 2026-04-06 4/200 2026-04-06 20:58 by lbsjt
[考研] 考研调剂一志愿北矿安全科学与工程学硕316 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 081700,311,求调剂 +16 冬十三 2026-04-04 17/850 2026-04-06 14:56 by 尚水阁主
[考研] 材料与化工371求调剂 +14 陪琳看海 2026-04-04 15/750 2026-04-06 06:59 by houyaoxu
[考研] 315求调剂 +14 欣喜777 2026-04-04 15/750 2026-04-06 06:25 by houyaoxu
[考研] 272求调剂 +4 电气李 2026-04-05 4/200 2026-04-05 10:41 by lbsjt
[考研] [调剂信息]085408光电信息 求调剂 总分291分数一英一 +3 iz11az 2026-04-02 3/150 2026-04-04 19:09 by 蓝云思雨
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 325分化学调剂 +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[考研] 材料调剂 +4 一样YWY 2026-04-03 4/200 2026-04-03 09:48 by 蓝云思雨
[考研] 一志愿大工学硕,求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 303分 0807学硕求调剂 +3 TYC3632 2026-04-01 3/150 2026-04-01 19:24 by lwk2004
信息提示
请填处理意见