24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sophy8069

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 455
散金: 10
帖子: 83
在线: 47.2小时
虫号: 577632
注册: 2008-06-27
性别: MM
专业: 原子与分子物理

[交流] 【求助】counterpoise中遇到的问题

我的分子体系是用DFT/B3LYP优化的，对不同的原子用的不同的基组，计算出来的能量还要BSSE吗？？用counterpoise后的输出文件中有很多行counterpoise:corrected energy=...........counterpoise:BSSE energy=...........那么在计算键能时应该用哪一行的修正呢？并且问下是E=E(AB)-E(A)-E(B)+corrected energy 还是E=E(AB)-E(A)-E(B)+BSSE energy呢????
请各位给小女子指点一下，多谢了！！！

[ Last edited by xuefei06 on 2009-2-21 at 23:17 ]

回复此楼

» 收录本帖的淘帖专辑推荐

高斯

» 猜你喜欢

书籍求助：汽车市场营销理论与实务（电子版）——章小平已经有0人回复
东方理-中科大联合博士生招聘已经有0人回复
物理化学论文润色/翻译怎么收费? 已经有50人回复
钴酸锂半电池小倍率容量上不去已经有1人回复
吉林大学材料物理本科生求问调剂信息已经有23人回复
2026第二届光电子与半导体器件前沿技术研讨会——光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有1人回复
光电子赋能·半导体创芯！已经有0人回复
Chemical Bonding at Surfaces and Interfaces,最经典的一本表面上化学相互作用教材已经有0人回复
QE计算电声耦合的时候报错Error in routine lambda (100)wrong or too many modes 已经有2人回复
分子动力学模拟合作需求已经有0人回复

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2009-02-18 17:10:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yuhuobuku

木虫 (正式写手)

应助: 2 (幼儿园)
贵宾: 0.09
金币: 1661.1
帖子: 332
在线: 6.4小时
虫号: 577652
注册: 2008-06-27
性别: MM
专业: 量子化学

★ ★ ★ ★ ★ ★
jghe(金币+3,VIP+0):感谢你参与讨论 2-19 20:24
gwdavid(金币+3,VIP+0):这种原创答案只是辛苦费也值6金币了！ 2-21 21:43

举个例子，比如计算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) 所以经过矫正的结合能由于加上了一个BSSE energy，并且这个值是大于零的，所以矫正后的结合能要比没矫正前的结合能大。

输出文件里面5个scf第一个是整体的能量(记为E(AB))，后面4个是单体能量，分别是2个部分在整体基组和各自基组下的能量。

我当时是用-Bq原子算的BSSE矫正，估计你是用CP方法吧。但原理都是一样的。

http://muchong.com/bbs/viewthread.php?tid=1099266