计算化学公社

标题: gamess 内存不足 [打印本页]

作者
Author: mozilla88 时间: 2016-1-16 11:03
标题: gamess 内存不足
本人刚刚接触GAMESS计算，遇到一个关于内存不足的问题，请教高手。谢谢先。
我计算的是一个小的团簇（共15个原子：13个Au,2个Cd），用的是网上找的LANL2DZ ECP势。按说这么小的体系，不应该需要太多的内存，但是计算运却给出下面的错误：
***** ERROR: MEMORY REQUEST EXCEEDS AVAILABLE MEMORY
PROCESS NO. 0 WORDS REQUIRED= 1066731 AVAILABLE= 1000000

我试着在$SYSTEM中加入MEMDDI=2000000，程序又提示我，没有权利改变机器的内存。
GAMESS 版本：5 DEC 2014 (R1)
希望得到高手的指点。
下面是我的inp文件：

$CONTRL SCFTYP=UHF RUNTYP=OPTIMIZE PP=READ
ICHARG= 1 MULT=1 DFTTYP=B3LYP MAXIT=200 $END
$SYSTEM TIMLIM=60 $END
$STATPT NSTEP=500 OPTTOL=1.0d-6  $END
$GUESS  GUESS=HUCKEL $END
$BASIS basnam(1)=aubs,aubs,aubs,aubs,aubs,aubs,aubs,aubs,aubs,aubs,aubs,aubs,
               aubs,cdbs,cdbs $end
$DATA
DFT geometry optimization
C1
Au  79.0 0.000000 0.000000 0.000000
Au  79.0 1.419860  -2.294120 0.514980
Au  79.0  -1.419850 2.294120  -0.514980
Au  79.0  -0.095180  -1.952200  -1.921100
Au  79.0 0.095180 1.952200 1.921090
Au  79.0 2.282270  -0.345800  -1.450620
Au  79.0  -2.282260 0.345800 1.450610
Au  79.0  -2.383420  -0.309530  -1.450210
Au  79.0 2.383430 0.309540 1.450210
Au  79.0  -0.051710 0.906160  -2.627790
Au  79.0 0.051710  -0.906150 2.627780
Au  79.0  -1.531320  -2.257770 0.570180
Au  79.0 1.531320 2.257770  -0.570190
Cd  48.0 2.759040  -3.699830  -2.328810
Cd  48.0  -2.759030 3.699830 2.328800
$END
$aubs
S 3
  1    2.8090000          -1.2021556
  2    1.5950000             1.6741578
  3    0.5327000             0.3526593
S 4
  1    2.8090000             1.1608481
  2    1.5950000          -1.8642846
  3    0.5327000          -1.0356230
  4    0.2826000             1.3064399
S 1
  1    0.0598000             1.0000000
P 3
  1    3.6840000          -0.2802681
  2    1.6660000             0.7818398
  3    0.5989000             0.4804776
P 2
  1    0.6838000          -0.0952078
  2    0.0977000             1.0299147
P 1
  1    0.0279000             1.0000000
D 2
  1    1.2870000             0.5844273
  2    0.4335000             0.5298161
D 1
  1    0.1396000             1.0000000

$end
$cdbs
S 2
  1    0.5095000          -1.2713002
  2    0.1924000             1.8002112
S 1
  1    0.0544000             1.0000000
P 2
  1    0.8270000          -0.1083020
  2    0.1287000             1.0367049
P 1
  1    0.0405000             1.0000000
D 3
  1    5.1480000             0.0703071
  2    1.9660000             0.5138427
  3    0.7360000             0.5416758
D 1
  1    0.2479000             1.0000000

$end
$ECP
  AU-ECP GEN    60    4
  5    ----- g-ul potential    -----
-60.0000000       1 622.6287956
-555.5292312       2 136.2843607
-168.0019785       2    33.1549781
-63.0399875       2    9.9894895
   -4.2516681       2    3.0481312
  6    ----- s-ul potential    -----
   3.0000000       0 194.7374304
   38.6020880       1 351.5327447
864.8370727       2 122.3270402
374.9935520       2    32.0914617
289.7910100       2    5.2451812
-152.4532773       2    4.4916223
  4    ----- p-ul potential    -----
   2.0000000       0 420.6158801
   73.8885625       1 109.4417815
326.6729872       2    34.1714280
126.5814591       2    5.9879750
  5    ----- d-ul potential    -----
   3.0000000       0 219.2666158
   55.6793149       1 122.7297786
449.1987335       2    63.1063369
215.0269091       2    18.3684520
   64.0840995       2    4.4972844
  5    ----- f-ul potential    -----
   4.0000000       0 108.5506037
   51.8065335       1    56.4795527
231.2183113       2    29.2069159
119.0047386       2    9.5440543
   15.3424188       2    2.8965118
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  AU-ECP
  CD-ECP GEN    36    3
  7    ----- f-ul potential    -----
   -0.0613214       0 350.0762490
-24.9355593       1    81.1757778
-87.7728765       2    24.7747818
-19.0665030       2    12.5760330
-16.1901254       2    4.9889712
   -3.9441254       2    1.8449910
   -0.3097476       2    0.5309494
  5    ----- s-ul potential    -----
   2.9050102       0    54.5786193
   39.6289198       1    16.0082316
   57.4435790       2    5.2767462
   90.2351166       2    0.9731070
-71.7139236       2    0.9305646
  4    ----- p-ul potential    -----
   1.9585463       0    41.3581766
   37.5669844       1    12.1214982
   44.0875104       2    3.4478806
   12.5678599       2    0.9362175
  6    ----- d-ul potential    -----
   2.9753499       0    85.3710499
   24.7054373       1    73.0500781
218.2163143       2    28.4262887
118.9945370       2    3.9450158
-82.4999671       2    3.5349379
   -0.2413758       2    0.6482194
  CD-ECP
$END

作者
Author: root 时间: 2016-1-16 13:47
默认只有1MWords即8兆内存。
MWORDS和MEMDDI的单位是1MW，你如果写2000000就是声明16TB的内存。。。

作者
Author: sobereva 时间: 2016-1-16 22:34
DFT计算不需要设MEMDDI，对后HF计算才需要
你应当设的是$SYSTEM里的MWORDS，MWORDS=x代表在每个节点上用x*8MB内存

PS：这种问题用Gaussian省事得多速度也更快

作者
Author: mozilla88 时间: 2016-1-18 12:04
非常感谢两位的解释，按sobereva的建议，在$SYSTEM中使用了MWORDS解决了问题。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)