Born-Oppenheimer approximation Gaussians

At this stage we are at the very beginning of development, implementation, and application of methods for quantum-mechanical calculations of molecular systems without assuming the Born-Oppenheimer approximation. So far we have done several calculations of ground and excited states of small diatomic molecules, extending them beyond two-electron systems and some preliminary calculations on triatomic systems. In the non-BO works, we have used three different correlated Gaussian basis sets. The simplest one without r,y premultipliers (4)j = exp[—r (A t (8> Is) "]) was used in atomic calculations the basis with premultipliers in the form of powers of rj exp[—r (Aj (8> /sjr])... 

The geminal ansatz still requires more effort than the standard one-electron approach of the independent particle model. It is therefore usually restricted to small molecules for feasibility reasons. As an example how the nonlinear optimization problem can be handled we refer to the stochastic variational approach [340]. However, the geminal ansatz as presented above has the useful feature that all elementary particles can be treated on the same footing. This means that we can actually use such an ansatz for total wave functions without employing the Born-Oppenheimer approximation, which exploits the fact that nuclei are much heavier than electrons. Hence, electrons and nuclei can be treated on the same footing [340-342] and even mixed approaches are possible, where protons and electrons are treated in the external field of heavier nuclei [343-346]. The integrals required for the matrix elements are hardly more complicated than those over one-electron Gaussians [338,339,347]. 

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. 

---