On the validity of the extended Koopmans’ theorem

ERRORS IN APPROXIMATE IONIZATION ENERGIES DUE TO THE ONE-ELECTRON SPACE TRUNCATION OF THE EKT EIGENPROBLEM

Unless the approximate wavefunction of the parent system is expressed in terms of explicitly correlated basis functions, the finite size of the generalized Fock matrix is unlikely to be the leading source of the truncation error in the ionization energy E produced by the EKT (extended Koopmans' theorem) formalism. This conclusion is drawn from a rigorous analysis that involves error partitioning into the parent- and ionized-system contributions, the former being governed by asymptotic power laws when the underlying wavefunction is assembled from a large number of spinorbitals and the latter arising from the truncation of the infinite-dimensional matrix V whose elements involve the 1-, 2- and 3-matrices of the parent system. Quite surprisingly, the decay of the second contribution with the number $n$ of the natural spinorbitals (NOs) employed in the construction of the truncated V turns out to be strongly system-dependent even in the simplest case of the 1S states of two-electron systems, following the n-5 power law for the helium atom while exhibiting an erratic behavior for the H- anion. This phenomenon, which stems from the presence of the so-called solitonic natural spinorbitals among the NOs, renders the extrapolation of the EKT approximates of E to the complete-basis-set limit generally unfeasible. However, attaining that limit is not contingent upon attempted reproduction of the ill-defined one-electron function known as 'the removal orbital', which does not have to be invoked in the derivation of EKT and whose expansion in terms of the NOs diverges.

State-Of-The-Art Computations of Vertical Electron Affinities with the Extended Koopmans' Theorem Integrated with the CCSD(T) Method

Accurate computation of electron affinities (EAs), within 0.10 eV, is one of the most challenging problems in modern computational quantum chemistry. The extended Koopmans' theorem (EKT) enables direct computations of electron affinities (EAs) from any level of the theory. In this research, the EKT approach based on the coupled-cluster singles and doubles with perturbative triples [CCSD(T)] method is applied to computations of EAs for the first time. For efficiency, the density-fitting (DF) technique is used for two-electron integrals. Further, the EKT-CCSD(T) method is applied to three test sets of atoms and closed- and open-shell molecules, denoted A16, C10, and O33, respectively, for comparison with the experimental electron affinities. For the A16, C10, and O33 sets, the EKT-CCSD(T) approach, along with the aug-cc-pV5Z basis set, provide mean absolute errors (MAEs) of 0.05, 0.08, and 0.09 eV, respectively. Hence, our results demonstrate that high-accuracy computations of EAs can be achieved with the EKT-CCSD(T) approach. Further, when the EKT-CCSD(T) approach is not computationally affordable, the EKT-MP2.5, EKT-LCCD, and EKT-CCSD methods can be considered, and their results are also reasonably accurate. The huge advantage of the EKT method for the computation of IPs is that it comes for free in an analytic gradient computation. Hence, one needs neither separate computations for neutral and ionic species, as in the case of common approaches, nor additional efforts to obtain IPs, as in the case of equation-of-motion approaches. Overall, we believe that the present research may open new avenues in EA computations.

Complete-active-space extended Koopmans theorem method

The complete-active-space (CAS) extended Koopmans theorem (EKT) method is defined as a special case of the EKT in which the reference state is a CAS configuration interaction (CI) expansion and the electron removal operator acts only on the active orbitals. With these restrictions, the EKT is equivalent to the CI procedure involving all hole-state configurations derived from the active space of the reference wavefunction and has properties analogous to those of the original Koopmans theorem. The equivalence is used to demonstrate in a transparent manner that the first ionization energy predicted by the EKT is in general not exact, i.e., not equal to the difference between the full CI energies of the neutral and the ion, but can approach the full CI result with arbitrary precision even within a finite basis set. The findings also reconcile various statements about the EKT found in the literature.

New constraints upon the electron-electron repulsion energy functional of the one-electron reduced density matrix

Three strict constraints upon the electron-electron repulsion energy functional of the one-electron reduced density matrix (the 1-matrix) are obtained by combining its invariance and stationary properties with the extended Koopmans' theorem. The constraints relate the quantities derived from the functional pertaining to an N-electron system with those of its (N-1)-electron congener. Together with the N-representability requirement for the 1-matrix of the congener, identities involving the electron-electron repulsion energies of the two systems and their derivatives with respect to the 1-matrices seriously narrow down the choices for potential approximate density-matrix functionals. This fact is well illustrated in the case of two-electron systems, where the validity of the new constraints is confirmed and found to originate from a nontrivial cancellation among different terms. Thus, the constraints provide a new tool for the construction and testing of new functionals that complements the previously known conditions such as the reproduction of the homogeneous gas energies and momentum distributions, convexity, and the N-representability of the associated 2-matrices.