利用Gaussian16在G4方法下计算D4D点群得S8分子能量与CCCBDB标准不符

alohomora100

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

-3184.644354Hartree，而自己计算得为

-2728.966969 Hartree

ldatea

你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不收敛报错的条件给改了

sobereva

应该是当前Gaussian版本的G4的bug。DE(HF)=                      455.588076太离谱了

paramecium86

试了试，G16 A03 结果也不太正常（以前帖子就有人提过）。 G16 C01 就跟CCCBDB上数据一致了。

sobereva

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

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

sobereva

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

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

alohomora100

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

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

ldatea

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

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

alohomora100

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

十分感谢

paramecium86

我找来G16 b01 SSE4.2的版本试了试。计算结果看起来是正常的。最后一步HF那步 也是正常收敛的。看起来这还不是B.01这个版本的general的bug。

量化小菜鸡

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

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