sobEDAw得到的ΔE_int与counterpoise得到未校正的ΔE_raw相等该如何处理？

chlorophyll · 发表于 Post on 7 day ago

本帖最后由 chlorophyll 于 2024-11-16 19:20 编辑

我打算探究六氟磷酸根与水分子之间的相互作用，在用sobEDAw进行能量分解的时候有些地方不太清楚，想在这里请教各位老师：
1.如图，左边是使用counterpoise方法计算得到的ΔE_raw(Eab-Ea-Eb)，以及校正后的ΔE_corr(Eab-Ea-Eb+E_BSSE)，右边是经过sobEDAw程序得到的ΔE_int，但为什么我这里ΔE_int的值未被校正呢？是不是这样的结果不能使用呢？因为在社长的博文提到：使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析 - 思想家公社的门口：量子化学·分子模拟·二次元

在sobEDA或sobEDAw能量分解计算中，如果计算每个单体都使用整个体系里所有原子的基函数，则给出的ΔE_int及其各个能量成份就都相当于是counterpoise方法校正后的。

2.如果出现这种情况的话，能不能用不较正的ΔE进行能量分解分析，最后计算总的结合能的时候再把E_BSSE加上呢？这样做会不会对能量分解的分析有影响呢？

下面是计算的log文件，泛函B3LYP，基组ma-TZVP

PF6-H2O-1.log (114.56 KB, 下载次数 Times of downloads: 1)

result.txt (2.82 KB, 下载次数 Times of downloads: 0)

sobereva · 发表于 Post on 6 day ago

博文说了，必须把sobEDA.sh脚本里iCP设成1，结果才相当于考虑counterpoise的情况

• iCP=0代表不用counterpoise校正，=1代表用（对应基组带w.CB后缀的情况）。分解弱相互作用的时候强烈建议用iCP=1以改进结果，分解化学键作用时则不建议iCP=1，不仅浪费很多时间，还可能令结果更差，见《计算化学键键能时考虑BSSE不仅是多余的甚至是有害的》（http://sobereva.com/381）。

chlorophyll · 发表于 Post on 3 day ago

sobereva 发表于 2024-11-17 02:28
博文说了，必须把sobEDA.sh脚本里iCP设成1，结果才相当于考虑counterpoise的情况

谢谢老师！是我疏忽漏看了这一处。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] sobEDAw得到的ΔE_int与counterpoise得到未校正的ΔE_raw相等该如何处理？

浏览过的版块