24小时热门版块排行榜    

查看: 1185  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天务委会开完了,明天出结果吗 +18 angus9576 2026-08-25 22/1100 2026-08-26 00:41 by merchancy
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 7K1CJE38xLG4 2026-08-25 3/150 2026-08-26 00:40 by cNXvBfCpiZOM
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +16 木水思豆 2026-08-25 19/950 2026-08-25 23:56 by dragonxp
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 2026国自然函评费到账 +20 羊腰板 2026-08-21 23/1150 2026-08-25 16:34 by zsna
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗?20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见