Nuclear-electronic orbital nonorthogonal configuration interaction approach

The nuclear-electronic orbital nonorthogonal configuration interaction (NEO-NOCI) approach is presented. In this framework, the hydrogen nuclei are treated quantum mechanically on the same level as the electrons, and a mixed nuclear-electronic time-independent Schrödinger equation is solved with molecular orbital techniques. For hydrogen transfer systems, the transferring hydrogen is represented by two basis function centers to allow delocalization of the nuclear wave function. In the two-state NEO-NOCI approach, the ground and excited state delocalized nuclear-electronic wave functions are expressed as linear combinations of two nonorthogonal localized nuclear-electronic wave functions obtained at the NEO-Hartree-Fock level. The advantages of the NEO-NOCI approach are the removal of the adiabatic separation between the electrons and the quantum nuclei, the computational efficiency, the potential for systematic improvement by enhancing the basis sets and number of configurations, and the applicability to a broad range of chemical systems. The tunneling splitting is determined by the energy difference between the two delocalized vibronic states. The hydrogen tunneling splittings calculated with the NEO-NOCI approach for the [He-H-He]+ model system with a range of fixed He-He distances are in excellent agreement with NEO-full CI and Fourier grid calculations. These benchmarking calculations indicate that NEO-NOCI is a promising approach for the calculation of delocalized, bilobal hydrogen wave functions and the corresponding hydrogen tunneling splittings.