24小时热门版块排行榜    

查看: 1057  |  回复: 4

zhimin1106

铁杆木虫 (正式写手)

[求助] 就真没人会吗?BSSE和ZPE同时校正问题

要获得经BSSE校正的相互作用能,直接用:E_interaction = E_AB - E_A,bAB - E_B,bAB,就可以了,就是使用counterpoise=2关键字计算5个体系能量后,第一个减去第二个再减去第三个,我现在的问题是,如果我的体系还要同时做ZPE校正的话,该如何处理?请详细给出过程,万分感激!

下面是我查到的说法,几个部分我也都算了,但就是不知道怎么才能得到最后经过ZPE和BSSE同时校正的作用能。也就是说具体从.out文件哪些数据算。
---------------------------------------------------------------------------------------------

零点能校正的计算是按通常的方法进行的,即:分别优化分子AB及其片段A和B,得到稳定构型后加FREQ计算零点能校正。这与BSSE的计算有两点不同:一是在BSSE的计算中,片段A和B的构型必须与它们在分子AB中的一样,而ZPE的计算中需要对A、B重新优化以得到稳定的构型才行,否则ZPE是无意义的;二是在BSSE的计算中,计算片段A、B所用的基集与计算分子AB一样(计算A时保留B的基集但把B中所有原子换位赝原子;同样计算B时则保留A的基集并把A中所有原子换为赝原子),而ZPE的计算中三者的基集一般不同,并且不能引入赝原子。

[ Last edited by zhimin1106 on 2013-10-6 at 21:50 ]
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

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

zhimin1106

铁杆木虫 (正式写手)

细节问题,没人能答吗?
2楼2013-10-19 10:51:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

showcat

金虫 (小有名气)

【答案】应助回帖

E=EAB-EA-EB+BSSE 其中E EAB EA EB都是在频率计算后用零点校正之后的能量就可以了。单位一致,都是Hartree
化工人,初入计算领域
3楼2013-11-28 16:12:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

marson

银虫 (正式写手)

【答案】应助回帖

使用BSSE校正的时候所用的所有的能量都用对于每个体系其分别经过zpe校正过能量即可,对应的都是log文件中的the sum of energy and zero point energy这个就可以了~
思想有多远,我们就能走多远~
4楼2013-11-28 20:47:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

有点钻牛角尖了。用bsse校正后的能量算zpe,与分别计算bsse和zpe的总能量相差不大。如果想提高精度,还不如算cbs,远比这个误差大得多。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2013-11-28 23:45:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhimin1106 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-17 04:28 by k06BKNblrGRH
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:14 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-17 02:38 by Tql5LQhh5rLK
[博后之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 mibUvS8DDCwf 2026-09-15 3/150 2026-09-16 23:28 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 vZfe6xYu34yj 2026-09-14 4/200 2026-09-16 22:40 by Tql5LQhh5rLK
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 21:14 by MHf40yA0Uhyb
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 QUjhNVAcOSff 2026-09-13 9/450 2026-09-16 18:36 by rRvxEj3zorCy
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:03 by i3KBRD19ERkG
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-16 15:48 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 0pYnUiPDfdnk 2026-09-14 3/150 2026-09-16 06:51 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:27 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:35 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 09:38 by VMDqgnlKimtZ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 07:02 by 5BDX0d0WFp7t
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见