A straightforward approach for evaluating the interaction energy is to calculate the difference between the objects' combined energy and all of their isolated energies. In the case of two objects, A and B, the interaction energy can be written as: 1 Δ E int = E ( A , B ) − ( E ( A ) + E ( B ) ) , {\displaystyle \Delta E_{\text{int}}=E(A,B)-\left(E(A)+E(B)\right),} where E ( A ) {\displaystyle E(A)} and E ( B ) {\displaystyle E(B)} are the energies of the isolated objects (monomers), and E ( A , B ) {\displaystyle E(A,B)} the energy of their interacting assembly (dimer).
For larger system, consisting of N objects, this procedure can be generalized to provide a total many-body interaction energy:
Δ E int = E ( A 1 , A 2 , … , A N ) − ∑ i = 1 N E ( A i ) . {\displaystyle \Delta E_{\text{int}}=E(A_{1},A_{2},\dots ,A_{N})-\sum _{i=1}^{N}E(A_{i}).}
By calculating the energies for monomers, dimers, trimers, etc., in an N-object system, a complete set of two-, three-, and up to N-body interaction energies can be derived.
The supermolecular approach has an important disadvantage in that the final interaction energy is usually much smaller than the total energies from which it is calculated, and therefore contains a much larger relative uncertainty. In the case where energies are derived from quantum chemical calculations using finite atom-centered basis functions, basis set superposition errors can also contribute some degree of artificial stabilization.
Theoretical and Computational Chemistry, 1999, Ideas of Quantum Chemistry, 2007 and Quantum Magnetic Resonance Imaging Diagnostics of Human Brain Disorders, 2010 ↩