To compute the vdW correction over DF energies, each fragment is replaced by an anisotropic QDO. The dipole-coupled dispersion energies between these oscillators (E ij QDO (r ij ) between i th and j th fragments) have been used as the vdW correction energy. The total dispersion energy is obtained by Eq.( ),