计算化学公社

 找回密码 Forget password
 注册 Register

求助高斯计算BSSE中的能量区别

查看数: 4065 | 评论数: 4 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-11-30 09:37

正文摘要:

我用高斯计算了BSSE,但是发现有两部分输出的能量,不懂有什么区别,求教各位大神。图一中的第一行不就是BSSE矫正过的E-AB吗?为何跟图二中的第一行不一致?还有就是图一中的倒数两行相加不应该就是图二中的第三行吗 ...

回复 Reply

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

谢谢回复,明白了
LittlePupil 发表于 Post on 2021-12-1 12:49:41
本帖最后由 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。

以上摘自前述博文,如果没完全理解,建议再认真读一遍。
这都差到小数点第九位了,已经完全超出量化计算的精度了,不明白你为什么死抓不放这种数值误差。
liuzl 发表于 Post on 2021-12-1 08:56:45
LittlePupil 发表于 2021-11-30 23:59
图一第一行是AB基函数下未经BSSE校正的E_AB。
原理上图一中的倒数两行相加就是图二中的第三行,不信你再敲 ...

不对啊,我看了这篇博文,图一第一行应该是BSSE矫正过后的E-AB,还有就是图一倒数两行相加不等于图二第三行,后面小数有偏差!
LittlePupil 发表于 Post on 2021-11-30 23:59:31
图一第一行是AB基函数下未经BSSE校正的E_AB。
原理上图一中的倒数两行相加就是图二中的第三行,不信你再敲敲计算器。
可以看看社长的这篇博文:《谈谈BSSE校正与Gaussian对它的处理

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-25 03:54 , Processed in 0.178098 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list