| 查看: 511 | 回复: 0 | ||
[求助]
NWchem 中QM/MM板块optimization问题
|
|
小弟是新手，用的NWchem,QM/MM板块。最近在做一个体系的optimization，pdb文件在附件中。选择结构中的ATP做量化处理。但总是不成功，都一个月很头疼，希望哪位大神可以施以援手，不胜感激！！ pre.nw文件如下： start reac permanent_dir ./perm scratch_dir ./data charge -2 prepare source atp.pdb new_top new_seq new_rst center orient modify segment 952 quantum update lists ignore write reac_ref.rst write reac_ref.pdb end task prepare 优化输入opt.nw start reac permanent_dir ./perm scratch_dir ./data md system reac_ref end basis * library "6-31G" end dft xc b3lyp end qmmm bqzone 20.0 region qm mm maxiter 10 1000 ncycles 2 density espfit xyz foo end task qmmm dft optimize prepare read reac_ref.rst write reac_opt.pdb end task prepare .out输出文件错误 Memory utilization after 1st SCF pass: Heap Space remaining (MW): 0.00 0 Stack Space remaining (MW): 12.99 12989676 convergence iter energy DeltaE RMS-Dens Diis-err time ---------------- ----- ----------------- --------- --------- --------- ------ d= 0,ls=0.0,diis 1 -2625.3970221245 -6.58D+03 1.64D+01 1.78D+01 216.4 Nuclear repulsion energy = 3949.30625746002 Bq nuclear interaction energy = 5.09231130369664 WARNING: error on integrated density = 0.11D-05 greater than required accuracy of 0.10D-05 d= 0,ls=0.5,diis 2 -2282.8510692315 3.43D+02 4.07D+00 9.39D+02 247.5 Nuclear repulsion energy = 3949.30625746002 Bq nuclear interaction energy = 5.09231130369664 d= 0,ls=0.5,diis 3 -2366.4917083419 -8.36D+01 2.05D+00 7.69D+02 278.5 、、、、、、、、、、、、 、、、、、、、、、、、、 、、、、、、、、、、、、 Nuclear repulsion energy = 3949.30625746002 Bq nuclear interaction energy = 5.09231130369664 d= 0,ls=0.5,diis 29 -192.8360685321 -2.07D+00 9.05D-02 7.95D+02 1084.8 Nuclear repulsion energy = 3949.30625746002 Bq nuclear interaction energy = 5.09231130369664 d= 0,ls=0.5,diis 30 -195.5449645017 -2.71D+00 1.25D+00 8.60D+02 1115.7 Nuclear repulsion energy = 3949.30625746002 Bq nuclear interaction energy = 5.09231130369664 Calculation failed to converge ------------------------------ 、、、、、、、、 、、、、、、、、 driver: task_gradient failed 0 ------------------------------------------------------------------------ ------------------------------------------------------------------------ current input line : 27: task qmmm dft optimize ------------------------------------------------------------------------ ------------------------------------------------------------------------ There is an error related to the specified geometry ------------------------------------------------------------------------ For more information see the NWChem manual at http://www.nwchem-sw.org/index.php/NWChem_Documentation For further details see manual section: 0:0:driver: task_gradient failed:: 0 (rank:0 hostname:liaojl pid:5750):ARMCI DASSERT fail. armci.c:ARMCI_Error():260 cond:0 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : reac_ref.pdb
- 附件 2 : atp.pdb
2014-11-10 18:55:02, 40.51 K
2014-11-10 18:57:57, 40.51 K
» 猜你喜欢
求国际会议网站
已经有1人回复
-
求取一些关于纳米材料和纳米技术相关的英文PPT。
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有230人回复
-
【复旦大学】二维材料方向招收2026年博士研究生1名
已经有0人回复
-
北京纳米能源与系统研究所 王中林院士/曹南颖研究员课题组2026级硕/博/博后招生
已经有10人回复
-
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生
已经有19人回复
-
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
