GAMESS做能量分解运行成功，但计算结果反常

Finesnow

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

我本是个娃娃

什么能量分解方法？如果是对于弱相互作用的能量分解，SAPT更严谨

zjxitcc

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

Finesnow

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

啊哈，是这个样子啊，我再好好改改让前辈见笑啦
又是您及时的指点我这个智障，再次谢谢您耐心的回答

Finesnow

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

首先谢谢您热心的回复
然后弱弱的我坦白的说......我并不知道SAPT是什么方法，我这就去找找资料

Finesnow

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

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

sobereva

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

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

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

Finesnow

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

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

zjxitcc

完全是没认真看手册就开始瞎算。至少有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下有很多例子，先自己练习。

Finesnow

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

恩恩，接受批评，确实是我急躁了，还不会走就想跑，就光想试着计算了。看来我这幼稚园水平是急不得的，还是得静下心好好看手册这次又暴露了我自己的无知，下次提问时我争取能够进步一些
才发现我每次发的帖子都有您不离不弃在指点和帮助，哈哈哈，真心的再次谢谢您