24小时热门版块排行榜    

查看: 883  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见