利用Gaussian16在PBEP86水平下计算所得的能量与ORCA的不匹配的原因

阿甘

各位老师，学生最近打算利用ORCA 5.0对研究体系在revDSD-PBEP86-D4/ma-def2-TZVPP水平下计算体系的电子能的时候，由于自己对ORCA的SCF收敛拿捏不准，于是打算先用Gaussian 16在PBEP86/ma-def2-TZVPP理论水平下计算得到fchk文件，并且通过stable=opt验证过波函数稳定性以后，再利用Multiwfn将fchk转化为mkl文件，再由orca_2mkl转化为gbw文件作为ORCA程序的初猜进行revDSD-PBEP86-D4的计算。但是结果发现，Gaussian16输出的结果为HF=-2775.7207617，而ORCA的SCF迭代时第一次能量却为-2766.8210806312，能量相差很大。不知道是自己哪一步没有做对，烦请各位老师帮忙解惑。
以下是Gaussian16的输入文件命令行：#p nosymm int=nobasistransform gen PBEP86 stable=opt（基组用的是社长发布的给def2系列基组加弥散的ma-TZVPP.txt文件）

以下是ORCA 5.0的输入文件命令行：
! UKS ma-def2-TZVPP def2-TZVPP/c def2/J tightscf RIJCOSX D4
%method
Exchange X_PBE
Correlation C_P86
ScalHFX 0.69
ScalDFX 0.31
ScalGGAC 0.4210
ScalLDAC 0.4210
ScalMP2C 1.0
D3S6 0.5132
D3A1 0.44
D3A2 3.6
D3S8 0.0
end
%mp2
DoSCS true
PS 0.5922
PT 0.0636
end

（以上参数是根据revDSD-PBEP86-D4原文写的）

hebrewsnabla

revDSD-PBEP86的自洽部分并不是PBEP86（这是一个纯泛函），而是一个杂化泛函，你得在高斯里面用iop修改设置

阿甘

hebrewsnabla 发表于 2022-9-16 11:30
revDSD-PBEP86的自洽部分并不是PBEP86（这是一个纯泛函），而是一个杂化泛函，你得在高斯里面用iop修改设置

啊那就是我理解错了revDSD原文的意思了，感谢老师告知。

阿甘

hebrewsnabla 发表于 2022-9-16 11:30
revDSD-PBEP86的自洽部分并不是PBEP86（这是一个纯泛函），而是一个杂化泛函，你得在高斯里面用iop修改设置

方才看到老师您写的python脚本，对学生这个问题的解决很有帮助~再次感谢