计算化学公社

标题: Gamess LMOEDA的DFT能量分解出现NaN [打印本页]

作者
Author:
brothers    时间: 2015-1-23 21:57
标题: Gamess LMOEDA的DFT能量分解出现NaN
本帖最后由 brothers 于 2015-1-24 01:40 编辑

这个问题在小木虫已经发过了,没人解答;把问题发给LMOEDA的苏培峰老师咨询,同样没有回复。烦请高人出现啊!!!)

在gamess提供的能量分解例题exam28.inp中,添加mplevl=2可以进行mp2的能量分解,但是用dfttyp=b3lyp就会有很多NaN。所以问题就是:Gamess是不是需要其他参数来控制在DFT水平的能量分解呢?
switch、soscf,swoff等关键字我也试了一下,结果都一样,NaN一大片,计算一直挂在那里不死不活的。

附上最开始出现NaN的部分

        DENSITY MATRIX CONVERGENCE THRESHOLD=  1.00E-05
     COARSE -> FINE DFT GRID SWITCH THRESHOLD=  3.00E-04 (SWITCH IN $DFT)
                   HF -> DFT SWITCH THRESHOLD=  5.00E-03 (SWOFF IN $DFT)
     MEMORY REQUIRED FOR RHF ITERS=    303447 WORDS.
DFT IS SWITCHED OFF, PERFORMING PURE SCF UNTIL SWOFF THRESHOLD IS REACHED.

ITER EX DEM     TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0  0      -56.0010149252   -56.0010149252   0.201191300   0.471902026
   2  1  0      -56.1766838888    -0.1756689637   0.086685749   0.045577797
   3  2  0      -56.1887253471    -0.0120414583   0.023893084   0.015751971
   4  3  0      -56.1902678935    -0.0015425464   0.016210481   0.009238234
   5  4  0      -56.1904867841    -0.0002188906   0.003782533   0.002451094
CONVERGED TO SWOFF, SO DFT CALCULATION IS NOW SWITCHED ON.
          * * *   INITIATING DIIS PROCEDURE   * * *
   6  5  0                 NaN              NaN   0.000000000           NaN
   7  6  0                 NaN              NaN   0.000000000           NaN
这个问题我已经连续在网上查了很多天,帮助文档也看了不少了,但是一直到找不到解决办法。
烦请各位老师指点迷津~~


作者
Author:
sobereva    时间: 2015-1-24 04:08
大抵是你的程序有问题,如果用的老版本,尝试最新版本。我这里dfttyp=b3lyp下做LMOEDA运行很正常,输入输出文件都附上了
作者
Author:
brothers    时间: 2015-1-24 09:48
还不能下载附件.....容我提升一下level~~
谢谢sob老师~这个论坛好给力啊,以后就驻扎这里了
作者
Author:
sobereva    时间: 2015-1-24 10:03
肯定能下载附件。只有资源分享板块下载附件才需要升到level 2
作者
Author:
brothers    时间: 2015-1-24 11:16
额(⊙o⊙)… sorry,低级错误 刚刚睡得迷迷糊糊的......手机上没有登陆
作者
Author:
brothers    时间: 2015-1-24 12:58
本帖最后由 brothers 于 2015-1-24 14:24 编辑
sobereva 发表于 2015-1-24 04:08
大抵是你的程序有问题,如果用的老版本,尝试最新版本。我这里dfttyp=b3lyp下做LMOEDA运行很正常,输入输出 ...

跑来实验室,用sob老师给的inp文件计算,同样出现很多NaN。看来真是gamess的问题。我用的是windows版本的。
被打败了.......
请问为什么只有DFT方法不奏效呢?用我的gamess,使用mp和ccsd方法均可以成功进行LMOEDA的能量分解~~。是编译的问题还是其他问题?

下载了一个64位的 晚上试试

作者
Author:
brothers    时间: 2015-1-24 23:16
成功在Gamess.13.64上运行了~~  果然是32位的gamess自身的问题




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3