计算化学公社

 找回密码 Forget password
 注册 Register
Views: 15604|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] Gamess LMOEDA的DFT能量分解出现NaN

[复制链接 Copy URL]

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 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
这个问题我已经连续在网上查了很多天,帮助文档也看了不少了,但是一直到找不到解决办法。
烦请各位老师指点迷津~~

十室之邑,必有忠信如丘焉,不如丘之好学也。

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

2#
发表于 Post on 2015-1-24 04:08:16 | 只看该作者 Only view this author
大抵是你的程序有问题,如果用的老版本,尝试最新版本。我这里dfttyp=b3lyp下做LMOEDA运行很正常,输入输出文件都附上了

exam28.inp

937 Bytes, 下载次数 Times of downloads: 96

out.txt

171.8 KB, 下载次数 Times of downloads: 48

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2015-1-24 09:48:48 | 只看该作者 Only view this author
还不能下载附件.....容我提升一下level~~
谢谢sob老师~这个论坛好给力啊,以后就驻扎这里了
十室之邑,必有忠信如丘焉,不如丘之好学也。

6万

帖子

99

威望

5万

eV
积分
120061

管理员

公社社长

4#
发表于 Post on 2015-1-24 10:03:24 | 只看该作者 Only view this author
肯定能下载附件。只有资源分享板块下载附件才需要升到level 2
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2015-1-24 11:16:54 来自手机 | 只看该作者 Only view this author
额(⊙o⊙)… sorry,低级错误 刚刚睡得迷迷糊糊的......手机上没有登陆

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

6#
 楼主 Author| 发表于 Post on 2015-1-24 12:58:48 | 只看该作者 Only view this author
本帖最后由 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位的 晚上试试
十室之邑,必有忠信如丘焉,不如丘之好学也。

239

帖子

0

威望

2340

eV
积分
2579

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2015-1-24 23:16:58 | 只看该作者 Only view this author
成功在Gamess.13.64上运行了~~  果然是32位的gamess自身的问题
十室之邑,必有忠信如丘焉,不如丘之好学也。

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-12 21:01 , Processed in 0.719735 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list