24小时热门版块排行榜    

查看: 991  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见