三分子中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)基于分子力场的能量贡献分析，绘制了能量贡献图。但目前我还不太清楚应该如何从这些能量贡献去解释体系之间的差异。老师是否有比较经典的参考文献或相关文章推荐?