计算化学公社

标题: GAMESS做能量分解运行成功，但计算结果反常 [打印本页]

作者
Author: Finesnow 时间: 2018-9-4 20:22
标题: GAMESS做能量分解运行成功，但计算结果反常
各位老师好，小的想把高斯优化后的结构用GAMESS进行能量分解，跌跌撞撞的能运行成功了，但是计算结果出错。初次接触计算的我自己纠结了半天不知道问题出在了哪里。再次诚挚的恳请各位前辈指点一二。

作者
Author: 我本是个娃娃 时间: 2018-9-4 20:34
什么能量分解方法？如果是对于弱相互作用的能量分解，SAPT更严谨

作者
Author: zjxitcc 时间: 2018-9-4 20:37
你之前提问过，当时我以为你写timlim=1是为了测试，现在还这么写

不知道从哪里学的

看ITER下的SCF迭代信息，问题在那里。timlim=1这句话意思是运行时间极限为1 min，你这体系在1分钟内SCF显然算不完收敛不了，1min到自动停掉了。
还有第一行你原来的DFT关键词也没了，生成的fort.7文件要自己改的，不是复制过来直接用。

作者
Author: Finesnow 时间: 2018-9-5 09:32

zjxitcc 发表于 2018-9-4 20:37
你之前提问过，当时我以为你写timlim=1是为了测试，现在还这么写不知道从哪里学的
看ITER下的SCF ...

啊哈，是这个样子啊，我再好好改改

让前辈见笑啦

又是您及时的指点我这个智障，再次谢谢您耐心的回答

作者
Author: Finesnow 时间: 2018-9-5 09:44

我本是个娃娃发表于 2018-9-4 20:34
什么能量分解方法？如果是对于弱相互作用的能量分解，SAPT更严谨

首先谢谢您热心的回复

然后弱弱的我坦白的说......我并不知道SAPT是什么方法，我这就去找找资料

作者
Author: Finesnow 时间: 2018-9-5 12:23

zjxitcc 发表于 2018-9-4 20:37
你之前提问过，当时我以为你写timlim=1是为了测试，现在还这么写不知道从哪里学的
看ITER下的SCF ...

前辈您好，我调整了timlim以后总能量仍然为零，REP和POL项还是显示NaN。另外加上DFT关键词以后计算会以“GRID DFT PROGRAM CAN'T DO UNSUPPORTED DFTTYP=M062X”错误退出。或许您能再指点一下吗？

作者
Author: sobereva 时间: 2018-9-5 12:39

Finesnow 发表于 2018-9-5 09:44
首先谢谢您热心的回复
然后弱弱的我坦白的说......我并不知道SAPT是什么方法，我这就去找找资料

这里第10节有简要说明
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

论坛首页google搜索框里搜SAPT也能找到一大票资料

作者
Author: Finesnow 时间: 2018-9-5 13:22

sobereva 发表于 2018-9-5 12:39
这里第10节有简要说明
Multiwfn支持的弱相互作用的分析方法概览
http://sobereva.com/252

好的好的，十分谢谢热心的您

作者
Author: zjxitcc 时间: 2018-9-5 14:07
完全是没认真看手册就开始瞎算。至少有5处错误。
DIRSCF一些自己不懂的关键词不要写；guess=huckel这行也不用写；mwords=200差不多才1.6GB，机子内存允许的话，这里稍微加大一些，算得会快一点；默认SCF圈数只有30圈，算大体系一般是不够用的，要加大一些，比如120.能量分解输出文件里有许多次ITER信息，包括总体系、每个单体在单体基组下、单体在总基组下等等，任何一个SCF不收敛都不行，不能只看一个。M062X在GAMESS里写作M06-2X，没看手册就臆想关键词。
GAMESS有个缺陷，就是SCF收敛性做得不好，经常导致能量分解过程中某一个SCF不收敛。
GAMESS是有PDF手册的，要自己看的啊，起码先把你写的$这几个下属关键词意思全看过一遍吧，你这MMULT值写的也是错的。还有目录test/eda下有很多例子，先自己练习。

作者
Author: Finesnow 时间: 2018-9-5 14:52

zjxitcc 发表于 2018-9-5 14:07
完全是没认真看手册就开始瞎算。至少有5处错误。
DIRSCF一些自己不懂的关键词不要写；guess=huckel这行也 ...

恩恩，接受批评，确实是我急躁了，还不会走就想跑，就光想试着计算了。看来我这幼稚园水平是急不得的，还是得静下心好好看手册

这次又暴露了我自己的无知，下次提问时我争取能够进步一些

才发现我每次发的帖子都有您不离不弃在指点和帮助，哈哈哈，真心的再次谢谢您

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