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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » 个人文集 » 【gchcaj个人文集】浅谈Orca计算BSSE校正
61/1返回列表
查看: 2896  |  回复: 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的回帖

lihb734

铁杆木虫 (职业作家)

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

我理解的是这个计算过程是包括了单体的准备能的，虽然这个能量很小啊。
赞一下 回复此楼
2楼2011-05-06 15:55:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ghcacj

荣誉版主 (著名写手)
引用回帖:
Originally posted by lihb734 at 2011-05-06 15:55:02:
我理解的是这个计算过程是包括了单体的准备能的，虽然这个能量很小啊。

不知道你指的是不是我贴出来的计算过程？BSSE能量是很小，而且我用的基组还算比较大，不过对于MP2，一般都校正的。
赞一下 回复此楼
3楼2011-05-06 18:00:43
已阅   回复此楼   关注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

荣誉版主 (著名写手)
引用回帖:
Originally posted by zpppanda111 at 2011-05-08 07:16:58:
呵呵，学习了

对于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 an ...

呵呵，学习了，But，Orca的人气真的比Gaussian差好多啊。
赞一下 回复此楼
5楼2011-05-09 10:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

green1986

木虫 (小有名气)

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

学习了
回复此楼
6楼2011-05-09 10:56:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ghcacj 的主题更新
61/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 可靠消息，8月21号(明天)放榜 +9 医学老男孩 2026-08-20 10/500 2026-08-20 10:21 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +20 yaoyewhu2008 2026-08-20 20/1000 2026-08-20 09:41 by 2026年国自然基�
[基金申请] 时间戳变了，能看出什么问题？ +12 基诺咪客 2026-08-17 16/800 2026-08-19 22:04 by kk 的小木虫
[基金申请] filecode，4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 重要消息，中午系统在维护 +10 yuleib84 2026-08-18 11/550 2026-08-19 07:30 by 羽毛枫f
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭