计算化学公社

标题: Gaussian-ORCA能量对比疑惑 [打印本页]

作者
Author: 柒月小鱼 时间: 2018-6-26 15:47
标题: Gaussian-ORCA能量对比疑惑
本帖最后由柒月小鱼于 2018-6-26 15:49 编辑

各位老师好！
我在计算一个绝热能的时候（T1-S0）开始用的Gaussian优化计算频率，而后用ORCA计算单点，发现二者差距达到5kcal，
而且Gaussian的计算结果明显更与实验值接近，是我的ORCA输入有问题么（我计算的是金属配合物）
还请各位老师指教一二
ORCA关键词
!RI-PWPB95 RIJK D3BJ def2-TZVP def2/JK Grid5 tightSCF nopop
%basis NewECP Ir "SDD" end
Gaussian关键词
#p freq nosymm b3lyp/genecp scf=(maxcyc=400,conver=8) EmpiricalDispersion=GD3BJ
同一原子坐标计算

作者
Author: Accelerator 时间: 2018-6-26 18:15
两个软件用的泛函和基组都不一样当然没有可比性了

作者
Author: 柒月小鱼 时间: 2018-6-26 19:30

Accelerator 发表于 2018-6-26 18:15
两个软件用的泛函和基组都不一样当然没有可比性了

是的
可是我纠结的是1.低级别的计算数据与实验更相近
2.在想我的ORCA是不是还是有些地方没有弄好
我重复文献的计算，用的是我Excel第三栏的基组，得到的结果与我高斯计算结果几乎一致

作者
Author: Graphite 时间: 2018-6-26 22:56

柒月小鱼发表于 2018-6-26 19:30
是的
可是我纠结的是1.低级别的计算数据与实验更相近
2.在想我的ORCA是不是还是有些地方没有弄好

泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应该固定泛函

作者
Author: sobereva 时间: 2018-6-26 23:06
问这个问题说明你对量化计算的基本思想还很糊涂

要么对比程序实现的差异，要么对比理论方法的差异，要么对比基组的差异，要么考察影响数值精度的设定的差异，要么对比理论方法对实验结果重现精度的差异...比其中一个就应当令其它因素保持不变。
你当前问题把乱七八糟因素全搅合在一起，根本什么都没法讨论

作者
Author: 柒月小鱼 时间: 2018-6-27 09:08

Graphite 发表于 2018-6-26 22:56
泛函之间没有固定的谁高谁低, 要看是否适合体系
比较泛函的适用性时应该固定基组, 比较基组的精度时也应 ...

是的，泛函之间没有固定的高低短长，这也是入门就该了解的
可以最近看到各种说双杂化来做单点计算的效果很好，所以想来试试
结果高级别的泛函基组与实验相差较远（这是做方法的一个分子，实验结果应该是没有问题的，多次被验证了）
当然文献中的能很好的重复实验结果，怀疑我的ORCA使用ing

作者
Author: 柒月小鱼 时间: 2018-6-27 09:12

sobereva 发表于 2018-6-26 23:06
问这个问题说明你对量化计算的基本思想还很糊涂

要么对比程序实现的差异，要么对比理论方法的差异，要么 ...

我是这样考虑的，
做高级别的计算，按理是比低级别的计算更理想，（尤其是实验已经非常成功的东西）
我没有去做程序对比也就是这个原因，考虑PWPB95是目前的金方法
\(^o^)/~

作者
Author: sobereva 时间: 2018-6-27 14:19

柒月小鱼发表于 2018-6-27 09:12
我是这样考虑的，
做高级别的计算，按理是比低级别的计算更理想，（尤其是实验已经非常成功的东西）
我 ...

那你就没必要去和Gaussian去对比，都用ORCA就完了
根本没必要用SDD，def2-TZVP本身对Ir就有定义，何况SDD本身比def2-TZVP还小，而过渡金属在整个体系里又占有最重要地位，刻意用SDD明显有害无益

倘若计算流程都合理，数据有可比性，B3LYP/6-31G**真的结果比PWPB95-D3BJ/def2-TZVP还好，那9成可能也只是对当前体系的巧合。

作者
Author: 柒月小鱼 时间: 2018-6-29 09:21

sobereva 发表于 2018-6-27 14:19
那你就没必要去和Gaussian去对比，都用ORCA就完了
根本没必要用SDD，def2-TZVP本身对Ir就有定义，何况 ...

老师早上好！
我现在将我的输出文件打包传上来了，希望您可以抽空看看，嘿嘿
其中ORCA使用的nohup命令，只有一个.out。第一个单点是单重态的，第二个单点是三重态的
麻烦老师了

作者
Author: 柒月小鱼 时间: 2018-7-3 16:37

sobereva 发表于 2018-6-27 14:19
那你就没必要去和Gaussian去对比，都用ORCA就完了
根本没必要用SDD，def2-TZVP本身对Ir就有定义，何况 ...

老师我在def2-QZVPP级别下计算了单点 PWPB95-D3BJ和B2PLYP-D3BJ
两者对比试验的效果确实不如用ORCA/B3LYP-D3BJ/def2-TZVP
那么我的单点计算是否使用高斯在B3LYP-D3BJ/def2-TZVP计算比较好，免得别人问在同种泛函还要换软件，虽然ORCA真的很快
测试了几天
还请老师指教

作者
Author: sobereva 时间: 2018-7-4 00:59

柒月小鱼发表于 2018-7-3 16:37
老师我在def2-QZVPP级别下计算了单点 PWPB95-D3BJ和B2PLYP-D3BJ
两者对比试验的效果确实不如用ORCA/B3LY ...

既然结果是这样，就用这个吧

作者
Author: 柒月小鱼 时间: 2018-7-4 09:42

sobereva 发表于 2018-7-4 00:59
既然结果是这样，就用这个吧

好的谢谢老师耐心指导

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