三分子中SAPT与EDA-FF能量分解差异及结果问题

Orange13

我构建了一个三明治型三体体系（ABC），其中A为苯酚，B为苯，C为不同含氮杂环化合物。所有体系的AB坐标相同，C在初始结构基础上进行替换，采用B3LYP-D3（BJ）/6-31G*方法并对重原子固定的约束优化，随后在优化结构上进行用SAPT进行二体SAPT2+（3）δMP2/jun-cc-pVDZ能量分解分析，以及三体MP2/jun-cc-pVDZ+VMFC校正计算。还用同样方法进行了EDA-FF能量分解方法进行能量分析。但在结果解释上遇到了一些困惑，想请教大家。目前得到的结果如下：
SAPT 能量分解：第二组中 2350 的相互作用能大于 1505，而第一组和第三组均表现为 1505 大于 2350。EDA-FF 能量分解：第二组和第三组表现为 2350 大于 1505，仅第一组表现为 1505 大于 2350。
针对这些结果，我主要有以下几个问题：
(1)为什么在SAPT能量分解里只有第二组表现为 2350 的相互作用能高于 1505，而第一、三组却相反？有什么方法去解释这个结果吗？
(2)SAPT与EDA-FF的能量分解趋势为什么不一样，是因为EDA-FF里无极化所以没有SAPT准吗？
(3)基于分子力场的能量贡献分析，绘制了能量贡献图。但目前我还不太清楚应该如何从这些能量贡献去解释体系之间的差异。老师是否有比较经典的参考文献或相关文章推荐?

sobereva

基于力场的EDA-FF原理上就不可能比诸如下文说的sobEDAw这样基于量子化学的严格的能量分解方法准，是否考虑可极化效应对于研究当前的弱相互作用都是相对次要的
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）
而且当前EDA-FF得到的色散作用项那么负，总相互作用能那么负，明显异常，大概率是用的方式不对。

若要了解不同复合物的相互作用能的差异本质，参考下文的讨论
谈谈分子间结合能的构成以及分解分析思想
http://sobereva.com/733（http://bbs.keinsci.com/thread-50890-1-1.html）

Orange13

sobereva 发表于 2026-7-14 09:02
基于力场的EDA-FF原理上就不可能比诸如下文说的sobEDAw这样基于量子化学的严格的能量分解方法准，是否考虑 ...

谢谢老师的解答，我想再请教您：那对于三体体系的 SAPT 能量分解分析，SAPT 是否足以解释该三体体系的相互作用差异？如果SAPT能量分解能解释我这体系，三组趋势不同（第二组2350相互作用能大于1505），在几何构型基本一致的前提下，我还不清楚怎么样解释各项贡献中趋势差异。请问我后续需要做那些分析吗？