24小时热门版块排行榜    

查看: 3956  |  回复: 9
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tailor

金虫 (小有名气)

[交流] 【求助】siesta计算中,如何进行BSSE修正?

即基组重叠误差,全称为Basis Set Superposition Error(BSSE)。对于双分子A和B的相互作用能应该表述为E_interaction = E_AB - E_A - E_B + E_BSSE。第三项即为修正项。对于弱相互作用,E_BSSE所占E_interaction的比例往往不小,甚至超过它,如果不进行校正,可能正负号都不对。
      这些也是从网上查到的,但是网上主要介绍的都是高斯算法中的BSSE修正,由于本人用siesta算,想麻烦各位高手多多赐教!!
回复此楼

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

atk&transiesta

» 猜你喜欢

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

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

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

夕阳西下

金虫 (著名写手)

★ ★
tailor(金币+2):谢谢参与
qasd(金币+1):谢谢 2010-09-11 18:57:35
即使基矢正交,A与B反应后的整个空间基矢也会发生改变。应该是原来两个空间的直积。如果A和B是同种原子应该不用修正。
7楼2010-09-11 12:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

ifmc1234

木虫 (著名写手)

★ ★ ★
tailor(金币+2):谢谢参与
qasd(金币+2):谢谢 2010-09-11 18:57:22
引用回帖:
Originally posted by tailor at 2010-09-06 21:06:56:
即基组重叠误差,全称为Basis Set Superposition Error(BSSE)。对于双分子A和B的相互作用能应该表述为E_interaction = E_AB - E_A - E_B + E_BSSE。第三项即为修正项。对于弱相互作用,E_BSSE所占E_interaction的 ...

不知道你算的是分子还是固体,建议不要用它来讨论。你可在userguide里搜一下,我没看有什么人用siesta讨论过BSSE,在userguide里的27页提到了BSSE。不过具体怎么做还不知道。说明书中只提到PAO.Basis里设置时,用Ghost原子可以处理BSSE。不过如果根据你上面的定义,我想可以这样E_interaction = E_AB - E_A - E_B + Sum_E_BSSE,其中Sum_E_BSSE为A,B单个的BSSE之和。也就是说用鬼原子来处理BSSE,然后根据上面求的
2楼2010-09-06 21:29:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)


tailor(金币+2):谢谢参与
引用回帖:
Originally posted by tailor at 2010-09-06 14:06:56:
即基组重叠误差,全称为Basis Set Superposition Error(BSSE)。对于双分子A和B的相互作用能应该表述为E_interaction = E_AB - E_A - E_B + E_BSSE。第三项即为修正项。对于弱相互作用,E_BSSE所占E_interaction的 ...

为什么要算这个?什么体系才学要考虑BSSE修正啊?
我不是很懂,但是对你的问题很感兴趣,想问问你。
弘德明志博学笃行
3楼2010-09-06 21:30:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tailor

金虫 (小有名气)

请教老外,得到的回复


zxzj05(金币+1):谢谢回帖交流 2010-09-25 22:34:07
引用回帖:
Originally posted by zhangguangping at 2010-09-06 21:30:41:

为什么要算这个?什么体系才学要考虑BSSE修正啊?
我不是很懂,但是对你的问题很感兴趣,想问问你。

What we calculated is the Basis Set Superposition Error (BSSE) that is an error occurring when you use a finite basis set to study, for example, a chemical reaction. In a tipical reaction you calculate the electronic structure of the reactants separately, let us say A and B, and then the electronic structure of the product, A+B. What happens in this case is that you are using different basis sets to represent the reactants in each step of the calculation. When you calculate the isolated reactant A, you use the basis set of A, and the electrons of A will be arranged spatially on these provided orbitals. However, when you put both systems together in the A+B, the electrons of A will have now a greater spatial freedom, since they can now occupy orbitals coming from B. The same happens with B. You have different basis sets to distribute the electrons of A in the two situations.
To correct this, you need to recalculate the energy of the isolated reactants (A and B) by considering the basis set of both reactants, i.e., the calculations of the isolated reactant A will now include the basis sets of B (without the system B itself). It needs to be done for each configuration that you take into account in the reaction of A+B. Now you will have the same basis set for the isolated reactants and the product.
You can look for more information searching for "counterpoise correction"
5楼2010-09-10 19:36:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见