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.     

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.     

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.     

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.     

Colle-Salvetti-type correction for electron-nucleus correlation in the nuclear orbital plus molecular orbital theory

A Colle-Salvetti (CS)-type electron-nucleus correction in the nuclear orbital plus molecular orbital theory is proposed. The CS-type correction is designed to satisfy the cusp condition for the electron-nucleus interaction. Since the CS-type correction is expressed in terms of the electron and nucleus densities, its evaluation is computationally feasible. Numerical assessment confirms that the CS-type correction performs well for the small G2 set.     

The dynamical correlation in spacer-mediated electron transfer couplings

The dynamical correlation effect in electron transfer (ET) coupling was studied in this work, for cases where electrons tunnel through a many-electron environment. The ET couplings for three different bridge-mediated model systems were calculated: (I) trans-alkyl chains [H2C-(CH2)n-CH2, n = 2-10], (II) two isomers of trans-1,4-dimethylenecyclohexane, and (III) two ethylenes spaced by a saturated ethane molecule. The couplings were calculated as half energy gaps of the two lowest adiabatic states. The dynamical correlation was included with spin-flip (SF) and ionization potential or electron affinity coupled-cluster singles and doubles (SF-CCSD and IP/EA-CCSD) and a DeltaCCSD scheme. The direct coupling (DC) scheme is also used as a way to obtain a solution with nondynamical correlation, since DC uses approximated eigenstates that are symmetry-restoring linear combinations of two symmetry-broken unrestricted Hartree-Fock configurations. For all cases tested except for one, results from the DC scheme closely follow the CCSD data, indicating that the dual-configuration solutions can be a good approximation of wave functions with nondynamical correlation included, but there exist exceptions. Comparing the DC results with SF-CCSD and IP or EA-CCSD data, we concluded that the dynamical correlation effect is small for most of the cases we tested.     

Scaling dynamical correlation energy from density functional theory correlation functionals

The scaling of dynamical correlation energy in molecules obtained by the correlation functionals of density functional theory (DFT) is examined. The approach taken is very similar to the scaled external correlation method of Brown and Truhlar but is based on the observation that DFT correlation functionals, especially the LYP, appear to represent the dynamical portion of the correlation energy in molecules. We examine whether higher accuracy in atomization energies can be gained by scaling without significant deterioration of the structural and spectroscopic properties of the molecules using four DFT functionals (BLYP, OLYP, B3LYP, and O3LYP) on 19 molecules including the six molecule AE6 database, the latter being representative of a much larger, 109 molecule training set. We show that, with molecule specific scale factors, nearly perfect agreement with experiment can be achieved in atomization energies without increasing the average errors in other molecular properties relative to the DFT calculation. We further show that it is possible to find optimal scale factors which reduce the mean unsigned error per bond to levels comparable to those of some multilevel multicoefficient methods.     

New functionals for correlation energy deduced in the framework of the correlation factor approach

The general characteristics of two-body density functionals (TBDF) are explored and two new correlation energy functionals are derived using the correlation factor approach. The optimization of the parameters entering the above functionals requires exact and accurate atomic correlation energies (ACE). We revised the ACE values in the literature and obtained a new set of “exact” ACE for atoms with 2 ≤ Z ≤ 10. Unfortunately, there exist some inaccuracies in the ACE values of the second-row atoms, which make unsuitable the inclusion of them in the optimization. The ACE calculated for the first period with the above functionals, using the optimized sets of parameters, are in excellent agreement with the exact ones, while the corresponding values calculated for the second-row atoms are between the precision margins estimated by us for the exact values. © 1997 John Wiley & Sons, Inc.     

Dyson-corrected orbital energies for the perturbative treatment of electron correlation

The effect of replacing the Hartree–Fock one-particle energies with ionization potentials obtained from inverse Dyson equation when calculating electron correlation energies perturbatively is investigated. Though the energy shifts vary from system to system, the slight decrease of the resulting excitation energies at around equilibrium geometries leads to a slight increase of the correlation energies in most cases. In the dissociation limit the inverse Dyson equation opens the gap, thus nondiverging potential curves emerge even at the restricted Hartree–Fock (RHF)+RS2 level. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 713–719, 1998     

A density functional for the correlation energy,deduced in the framework of the correlation factor approach

An approximate density functional is deduced from a wave function within the correlation factor method. The new functional does not include terms depending on the gradient of the density, but shows the simplicity of local density functionals without spin polarization. However, it includes correctly the inhomogeneity effects and, also, the nonlocal nature of an electronic system. The approach adopted here stresses the goodness of the expression taken by Colle and Salvetti for building a correlation factor and, at the same time, allows us to gain light on the nature of the deficiencies of those functionals obtained, up to now, from the perspective of the Hohenberg and Kohn theorem.     

Electron correlation: Keeping close to an orbital description

This article deals with excited-state properties and spectroscopies which probe them. We argue that the orbitals obtained in local density theory form an excellent basis for discussing excited states and that the self-energy operator (or optical potential) associated with one-electron-type excited states is closely diagonal in such a basis. This diagonality often leads to considerable simplifications. We also discuss the reasons for the successes of the “GW -approximation” for the optical potential in predicting excited states. We further discuss a recent approach to calculate spectroscopic data. In particular, we present explicit results of this approach to photoelectron emission. The approach treats correlations within the target exactly but neglects some terms in the coupling between target and probe electron (= the photoemitted electron). This leads to closed expressions which exhibit a number of qualitative features. One is the cancellation of intrinsic and extrinsic losses at threshold excitation. Another is the broadening of the quasi-particles; we have energy-broadening for the occupied states but broadening through decaying wave functions for unoccupied states. © 1995 John Wiley & Sons, Inc.     

Domain-based local pair natural orbital methods within the correlation consistent composite approach

Ab initio composite approaches have been utilized to model and predict main group thermochemistry within 1 kcal mol⁻¹, on average, from well-established reliable experiments, primarily for molecules with less than 30 atoms. For molecules of increasing size and complexity, such as biomolecular complexes, composite methodologies have been limited in their application. Therefore, the domain-based local pair natural orbital (DLPNO) methods have been implemented within the correlation consistent composite approach (ccCA) framework, namely DLPNO-ccCA, to reduce the computational cost (disk space, CPU (central processing unit) time, memory) and predict energetic properties such as enthalpies of formation, noncovalent interactions, and conformation energies for organic biomolecular complexes including one of the largest molecules examined via composite strategies, within 1 kcal mol⁻¹, after calibration with 119 molecules and a set of linear alkanes. © 2019 Wiley Periodicals, Inc.     

Electron correlation effects on the shape of the Kohn–Sham molecular orbital

Even systems in which strong electron correlation effects are present, such as the large near-degeneracy correlation in a dissociating electron pair bond exemplified by stretched H₂, are represented in the Kohn–Sham (KS) model of non-interacting electrons by a determinantal wavefunction built from the KS molecular orbitals. As a contribution to the discussion on the status and meaning of the KS orbitals we investigate, for the prototype system of H₂ at large bond distance, and also for a one-dimensional molecular model, how the electron correlation effects show up in the shape of the KS σ _g orbital. KS orbitals φ^HL and φ^FCI obtained from the correlated Heitler-London and full configuration interaction wavefunctions are compared to the orbital φ^LCAO, the traditional linear combination of atomic orbitals (LCAO) form of the (approximate) Hartree-Fock orbital. Electron correlation manifests itself in an essentially non-LCAO structure of the KS orbitals φ^HL and φ^FCI around the bond midpoint, which shows up particularly clearly in the Laplacian of the KS orbital. There are corresponding features in the kinetic energy density t _s of the KS system (a well around the bond midpoint) and in the one-electron KS potential v _s (a peak). The KS features are lacking in the Hartree-Fock orbital, in a minimal LCAO approximation as well as in the exact one. Received: 11 December 1996 / Accepted: 10 January 1997     

Tempered orbital energies in SCF MO calculations and their relation to the ordinate in Mulliken-Walsh correlation diagrams and extended Hückel orbital energies

An “average state” of a molecule is defined by distributing the electrons equally among the valence orbitals of a minimal basis set Hartree-Fock calculation. The resulting eigenvalues, called tempered orbital energies, behave much more like the Mulliken-Walsh diagram energies or extended Hückel eigenvalues than do the Hartree-Fock canonical orbital energies.