Gaussian16里用PBE0关键词计算的实际上是PBE0-DH双杂化泛函

sobereva

这个帖子目的是给Gaussian初学者们提个醒。

Gaussian里，知名的PBE0泛函应当写为PBE1PBE，这点手册里已经写明了。在G09里，直接写PBE0程序会报错，因为不认识
在Gaussian里有一个规则是只要不引起歧义，关键词可以最大程度简化，比如empiricaldispersion=GD3BJ可以简写为em=GD3BJ
从G16开始，支持了一个精度很差的双杂化泛函PBE0-DH。于是乎，根据关键词可简化规则，在G16里面写PBE0的话程序并不会报错，而给你算的是PBE0-DH。倘若你读取的是这个级别对应的电子能量那还好（比如读的是输出文件末尾MP2=后面的值，见《谈谈该从Gaussian输出文件中的什么地方读电子能量》http://sobereva.com/488），虽然精度不怎么样但还算勉强能用，而如果读取的是SCF Done或者末尾的HF=后面的值，那得到的根本什么都不是（实际上对应于PBE0-DH里杂化泛函部分，并不等于PBE0泛函的结果），毫无实际意义，完全不能用来发表。

所以，如果你之前用G16，取的是写了PBE0关键词得到的数据，并且当成是PBE0泛函的结果，需要把所有数据重算；如果数据已经发表了，对数据有责任心的话，应当写个correction。

Gaussian的这个关键词可简化规则在PBE0上很坑爹，Gaussian开发者应该将这种情况作为一个特例考虑，不允许只写PBE0。

为了让初学者更清楚关系，这里总结一下，对于G16（至少截止到B.01）而言：
PBE关键词=报错
PBEPBE关键词=PBE
PBE1PBE关键词=PBE0
PBE0关键词=PBE0-DH

Acceleralater

感谢sob老师，跑了一天看到收敛了但几何优化没做，还在纳闷是为啥，原来算成双杂化泛函了，难怪跑不动。

风起时想你

刚刚犯错，特意在公社检索一下。果然社长未雨绸缪，赞！

luzujia

同楼上，刚从坑里爬出来

lao7

刚刚犯错，特来看看！谢谢！

granvia

G16有时SCF没收敛的计算也给Normal termination，曾被坑过，当时处理数据还以为有重大发现了。还有个bug就是，开壳层体系算NBO，如果输出chk文件，那么输出的只有alpha轨道信息