计算化学公社

标题: sobEDAw得到的ΔE_int与counterpoise得到未校正的ΔE_raw相等该如何处理？ [打印本页]

作者
Author: chlorophyll 时间: 2024-11-16 19:17
标题: sobEDAw得到的ΔE_int与counterpoise得到未校正的ΔE_raw相等该如何处理？
本帖最后由 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方法校正后的。

(, 下载次数 Times of downloads: 5)

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

下面是计算的log文件，泛函B3LYP，基组ma-TZVP
(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 0)

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

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

作者
Author: chlorophyll 时间: 2024-11-20 15:02

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)