计算化学公社

标题: 利用Gaussian16在G4方法下计算D4D点群得S8分子能量与CCCBDB标准不符 [打印本页]

作者
Author: alohomora100 时间: 2020-9-22 14:45
标题: 利用Gaussian16在G4方法下计算D4D点群得S8分子能量与CCCBDB标准不符

今天利用Gaussian16，在G4水平下计算D4D点群得S8，从输出文件中得到得能量与CCCBDB数据库中给出的标准数值不符，不知原因在哪，还望各位帮助。谢谢。CCCBDB给出的数值为

-3184.644354Hartree，而自己计算得为

-2728.966969 Hartree

作者
Author: ldatea 时间: 2020-9-22 17:33
本帖最后由 ldatea 于 2020-9-22 17:43 编辑

你G4最后一步的HF在SCF部分没有收敛，相差的能量非常多（能相差这么多我也没想到），但是接着往下算了。
你看输出文件里E(CCSD(T))=          -3181.186326  其实到这里已经可以下结论，最终的结果肯定在-3181附近了。
回到CCSD(T)这一步的SCF，发现
SCF Done:  E(RHF) =  -3180.07408979    A.U. after 11 cycles
这里的HF用的是一个小基组，而最后一步的HF用的是非常大的基组，不可能出现（HF/小基组）能量小于（HF/大基组）能量
CoulSu:  requested number of processors reduced to: 5 ShMem 1 Linda.
Density matrix breaks symmetry, PCut= 1.00D-07
Density has only Abelian symmetry.
CoulSu:  requested number of processors reduced to: 5 ShMem 1 Linda.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RHF) =  -2814.06847136    A.U. after  129 cycles
         NFock=128  Conv=0.43D-02    -V/T= 1.9675
SCF density gives NOpUse=  4 NOpAll= 16.

这个能量相差这么多，导致你得到离谱的结果。

http://bbs.keinsci.com/thread-1132-1-1.html里面有详细介绍G4的流程。

很怪的是，SCF不收敛，居然不报错，我看最后一步计算路径里也没有iop(5/13=1)。
好像哪里看到g16把scf不收敛报错的条件给改了

作者
Author: sobereva 时间: 2020-9-22 17:36
应该是当前Gaussian版本的G4的bug。DE(HF)= 455.588076太离谱了

作者
Author: paramecium86 时间: 2020-9-22 18:28
试了试，G16 A03 结果也不太正常（以前帖子就有人提过）。 G16 C01 就跟CCCBDB上数据一致了。

作者
Author: sobereva 时间: 2020-9-22 18:42

paramecium86 发表于 2020-9-22 18:28
试了试，G16 A03 结果也不太正常（以前帖子就有人提过）。 G16 C01 就跟CCCBDB上数据一致了。

A.03不止是不正常，输出信息形式都是错的

作者
Author: sobereva 时间: 2020-9-22 18:43

ldatea 发表于 2020-9-22 17:33
你G4最后一步的HF在SCF部分没有收敛，相差的能量非常多（能相差这么多我也没想到），但是接着往下算了。
...

后来的G16即便SCF没收敛也不会报错终止。这是很讨厌的设定，会坑不少人。

作者
Author: alohomora100 时间: 2020-9-22 19:37

ldatea 发表于 2020-9-22 17:33
你G4最后一步的HF在SCF部分没有收敛，相差的能量非常多（能相差这么多我也没想到），但是接着往下算了。
...

那要把最后一步单独出来计算吗？把G4每一步独立计算得那种？谢谢

作者
Author: ldatea 时间: 2020-9-22 20:06
本帖最后由 ldatea 于 2020-9-22 20:12 编辑

alohomora100 发表于 2020-9-22 19:37
那要把最后一步单独出来计算吗？把G4每一步独立计算得那种？谢谢

可以的，只要算最后一步就可以了，你直接找输出文件里最后一步的关键词，就是总能量不能自动算了。
最好用#p 并且可能要用额外的关键词减少scf不收敛的概率。
你可以检查一下几何优化的结构和网页上提供的G4差别大不大，我看你的核排斥能和它的有点不一样。几何优化如果能对上，一般来说，出问题的概率就比较低了。
不考虑时间的话，全重算一遍也可以。前面的CCSD(T)如果保存了收敛的波幅，后面可以用ccsd=readam减少花的时间，我印象中G4默认保存收敛的波幅。已经优化完成的结构可以用G4=noopt省略第一步几何优化。如果你知道怎么手动把这些能量得到最终的能量，就没必要重算。

作者
Author: alohomora100 时间: 2020-9-22 21:48

ldatea 发表于 2020-9-22 20:06
可以的，只要算最后一步就可以了，你直接找输出文件里最后一步的关键词，就是总能量不能自动算了。
最好 ...

十分感谢

作者
Author: paramecium86 时间: 2020-9-23 02:07
我找来G16 b01 SSE4.2的版本试了试。计算结果看起来是正常的。最后一步HF那步也是正常收敛的。看起来这还不是B.01这个版本的general的bug。

作者
Author: 量化小菜鸡 时间: 2022-1-28 17:14
本帖最后由量化小菜鸡于 2022-1-28 23:03 编辑

sobereva 发表于 2020-9-22 18:43
后来的G16即便SCF没收敛也不会报错终止。这是很讨厌的设定，会坑不少人。

G16 A.03那个最后ARCHIVE对应的G4能量还是对的。

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