Perturbed matrix method surfaces

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. 

Theory. Usually we do not solve the fundamental equations directly. We use a theory, for example, Har-tree-Fock theory [3], Moller-Plesset perturbation theory [4], coupled-cluster theory [5], Kohn s [6, 7], Newton s [8], or Schlessinger s [9] variational principle for scattering amplitudes, the quasiclassical trajectory method [10], the trajectory surface hopping method [11], classical S-matrix theory [12], the close-coupling approximation... 

The parameterization of the exo-anomeric terms was based on the difference between the energy surfaces for 2-methoxytetrahydropyran (2) and 2-ethyltetrahydropyran (3) (see Figure 1), as calculated with the semi-empirical PCILO (perturbation of the configuration interaction matrix from localized bond orbitals) method (equations 7 and 8). ... 

---