相互作用能（Interaction Energy） 本质上是两个（或多个）分子从无限远到结合状态时能量的变化量，用以衡量分子间吸引或排斥的程度。
在理论化学中，它通常表示为：

若结果为负值，说明体系稳定化（结合形成），为正值则说明排斥。

从能量构成角度看，相互作用能并非单一来源，而是多种量子效应的综合结果：

• 静电作用（Electrostatics）：由分子电荷分布决定，是相互作用的长程主导力。

• 交换–排斥（Exchange–Repulsion）：源于泡利不相容原理，是短程排斥项。

• 极化与轨道相互作用（Polarization / Orbital Interaction）：反映电子云重叠与电荷转移，是共价性成分的体现。

• 色散作用（Dispersion）：源于瞬时偶极–诱导偶极，是范德华吸引的量子基础。

在弱相互作用体系中，静电与色散贡献常相当；而在配位或氢键主导的体系中，轨道极化占据更高比例。

为了理解这些相互作用的“能量组成”，研究者通常采用 能量分解分析（Energy Decomposition Analysis, EDA），
将总相互作用能分拆为以下几项：

这种分解可揭示结合的主导能量来源：

• 若静电项占主导，则是极性吸引型结合（如氢键）；

• 若色散贡献突出，说明范德华或π–π 相互作用为主；

• 若极化项明显，则体系具有部分共价性或电荷转移特征。

相互作用不仅可由能量表示，也可通过电子密度拓扑直接可视化。
在 Multiwfn 或 TopMod 等程序中，可以基于 Gaussian 输出波函数文件进行以下分析：

✴️ AIM（Atoms in Molecules）分析

由 Bader 提出，通过识别电子密度 ρ(r) 的键临界点 (BCP) 判断是否存在相互作用：

• 若 ρ(r) 高且 ∇²ρ(r) < 0 → 共价键性质；

• 若 ρ(r) 较低且 ∇²ρ(r) > 0 → 弱相互作用（如氢键、范德华力）。

AIM 可提供定量指标：ρ(BCP)、H(r)、∇²ρ(r)，用于区分键强与性质。

✴️ NCI 与 IGMH 分析

利用sign(λ₂)ρ与RDG构建三维等值面，
可以直接显示氢键、范德华与π–π 作用区域。

• 蓝色区域代表吸引作用（氢键、静电）；

• 绿色区域对应范德华作用；

• 红色区域为排斥区域。

这些分析常与能量项结合，用于解释配体–底物复合物或催化过渡态的稳定来源。

要真正“看到”这些相互作用，需要三者协作：

结合流程为：
1️⃣ Gaussian 优化结构与计算能量；
2️⃣ 将 .fchk 文件导入 Multiwfn，进行 AIM / NCI / EDA 分析；
3️⃣ 结果输出 .cube 文件导入 VMD，生成电子密度或相互作用等值面图。