24小时热门版块排行榜    

查看: 984  |  回复: 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能看出啥? +7 要乐观耀哥 2026-08-10 19/950 2026-08-12 10:55 by 要乐观耀哥
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 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-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 化学口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 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见