博文“谈谈BSSE校正与Gaussian对它的处理”中遇到的一些困惑

xujian · 发表于 Post on 2016-10-20 20:07:22

看了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两个基组，谢谢！

冰释之川 · 发表于 Post on 2016-10-20 20:17:58

本帖最后由冰释之川于 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）这两种不同算法的差异！

dreamyeye · 发表于 Post on 2016-10-20 22:20:14

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

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

冰释之川 · 发表于 Post on 2016-10-20 23:16:43

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

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

dreamyeye · 发表于 Post on 2016-10-21 09:08:47

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

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

liyuanhe211 · 发表于 Post on 2016-10-21 09:13:44

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

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

limen · 发表于 Post on 2023-7-9 15:42:24

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

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

wzkchem5 · 发表于 Post on 2023-7-9 16:09:55

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 博文“谈谈BSSE校正与Gaussian对它的处理”中遇到的一些困惑

评分 Rate