应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Gromacs/Amber/NAMD » 【求助】gromacs出错
81/1返回列表
查看: 699  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

» 猜你喜欢
1楼 2009-03-19 18:52:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuefei06

木虫 (正式写手)
★ ★
lei0736(金币+2,VIP+0):谢谢 3-20 09:43
你使用的是什么力场，有添加amber力场吗，可能是你的include的次序有问题，能不能看下你的top文件呢
赞一下(10人) 回复此楼
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


你调下试试
赞一下(10人) 回复此楼
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力场文件看看里面，水的原子类型；或者直接找到水的各种类型文件看看。
赞一下(10人) 回复此楼
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/里看到很多力场文件
赞一下(10人) 回复此楼
8楼2009-03-25 00:25:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wangyibo20 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:40 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 FChKX8meu1ds 2026-11-09 5/250 2026-11-11 04:23 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 cqv1zKnzd7hr 2026-11-09 7/350 2026-11-11 04:21 by Z7qfCsdIiD1T
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 04:12 by Z7qfCsdIiD1T
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 04:01 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一区文章，我: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 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[考博] 申请环境，微生物方向博士 +3 浅笑09 2026-11-08 3/150 2026-11-10 12:44 by 可淡不可忘
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 zSbvVX9JFl9J 2026-11-09 4/200 2026-11-10 01:29 by FChKX8meu1ds
[考博] 27博士申请 +3 白白白ly 2026-11-09 3/150 2026-11-09 21:39 by szq0819
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 15/750 2026-11-09 13:19 by Equinoxhua
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +5 5sJFDKp2kwuI 2026-11-06 16/800 2026-11-09 11:47 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售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,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
信息提示
关闭
请填处理意见
关闭