计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 求助：用sobEDA计算出的静电项是正值，是哪里出问题了呢

求助：用sobEDA计算出的静电项是正值，是哪里出问题了呢

查看数: 357 | 评论数: 3 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

liruixue

发布时间: 2025-8-27 11:14

正文摘要:

本帖最后由 liruixue 于 2025-8-27 11:20 编辑各位老师好，我在计算sobEDA能量分解时，得出的静电项是正值，是哪里出问题了呢?用GKS-EDA算出的静电项是负值

回复 Reply

sobereva 发表于 Post on 2025-8-29 05:14:45

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

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

liruixue 发表于 Post on 2025-8-28 21:08:36

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

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

sobereva 发表于 Post on 2025-8-28 04:43:08

不要用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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register