计算化学公社

标题: Gaussian优化结构得到结合能和PSI4能量分解Total sSAPT0能量的区别 [打印本页]

作者
Author: Jerryluo 时间: 2023-4-16 01:49
标题: Gaussian优化结构得到结合能和PSI4能量分解Total sSAPT0能量的区别
请问 Gaussian优化结构得到A和B的结合能（EAB-EA-EB），在PSI4中A和B作用能量分解也有一个Total SAPT0能量，此时要选择哪一个能量作为A和B的结合能。还是说Gaussian和PSI4计算结果都应当在文章列出而不冲突（两者结果有较大差别，分别为-6kal/mol和-10kcal/mol）。

作者
Author: sobereva 时间: 2023-4-16 13:59
Gaussian优化结构根本得不到结合能。表述必须严谨准确

计算级别主动交代清楚，不要有任何歧义和隐藏的关键性细节
仔细看此文
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

作者
Author: Jerryluo 时间: 2023-4-16 21:52

sobereva 发表于 2023-4-16 13:59
Gaussian优化结构根本得不到结合能。表述必须严谨准确

计算级别主动交代清楚，不要有任何歧义和隐藏的关 ...

A,B,AB复合物分别在 b3lyp-D3/6-31g* 水平优化结构，能量在b3lyp-D3/def2-tzvp水平计算，EAB-EA-EB得到结合能。然后AB复合物结构用PSI4在sSAPT0/jun-cc-pVDZ水平进行能量分解，A，B分别为片段1和片段2

作者
Author: sobereva 时间: 2023-4-16 23:32

Jerryluo 发表于 2023-4-16 21:52
A,B,AB复合物分别在 b3lyp-D3/6-31g* 水平优化结构，能量在b3lyp-D3/def2-tzvp水平计算，EAB-EA-EB得到结 ...

b3lyp-D3/def2-tzvp直接算的结合能没法用，BSSE太大（假设你研究的是弱相互作用），至少也得考虑counterpoise或gCP做BSSE校正。用B3LYP-D3(BJ)/ma-def2-TZVP结合counterpoise能得到质量好的多的结果。但是否精确，仍然得看具体体系。
如果对BSSE都没概念，看此文
谈谈BSSE校正与Gaussian对它的处理
http://sobereva.com/46

只有两个都是高精度的级别，算的相互作用能才可能相差不大。

作者
Author: Jerryluo 时间: 2023-4-16 23:56

sobereva 发表于 2023-4-16 23:32
b3lyp-D3/def2-tzvp直接算的结合能没法用，BSSE太大（假设你研究的是弱相互作用），至少也得考虑counterp ...

好的谢谢老师解答

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