计算化学公社

标题: 求助高斯计算BSSE中的能量区别 [打印本页]

作者
Author: liuzl 时间: 2021-11-30 09:37
标题: 求助高斯计算BSSE中的能量区别
我用高斯计算了BSSE,但是发现有两部分输出的能量，不懂有什么区别，求教各位大神。图一中的第一行不就是BSSE矫正过的E-AB吗？为何跟图二中的第一行不一致？还有就是图一中的倒数两行相加不应该就是图二中的第三行吗，为什么值也不一致？

作者
Author: LittlePupil 时间: 2021-11-30 23:59
图一第一行是AB基函数下未经BSSE校正的E_AB。
原理上图一中的倒数两行相加就是图二中的第三行，不信你再敲敲计算器。
可以看看社长的这篇博文：《谈谈BSSE校正与Gaussian对它的处理》

作者
Author: liuzl 时间: 2021-12-1 08:56

LittlePupil 发表于 2021-11-30 23:59
图一第一行是AB基函数下未经BSSE校正的E_AB。
原理上图一中的倒数两行相加就是图二中的第三行，不信你再敲 ...

不对啊，我看了这篇博文，图一第一行应该是BSSE矫正过后的E-AB，还有就是图一倒数两行相加不等于图二第三行，后面小数有偏差！

作者
Author: LittlePupil 时间: 2021-12-1 12:49
本帖最后由 LittlePupil 于 2021-12-1 13:03 编辑

liuzl 发表于 2021-12-1 08:56
不对啊，我看了这篇博文，图一第一行应该是BSSE矫正过后的E-AB，还有就是图一倒数两行相加不等于图二第三 ...

若要计算A、B两个分子的E_BSSE，并且顺带计算校正后的相互作用能，在Gaussian中使用counterpoise=2关键字（可简写为counter=2）。会自动依次计算5个体系，依次输出以下能量：
E_AB：A、B基函数下AB复合物的能量

...

如果你用的是G09早期版本或者更早的Gaussian，任务最后会输出"Counterpoise: BSSE energy"，这即是E_BSSE，即(E_A - E_A,bAB) + (E_B - E_B,bAB)。还会看到"Counterpoise: corrected energy"，记为E_corrected，这是消除了因单体基函数重叠造成的虚假能量降低后的AB复合物能量，E_corrected = E_AB + E_BSSE。

以上摘自前述博文，如果没完全理解，建议再认真读一遍。
这都差到小数点第九位了，已经完全超出量化计算的精度了，不明白你为什么死抓不放这种数值误差。

作者
Author: liuzl 时间: 2021-12-1 14:50

LittlePupil 发表于 2021-12-1 12:49
...
以上摘自前述博文，如果没完全理解，建议再认真读一遍。
这都差到小数点第九位了，已经完全超出 ...

谢谢回复，明白了

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