LittlePupil 发表于 2021-12-1 12:49 谢谢回复,明白了 |
本帖最后由 LittlePupil 于 2021-12-1 13:03 编辑 liuzl 发表于 2021-12-1 08:56 若要计算A、B两个分子的E_BSSE,并且顺带计算校正后的相互作用能,在Gaussian中使用counterpoise=2关键字(可简写为counter=2)。会自动依次计算5个体系,依次输出以下能量: ... 如果你用的是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。 以上摘自前述博文,如果没完全理解,建议再认真读一遍。 这都差到小数点第九位了,已经完全超出量化计算的精度了,不明白你为什么死抓不放这种数值误差。 |
LittlePupil 发表于 2021-11-30 23:59 不对啊,我看了这篇博文,图一第一行应该是BSSE矫正过后的E-AB,还有就是图一倒数两行相加不等于图二第三行,后面小数有偏差! |
|
图一第一行是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.