计算化学公社

标题: 能量分解得到的相互作用能和单点能不一致 [打印本页]

作者
Author: 七尺贱 时间: 2023-1-4 16:04
标题: 能量分解得到的相互作用能和单点能不一致
各位老师好，我研究的体系是富勒烯与碳纳米环CPP复合的弱相互作用体系。前三个是一种体系，后三个是一种体系。
在计算结合能时，在W-B97MV/def2-TZVPP级别下，采用Counterpoise方式考虑了BSSE的结合能，可以看到在该级别下，后三个体系的结合能更负。
但由于体系太大，只能使用SAPT0/def2-SVP去进行能量分解，得到的相互作用能和单点能计算是不一样的,变成后三个体系更正了，我该怎么解释呢？可以从基组太小无法正确描述体系的相互作用能来解释吗？或者说能量分解无法考虑BSSE校正，导致定性错误？

作者
Author: yflchx 时间: 2023-1-4 16:28
因为误差相抵，sSAPT0和jun-cc-pVDZ配合精度要高一些。
见：JCP 140, 094106(2014)。不能肯定，可以试一下，【文献里解释了什么是sSAPT0】。

作者
Author: 七尺贱 时间: 2023-1-4 16:34

yflchx 发表于 2023-1-4 16:28
因为误差相抵，sSAPT0和jun-cc-pVDZ配合精度要高一些。
见：JCP 140, 094106(2014)。不能肯定，可以试一下 ...

谢谢我去看下

作者
Author: zjxitcc 时间: 2023-1-4 16:56
GKS-EDA可以算，见《GKS-EDA计算简介》

作者
Author: 七尺贱 时间: 2023-1-4 17:09

zjxitcc 发表于 2023-1-4 16:56
GKS-EDA可以算，见《GKS-EDA计算简介》

谢谢老师，这个是审稿人问的，暂时换不了软件

作者
Author: sobereva 时间: 2023-1-4 22:38
高阶SAPT算不动而不得已用SAPT0时，对普通体系要用就用sSAPT0/aug-cc-VDZ，这是为什么我的Carbon 171 (2021) 514-523一文里用sSAPT0/aug-cc-VDZ：

we employed SAPT method
again to characterize the nature of intermolecular interaction. Since
the calculation cost of the SAPT2t(3)dMP2 method is prohibitively
high for the cyclo[18]carbon dimer, the scaled SAPT0 method with
jun-cc-pVDZ basis set was adapted instead.

下文也都明确说过此问题
使用PSI4做对称匹配微扰理论(SAPT)能量分解计算
http://sobereva.com/526（http://bbs.keinsci.com/thread-15902-1-1.html）

就算不说误差低消原因，def2-SVP本身也忒烂了，弥散函数都没有

SAPT本来就不存在BSSE问题，稍微了解方法的原理便知。

当前的“单点能计算”是严重不当用词，要说单点能求差，或者supramolecule方式计算相互作用能。

作者
Author: 七尺贱 时间: 2023-1-5 01:38

sobereva 发表于 2023-1-4 22:38
高阶SAPT算不动而不得已用SAPT0时，对普通体系要用就用sSAPT0/aug-cc-VDZ，这是为什么我的Carbon 171 (2021 ...

谢谢老师回复，我的体系是金属富勒烯与CPP复合，因为有Sc原子，所以用不了aug-cc-PVDZ，我用def2-svp也是想定性的说明体系主要是色散作用为主的弱相互作用

作者
Author: sobereva 时间: 2023-1-5 02:46

七尺贱发表于 2023-1-5 01:38
谢谢老师回复，我的体系是金属富勒烯与CPP复合，因为有Sc原子，所以用不了aug-cc-PVDZ，我用def2-svp也是 ...

aug-cc-pVDZ对Sc有定义

作者
Author: 七尺贱 时间: 2023-1-5 10:47

sobereva 发表于 2023-1-5 02:46
aug-cc-pVDZ对Sc有定义

谢谢老师

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