计算化学公社

标题: gCP和DFT在RG18数据集表现极差 [打印本页]

作者
Author: 惠成功QUB 时间: 2021-10-13 21:23
标题: gCP和DFT在RG18数据集表现极差
本帖最后由惠成功QUB 于 2021-10-13 21:31 编辑

我最近把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

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

作者
Author: 惠成功QUB 时间: 2021-10-14 01:16

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

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

作者
Author: scf 时间: 2021-10-14 01:24
我的经验是加了gcp几何优化不好收敛了。。

作者
Author: 惠成功QUB 时间: 2021-10-14 01:27

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

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

作者
Author: scf 时间: 2021-10-14 01:28
当时拿orca4测试的，没试过orca5

作者
Author: 冰释之川 时间: 2021-10-14 08:10
那你使用def2-TZVPP基组，不加gCP 误差是多少？

作者
Author: wzkchem5 时间: 2021-10-14 15:20

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

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

作者
Author: 惠成功QUB 时间: 2021-10-14 16:39

冰释之川发表于 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

作者
Author: 冰释之川 时间: 2021-10-14 16:57
本帖最后由冰释之川于 2021-10-14 17:00 编辑

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

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

作者
Author: 惠成功QUB 时间: 2021-10-15 04:22

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

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

作者
Author: wzkchem5 时间: 2021-10-15 14:42

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

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

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