计算化学公社

标题: 博文“谈谈BSSE校正与Gaussian对它的处理”中遇到的一些困惑 [打印本页]

作者
Author: xujian 时间: 2016-10-20 20:07
标题: 博文“谈谈BSSE校正与Gaussian对它的处理”中遇到的一些困惑
看了sob老师的关于BSSE校正与Gaussian对它的处理博文，文中有这么一段要计算A、B两个分子的相互作用能，在Gaussian中使用counterpoise=2关键字（可简写为counter=2），会计算5个体系，输出的能量按照如下顺序：
E_AB：A、B基组下AB复合物的能量
E_A,bAB：A、B基组下A的能量
E_B,bAB：A、B基组下B的能量
E_A：A基组下A的能量
E_B：B基组下B的能量请问老师，为什么这里有A、B两个基组，谢谢！

作者
Author: 冰释之川 时间: 2016-10-20 20:17
本帖最后由冰释之川于 2016-10-20 20:33 编辑

按道理说，精确的分子相互作用能就应该是做差的E_A, E_B 要与 E_AB 在同一个基组（A、B）下的，这里的A、B被你理解成两个基组有点误导人，其实应该说是A分子在 A-B 整个基函数集合环境下（B相当于是虚的分子，只提供基函数）。说白了，你用同一个基组比如6-31G(d), 对于A分子所展开的基函数的总数，和对于B分子所展开的基函数的总数，是不一样的。大家都是同一个基组，基函数总数目不一样。 BSSE 其实就是为了校正（1）大家都是用（A+B）的基函数数目下的计算值E_A(A+B);E_B(A+B) 与（2）分别采用A的基函数数目下的E_A(A)与B下的E_B(B)，（1）与（2）这两种不同算法的差异！

作者
Author: dreamyeye 时间: 2016-10-20 22:20

冰释之川发表于 2016-10-20 20:17
按道理说，精确的分子相互作用能就应该是做差的E_A, E_B 要与 E_AB 在同一个基组（A、B）下的，这里的A、B ...

就是不同的基组，基函数数目不一定相同，内容也不见得完全相同，说是不同的基组是对的，只是基组的名字相同而已。

作者
Author: 冰释之川 时间: 2016-10-20 23:16

dreamyeye 发表于 2016-10-20 22:20
就是不同的基组，基函数数目不一定相同，内容也不见得完全相同，说是不同的基组是对的，只是基组的名字相 ...

所以说，还是从基函数这个角度出发理解会更好，因为如果说“A与B”两个基组，有可能会认为是名字上不同的基组，其实这两个都是同名基组，只不过里面包含的基函数数目和信息不一样而已

作者
Author: dreamyeye 时间: 2016-10-21 09:08

冰释之川发表于 2016-10-20 23:16
所以说，还是从基函数这个角度出发理解会更好，因为如果说“A与B”两个基组，有可能会认为是名字上不同的 ...

感觉不应从名字上理解，因为名字是人为定义的，而内容是绝对的。

作者
Author: liyuanhe211 时间: 2016-10-21 09:13

dreamyeye 发表于 2016-10-21 09:08
感觉不应从名字上理解，因为名字是人为定义的，而内容是绝对的。

你完全没有理解冰释之川的意思。。。

作者
Author: limen 时间: 2023-7-9 15:42

冰释之川发表于 2016-10-20 20:17
按道理说，精确的分子相互作用能就应该是做差的E_A, E_B 要与 E_AB 在同一个基组（A、B）下的，这里的A、B ...

你好，请问为什么大基组可以减小BSSE误差？基函数总数目由什么决定？

作者
Author: wzkchem5 时间: 2023-7-9 16:09

limen 发表于 2023-7-9 08:42
你好，请问为什么大基组可以减小BSSE误差？基函数总数目由什么决定？

把这个以及它里面的链接看了自然就懂了http://sobereva.com/46

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