Potential energy surfaces, calculation perturbation methods

In Chapter IX, Liang et al. present an approach, termed as the crude Bom-Oppenheimer approximation, which is based on the Born-Oppen-heimer approximation but employs the straightforward perturbation method. Within their chapter they develop this approximation to become a practical method for computing potential energy surfaces. They show that to carry out different orders of perturbation, the ability to calculate the matrix elements of the derivatives of the Coulomb interaction with respect to nuclear coordinates is essential. For this purpose, they study a diatomic molecule, and by doing that demonstrate the basic skill to compute the relevant matrix elements for the Gaussian basis sets. Finally, they apply this approach to the H2 molecule and show that the calculated equilibrium position and foree constant fit reasonable well those obtained by other approaches. 

On the theoretical side the H20-He systems has a sufficiently small number of electrons to be tackled by the most sophisticated quantum-chemical techniques, and in the last two decades several calculations by various methods of electronic structure theory have been attempted [77-80]. More recently, new sophisticated calculations appeared in the literature they exploited combined symmetry - adapted perturbation theory SAPT and CCSD(T), purely ab initio SAPT [81,82], and valence bond methods [83]. A thorough comparison of the topology, the properties of the stationary points, and the anisotropy of potential energy surfaces obtained with coupled cluster, Moller-Plesset, and valence bond methods has been recently presented [83]. 

---