Search

Multicomponent density functional theory (DFT) has many practical advantages for incorporating nuclear quantum effects into quantum chemistry calculations. Within the nuclear-electronic orbital (NEO) framework, specified nuclei, typically protons, are treated quantum mechanically on the same level as the electrons. Previously, electron-proton correlation functionals based on the local density approximation (LDA), denoted epc17 and epc18, were developed and shown to provide more accurate proton densities and energies compared to the neglect of electron-proton correlation, but a quantitatively accurate description of both densities and energies simultaneously has remained elusive. Herein, an electron-proton correlation functional that depends on the electron and proton density gradients, as well as the densities, is derived and implemented. Compared to the LDA functionals, the resulting generalized gradient approximation functional, denoted epc19, is able to simultaneously provide accurate proton densities and energies, as well as reproduce the impact of nuclear quantum effects on optimized geometries. In addition, without further parameterization, the NEO-DFT/epc19 method provides accurate densities and energies for deuterium as well as hydrogen. These results demonstrate that the form of the epc19 functional is able to capture the essential aspects of electron-proton correlation and highlights the importance of including gradient terms. This approach will enable the exploration of nuclear quantum effects and isotope effects in a wide range of systems. [ABSTRACT FROM AUTHOR]

Multicomponent quantum chemistry allows the quantum mechanical treatment of electrons and specified protons on the same level. Typically the goal is to identify a self-consistent-field (SCF) solution that is the global minimum associated with the molecular orbital coefficients of the underlying Hartree-Fock (HF) or density functional theory (DFT) calculation. To determine whether the solution is a minimum or a saddle point, herein we derive the stability conditions for multicomponent HF and DFT in the nuclear-electronic orbital (NEO) framework. The gradient is always zero for an SCF solution, whereas the Hessian must be positive semi-definite for the solution to be a minimum rather than a saddle point. The stability matrices for NEO-HF and NEO-DFT have the same matrix structures, which are identical to the working matrices of their corresponding linear response time-dependent theories (NEO-TDHF and NEO-TDDFT) but with a different metric. A negative eigenvalue of the stability matrix is a necessary but not sufficient condition for the corresponding NEO-TDHF or NEO-TDDFT working equation to have an imaginary eigenvalue solution. Electron-proton systems could potentially exhibit three types of instabilities: electronic, protonic, and electron-proton vibronic instabilities. The internal and external stabilities for theories with different constraints on the spin and spatial orbitals can be analyzed. This stability analysis is a useful tool for characterizing SCF solutions and is helpful when searching for lower-energy solutions. Initial applications to HCN, HNC, and 2-cyanomalonaldehyde, in conjunction with NEO ∆SCF calculations, highlight possible connections between stationary points in nuclear coordinate space for conventional electronic structure calculations and stationary points in orbital space for NEO calculations. [ABSTRACT FROM AUTHOR]