计算化学公社

标题: 求助sobEDAw进行EDA分析结果出现数值很大，三百多kcal/mol这合理吗 [打印本页]

作者
Author: 悠哉2020 时间: 2026-2-23 13:05
标题: 求助sobEDAw进行EDA分析结果出现数值很大，三百多kcal/mol这合理吗
各位老师好，在用sobEDAw进行能量分解分析时，这几个值都很大 Exchange (E_x): -107.53 kcal/mol
Pauli repulsion (E_rep): 323.80 kcal/mol Exchange-repulsion (E_xrep = E_x + E_rep): 216.26 kcal/mol
Orbital (E_orb): -121.21 kcal/mol请问这是什么原因，这数据合理吗。

作者
Author: kirSTY 时间: 2026-2-23 13:33

1. 检查fragments的电荷、自旋度设置是否正确
2. 另外还可以用ELF、RDG分析等从另一个角度核验当前结果

作者
Author: wzkchem5 时间: 2026-2-24 09:39
切断共价键、强配位键的结果就是会这么大。可以看一下sobEDA原文献里切断共价键的sobEDA算例。
说实话，对于这么强的相互作用，用sobEDAw的结果可能未必有意义

作者
Author: sobereva 时间: 2026-2-25 04:01
分清楚sobEDA和sobEDAw
后者仅限于分解弱相互作用。而当前你做的是sobEDA（sh脚本里可见sobEDAw=0），并非你说的sobEDAw。对于化学键强度的能量分解，当前结果并没什么必然问题。sobEDA原文里有诸多例子可以参考

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