24小时热门版块排行榜    

查看: 3968  |  回复: 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的回帖

zhangguangping

木虫 (著名写手)


zxzj05(金币+1):谢谢回帖交流 2010-09-25 22:34:19
引用回帖:
Originally posted by tailor at 2010-09-10 12:36:33:



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 ...

意思就是说因为计算一个体系的时候前后,使用的基组不一样大造成了一定的误差,这个误差需要修正。如果第二次引入的那部分基函数和第一次的基函数都正交的话,也就是各自是一个子空间的话,这个修正就不用做了吧?
弘德明志博学笃行
6楼2010-09-11 10:02:01
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +12 爱看书的可乐 2026-08-26 14/700 2026-08-26 12:49 by coolysnow
[基金申请] 为什么国自然不能直接公布 +3 bjdxyxy 2026-08-26 3/150 2026-08-26 12:48 by 小丑鱼爱读书
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +3 yuleib84 2026-08-26 3/150 2026-08-26 12:34 by xiacongxu
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +24 木水思豆 2026-08-25 27/1350 2026-08-26 08:26 by yanchen918
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 19/950 2026-08-25 17:57 by chick875
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天,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
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见