计算化学公社

标题: 求助：用sobEDA计算出的静电项是正值，是哪里出问题了呢 [打印本页]

作者
Author: liruixue 时间: 2025-8-27 11:14
标题: 求助：用sobEDA计算出的静电项是正值，是哪里出问题了呢
本帖最后由 liruixue 于 2025-8-27 11:20 编辑

各位老师好，我在计算sobEDA能量分解时，得出的静电项是正值，是哪里出问题了呢?用GKS-EDA算出的静电项是负值

作者
Author: sobereva 时间: 2025-8-28 04:43
不要用M06L，没任何好处
我用B3LYP-D3(BJ)/def2-TZVP重新做了sobEDA计算，结果非常合理，静电部分也不是负值
Total interaction energy: -22.83 kcal/mol

Physical components of interaction energy derived by sobEDA:
Electrostatic (E_els): -27.70 kcal/mol
Exchange (E_x): -22.20 kcal/mol
Pauli repulsion (E_rep):    73.81 kcal/mol
Exchange-repulsion (E_xrep = E_x + E_rep):    51.61 kcal/mol
Orbital (E_orb): -30.39 kcal/mol
DFT correlation (E_DFTc):    -7.65 kcal/mol
Dispersion correction (E_dc):    -8.70 kcal/mol
Coulomb correlation (E_c = E_DFTc + E_dc): -16.35 kcal/mol

作者
Author: liruixue 时间: 2025-8-28 21:08

sobereva 发表于 2025-8-28 04:43
不要用M06L，没任何好处
我用B3LYP-D3(BJ)/def2-TZVP重新做了sobEDA计算，结果非常合理，静电部分也不是负 ...

谢谢老师，老师我想请问一下为什么M06L不适合做sobEDA呢，审稿人提出用SAPT等方法证明GKS-EDA数据的可靠性，但是因为条件有限算不动SAPT，就想着用sobEDA代替SAPT

作者
Author: sobereva 时间: 2025-8-29 05:14

liruixue 发表于 2025-8-28 21:08
谢谢老师，老师我想请问一下为什么M06L不适合做sobEDA呢，审稿人提出用SAPT等方法证明GKS-EDA数据的可靠 ...

干嘛用这个，又没任何优点，本身又不适合算当前体系

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