应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 个人文集 » 【gchcaj个人文集】浅谈Orca计算BSSE校正
21/1返回列表
查看: 2375  |  回复: 5
查看全部回帖 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价3次，作者ghcacj增加金币 2.4 个
本帖产生 1 个 QC强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ghcacj

荣誉版主 (著名写手)

[资源] 【gchcaj个人文集】浅谈Orca计算BSSE校正

今天看到lihb734兄写了个关于BSSE校正的帖子，颇有启发，所以我也来写一个关于使用Orca校正BSSE的，大家看到的Gaussian的BSSE校正比较多，可以说，Gaussian在计算BSSE的input文件上比较方便，给出两个分子的坐标以及Counterpoise的关键词即可，上手容易，但是细究中间的原理过程，我个人orca在这里好些，废话不多说，先贴出BSSE的定义是，这个是Boys和Bernardi二位搞出来的，如公式1：
（1）
     初看起来比较复杂，其实不难，表达式中的下表表示计算的片段，AB表示二聚体，A或B表示单聚体；上标表示分子所处的优化几何坐标，AB表示处于二聚体的优化几何坐标，A或B表示单聚体优化几何坐标；括号里表示使用的基组。这样就可以把这个分成七项四个大项，第一项是二聚体的能量，第二项是两个单聚体能量之和，第三项是处于二聚体优化几何坐标下的使用单聚体部分的基组计算两个单聚体的能量之和，第四项是处于二聚体优化几何坐标下的使用二聚体部分的基组计算两个单聚体的能量之和。所以公式2可以改写为：
（2）
    这里dimer表示二聚体，monomer2表示原公式的第四项，monomer1表示原公式的第三项，monomer0表示原公式的第二项。如果只计算前三项，就是被校正后的能量。
        现在用一个例子来说明下，也是两个水分子的作用，在Orca手册91-93页：
CODE:
        #
        # BSSE test
        #
        # --------------------------------------------
        # First the monomer. It is a waste of course
        # to run the monomer twice ...
        # --------------------------------------------
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_0"
        * xyz 0 1
         O     7.405639     6.725069     7.710504
         H     7.029206     6.234628     8.442160
         H     8.247948     6.296600     7.554030
        *
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_0"
        * xyz 0 1
         O     7.405639     6.725069     7.710504
         H     7.029206     6.234628     8.442160
         H     8.247948     6.296600     7.554030
        *
        # --------------------------------------------
        # now the dimer
        # --------------------------------------------
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "dimer"
        * xyz 0 1
        O   7.439917   6.726792   7.762120
        O   5.752050   6.489306   5.407671
        H   7.025510   6.226170   8.467436
        H   8.274883   6.280259   7.609894
        H   6.313507   6.644667   6.176902
        H   5.522285   7.367132   5.103852
        *
        # --------------------------------------------
        # Now the calculations of the monomer at the  
        # dimer geometry
        # --------------------------------------------
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_1"
        * xyz 0 1
        O   7.439917   6.726792   7.762120
        H   7.025510   6.226170   8.467436
        H   8.274883   6.280259   7.609894
        *
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_1"
        * xyz 0 1
        O   5.752050   6.489306   5.407671
        H   6.313507   6.644667   6.176902
        H   5.522285   7.367132   5.103852
        *
        # --------------------------------------------
        # Now the calculation of the monomer at the  
        # dimer geometry but with the dimer basis set
        # --------------------------------------------
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_2"
        * xyz 0 1
        O   7.439917   6.726792   7.762120
        O :  5.752050   6.489306   5.407671
        H   7.025510   6.226170   8.467436
        H   8.274883   6.280259   7.609894
        H :  6.313507   6.644667   6.176902
        H :  5.522285   7.367132   5.103852
        *
        $new_job
        ! RHF MP2 TZVPP VeryTightSCF XYZFile PModel
        %id "monomer_2"
        * xyz 0 1
        O :  7.439917   6.726792   7.762120
        O   5.752050   6.489306   5.407671
        H :  7.025510   6.226170   8.467436
        H :  8.274883   6.280259   7.609894
        H   6.313507   6.644667   6.176902
        H   5.522285   7.367132   5.103852
        *

Orca的常规语法这里不再解释，只是介绍下输入文件中的冒号作用，这个冒号是专门为这个BSSE计算考虑的，把冒号跟在原子后面，该原子的电子和核电荷均被删除，但是保留原子以被分配基组。
    这个任务其实是通过建立多个任务，再用关键词%id将每大项的两个结果能量联系起来，部分输出文件如下：
Monomer      :   -152.647062118 Eh
Dimer        :    -152.655623625 Eh   -5.372 kcal/mol
Monomer at dimer geometry:   -152.647006948 Eh    0.035 kcal/mol
Same with AB Basis set  :     -152.648364970 Eh   -0.818 kcal/mol
       按照公式2，这时候就不难计算了，ΔE=-152.655623625 -(-152.648364970+152.647006948) +152.647062118=-0.007203485 Eh
       计算到了这里，应该说Orca不如Gaussian那么智能，因为Gaussian一个关键词，Orca要写这么多，但是做计算的经验告诉我，越是书写麻烦，说明自定义的空间就越大，这里就体现如果我们要做一个Scan的话，Orca不用一个点一个点的计算，而且可以将BSSE校正中的每项拆开，分别做Scan或者SP，最后通过公式进行处理即可。这点Gaussian我好像发现不行，因为Gaussian的Counterpoise包含了太多的东西，中间的东西无法单独调用（本人用Gaussian较少，如果说错了，请指出，不过别骂娘就好）。这样在计算时方便些（个人感觉）
我也来个免责声明：
本人主业是Monte Carlo，对于QC是外行，所以这个基本是本人的一家之言，如果有问题请指出，不承担任何学术上的责任。

[ Last edited by 375642546 on 2011-5-6 at 17:25 ]
回复此楼

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

量化软件学习

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-05-06 12:45:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zpppanda111

金虫 (小有名气)

★★★★★ 五星级,优秀推荐

引用回帖:
Originally posted by ghcacj at 2011-05-06 18:00:43:
不知道你指的是不是我贴出来的计算过程？BSSE能量是很小，而且我用的基组还算比较大，不过对于MP2，一般都校正的。

呵呵，学习了

对于BSSE，也会导致大的误差，(almost 10 kcal mol-1 for benzene–Na+ complex at the MP2(full)/aug-cc-pVTZ level.)
C. Estarellas, X. Lucas, A. Frontera , D. Quiñonero and P. M. Deyà, Chem. Phys. Lett., 2010, 489, 254.),
赞一下 回复此楼
4楼2011-05-08 07:16:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghcacj 的主题更新
21/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭