Alternative forms and transferability of electron-proton correlation functionals in nuclear-electronic orbital density functional theory

Capturing Correlation Effects in Positron Binding to Atoms and Molecules

A major challenge in contemporary electronic structure theory involves the development of methods to describe in a balanced manner the contribution of correlation effects to energy differences. This challenge can be even greater for multicomponent systems containing more than one type of quantum particle. In the present work, we describe a flexible code for carrying out self-consistent field and configuration interaction (CI) calculations on multicomponent systems and use it to generate trial wave functions for use in diffusion Monte Carlo (DMC) calculations of the positron affinity of Be, Be2, Be4, Mg, CS2, and benzene. The resulting positron affinities (PAs) are in good agreement with the best values from the literature.

Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory

Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.

Symmetry-Projected Nuclear-Electronic Hartree-Fock: Eliminating Rotational Energy Contamination

We present a symmetry projection technique for enforcing rotational and parity symmetries in nuclear-electronic Hartree-Fock wave functions, which treat electrons and nuclei on equal footing. The molecular Hamiltonian obeys rotational and parity inversion symmetries, which are, however, broken by expanding in Gaussian basis sets that are fixed in space. We generate a trial wave function with the correct symmetry properties by projecting the wave function onto representations of the three-dimensional rotation group, i.e., the special orthogonal group in three dimensions SO(3). As a consequence, the wave function becomes an eigenfunction of the angular momentum operator which (i) eliminates the contamination of the ground-state wave function by highly excited rotational states arising from the broken rotational symmetry and (ii) enables the targeting of specific rotational states of the molecule. We demonstrate the efficiency of the symmetry projection technique by calculating the energies of the low-lying rotational states of the H2 and H3+ molecules.

MC-QTAIM analysis reveals an exotic bond in coherently quantum superposed malonaldehyde

The proton between the two oxygen atoms of the malonaldehyde molecule experiences an effective double-well potential in which the proton's wavefunction is delocalized between the two wells. Herein we employ a state-of-the-art multi-component quantum theory of atoms in molecules partitioning scheme to obtain the molecular structure, i.e. atoms in molecules and bonding network, from the superposed ab initio wavefunctions of malonaldehyde. In contrast to the familiar clamped-proton portrayal of malonaldehyde, in which the proton forms a hydrogen basin, for the superposed states the hydrogen basin disappears and two novel hybrid oxygen-hydrogen basins appear instead, with an even distribution of the proton population between the two basins. The interaction between the hybrid basins is stabilizing thanks to an unprecedented mechanism. This involves the stabilizing classical Coulomb interaction of the one-proton density in one of the basins with one-electron density in the other basin. This stabilizing mechanism yields a bond foreign to the known bonding modes in chemistry.

Second-Order Self-Consistent Field Algorithms: From Classical to Quantum Nuclei

This work presents a general framework for deriving exact and approximate Newton self-consistent field (SCF) orbital optimization algorithms by leveraging concepts borrowed from differential geometry. Within this framework, we extend the augmented Roothaan-Hall (ARH) algorithm to unrestricted electronic and nuclear-electronic calculations. We demonstrate that ARH yields an excellent compromise between stability and computational cost for SCF problems that are hard to converge with conventional first-order optimization strategies. In the electronic case, we show that ARH overcomes the slow convergence of orbitals in strongly correlated molecules with the example of several iron-sulfur clusters. For nuclear-electronic calculations, ARH significantly enhances the convergence already for small molecules, as demonstrated for a series of protonated water clusters.

Hydrogen-Atom Electronic Basis Sets for Multicomponent Quantum Chemistry

Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.

Nuclear-Electronic Orbital Approach to Quantization of Protons in Periodic Electronic Structure Calculations

The nuclear-electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born-Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems such as two-dimensional materials and liquid-semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide-water interface. The zero-point energy effects of the protons, as well as electron-proton correlation, are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

Quantum Proton Effects from Density Matrix Renormalization Group Calculations

We recently introduced [J. Chem. Phys. 2020, 152, 204103] the nuclear-electronic all-particle density matrix renormalization group (NEAP-DMRG) method to solve the molecular Schrödinger equation, based on a stochastically optimized orbital basis, without invoking the Born-Oppenheimer approximation. In this work, we combine the DMRG method with the nuclear-electronic Hartree-Fock (NEHF-DMRG) approach, treating nuclei and electrons on the same footing. Inter- and intraspecies correlations are described within the DMRG method without truncating the excitation degree of the full configuration interaction wave function. We extend the concept of orbital entanglement and mutual information to nuclear-electronic wave functions and demonstrate that they are reliable metrics to detect strong correlation effects. We apply the NEHF-DMRG method to the HeHHe⁺ molecular ion, to obtain accurate proton densities, ground-state total energies, and vibrational transition frequencies by comparison with state-of-the-art data obtained with grid-based approaches and modern configuration interaction methods. For HCN, we improve on the accuracy of the latter approaches with respect to both the ground-state absolute energy and proton density, which is a major challenge for multireference nuclear-electronic state-of-the-art methods.

Multicomponent MP4 and the inclusion of triple excitations in multicomponent many-body methods

This study implements the full multicomponent third-order (MP3) and fourth-order (MP4) many-body perturbation theory methods for the first time. Previous multicomponent studies have only implemented a subset of the full contributions, and the present implementation is the first multicomponent many-body method to include any connected triples contribution to the electron-proton correlation energy. The multicomponent MP3 method is shown to be comparable in accuracy to the multicomponent coupled-cluster doubles method for the calculation of proton affinities, while the multicomponent MP4 method is of similar accuracy as the multicomponent coupled-cluster singles and doubles method. From the results in this study, it is hypothesized that the relative accuracy of multicomponent methods is more similar to their single-component counterparts than previously assumed. It is demonstrated that for multicomponent MP4, the fourth-order triple-excitation contributions can be split into electron-electron and electron-proton contributions and the electron-electron contributions ignored with very little loss of accuracy of protonic properties.

Direct Dynamics with Nuclear-Electronic Orbital Density Functional Theory

Direct dynamics simulations of chemical reactions typically require the selection of a method for generating the potential energy surfaces and a method for the dynamical propagation of the nuclei on these surfaces. The nuclear-electronic orbital (NEO) framework avoids this Born-Oppenheimer separation by treating specified nuclei on the same level as the electrons with wave function methods or density functional theory (DFT). The NEO approach is particularly applicable to proton, hydride, and proton-coupled electron transfer reactions, where the transferring proton(s) and all electrons are treated quantum mechanically. In this manner, the zero-point energy, density delocalization, and anharmonicity of the transferring protons are inherently and efficiently included in the energies, optimized geometries, and dynamics.This Account describes how various NEO methods can be used for direct dynamics simulations on electron-proton vibronic surfaces. The strengths and limitations of these approaches are discussed, and illustrative examples are presented. The NEO-DFT method can be used to simulate chemical reactions on the ground state vibronic surface, as illustrated by the application to hydride transfer in C₄H₉⁺. The NEO multistate DFT (NEO-MSDFT) method is useful for simulating ground state reactions in which the proton density becomes bilobal during the dynamics, a characteristic of hydrogen tunneling, as illustrated by proton transfer in malonaldehyde. The NEO time-dependent DFT (NEO-TDDFT) method produces excited electronic, vibrational, and vibronic surfaces. The application of linear-response NEO-TDDFT to H₂ and H₃⁺, as well as the partially and fully deuterated counterparts, shows that this approach produces accurate fundamental vibrational excitation energies when all nuclei and all electrons are treated quantum mechanically. Moreover, when only specified nuclei are treated quantum mechanically, this approach can be used to optimize geometries on excited state vibronic surfaces, as illustrated by photoinduced single and double proton transfer systems, and to conduct adiabatic dynamics on these surfaces. The real-time NEO-TDDFT method provides an alternative approach for simulating nonequilibrium nuclear-electronic dynamics of such systems. These various NEO methods can be combined with nonadiabatic dynamics methods such as Ehrenfest and surface hopping dynamics to include the nonadiabatic effects between the quantum and classical subsystems. The real-time NEO-TDDFT Ehrenfest dynamics simulation of excited state intramolecular proton transfer in o-hydroxybenzaldehyde illustrates the power of this type of combined approach. The field of multicomponent quantum chemistry is in the early stages, and the methods discussed herein provide the foundation for a wide range of promising future directions to be explored. An appealing future direction is the expansion of the real-time NEO-TDDFT method to describe the dynamics of all nuclei and electrons on the same level. Direct dynamics simulations using NEO wave function methods such as equation-of-motion coupled cluster or multiconfigurational approaches are also attractive but computationally expensive options. The further development of NEO direct dynamics methods will enable the simulation of the nuclear-electronic dynamics for a vast array of chemical and biological processes that extend beyond the Born-Oppenheimer approximation.

Nucleus-electron correlation revising molecular bonding fingerprints from the exact wavefunction factorization

We present a novel theory and implementation for computing coupled electronic and quantal nuclear subsystems on a single potential energy surface, moving beyond the standard Born-Oppenheimer (BO) separation of nuclei and electrons. We formulate an exact self-consistent nucleus-electron embedding potential from the single product molecular wavefunction and demonstrate that the fundamental behavior of the correlated nucleus-electron can be computed for mean-field electrons that are responsive to a quantal anharmonic vibration of selected nuclei in a discrete variable representation. Geometric gauge choices are discussed and necessary for formulating energy invariant biorthogonal electronic equations. Our method is further applied to characterize vibrationally averaged molecular bonding properties of molecular energetics, bond lengths, and protonic and electron densities. Moreover, post-Hartree-Fock electron correlation can be conveniently computed on the basis of nucleus-electron coupled molecular orbitals, as demonstrated for correlated models of second-order Møllet-Plesset perturbation and full configuration interaction theories. Our approach not only accurately quantifies non-classical nucleus-electron couplings for revising molecular bonding properties but also provides an alternative time-independent approach for deploying non-BO molecular quantum chemistry.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Nuclear-electronic orbital methods: Foundations and prospects

The incorporation of nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations and molecular dynamics simulations is a longstanding challenge. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons with wave function and density functional theory methods. This approach inherently includes nuclear delocalization and zero-point energy in molecular energy calculations, geometry optimizations, reaction paths, and dynamics. It can also provide accurate descriptions of excited electronic, vibrational, and vibronic states as well as nuclear tunneling and nonadiabatic dynamics. Nonequilibrium nuclear-electronic dynamics simulations beyond the Born-Oppenheimer approximation can be used to investigate a wide range of excited state processes. This Perspective provides an overview of the foundational NEO methods and enumerates the prospects for using these methods as building blocks for future developments. The conceptual simplicity and computational efficiency of the NEO approach will enhance its accessibility and applicability to diverse chemical and biological systems.

Analytical Gradients for Nuclear-Electronic Orbital Time-Dependent Density Functional Theory: Excited-State Geometry Optimizations and Adiabatic Excitation Energies

The computational investigation of photochemical processes often entails the calculation of excited-state geometries, energies, and energy gradients. The nuclear-electronic orbital (NEO) approach treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons, thereby including the associated nuclear quantum effects and non-Born-Oppenheimer behavior into quantum chemistry calculations. The multicomponent density functional theory (NEO-DFT) and time-dependent DFT (NEO-TDDFT) methods allow efficient calculations of ground and excited states, respectively. Herein, the analytical gradients are derived and implemented for the NEO-TDDFT method and the associated Tamm-Dancoff approximation (NEO-TDA). The programmable equations for these analytical gradients as well as the NEO-DFT analytical Hessian are provided. The NEO approach includes the anharmonic zero-point energy (ZPE) and density delocalization associated with the quantum protons as well as vibronic mixing in geometry optimizations and energy calculations of ground and excited states. The harmonic ZPE associated with the other nuclei can be computed via the NEO Hessian. This approach is used to compute the 0-0 adiabatic excitation energies for a set of nine small molecules with all protons quantized, exhibiting slight improvement over the conventional electronic approach. Geometry optimizations of two excited-state intramolecular proton-transfer systems, [2,2'-bipyridyl]-3-ol and [2,2'-bipyridyl]-3,3'-diol, are performed with one and two quantized protons, respectively. The NEO calculations for these systems produce electronically excited-state geometries with stronger intramolecular hydrogen bonds and similar relative stabilities compared to conventional electronic methods. This work provides the foundation for nonadiabatic dynamics simulations of fundamental processes such as photoinduced proton transfer and proton-coupled electron transfer.

Multicomponent CASSCF Revisited: Large Active Spaces Are Needed for Qualitatively Accurate Protonic Densities

Multicomponent methods seek to treat select nuclei, typically protons, fully quantum mechanically and equivalent to the electrons of a chemical system. In such methods, it is well-known that due to the neglect of electron-proton correlation, a Hartree-Fock (HF) description of the electron-proton interaction catastrophically fails leading to qualitatively incorrect protonic properties. In single-component quantum chemistry, the qualitative failure of HF is normally indicative of the need for multireference methods such as complete active space self-consistent field (CASSCF). While a multicomponent CASSCF method was implemented nearly 20 years ago, it is only able to perform calculations with very small active spaces (∼10⁵ multicomponent configurations). Therefore, in order to extend the realm of applicability of the multicomponent CASSCF method, this study derives and implements a new two-step multicomponent CASSCF method that uses multicomponent heat-bath configuration interaction for the configuration interaction step, enabling calculations with very large active spaces (up to 16 electrons in 48 orbitals). We find that large electronic active spaces are needed to obtain qualitatively accurate protonic densities for the HCN and FHF^- molecules. Additionally, the multicomponent CASSCF method implemented here should have further applications for double-well protonic potentials and systems that are inherently electronically multireference.

Nuclear-electronic all-particle density matrix renormalization group

We introduce the Nuclear-Electronic All-Particle Density Matrix Renormalization Group (NEAP-DMRG) method for solving the time-independent Schrödinger equation simultaneously for electrons and other quantum species. In contrast to the already existing multicomponent approaches, in this work, we construct from the outset a multi-reference trial wave function with stochastically optimized non-orthogonal Gaussian orbitals. By iterative refining of the Gaussians' positions and widths, we obtain a compact multi-reference expansion for the multicomponent wave function. We extend the DMRG algorithm to multicomponent wave functions to take into account inter- and intra-species correlation effects. The efficient parameterization of the total wave function as a matrix product state allows NEAP-DMRG to accurately approximate the full configuration interaction energies of molecular systems with more than three nuclei and 12 particles in total, which is currently a major challenge for other multicomponent approaches. We present the NEAP-DMRG results for two few-body systems, i.e., H₂ and H₃ ⁺, and one larger system, namely, BH₃.

Molecular Vibrational Frequencies with Multiple Quantum Protons within the Nuclear-Electronic Orbital Framework

The nuclear-electronic orbital (NEO) approach treats all electrons and specified nuclei, typically protons, on the same quantum mechanical level. Proton vibrational excitations can be calculated using multicomponent time-dependent density functional theory (NEO-TDDFT) for fixed classical nuclei. Recently the NEO-DFT(V) approach was developed to enable the calculation of molecular vibrational frequencies for modes composed of both classical and quantum nuclei. This approach uses input from NEO-TDDFT to construct an extended NEO Hessian that depends on the expectation values of the quantum protons as well as the classical nuclear coordinates. Herein strategies are devised for extending these approaches to molecules with multiple quantum protons in a self-contained, effective, and computationally practical manner. The NEO-TDDFT method is shown to describe vibrational excitations corresponding to collective nuclear motions, such as linear combinations of proton vibrational excitations. The NEO-DFT(V) approach is shown to incorporate the most significant anharmonic effects in the molecular vibrations, particularly for the hydrogen stretching modes. These theoretical strategies pave the way for a wide range of multicomponent quantum chemistry applications.

Multicomponent density functional theory: Including the density gradient in the electron-proton correlation functional for hydrogen and deuterium

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.

Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies

The nuclear-electronic orbital (NEO) framework enables computationally practical coupled cluster calculations of multicomponent molecular systems, in which all electrons and specified nuclei, typically protons, are treated quantum mechanically. In addition to energies, computing accurate proton densities is essential for the calculation of reliable molecular properties, including vibrationally averaged geometries and vibrational frequencies. Herein, the Lagrangian formalism for the multicomponent coupled cluster with single and double excitations (NEO-CCSD) method is derived and implemented. The multicomponent coupled cluster with double excitations method using optimized Brueckner orbitals, denoted as NEO-BCCD, is also developed. Both of these methods are used to compute the proton densities for two molecular systems. The results illustrate that orbital relaxation effects, which can be included either indirectly with the NEO-CCSD method or directly with the NEO-BCCD method, are critical for computing even qualitatively accurate proton densities. Both methods are also able to provide accurate proton affinities and vibrationally averaged optimized geometries. This Lagrangian formalism will enable the calculation of other properties such as analytical nuclear gradients and Hessians with NEO coupled cluster methods. Moreover, the accuracy of these methods may be improved systematically by the inclusion of higher-order excitations. Thus, this work provides the foundation for a wide range of future methodological developments and applications within the NEO framework.

Diagonal Born-Oppenheimer Corrections within the Nuclear-Electronic Orbital Framework

The nuclear-electronic orbital (NEO) method treats specified nuclei, typically protons, quantum mechanically on the same level as the electrons. This approach invokes the Born-Oppenheimer separation between the quantum and classical nuclei, as well as the conventional separation between the electrons and classical nuclei. To test the validity of this additional adiabatic approximation, herein the diagonal Born-Oppenheimer correction (DBOC) within the NEO framework is derived, analyzed, and calculated numerically for a set of eight molecules. Inclusion of the NEO DBOC is found to change the equilibrium bond lengths by only ∼10^-4 Å and the heavy atom vibrational stretching frequencies by ∼1-2 cm^-1 per quantum proton bonded to an atom participating in the vibrational mode. These results imply that the DBOC does not significantly impact molecular properties computed with the NEO approach, although it can be included when necessary. Understanding the physical characteristics and quantitative contributions of the DBOC has broad implications for applications of multicomponent density functional theory and wave function methods.