24小时热门版块排行榜    

查看: 986  |  回复: 2

sdnuzll

金虫 (初入文坛)

[求助] nwchem:task scf optimize 已有1人参与

各位大虾,小妹最近刚开始学习nwchem,运行了examples里给的qmd中h2o-scf.nw改task为task scf optimize,出现了下面的错误
grad: scratch alloc failed 48333597
current input line :
    41: task scf optimize
------------------------------------------------------------------------
------------------------------------------------------------------------
------------------------------------------------------------------------
For more information see the NWChem manual at http://www.nwchem-sw.org/index.php/NWChem_Documentation


For further details see manual section:                                                                                                                                          
Last System Error Message from Task 0:: Inappropriate ioctl for device
0:0:grad: scratch alloc failed:: 48333597
(rank:0 hostname:node1 pid:14826):ARMCI DASSERT fail. armci.c:ARMCI_Error():260 cond:0
  0: ARMCI aborting 0 (0).
  0: ARMCI aborting 0 (0).
system error message: Invalid argument
有没有人知道是什么错误,怎么改
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

carbonNT

新虫 (初入文坛)

【答案】应助回帖

★ ★
sdnuzll: 金币+2 2014-08-25 12:47:10
找到答案了吗?如果没有,请提供*.nw文件。我可看一下。
2楼2014-08-16 11:22:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sdnuzll

金虫 (初入文坛)

引用回帖:
2楼: Originally posted by carbonNT at 2014-08-16 09:22:34
找到答案了吗?如果没有,请提供*.nw文件。我可看一下。

使用nwchem困难多多,能帮我看下这个问题吗?
用qm/mm做动力学,输出文件由开始的

    atom               coordinates                        gradient
                 x          y          z           x          y          z
   1 Li      0.000000   0.000000  18.636124    0.005466   0.000949  -0.007030
   2 O       2.342822  -1.771483 -15.354033    0.011121  -0.009646   0.005947
   3 H       1.106456  -0.444499 -14.740569   -0.013872   0.022380   0.002479
   4 H       1.952188  -3.312149 -14.464895   -0.002526  -0.004838  -0.006080
   5 O      -0.612316   1.825524 -14.037520    0.003461  -0.014901  -0.006541
   6 H      -0.507252   2.727578 -12.408485    0.003863   0.007021   0.010092
   7 H      -2.345699   1.933562 -14.591228    0.000798  -0.001067  -0.000851
   8 O      -0.164542   4.313244  -9.592731   -0.008626  -0.003181  -0.005363
   9 H      -0.810780   3.755016  -7.981168    0.004446  -0.002104   0.000130
  10 H       0.000000   6.125855  -9.465091    0.005700   0.002564   0.000586

变成
atom               coordinates                        gradient
                 x          y          z           x          y          z
   1 Li     -1.122517  -0.472775 139.919967   -0.000004  -0.000032   0.000002
   2 O      16.759796 -13.363799 **********    0.004713  -0.003958   0.012128
   3 H      16.042379 -11.805762 **********   -0.004820   0.008183  -0.007343
   4 H      17.139576 -13.055927 **********    0.000107  -0.004228  -0.004784
   5 O      -5.176383  14.046525 **********   -0.015292  -0.001156  -0.007835
   6 H      -5.130719  15.343895 **********    0.002457   0.003242   0.004273
   7 H      -6.835613  14.076591 **********    0.012833  -0.002084   0.003563
   8 O      -1.471080  32.603361 -70.743400   -0.003689   0.001025   0.000028
   9 H      -3.132862  32.317488 -70.059433    0.002053  -0.001000   0.000386
  10 H      -1.656437  33.959569 -71.946970    0.001637  -0.000025  -0.000414

nw"如下:
start Li-3w
  memory 1 gb
  permanent_dir /tmp/perm
  scratch_dir /tmp/data
  print medium
  ECHO

  charge 1
  prepare
  source Li-3w.pdb
  new_top new_seq
  new_rst
  modify segment 1  quantum
  modify segment 2  quantum
  modify segment 3  quantum
  modify segment 4  quantum
  center
  orient
  solvate box 3.0
  update lists
  ignore
  write Li-3w_qmmm.pdb
  write Li-3w_qmmm.rst
  end
  task prepare

  md
   system Li-3w_qmmm
   step 0.001 equil 20000 data 10000
   cutoff 0.9
   isotherm
   isobar
   print step 10 stat 100
   print step 10 stat 100

update pairs 10
   record rest 100
   record prop 100 sync 10 times 0
   record coord 100
   profile
   load pairs 10 size 0.75
   test 10
  end

  basis
   Li library aug-cc-pvdz
   O  library aug-cc-pvdz
   H  library aug-cc-pvdz
   end
   dft
   xc b3lyp
   mult 1
   direct
   iterations 600
   print medium
  end

  qmmm
   bq_zone  15.0
   region    qm solvent
   maxiter   25    300
   ncycles   50
   density   espfit
   xyz foo
   link_atoms hydrogen
  end

  task qmmm dft energy
  task qmmm dft dynamics

初始结构和基组都改过,但不起作用。到底是什么原因?求指导!!!
3楼2014-09-02 16:44:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sdnuzll 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] filecode +13 documentary 2026-08-10 15/750 2026-08-12 15:56 by 云上清扬
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 35/1750 2026-08-12 14:38 by sdfapple719
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
信息提示
请填处理意见