计算化学公社

标题: 计算BSSE数据结果与高斯输出不一致 [打印本页]

作者
Author: troll 时间: 2019-2-26 15:27
标题: 计算BSSE数据结果与高斯输出不一致
本帖最后由 troll 于 2019-2-26 15:29 编辑

各位老师好，我使用Gaussian计算苯分子和Ar之间的相互作用能，发现输出文件中的BSSE enrergy与我手算的结果不一致，请各位帮忙看看哪里有问题。计算设置为 ccsd(t)/cc-pvtz counterpoise=2 ，软件版本为G09 C。输出文件的结果为

Counterpoise: BSSE energy = 0.001106516798，

复制代码

然而如果按照sobereva老师的文章[谈谈BSSE校正与Gaussian对它的处理](http://sobereva.com/46)中公式进行计算，

E_BSSE = (E_A - E_A,bAB) + (E_B - E_B,bAB)

我提取了输出文件中的下列五行结果，分别为E_AB, E_A,bAB,  E_B,bAB, E_A, E_B。根据上面公式计算得到的E_BSSE为 0.000138676，与高斯输出的BSSE能量不一致。

SCF Done:  E(RHF) =  -757.575422260    A.U. after 12 cycles
SCF Done:  E(RHF) =  -230.778537839    A.U. after 9 cycles
SCF Done:  E(RHF) =  -526.813236002    A.U. after 8 cycles
SCF Done:  E(RHF) =  -230.778501460    A.U. after 9 cycles
SCF Done:  E(RHF) =  -526.813133705    A.U. after 7 cycles

作者
Author: sobereva 时间: 2019-2-26 16:23
你当前是做CCSD(T)计算，显然不能读E(RHF)，得读每次产生的CCSD(T)电子能量

作者
Author: troll 时间: 2019-2-27 08:33
本帖最后由 troll 于 2019-2-27 09:39 编辑

sobereva 发表于 2019-2-26 16:23
你当前是做CCSD(T)计算，显然不能读E(RHF)，得读每次产生的CCSD(T)电子能量

请问是 CCSD(T)= -0.52704306745D+03是电子能量吗

作者
Author: sobereva 时间: 2019-2-27 14:52

troll 发表于 2019-2-27 08:33
请问是 CCSD(T)= -0.52704306745D+03是电子能量吗

是

作者
Author: dgwfly 时间: 2022-1-31 10:45
本帖最后由 dgwfly 于 2022-1-31 10:49 编辑

想问一下，counterpoise=2 最后得到的 complexation energy = XXX kcal/mole (corrected) 可以直接拿来做为binding energy用吗，又看到了一篇博文" Gaussian教程 | 氢键与卤键的计算以及BSSE校正
"，说还是算两个单体的能量，然后自己减一下.

作者
Author: sobereva 时间: 2022-2-5 05:01

dgwfly 发表于 2022-1-31 10:45
想问一下，counterpoise=2 最后得到的 complexation energy = XXX kcal/mole (corrected) 可以直接拿 ...

看
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

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