Gaussian16计算求助，按照文献一模一样方法基组却得不到文献结果，不知道哪儿算错了

baicao

各位老师好，好不容易解决了电脑算CCSD的问题然后按文章一模一样的方法(CCSD(T)/cc-pVDZ//B3LYP-D3/def2-TZVP) 计算单线态/三线态卡宾能量差结果又出了问题。附件文献中Table3中DPC的Es/t能量差是3.3 kcal/mol。然后我比较我计算结果CCSD(T)的两种卡宾单点能差是1.35 kcal/mol, 如果加上B3LYP-D3/def2-TZVP零点能矫正值后是1.96 kcal/mol。 怎么也算不出文献的3.3 kcal/mol。文章还特地说明了没有用矫正因子。各位老师能帮忙看看哪里出了问题吗？
有个输出文件上传不了就截个图了

liyuanhe211

看到“CCSD(T)/cc-pVDZ//B3LYP-D3/def2-TZVP”这种计算级别的选择，这文章的计算结果就别看了。（当然，这不代表它的结果是无法重复出来的。）

naonao5205

是G（T） 不是U（0） 我用手機加出來是3.01kcal/mol 当然 这个级别搭配很扯

sobereva

作者没有丝毫的基组选用常识。作者用DLPNO-CCSD(T)居然还用这么点基组，太荒诞了。

仔细看
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

baicao

sobereva 发表于 2020-9-30 19:00
作者没有丝毫的基组选用常识。作者用DLPNO-CCSD(T)居然还用这么点基组，太荒诞了。

仔细看

谢谢老师，我当时也很奇怪文章基组的选用，虽然他的结果和实验还挺吻合的，我会尝试换个更高精度的基组计算的。老师我还有个问题文章里面的单线态/三线态卡宾deltaE是什么，是CCSD(T)计算结果下的单点能差吗？以前主要算一些简单的势能面，没有用到过这么高精度的计算。然后有一篇类似文献（基组选用也不靠谱），文章倒数第二页这个能量的计算看不太懂，这个是不是我没有得到与文献相符结果的原因？

baicao

liyuanhe211 发表于 2020-9-30 02:27
看到“CCSD(T)/cc-pVDZ//B3LYP-D3/def2-TZVP”这种计算级别的选择，这文章的计算结果就别看了。（当然，这 ...

谢谢老师， 我看了sob老师的帖子，确实不太靠谱，我会用高级别的基组算的，只是重复文献的计算不知道怎么计算文章报道的这个能量差，是仅仅CCSD（T）下的单点能差吗？

wzkchem5

baicao 发表于 2020-9-30 20:49
谢谢老师，我当时也很奇怪文章基组的选用，虽然他的结果和实验还挺吻合的，我会尝试换个更高精度的基组计 ...

这个应该说的是用小体系计算出B3LYP对于这种问题相对于CCSD(T)的平均误差，然后拿这个值来校正大体系的B3LYP计算值，以获得接近CCSD(T)的精度。
这个应该和你的计算无关，因为如果用了这种方法，文章必然会说，不会默认使用。但是有一点是值得注意的，就是算能量差的时候加了B3LYP水平下的ZPE校正

baicao

wzkchem5 发表于 2020-9-30 22:05
这个应该说的是用小体系计算出B3LYP对于这种问题相对于CCSD(T)的平均误差，然后拿这个值来校正大体系的B3 ...

十分感谢，我算上ZPE矫正时得到的结果是1.96 kcal/mol,与文章的3.3 kcal/mol还是有出入。但是我一眼看到单点能的差是1.35kcal/mol，然后这个加上ZPE矫正的差值是1.96 kcal/mol，两者相加刚好3.3 kcal/mol. 但觉得这没有任何依据。