计算化学公社

标题: 关于使用不同软件/方法进行EDA计算的结果问题 [打印本页]

作者
Author: wangzhe 时间: yesterday 08:54
标题: 关于使用不同软件/方法进行EDA计算的结果问题
各位老师好，请教关于使用Q-Chem (ALMO-EDA), sobEDA以及 Psi4的SAPT计算时结果的差异问题。
参考了某篇论文中的例子，论文中使用了Q-Chem在M06/6-311+G(d,p)方法下进行了ALMO-EDA计算，算出来的E(elst) = -3.6 kcal/mol
本人用Q-Chem在同样的方法得出来的结果也是E(elst) = -3.6 kcal/mol，用Psi4的SAPT0/jun-cc-pvdz方法下算出来的E(elst) = -4.6 kcal/mol
另外用sob老师的sobEDA在M06/6-311+G(d,p)方法下计算发现E(elst) = -88.0 kcal/mol，但是total interaction energy和Q-Chem的ALMO-EDA结果一致，都为-5.2 kcal/mol
请问sobEDA的E(elst)项和其他方法有什么区别吗？

作者
Author: sobereva 时间: yesterday 09:43
我不知道你算的什么体系，没法评论
对于化学键，那种大小E(elst)很正常，对于弱相互作用则不正常（可能操作不当）
对于弱相互作用，一般情况强烈建议使用B3LYP-D3(BJ)做sobEDAw能量分解，而不要用明尼苏达系列泛函，本身sobEDAw也没有给明尼苏达泛函拟合的参数。
sobEDA/sobEDAw原文里https://doi.org/10.1021/acs.jpca.3c04374给了弱相互作用分解的例子，可以参考其结果，可见sobEDAw和高阶SAPT在各项上都吻合很好

作者
Author: wangzhe 时间: yesterday 17:43

sobereva 发表于 2026-5-22 10:43
我不知道你算的什么体系，没法评论
对于化学键，那种大小E(elst)很正常，对于弱相互作用则不正常（可能操 ...

感谢老师回复，我再检查一下

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