Impact of nuclear quantum effects on the molecular structure of bihalides and the hydrogen fluoride dimer

The structural impact of nuclear quantum effects is investigated for a set of bihalides, [XHX](-), X = F, Cl, and Br, and the hydrogen fluoride dimer. Structures are calculated with the vibrational self-consistent-field (VSCF) method, the second-order vibrational perturbation theory method (VPT2), and the nuclear-electronic orbital (NEO) approach. In the VSCF and VPT2 methods, the vibrationally averaged geometries are calculated for the Born-Oppenheimer electronic potential energy surface. In the NEO approach, the hydrogen nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wave functions are calculated variationally with molecular orbital methods. Electron-electron and electron-proton dynamical correlation effects are included in the NEO approach using second-order perturbation theory (NEO-MP2). The nuclear quantum effects are found to alter the distances between the heavy atoms by 0.02-0.05 A for the systems studied. These effects are of similar magnitude as the electron correlation effects. For the bihalides, inclusion of the nuclear quantum effects with the NEO-MP2 or the VSCF method increases the X-X distance. The bihalide X-X distances are similar for both methods and are consistent with two-dimensional grid calculations and experimental values, thereby validating the use of the computationally efficient NEO-MP2 method for these types of systems. For the hydrogen fluoride dimer, inclusion of nuclear quantum effects decreases the F-F distance with the NEO-MP2 method and increases the F-F distance with the VSCF and VPT2 methods. The VPT2 F-F distances for the hydrogen fluoride dimer and the deuterated form are consistent with the experimentally determined values. The NEO-MP2 F-F distance is in excellent agreement with the distance obtained experimentally for a model that removes the large amplitude bending motions. The analysis of these calculations provides insight into the significance of electron-electron and electron-proton correlation, anharmonicity of the vibrational modes, and nonadiabatic effects for hydrogen-bonded systems.     

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

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.     

Analysis of electron-positron wavefunctions in the nuclear-electronic orbital framework

The nuclear-electronic orbital explicitly correlated Hartree-Fock (NEO-XCHF) approach is extended and applied to the positronic systems PsH, LiPs, and e(+)LiH. In this implementation, all electrons and positrons are treated quantum mechanically, and all nuclei are treated classically. This approach utilizes molecular orbital techniques with Gaussian basis sets for the electrons and positrons and includes electron-positron correlation with explicitly correlated Gaussian-type geminal functions. An efficient strategy is developed to reduce the number of variational parameters in the NEO-XCHF calculations. The annihilation rates, electron and positron densities, and electron-positron contact densities are compared to available results from higher-level calculations. Our analysis illustrates that the NEO-XCHF method produces qualitative to semi-quantitative results for these properties at a relatively low computational cost by treating only the essential electron-positron correlation explicitly. The NEO-HF method, which does not include explicit correlation and therefore is extremely efficient, is found to provide qualitatively accurate electron-positron contact densities for the e(+)LiH system but not for the LiPs system. Thus, the utility of the NEO-HF method for determining where annihilation occurs is system dependent and not generally reliable. The NEO-XCHF method, however, provides a computationally practical and reliable approach for determining where annihilation will occur in positronic systems.     

Density functional theory treatment of electron correlation in the nuclear-electronic orbital approach

This paper presents the nuclear-electronic orbital density functional theory [NEO-DFT(ee)] method for including electron-electron correlation and nuclear quantum effects self-consistently in quantum chemical calculations. The NEO approach is designed to treat a relatively small number of nuclei quantum mechanically, while the remaining nuclei are treated classically. In the NEO-DFT(ee) approach, the correlated electron density is used to obtain the nuclear molecular orbitals, and the resulting nuclear density is used to obtain the correlated electron density during an iterative procedure that continues until convergence of both the nuclear and electronic densities. This approach includes feedback between the correlated electron density and the nuclear wavefunction. The application of this approach to bihalides and acetylene indicates that the nuclear quantum effects do not significantly impact the electron correlation energy, but the quantum nuclear energy is enhanced in the NEO-DFT(ee) B3LYP method. The excellent agreement of the NEO-DFT(ee)-optimized bihalide structures with the vibrationally averaged geometries from grid-based quantum dynamical methods provides validation for the NEO-DFT(ee) approach. Electron-proton correlation could be included by the development of an electron-nucleus correlation functional. Alternatively, explicit electron-proton correlation could be included directly into the NEO self-consistent-field framework with Gaussian-type geminal functions.     

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 Schrodinger 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.     

Multicomponent density functional theory study of the interplay between electron-electron and electron-proton correlation

The interplay between electron-electron and electron-proton correlation is investigated within the framework of the nuclear-electronic orbital density functional theory (NEO-DFT) approach, which treats electrons and select protons quantum mechanically on the same level. Recently two electron-proton correlation functionals were developed from the electron-proton pair densities obtained from explicitly correlated wavefunctions. In these previous derivations, the kinetic energy contribution arising from electron-proton correlation was neglected. In this paper, an electron-proton correlation functional that includes this kinetic energy contribution is derived using the adiabatic connection formula in multicomponent DFT. The performance of the NEO-DFT approach using all three electron-proton correlation functionals in conjunction with three well-established electronic exchange-correlation functionals is assessed. NEO-DFT calculations with these electron-proton correlation functionals capture the increase in the hydrogen vibrational stretching frequencies arising from the inclusion of electron-electron correlation in model systems. Electron-proton and electron-electron correlation are found to be uncoupled and predominantly additive effects to the total energy for the model systems studied. Thus, electron-proton correlation functionals and electronic exchange-correlation functionals can be developed independently and subsequently combined together without re-parameterization.     

Analysis of the nuclear-electronic orbital method for model hydrogen transfer systems

Fundamental issues associated with the application of the nuclear-electronic orbital (NEO) approach to hydrogen transfer systems are addressed. In the NEO approach, specified nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wavefunctions are calculated with molecular orbital methods. The positions of the nuclear basis function centers are optimized variationally. In the application of the NEO approach to hydrogen transfer systems, the hydrogen nuclei and all electrons are treated quantum mechanically. Within the NEO framework, the transferring hydrogen atom can be represented by two basis function centers to allow delocalization of the proton vibrational wavefunction. In this paper, the NEO approach is applied to the [He-H-He]+ and [He-H-He]++ model systems. Analyses of technical issues pertaining to flexibility of the basis set to describe both single and double well proton potential energy surfaces, linear dependency of the hydrogen basis functions, multiple minima in the basis function center optimization, convergence of the number of hydrogen basis function centers, and basis set superposition error are presented. The accuracy of the NEO approach is tested by comparison to grid calculations for these model systems.     

Multi-component molecular orbital theory for electrons and nuclei including many-body effect with full configuration interaction treatment: isotope effects on hydrogen molecules

A multi-component molecular orbital (MC_MO) theory is developed for a combined quantum system of electrons and nuclei with the full configuration interaction (CI) scheme of Cartesian Gaussian-type functions. The technique of graphical unitary group approach (GUGA) is modified to obtain the CI matrix elements for many kinds of quantum particles efficiently. The optimum basis sets for quantum nuclei are proposed with the fully variational procedure. The average internuclear distances, dipole polarizabilities, and nuclear vibrational excitation energies for isotopic hydrogen molecules calculated with optimized basis sets are found to adequately reproduce the corresponding experimental values.     

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.     

Quantum chemistry beyond Born–Oppenheimer approximation on a quantum computer: A simulated phase estimation study

We present an efficient quantum algorithm for beyond‐Born–Oppenheimer molecular energy computations. Our approach combines the quantum full configuration interaction method with the nuclear orbital plus molecular orbital method. We give the details of the algorithm and demonstrate its performance by classical simulations. Two isotopomers of the hydrogen molecule (H₂, HT) were chosen as representative examples and calculations of the lowest rotationless vibrational transition energies were simulated. © 2016 Wiley Periodicals, Inc.     

On the extent of intramolecular hydrogen bonding in gas-phase and hydrated 1,2-ethanediol

We investigate the quantum dynamical nature of hydrogen bonding in 1,2-ethanediol and monohydrated 1,2-ethanediol using different levels of ab initio theory. Global full-dimensional potential energy surfaces were constructed from PW91/cc-pVDZ, B3LYP/cc-pVDZ, and MP2/cc-pVDZ ab initio data for gas-phase and monohydrated 1,2-ethanediol, using a modified Shepard interpolation scheme. Zero-point energies and nuclear vibrational wave functions were calculated on these surfaces using the quantum diffusion Monte Carlo algorithm. The nature of intra- and intermolecular hydrogen bonding in these molecules was investigated by considering a ground-state nuclear vibrational wavefunction with reduced complete nuclear permutation and inversion (CNPI) symmetry. Separate wavefunction histograms were determined from the ground-state nuclear vibrational wavefunction by projection into bondlength coordinates. The O-H and O-O wavefunction histograms and vibrationally averaged distances were then used to probe the extent of intra- and intermolecular hydrogen bonding. From these data, we conclude that gas-phase ethanediol may possess a weak hydrogen bond, with a relatively short O-O distance but no detectable proton delocalization. Monohydrated ethanediol was found to exhibit no intramolecular hydrogen bonding but instead possessed two intermolecular hydrogen bonds, indicated by both shortening of the O-O distance and significant proton delocalization. The degree of proton delocalization and shortening of the vibrationally averaged O-O distance was found to be dependent on the ab initio method used to generate the potential energy surface (PES) data set.     

Modeling positrons in molecular electronic structure calculations with the nuclear-electronic orbital method

The nuclear-electronic orbital (NEO) method was modified and extended to positron systems for studying mixed positronic-electronic wavefunctions, replacing the mass of the proton with the mass of the positron. Within the modified NEO framework, the NEO-HF (Hartree-Fock) method provides the energy corresponding to the single-configuration mixed positronic-electronic wavefunction, minimized with respect to the molecular orbitals expressed as linear combinations of Gaussian basis functions. The electron-electron and electron-positron correlation can be treated in the NEO framework with second-order perturbation theory (NEO-MP2) or multiconfigurational methods such as the full configuration interaction (NEO-FCI) and complete active space self-consistent-field (NEO-CASSCF) methods. In addition to implementing these methods for positronic systems, strategies for calculating electron-positron annihilation rates using NEO-HF, NEO-MP2, and NEO-FCI wavefunctions were also developed. To apply the NEO method to the positronium hydride (PsH) system, positronic and electronic basis sets were optimized at the NEO-FCI level and used to compute NEO-MP2 and NEO-FCI energies and annihilation rates. The effects of basis set size on NEO-MP2 and NEO-FCI correlation energies and annihilation rates were compared. Even-tempered electronic and positronic basis sets were also optimized for the e+LiH molecule at the NEO-MP2 level and used to compute the equilibrium bond length and vibrational energy.     

Simultaneous determination of nuclear and electronic wave functions without Born–Oppenheimer approximation: Ab initio NO+MO/HF theory

We develop a simultaneous determination method of nuclear and electronic wave functions without the Born–Oppenheimer approximation. We examine two expanding methods, namely, molecular orbital (MO)-type and valence bond (VB)-type expansions for a nuclear orbital, which is a one-particle wave function of a nucleus. The VB-type expansion is shown to be more accurate than the MO-type one because of the local nature of the nuclei. We also investigate the basis function expansion of the nuclear orbital and propose a scheme to determine the orbital exponent for the nuclear basis function. Numerical calculations confirm the accuracy and feasibility of the present method. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Atoms in molecules: beyond Born–Oppenheimer paradigm

This contribution presents the first atoms in molecules study that goes beyond the Born–Oppenheimer paradigm employing the newly developed two-component quantum theory of atoms in molecules (TC-QTAIM). The LiH, LiD, and LiT systems containing quantum instead of clamped hydrogen nuclei are used as typical examples. The computational analysis that is done on non-adiabatic wavefunctions derived from the fully variational multicomponent molecular orbital approach (FV-MC-MO) results in hydrogen atomic basins without any clamped nucleus. The topological analysis of the Γ-field, the field that replaces the usual one-electron density used in the orthodox topological analysis, reveals delicate differences among the considered systems. The calculation of basin properties also demonstrates that the TC-QTAIM differentiates among atomic basins containing isotopes. Since the nuclear dynamics is contained intrinsically in non-adiabatic wavefunctions, the nuclear contribution to both topological analysis and basin properties naturally emerges from the TC-QTAIM analysis resolving the long-standing obstacle of consistent incorporation of nuclear dynamics within the context of the orthodox QTAIM. Also, a similar analysis is done on non-adiabatic wavefunctions describing excited instead of ground nuclear vibrations of the considered systems demonstrating the fact that TC-QTAIM is capable of being employed for both ground and excited nuclear vibrational states.     

Quantum confinement in hydrogen bond

In this work, the quantum confinement effect is proposed as the cause of the displacement of the vibrational spectrum of molecular groups that involve hydrogen bonds. In this approach, the hydrogen bond imposes a space barrier to hydrogen and constrains its oscillatory motion. We studied the vibrational transitions through the Morse potential, for the NH and OH molecular groups inside macromolecules in situation of confinement (when hydrogen bonding is formed) and nonconfinement (when there is no hydrogen bonding). The energies were obtained through the variational method with the trial wave functions obtained from supersymmetric quantum mechanics formalism. The results indicate that it is possible to distinguish the emission peaks related to the existence of the hydrogen bonds. These analytical results were satisfactorily compared with experimental results obtained from infrared spectroscopy. © 2015 Wiley Periodicals, Inc.     

A second quantization formulation of multimode dynamics

A new formalism for calculating and analyzing many-mode quantum dynamics is presented. The formalism is similar in spirit to the second quantization formulation of electronic structure theory. The similarity means that similar techniques can be employed for calculating the many-mode nuclear wave function. As a consequence a new formulation of the vibrational self-consistent-field (VSCF) method can be developed. Another result is that the formalism opens up for the construction of new methods that go beyond the VSCF level. A vibrational coupled cluster (VCC) theory is constructed using the new formalism. The size-extensivity concept is introduced in the context of multimode wave functions and the size extensivity of approximate VCC methods is illustrated in comparison with the non-size-extensive vibrational configuration interaction method.     

A generalized any particle propagator theory: Assessment of nuclear quantum effects on electron propagator calculations

In this work we propose an extended propagator theory for electrons and other types of quantum particles. This new approach has been implemented in the LOWDIN package and applied to sample calculations of atomic and small molecular systems to determine its accuracy and performance. As a first application of the method we have studied the nuclear quantum effects on electron ionization energies. We have observed that ionization energies of atoms are similar to those obtained with the electron propagator approach. However, for molecular systems containing hydrogen atoms there are improvements in the quality of the results with the inclusion of nuclear quantum effects. An energy term analysis has allowed us to conclude that nuclear quantum effects are important for zero order energies whereas propagator results correct the electron and electron-nuclear correlation terms. Results presented for a series of n-alkanes have revealed the potential of this method for the accurate calculation of ionization energies of a wide variety of molecular systems containing hydrogen nuclei. The proposed methodology will also be applicable to exotic molecular systems containing positrons or muons.