| 查看: 459 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
[求助]如何使用完全均衡校正法消除BSSE
|
|||
|
计算复合物结构的相互作用能时，对由基函数引起的基组叠加误差（BSSE），采用Boys和Bernardi的完全均衡校正法可以消除，但是我不知道如何具体操作？ 例如a+b生成ab，已经由Gauss计算出Ea, Eb和Eab，一般来说，作用能可以用Ea+Eb-Eab计算，但是不太准确。请问大牛们，如何消除BSSE？  |
回复此楼» 猜你喜欢
河南理工大学化工学院招收调剂生啦！
已经有1人回复
-
硕士研究生调剂
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有189人回复
-
投稿期刊求助
已经有1人回复
-
江西师大化学工程学院调剂今晚0点开放！ 缺额100个，速冲！
已经有0人回复
-
图书求助：Fluorescent Probes for Bioactive Species: Design
已经有0人回复
-
同系列相似结构的单层二维COF真空层厚度该如何设置？
已经有0人回复
spkeey
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1244.9
- 散金: 3
- 帖子: 308
- 在线: 14小时
- 虫号: 320174
- 注册: 2007-03-10
- 性别: GG
- 专业: 理论和计算化学
2楼2008-12-19 12:18:15
yuhuobuku
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 贵宾: 0.09
- 金币: 1661.1
- 帖子: 332
- 在线: 6.4小时
- 虫号: 577652
- 注册: 2008-06-27
- 性别: MM
- 专业: 量子化学
★ ★ ★ ★ ★ ★
linlin1011(金币+2,VIP+0):说明得太详细了，非常感谢！以后有问题还得向您多多请教。
lei0736(金币+4,VIP+0):谢谢
linlin1011(金币+2,VIP+0):说明得太详细了，非常感谢！以后有问题还得向您多多请教。
lei0736(金币+4,VIP+0):谢谢
|
举个例子，比如计算AB体系的BSSE矫正： 需要得出的五个能量分别是： 第一个是整体的能量E(AB) 第二个是单独计算A的能量E(A) 第三个是单独计算B的能量E(B) 第四个是计算在AB基组下，A的能量E(A/AB) 第五个是计算在AB基组下，B的能量E(B/AB) 那么，BSSE energy =E(A)+E(B)-E(A/AB)-E(B/AB)。计算E(A/AB)的基函数是计算AB的基函数，而计算E(A)的基函数只是计算AB基函数的一部分，基函数越多，能量计算的越准确，能量计算越准确，就越低，所以E(A/AB) 输出文件里面5个scf第一个是整体的能量(记为E(AB))，后面4个是单体能量，分别是2个部分在整体基组和各自基组下的能量。 我当时是用-Bq原子算的BSSE矫正， 估计你是用CP方法吧。但原理都是一样的。 http://muchong.com/bbs/viewthread.php?tid=1099266 [ Last edited by yuhuobuku on 2008-12-19 at 12:23 ] |
3楼2008-12-19 12:22:33
|
ONIOM模型 （DFT/genecp: DFT/3-21G) 一个钌的五配位化合物 能量：ONIOM: extrapolated energy = -1813.606312467803 该配合物加上一个三苯基膦后 后能量反而升高 ，请问这是为什么，是我的计算方法不对么？ 三苯基膦 同样用（DFT/3-21g**: DFT/3-21G) 与五配位 六配位中的化合物都采用相同的机组和方法，三苯基膦的能量是：ONIOM: extrapolated energy = -1030.872640175848 六配位的能量是： ONIOM: extrapolated energy = -2844.450299687489 这样求出的6配位与5配位的能量之差为：17.98kcal/mol，即六配位的能量更高，更不稳定，请问这是为什么？ 知道的朋友请帮忙解释一下！谢谢！ 三苯基膦与5，6配位化合物中的三苯基膦所用基组完全相同，请问楼上的我需要bsse校正么？ |
4楼2008-12-19 15:54:16
5