Koopmans' theorem for large molecular systems within density functional theory

It is shown that in density functional theory (DFT), Koopmans' theorem for a large molecular system can be stated as follows: The ionization energy of the system equals the negative of the highest occupied molecular orbital (HOMO) energy plus the Coulomb electrostatic energy of removing an electron from the system, or equivalently, the ionization energy of an N-electron system is the negative of the arithmetic average of the HOMO energy of this system and the lowest unoccupied molecular orbital (LUMO) energy of the (N - 1)-electron system. Relations between this DFT Koopmans' theorem and its existing counterparts in the literature are discussed. Some of the previous results are generalized and some are simplified. DFT calculation results of a fullerene molecule, a finite single-walled carbon nanotube and a finite boron nitride nanotube are presented, indicating that this Koopmans' theorem approximately holds, even if the orbital relaxation is taken into consideration.     

Characterization of the temporary anion states on perfluoroalkanes via stabilized Koopmans' theorem in long-range corrected density functional theory

The stabilized Koopmans' theorem (SKT) in long-range corrected density functional theory is used to characterize the temporary anion states of perfluoro-n-alkanes (n-PFAs) from C(2) to C(5), and perfluorocycloalkanes (c-PFAs) from C(3) to C(4). In this approach, stabilization is accomplished by varying the exponents of appropriate diffuse functions. The energies of temporary anion states are then identified by investigating the relationship between the resultant eigenvalues and scale parameter. The characteristics of resonance orbitals are also examined. For the lowest unfilled orbitals of perfluoroalkanes, results indicate that they are mainly from the π-bonding interactions between all neighboring C atoms. In addition, their energies decrease as the sizes of the perfluoroalkanes increase. Moreover, the energies of the c-C(3)F(6)/c-C(4)F(8) are lower than those of the corresponding n-C(3)F(8)/n-C(4)F(10). When compared with experimental data, our SKT calculations can yield conformable results. Thus, this SKT approach can provide more information on the resonance states of perfluoroalkanes.     

The extended Koopmans' theorem: vertical ionization potentials from natural orbital functional theory

The Piris natural orbital functional, PNOF5, has been used to predict vertical ionization potentials of a selected set of 30 organic and inorganic spin-compensated molecules by means of the extended Koopmans' theorem. Electron affinities of 10 selected radicals have also been estimated as the inverse of the ionization potentials of the anionic species, calculated at the experimental geometries of the neutral radicals. The basis set limit effects have been assessed by inspecting the data obtained for the Dunning's basis set series cc-pVXZ and aug-cc-pVXZ (X = D, T, Q, 5). The performance of the PNOF5 is established by carrying out a statistical analysis of the mean absolute errors (MAEs) with respect to the experiment values. The calculated PNOF5 ionization potentials and electron affinities agree satisfactorily with the corresponding experimental data, with MAEs smaller than 0.5 eV.     

Extended Koopmans' theorem at the second-order perturbation theory

The extended Koopmans' theorem (EKT), when combined with the second-order Møller−Plesset (MP2) perturbation theory through the relaxed density matrix approach [J. Cioslowski, P. Piskorz, and G. Liu, J. Chem. Phys. 1997, 107, 6,804], provides a straightforward way to calculate the ionization potentials (IPs) as an one electron quantity. However, such an EKT-MP2 method often suffers from the negative occupation problem, failing to provide the complete IP spectra for a system of interest. Here a small positive number scheme is proposed to cure this problem so as to remove the associated unphysical results. In order to obtain an in-depth physical interpretation of the EKT-MP2 method, we introduce a Koopmans-type quantity, named KT-MP2, based on which the respective contribution from the relaxation and the correlation parts in the EKT-MP2 results are recognized. Furthermore, the close relationship between the EKT-MP2 method and the derivative approach of the MP2 energy with respect to the orbital occupation numbers [N. Q. Su and X. Xu, J. Chem. Theory Comput. 2015, 11, 4,677] is revealed. When these MP2-based methods are applied to a set of atoms and molecules, new insights are gained on the role played by the relaxation and the correlation effects in the electron ionization processes.     

Spin virial theorem in the density functional theory

The spin virial theorem is derived in the density functional theory. The theorem establishes a relation between the differences of spin-up and -down kinetic and potential energies. The theorem is useful for checking the accuracy of spin orbitals. As an illustration, the example of the Xα method is studied. © 1994 John Wiley & Sons, Inc.     

Koopmans' theorem and virtual orbital energies in the general SCF theory

The equality of the ionization potential and the orbital energy of the electron being removed is investigated using a general SCF theory for open-shell configurations. The significance of virtual orbital energies is investigated in the same context.

---

Zusammenfassung Die Gleichsetzung von Ionisierungsenergie und entsprechender Einelektronenenergie wird für Konfigurationen mit unabgeschlossenen Schalen mit Hilfe einer allgemeinen SCE-Theorie geprüft. In diesem Rahmen wird auch die Bedeutung von virtuellen Einelektronenenergien untersucht.

---

Résumé L'équation entre le potentiel d'ionisation et l'énergie orbitale de l'électron correspondant est examinée à l'aide d'une théorie générale SCE pour les configurations à couches ouvertes. Au cadre de cette théorie, la signification des énergies d'orbitales inoccupées est étudiée.

     

The extended Koopmans' theorem Fock operator

By expanding the wave function of a system of N particles in terms of products of functions of one and (N-1) particles, the one-particle, nonlocal operator F?_EKT (extended Koopmans' theorem) is determined. It is shown that although this operator is nonhermitian, its eigenvalues and eigenfunctions represent the ionization energies and occupied orbitals, respectively. The eigenfunctions of F?_EKT are the one-particle functions that enter into the expansion of the wave function of the system as partners of the (N-1)-particle wave functions. The eingenvalues are also one-particle energies that, multipled by the orbital occupancy probalities, enter the expression for the total N-particle energy of the system.     

The virial theorem as an entail for Koopmans' theorem in atomic calculations

The frozen approximation in the Koopmans' theorem suggests that the virial theorem is not obeyed. By imposing the virial theorem to Koopmans' theorem, we observe that the ionization potential of an atomic orbital is directly related to the respective kinetic energy and that the virial theorem introduces some reorganization effect on the electronic cloud. The quantity of reorganization introduced is not hazard, depending on the type of atom as well as the atomic orbital.     

On the quantum chemical origin for the nonvalidity of Koopmans' theorem in transitionmetal compounds

The quantum chemical origin for the nonvalidity of Koopmans' theorem in transitionmetal compounds of the 3d series is analyzed by means of the Green's function formalism applied in the framework of a semiempirical INDO Hamiltonian. In the case of ferrocene (1), cyclobutadiene iron tricarbonyl (2) and irontetracarbonyl dihydride (3) the self-energy part of a geometric approximation has been partitioned into relaxation and correlation (pair removal, pair relaxation) increments. The breakdown of Koopmans' theorem for strongly localized MOs with large Fe 3d amplitudes is predominantly the result of electronic relaxation lowering the calculated ionization potentials. On the other hand the variation of the pair correlation energy in the cationic hole-state is by no means negligible and acts into the opposite direction as the relaxation increment. These significant pair relaxation contributions explain the wellknown failtures of the ΔSCF approach in combination with large scaleab initio bases. The loss of ground state pair correlation in the outer valence region is small in comparison to relaxation and pair relaxation. The magnitude of the aforementioned reorganization increments has been studied as a function of the localization properties of the MOs and as a function of the one-electron energies of the available particle- and hole-states. The computational findings derived with the INDO model are compared with recentab initio studies.     

A modification of Koopmans' theorem by imposition of the virial theorem in molecules

Imposition of the virial theorem on Koopmans' theorem permits the introduction of some relaxation effect in the electronic cloud of atomic (less than 5%) or molecular (less than 1.3% for the systems studied) systems and a partitioning of the ionization energy. The method is applied in some diatomic hydrides. It is observed that the imposition of the virial theorem improves the ionization of the innermost molecular orbitals significantly, while the improvement is negligible for the outermost orbitals. The ionization energy is divided among three different terms that elucidate some aspects of the nature of the ionization process.     

On correlated electron-nuclear dynamics using time-dependent density functional theory

We discuss possibilities and challenges for describing correlated electron and nuclear dynamics within a surface-hopping framework using time-dependent density functional theory (TDDFT) for the electron dynamics. We discuss the recent surface-hopping method proposed by Craig et al. [Phys. Rev. Lett. 95, 163001 (2005)] that is based on Kohn-Sham potential energy surfaces. Limitations of this approach arise due to the Kohn-Sham surfaces generally having different gradients than the true TDDFT-corrected ones. Two mechanisms of the linear response procedure cause this effect: we illustrate these with examples.     

Embedding wave function theory in density functional theory

We present a framework for embedding a highly accurate coupled-cluster calculation within a larger density functional calculation. We use a perturbative buffer to help insulate the coupled-cluster region from the rest of the system. Regions are defined, not in real space, but in Hilbert space, though connection between the two can be made by spatial localization of single-particle orbitals. Relations between our embedding approach and some similar techniques are discussed. We present results for small sample systems for which we can extract essentially exact results, demonstrating that our approach seems to work quite well and is generally more reliable than some of the related approaches due to the introduction of additional interaction terms.     

Molecular fragments in density functional theory

The second-order density functional approach to the partitioning of the molecular density of Cedillo, Chattaraj, and Parr (Int. J. Quantum Chem. 2000, 77, 403-407) is used, together with a local assumption for the function that projects the total density into its components, to show that the distribution function adopts a stockholders form, in terms of the local softness of the isolated fragments, and that the molecular Fukui function is distributed in the molecular fragments in the same proportion as the electronic density.     

Extended Koopmans' theorem ionization potentials for beryllium atom shake-up transitions

The ionization potentials were calculated for Be using the extended Koopmans' theorem (EKT ) using several full configuration interaction (CI ) and multiconfigurational-self-consistent-field (MCSCF ) wave functions as reference wave functions. The wave functions used account for 89.7–96.7% of the correlation energy. Comparisons are made with experimental values and with δCI values calculated as the difference in energy obtained from CI wave functions for Be and Be⁺. The best EKT IP differed from the δCI value by 0.0003 eV for the lowest IP and by 0.0006 eV for ionization into the lowest ²P state of Be⁺. A calculation of ionization into the second ²P state of Be⁺ requires diffuse orbitals that are unimportant in the wave function for the ground state of Be. This results in small natural orbital occupation numbers for natural orbitals needed in the EKT calculation. © 1994 John Wiley & Sons, Inc.     

Multigrid methods in density functional theory

A multigrid method for real-space solution of the Kohr-Sham equations is presented. By using this multiscale approach, the problem of critical slowing down typical of iterative real-space solvers is overcome. The method scales linearly in computer time with the number of electrons if the orbitals are localized. Here, we describe details of our multigrid method, present preliminary many-electron numerical results illustrating the efficiency of the solver, and discuss its strengths and limitations. © 1997 John Wiley & Sons, Inc.     

Orbital-corrected orbital-free density functional theory

A new implementation of density functional theory (DFT), namely orbital-corrected orbital-free (OO) DFT, has been developed. With at most two non-self-consistent iterations, OO-DFT accomplishes the accuracy comparable to fully self-consistent Kohn-Sham DFT as demonstrated by its application on the cubic-diamond Si and the face-centered-cubic Ag systems. Our work provides a new impetus to further improve orbital-free DFT method and presents a robust means to significantly lower the cost associated with general applications of linear-scaling Kohn-Sham DFT methods on large systems of thousands of atoms within different chemical bonding environment.