| 查看: 1044 | 回复: 7 | ||||
| 当前主题已经存档。 | ||||
sophy8069银虫 (小有名气)
|
[交流]
【求助】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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯 |
» 猜你喜欢
博士生招生, 材料模拟方向，University of South Carolina
已经有1人回复
-
瑞典林雪平大学范立洲课题组现招聘博士后、岗位制博士生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有194人回复
-
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚聘博士后研究员
已经有4人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有1人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
中国科大仿生界面材料科学全国重点实验室戴懿涛研究员/江雷院士团队招聘科研助理
已经有0人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有1人回复
-
寻求物理化学学报最新的投稿模板！
已经有4人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2009-02-18 20:24:07
sophy8069
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 455
- 散金: 10
- 帖子: 83
- 在线: 47.2小时
- 虫号: 577632
- 注册: 2008-06-27
- 性别: MM
- 专业: 原子与分子物理
3楼2009-02-18 21:20:30
sophy8069
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 455
- 散金: 10
- 帖子: 83
- 在线: 47.2小时
- 虫号: 577632
- 注册: 2008-06-27
- 性别: MM
- 专业: 原子与分子物理
4楼2009-02-18 21:37:54
5楼2009-02-19 09:18:01
sophy8069
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 455
- 散金: 10
- 帖子: 83
- 在线: 47.2小时
- 虫号: 577632
- 注册: 2008-06-27
- 性别: MM
- 专业: 原子与分子物理
6楼2009-02-19 11:04:58
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
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) 输出文件里面5个scf第一个是整体的能量(记为E(AB))，后面4个是单体能量，分别是2个部分在整体基组和各自基组下的能量。 我当时是用-Bq原子算的BSSE矫正， 估计你是用CP方法吧。但原理都是一样的。 http://muchong.com/bbs/viewthread.php?tid=1099266 |
7楼2009-02-19 15:22:19
sophy8069
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 455
- 散金: 10
- 帖子: 83
- 在线: 47.2小时
- 虫号: 577632
- 注册: 2008-06-27
- 性别: MM
- 专业: 原子与分子物理
|
可我得到的 SCF Done: E(RB+HF-LYP) = -0.927037093185 A.U. after 5 cycles Convg = 0.5768D-06 -V/T = 2.0727 S**2 = 0.0000 Cartesian Forces: Max 0.016383271 RMS 0.005381735 Counterpoise: corrected energy = -929.170121753971 Counterpoise: BSSE energy = 0.008432377925 不只5个啊 这样的输出有好多个啊 我选的就是最后的一个，也就是 Maximum Force 0.000372 0.000450 YES RMS Force 0.000131 0.000300 YES Maximum Displacement 0.001575 0.001800 YES RMS Displacement 0.000648 0.001200 YES 这个上面的输出 |
8楼2009-02-21 15:45:33
