计算化学公社

标题: 求助：BSSE校正后出现较大的正值是什么原因？ [打印本页]

作者
Author: suosuosky 时间: 2019-4-26 11:59
标题: 求助：BSSE校正后出现较大的正值是什么原因？
请问老师，在对复合物进行BSSE校正时，出现较大的正数，是什么原因？
谢谢！

输入文件：
#p B2PLYPD3/jul-cc-pVTZ counterpoise=2

输出结果：
Counterpoise corrected energy = -1957.912284745451
               BSSE energy =    0.015461552645
            sum of fragments = -1957.924790627418
         complexation energy =    -1.85 kcal/mole (raw)
         complexation energy =    7.85 kcal/mole (corrected)

作者
Author: sobereva 时间: 2019-4-26 12:39
说不定是算某些能量时碰巧没收敛到稳定波函数

作者
Author: suosuosky 时间: 2019-4-26 13:23

sobereva 发表于 2019-4-26 12:39
说不定是算某些能量时碰巧没收敛到稳定波函数

谢谢老师。请问这种情况下应该怎么办？其他系列的复合物都是使这个级别计算出的结合能，那应该换一种方法或基组试试吗？

作者
Author: sobereva 时间: 2019-4-27 00:51
利用鬼原子Bq手动实现counterpoise任务中的计算二聚体基函数下各自单体能量的任务，所有计算任务都加stable关键词检测波函数稳定性，确定了波函数稳定，再按照counterpoise公式计算BSSE校正能加到常规方式算的结合能上

作者
Author: suosuosky 时间: 2019-4-28 10:28
本帖最后由 suosuosky 于 2019-4-28 10:32 编辑

sobereva 发表于 2019-4-27 00:51
利用鬼原子Bq手动实现counterpoise任务中的计算二聚体基函数下各自单体能量的任务，所有计算任务都加stable ...

谢谢sob老师的回复。

除了出现上面BSSE校正后的结合能是正数之外，我在计算复合物中还遇到一个问题，我用两种方法计算了一个复合物体系的结合能：

方法（1）BE-Eab-Ea-Eb，其中E均为优化后的单体和复合物的经零点能校正后的能量，得到的数值为 -1.99 kcal/mol

方法（2）用优化好的AB构型，进行BSSE校正，最后输出的：
Counterpoise corrected energy = -877.197861806622
               BSSE energy =    0.001494458230
            sum of fragments = -877.163834122510
         complexation energy =    -22.29 kcal/mole (raw)
         complexation energy = -21.35 kcal/mole (corrected)

为什么两种方法计算出的结合能相差这么多，我已经反复检查了好几遍，能量没有找错。此外，在进行BSSE校正的输入文件中，我使用了不同类型的方法，包括m062x, mp2, b3plyp-d3(附件)，结果都是 -20多 kcal/mol。

我计算了10个体系中，有3个复合物体系出现这种问题，其他复合物体系计算结果中这两种方法得到的结合能相差都比较小。

请问老师是什么原因呢？应该怎样解决呢?  谢谢老师

作者
Author: suosuosky 时间: 2019-4-28 11:00
本帖最后由 suosuosky 于 2019-4-29 18:38 编辑

suosuosky 发表于 2019-4-28 10:28
谢谢sob老师的回复。

除了出现上面BSSE校正后的结合能是正数之外，我在计算复合物中还遇到一个问题， ...

附件这是另外两个体系的，把Si换成Ge之后，结果Si是正常的，Ge是异常大。

按照理论上讲，应该不会出现这么大的变化。

请问sob老师，这种情况应该如何解决？

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