TY - JOUR

T1 - Alternative wavefunction ansatz for including explicit electron-proton correlation in the nuclear-electronic orbital approach

AU - Ko, Chaehyuk

AU - Pak, Michael V.

AU - Swalina, Chet

AU - Hammes-Schiffer, Sharon

N1 - Funding Information:
We thank Dr. Ari Chakraborty for helpful input in the early stages of this project and Andrew Sirjoosingh for helpful comments on the manuscript. We gratefully acknowledge funding from AFOSR Grant No. FA9550-10-1-0081.

PY - 2011/8/7

Y1 - 2011/8/7

N2 - The nuclear-electronic orbital (NEO) approach treats specified nuclei quantum mechanically on the same level as the electrons with molecular orbital techniques. The explicitly correlated Hartree-Fock (NEO-XCHF) approach was developed to incorporate electron-nucleus dynamical correlation directly into the variational optimization of the nuclear-electronic wavefunction. In the original version of this approach, the Hartree-Fock wavefunction is multiplied by (1+G), where G is a geminal operator expressed as a sum of Gaussian type geminal functions that depend on the electron-proton distance. Herein, a new wavefunction ansatz is proposed to avoid the computation of five- and six-particle integrals and to simplify the computation of the lower dimensional integrals involving the geminal functions. In the new ansatz, denoted NEO-XCHF2, the Hartree-Fock wavefunction is multiplied by 1+G rather than (1+G). Although the NEO-XCHF2 ansatz eliminates the integrals that are quadratic in the geminal functions, it introduces terms in the kinetic energy integrals with no known analytical solution. A truncated expansion scheme is devised to approximate these problematic terms. An alternative hybrid approach, in which the kinetic energy terms are calculated with the original NEO-XCHF ansatz and the potential energy terms are calculated with the NEO-XCHF2 ansatz, is also implemented. Applications to a series of model systems with up to four electrons provide validation for the NEO-XCHF2 approach and the treatments of the kinetic energy terms.

AB - The nuclear-electronic orbital (NEO) approach treats specified nuclei quantum mechanically on the same level as the electrons with molecular orbital techniques. The explicitly correlated Hartree-Fock (NEO-XCHF) approach was developed to incorporate electron-nucleus dynamical correlation directly into the variational optimization of the nuclear-electronic wavefunction. In the original version of this approach, the Hartree-Fock wavefunction is multiplied by (1+G), where G is a geminal operator expressed as a sum of Gaussian type geminal functions that depend on the electron-proton distance. Herein, a new wavefunction ansatz is proposed to avoid the computation of five- and six-particle integrals and to simplify the computation of the lower dimensional integrals involving the geminal functions. In the new ansatz, denoted NEO-XCHF2, the Hartree-Fock wavefunction is multiplied by 1+G rather than (1+G). Although the NEO-XCHF2 ansatz eliminates the integrals that are quadratic in the geminal functions, it introduces terms in the kinetic energy integrals with no known analytical solution. A truncated expansion scheme is devised to approximate these problematic terms. An alternative hybrid approach, in which the kinetic energy terms are calculated with the original NEO-XCHF ansatz and the potential energy terms are calculated with the NEO-XCHF2 ansatz, is also implemented. Applications to a series of model systems with up to four electrons provide validation for the NEO-XCHF2 approach and the treatments of the kinetic energy terms.

UR - http://www.scopus.com/inward/record.url?scp=80051702361&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80051702361&partnerID=8YFLogxK

U2 - 10.1063/1.3611054

DO - 10.1063/1.3611054

M3 - Article

C2 - 21823689

AN - SCOPUS:80051702361

VL - 135

JO - Journal of Chemical Physics

JF - Journal of Chemical Physics

SN - 0021-9606

IS - 5

M1 - 054106

ER -