24小时热门版块排行榜    

查看: 1256  |  回复: 4

xianglei1118

新虫 (初入文坛)

[交流] gromacs运行问题 已有3人参与

有没有大神可以指点一下小弟:
       我建立了一个水团,然后给水中加入了NACL离子(在加入氯化钠时,我把top文件中的ions.itp给删了,加入氯化钠后又重新加入了#include “ions.itp”),然后用grompp_mpi运行能量最小化的时候,系统提示了一个问题:
Program grompp_mpi, VERSION 4.6.5
Source code file: /home/dell/gromacs-4.6.5/src/kernel/toppush.c, line: 1336

Fatal error:
Atomtype Zn2+ not found
For more information and tips for troubleshooting, please check the GROMACS
website at http://www.gromacs.org/Documentation/Errors
我自己的想法是:可能是top文件中包含的itp文件不匹配,发现gromacs自带的ions.itp文件里定义的是ZN2+,这个问题不知道怎么样来解决,谢谢前辈们的指点。这是我的top文件内容:#include "ffoplsaa.itp"

;[ defaults ]
; nbfunc        comb-rule        gen-pairs        fudgeLJ        fudgeQQ
;  1                2                  no                1.0        1.0

[atomtypes]
;name     mass      charge   ptype    sigma        epsilon
IW     0             0.000       D   0.0           0.0
OWT4   15.99940      0.000       A   0.31589       0.77490
HW     1.00800       0.000       A   0.00000E+00   0.00000E+00

# include "ions.itp"



[moleculetype]
; name nrexcl
water  1

[atoms]
; nr type resnr residu atom cgnr charge
1     OWT4 1     water  OW1  1     0          15.994
2     HW   1     water  HW2  1     0.5564    1.008
3     HW   1     water  HW3  1     0.5564    1.008
4     IW   1     water  MW4  1    -1.1128    0.0

[constraints]
;i j funct doh  dhh
1       2       1       0.09572
1       3       1       0.09572
2       3       1       0.15139

[exclusions]
1       2       3       4
2       1       3       4
3       1       2       4
4       1       2       3



; The position of the dummy is computed as follows:
;
;                O
;               
;                    D
;          
;        H                H
;
; const = distance (OD) / [ cos (angle(DOH))         * distance (OH) ]
;          0.01546 nm        / [ cos (52.26 deg)        * 0.09572 nm        ]

; Dummy pos x4 = x1 + a*(x2-x1) + b*(x3-x1)


[dummies3]
; Dummy from                        funct        a                b
4       1       2       3       1       0.13193828      0.13193828



[system]
water tip4p/2005

[molecules]
water       2116
NA+              32
CL-              32
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
oplsAA ions.itp里面没发现有Zn的参数。不知道你的力场包是哪里来的。。
2楼2016-11-05 02:53:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awaken2013

铜虫 (小有名气)

而且你的molecules里面没有Zn的数目..
3楼2016-11-05 02:54:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)

不要再问DS了。

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖仅楼主可见
4楼2016-11-07 09:15:43
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

青青白白123

新虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
您好 请问最后怎么解决的?
5楼2018-03-20 21:17:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xianglei1118 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 PYZmcHr5kSaq 2026-11-05 5/250 2026-11-05 21:27 by qL1cU46YIdK0
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-05 20:54 by qL1cU46YIdK0
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-05 20:43 by qL1cU46YIdK0
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 9/450 2026-11-05 20:34 by qL1cU46YIdK0
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 20:23 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-05 20:21 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 11/550 2026-11-05 20:21 by qL1cU46YIdK0
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ttH3f1oA6ace 2026-11-04 8/400 2026-11-05 20:21 by qL1cU46YIdK0
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-05 20:10 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 20:10 by qL1cU46YIdK0
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 kvewwGtKktqB 2026-11-04 8/400 2026-11-05 17:47 by My20PsvIExYO
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 kvewwGtKktqB 2026-11-04 13/650 2026-11-05 17:47 by My20PsvIExYO
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 ncTVBM6AOJHF 2026-11-04 10/500 2026-11-05 17:36 by My20PsvIExYO
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 17:26 by My20PsvIExYO
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 16/800 2026-11-05 17:25 by My20PsvIExYO
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 19/950 2026-11-05 17:25 by My20PsvIExYO
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 ncTVBM6AOJHF 2026-11-04 9/450 2026-11-05 17:25 by My20PsvIExYO
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ncTVBM6AOJHF 2026-11-04 13/650 2026-11-05 17:15 by My20PsvIExYO
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 ttH3f1oA6ace 2026-11-04 9/450 2026-11-05 15:46 by My20PsvIExYO
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见