gCP和DFT在RG18数据集表现极差

惠成功QUB

我最近把GMTKN55中的非键相互作用部分下载下来自己测了一下。
发现如果加上gCP修正在RG18数据集上BLYP和B3LYP就表现极差。Grimme当时用def2-QZVP不到1kcal/mol的误差加上gCP我用def2-TZVPP就10+kcal/mol了

我测试了ORCA4.2.1和ORCA5.0.1，结果几乎一样。
我还测试了B97-3c，毕竟B97-3c自带gCP，B97-3c就没问题。

我还测试了BLYP/def2-SVP GCP(DFT/SVP)，误差还是很大。

我的测试文件在我的github上可以找到
https://github.com/huichenggong/GMTKN55

可能gCP不适合稀有气体（RG18数据集主要是稀有气体）。也可能gCP(DFT/TZ)的bug

量化小菜鸡

很有趣的结果，你有试过D4吗？B97-3c应该是参数拟合的多，所以没这个问题。

惠成功QUB

量化小菜鸡 发表于 2021-10-14 01:12
很有趣的结果，你有试过D4吗？B97-3c应该是参数拟合的多，所以没这个问题。

以上的结果都是D3BJ，我没试过D4。但是我觉得这不应该呀，没理由D3+gCP不行，D4+gCP就行

scf

我的经验是加了gcp几何优化不好收敛了。。

惠成功QUB

scf 发表于 2021-10-14 01:24
我的经验是加了gcp几何优化不好收敛了。。

你这样一说让我不禁怀疑gCP是不是有bug

scf

当时拿orca4测试的，没试过orca5

冰释之川

那你使用def2-TZVPP基组，不加gCP 误差是多少？

wzkchem5

惠成功QUB 发表于 2021-10-13 18:27
你这样一说让我不禁怀疑gCP是不是有bug

可能拟合参数的时候没注意势能面平滑性问题？我记得比较早版本的DFT-D也容易导致几何优化不容易收敛

惠成功QUB

冰释之川 发表于 2021-10-14 08:10
那你使用def2-TZVPP基组，不加gCP 误差是多少？

我的BLYP/def2-TZVPP RG18上的AUE=0.109kcal/mol，图中的紫色柱子
Grimme当时BLYP/def2-QZVP在RG18上是0.06kcal/mol
http://www.thch.uni-bonn.de/tc.o ... 18/BLYP/result.html

冰释之川

惠成功QUB 发表于 2021-10-14 16:39
我的BLYP/def2-TZVPP RG18上的AUE=0.109kcal/mol，图中的紫色柱子
Grimme当时BLYP/def2-QZVP在RG18上是0 ...

那 加了 gCP之后确实有问题，可以反馈给Grimme

惠成功QUB

冰释之川 发表于 2021-10-14 16:57
那 加了 gCP之后确实有问题，可以反馈给Grimme

Grimme回复了邮件，怀疑是参数抄错了。具体是什么还有待检查。等正式结果出来我在来回复

wzkchem5

惠成功QUB 发表于 2021-10-14 21:22
Grimme回复了邮件，怀疑是参数抄错了。具体是什么还有待检查。等正式结果出来我在来回复

我去这都行。。。
不过这样的话，有希望在ORCA 5.0.2就修好，最近ORCA打算推出5.0.2，但是还在等几个比较棘手的bugfix，gCP的bugfix估计能搭上这趟车