计算化学公社

标题: GD3（BJ）导致吸附能近20千卡的能量差？ [打印本页]

作者
Author: gauss98 时间: 2018-4-19 16:48
标题: GD3（BJ）导致吸附能近20千卡的能量差？
感谢sob大神介绍，我近来计算催化剂体系一直用 B3LYP-D3，最近又由gaussian16改到B3LYP-D3（bj），感觉效果都还不错，
但是这两天计算 MgCl2团簇吸附 EB（苯甲酸乙酯），发现在def2svp基组下用 B3LYP-D3(BJ) 吸附能高达 40kcal/mol，感觉明显偏高了，后来取消GD3，吸附能为21kcal/mol
两者相差达 19kcal/mol

原来觉得也就是两三个千卡的修正，这下导致了近20千卡的能量差，感觉很不正常，都不知道该相信哪个数据了。大神们有什么建议吗？谢谢！！

作者
Author: yflchx 时间: 2018-4-19 19:53
D3BJ参数拟合时直接用的def2-QZVP基组（我印象中拟合参数时由于def2-QZVP够大而没有考虑BSSE）。

如果用def2-SVP基组（非D3BJ参数标配基组），BSSE会很显著，考虑gCP测试一下？

作者
Author: gauss98 时间: 2018-4-20 00:16

yflchx 发表于 2018-4-19 19:53
D3BJ参数拟合时直接用的def2-QZVP基组（我印象中拟合参数时由于def2-QZVP够大而没有考虑BSSE）。

如果用 ...

gcp能在gaussian16里直接用么？

我用了6-311+G**基组试了试，加与不加GD3(BJ)吸附能仍然相差16千卡

作者
Author: sobereva 时间: 2018-4-20 02:51

gauss98 发表于 2018-4-20 00:16
gcp能在gaussian16里直接用么？

我用了6-311+G**基组试了试，加与不加GD3(BJ)吸附能仍然相差16千卡

没法直接用。你做个counterpoise计算来评估就行了。gCP在Grimme网站上也有现成的工具
你也可以试试零阻尼看什么情况
用独立的DFT-D3程序可以把色散校正能分解为片段对的贡献，也可以由此分析一下到底是哪些原子间的D3校正产生那么大的影响，从而评估合理性

对于离子体系，由于DFT-D3只考虑几何结构，没有考虑电子结构，因此加了之后结果未必绝对合理。可以和M06-2X对比一下，或者尝试和其它色散校正方法（依赖于电子结构的）做个对比，比如XDM、vdW-DF。

作者
Author: Jasminer 时间: 2018-4-20 10:38
参阅一下这个？
Statistical Analysis of Semiclassical Dispersion Corrections
Journal of Chemical Theory and Computation
DOI : 10.1021/acs.jctc.8b00078

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