24小时热门版块排行榜    

查看: 896  |  回复: 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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophy8069

银虫 (小有名气)

可我得到的
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

jghe

金虫 (著名写手)

★ ★
xuefei06(金币+2,VIP+0):thanks! 2-18 23:31
要用这个公式E=E(AB)-E(A)-E(B) 就要考虑BSSE校正了,一般如果不是若相互作体系不用考虑
输出文件中第一个是总能量,第二个是Part1(校正过)的,第三个是Part2(校正过)的 ,第四个是Part1(未校正过)的,第五个是Part2(未校正过)的。
自己把能量列出来,比一比就知道了
2楼2009-02-18 20:24:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophy8069

银虫 (小有名气)

你说的是没用counterpoise时的输出文件还是用了的啊??再有我用counterpoise算得的输出文件中有好多行counterpoise输出结果,我不知道用哪一次的,是不是用最后收敛前的counterpoise????
3楼2009-02-18 21:20:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophy8069

银虫 (小有名气)

请问在counterpoise输出文件从哪可以找到最终的BSSE energy的值啊???哪个好心人给解答一下啊!!!
4楼2009-02-18 21:37:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见