| 查看: 3953 | 回复: 9 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】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 |
» 猜你喜欢
人气不行了
已经有10人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有6人回复
-
我面上完蛋了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
-
filecode，4个jtjc了
已经有17人回复
-
没有任何消息-是不是就凉了
已经有6人回复
-
2026国自然函评费到账
已经有19人回复
-
建议基金发布提前给出明确的时间点
已经有16人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有21人回复
-
估计是周四
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何萃取4-氯丁醇
已经有6人回复
- gaussian计算
已经有4人回复
- 急！ 关于Ms从Windows提交到Linux计算的问题？
已经有3人回复
- 请教：orca 中用CCSD(T)计算单点能出错--提示内存不足 该怎么办？
已经有10人回复
- 难道每次计算前都咬进行收敛性测试吗？
已经有14人回复
- siesta 计算电子占据数的问题
已经有3人回复
- siesta计算dipole moment（电偶极矩）
已经有7人回复
- 关于SIESTA中如何计算功函数的一种做法探讨
已经有8人回复
- siesta计算弱相互作用VDW
已经有17人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】SIESTA 中ATOM生成赝势测试后如何看赝势的好坏
已经有4人回复
- 【求助】如何在siesta中设置ghost state原子的计算
已经有16人回复
- 【求助】BSSE计算后如何进行零点能校正？
已经有8人回复
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
3楼2010-09-06 21:30:41
ifmc1234
木虫 (著名写手)
- 应助: 18 (小学生)
- 金币: 2797.9
- 散金: 6001
- 红花: 14
- 帖子: 1761
- 在线: 572.6小时
- 虫号: 825500
- 注册: 2009-08-10
- 专业: 理论和计算化学
★ ★ ★
tailor(金币+2):谢谢参与
qasd(金币+2):谢谢 2010-09-11 18:57:22
tailor(金币+2):谢谢参与
qasd(金币+2):谢谢 2010-09-11 18:57:22
|
不知道你算的是分子还是固体，建议不要用它来讨论。你可在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
请教老外，得到的回复
★
zxzj05(金币+1):谢谢回帖交流 2010-09-25 22:34:07
zxzj05(金币+1):谢谢回帖交流 2010-09-25 22:34:07
|
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
zhangguangping
木虫 (著名写手)
- 1ST强帖: 23
- 应助: 71 (初中生)
- 贵宾: 0.031
- 金币: 1974.1
- 散金: 9059
- 红花: 79
- 帖子: 2681
- 在线: 2128.8小时
- 虫号: 529624
- 注册: 2008-03-20
- 性别: GG
- 专业: 原子和分子物理
6楼2010-09-11 10:02:01
20