Gaussian与orca单点能量差不一致

rookiey

大家好，我在用Gaussian和orca对2个体系进行优化时发现能量差差很多。下面是我的计算过程：
gaussian:
化合物1与化合物2相同操作
优化：
%nprocs=16
%mem=8GB
# b3lyp/def2svp opt freq

0 1
坐标

化合物1得到自由能-3915.686515
化合物1得到自由能-3915.717136
能量差为：2-1=-19.2kJ/mol
随后用优化的结构算单点：                                                                                                                                                      
%nprocs=16
%mem=8GB
# b3lyp/def2tzvpp geom=check guess=read EmpiricalDispersion=GD3BJ


0 1

得到化合物1的电子能量：-3919.23918731
得到化合物2的电子能量：--3919.25874197
能量差为：2-1=-12.3kJ/mol
orca：
对化合物1和2，初始坐标都相同
!rks B3LYP def2-svp def2-svp/J RIJCOSX Grid4 NoFinalGrid D3BJ
!opt freq uno uco
%maxcore 3500
%pal nprocs 16 end
%scf maxiter 1000
     end
%geom maxiter 300
      end
*xyz 0 1
坐标

得到化合物1的电子能量-3914.97954167466
得到化合物2的电子能量-3914.99447139757
能量差为：2-1=-9.6kJ/mol

请问大家这种差异是合理的吗？
问题2:
以化合物1为例，用Gaussian计算时，发现能量高低为：s=1<s=0（rks）<s=0（uks）
用orca计算时，能量高低为：s=1<s=0（uks）<s=0（rks）
请问为什么会出现这种不一致的情况？应该以哪个为准或是应该怎么调节才能一致？

GoldenBaby

我记得社长好像说过orca和Gaussian的B3LYP泛函用的不是一个定义，要使用相同泛函比较的

ymygca

1. orca b3lyp默认使用vwn5, gaussian b3lyp默认使用vwn3。orca 中的b3lyp/g 才是gaussian中的版本
https://chemistry.stackexchange. ... -common-qc-programs
2. orca 使用了density fitting, gaussian 没有
3. 如果你用的gaussian 16,默认grid=ultrafine, 而orca(grid=4)相当于gaussian(grid=fine)

sobereva

此类问题的基本逻辑：比较的时候尽量把因素简化，别同时牵扯一大堆对结果会产生影响的因素，否则很难搞清楚。从最简单的HF/STO-3G单点开始对比，一点一点把影响因素加上去，不断对比，自然什么因素导致的很容易就搞清楚了

rookiey

sobereva 发表于 2019-3-19 23:14
此类问题的基本逻辑：比较的时候尽量把因素简化，别同时牵扯一大堆对结果会产生影响的因素，否则很难搞 ...

好的，感谢sob老师！

rookiey

GoldenBaby 发表于 2019-3-19 23:02
我记得社长好像说过orca和Gaussian的B3LYP泛函用的不是一个定义，要使用相同泛函比较的

对对，我刚刚看到了，想要比较orca要使用B3LYP/G

rookiey

ymygca 发表于 2019-3-19 23:03
1. orca b3lyp默认使用vwn5, gaussian b3lyp默认使用vwn3。orca 中的b3lyp/g 才是gaussian中的版本
https: ...

明白了！非常感谢！