On Koopmans' theorem in density functional theory

This paper clarifies why long-range corrected (LC) density functional theory gives orbital energies quantitatively. First, the highest occupied molecular orbital and the lowest unoccupied molecular orbital energies of typical molecules are compared with the minus vertical ionization potentials (IPs) and electron affinities (EAs), respectively. Consequently, only LC exchange functionals are found to give the orbital energies close to the minus IPs and EAs, while other functionals considerably underestimate them. The reproducibility of orbital energies is hardly affected by the difference in the short-range part of LC functionals. Fractional occupation calculations are then carried out to clarify the reason for the accurate orbital energies of LC functionals. As a result, only LC functionals are found to keep the orbital energies almost constant for fractional occupied orbitals. The direct orbital energy dependence on the fractional occupation is expressed by the exchange self-interaction (SI) energy through the potential derivative of the exchange functional plus the Coulomb SI energy. On the basis of this, the exchange SI energies through the potential derivatives are compared with the minus Coulomb SI energy. Consequently, these are revealed to be cancelled out only by LC functionals except for H, He, and Ne atoms.     

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.     

Water cluster anions studied by the long-range corrected density functional theory

Long-range corrected density functional theory (LC-DFT) is applied to a series of small water cluster anions(n= 2-6) to compute their vertical detachment energies (VDEs). The LC scheme is shown to eliminate an unphysical overestimation of the electron-water attraction in the hybrid functional by properly accounting for the long-range exchange repulsions. It is shown that a correct correlation energy behavior for a rapidly varying density is also important for describing a spatially extent, excess electron. The one-parameter progressive (OP) correlation functional, which satisfies this condition, leads to a remarkable improvement in the calculated VDE over the conventional one. The LC-BOP method produces highly accurate VDEs with a mean absolute deviation of 13.8 meV from the reference CCSD(T) results, reducing the error of B3LYP by more than 15 times. LC-BOP is found to be more accurate than MP2 which yields an excess electron underbound by 43.6 meV. The effect of basis sets on the calculated VDE is also examined. The aug-cc-pVDZ basis set with an extra diffuse function is found to be more accurate and reliable than the extended Pople-type basis sets used in the previous works. The extrapolation of the calculated VDE of different electron binding motifs is compared with the VDEs of experimentally observed three isomers (Verlet, J. R. R.; Bragg,A. E.; Kammrath, A.; Cheshnovsky, O.; Neumark, D. M. Science 2005, 307, 93).     

Spin-orbit relativistic long-range corrected time-dependent density functional theory for investigating spin-forbidden transitions in photochemical reactions

A long-range corrected (LC) time-dependent density functional theory (TDDFT) incorporating relativistic effects with spin-orbit couplings is presented. The relativistic effects are based on the two-component zeroth-order regular approximation Hamiltonian. Before calculating the electronic excitations, we calculated the ionization potentials (IPs) of alkaline metal, alkaline-earth metal, group 12 transition metal, and rare gas atoms as the minus orbital (spinor) energies on the basis of Koopmans' theorem. We found that both long-range exchange and spin-orbit coupling effects are required to obtain Koopmans' IPs, i.e., the orbital (spinor) energies, quantitatively in DFT calculations even for first-row transition metals and systems containing large short-range exchange effects. We then calculated the valence excitations of group 12 transition metal atoms and the Rydberg excitations of rare gas atoms using spin-orbit relativistic LC-TDDFT. We found that the long-range exchange and spin-orbit coupling effects significantly contribute to the electronic spectra of even light atoms if the atoms have low-lying excitations between orbital spinors of quite different electron distributions.     

Excited state geometry optimizations by analytical energy gradient of long-range corrected time-dependent density functional theory

An analytical excitation energy gradient of long-range corrected time-dependent density functional theory (LC-TDDFT) is presented. This is based on a previous analytical TDDFT gradient formalism, which avoids solving the coupled-perturbed Kohn-Sham equation for each nuclear degree of freedom. In LC-TDDFT, exchange interactions are evaluated by combining the short-range part of a DFT exchange functional with the long-range part of the Hartree-Fock exchange integral. This LC-TDDFT gradient was first examined by calculating the excited state geometries and adiabatic excitation energies of small typical molecules and a small protonated Schiff base. As a result, we found that long-range interactions play a significant role even in valence excited states of small systems. This analytical LC-TDDFT gradient was also applied to the investigations of small twisted intramolecular charge transfer (TICT) systems. By comparing with calculated ab initio multireference perturbation theory and experimental results, we found that LC-TDDFT gave much more accurate absorption and fluorescence energies of these systems than those of conventional TDDFTs using pure and hybrid functionals. For optimized excited state geometries, LC-TDDFT provided fairly different twisting and wagging angles of these small TICT systems in comparison with conventional TDDFT results.     

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.     

Application of the stabilization method to temporary anion states of CH3CN, CH3NC, CH3SCN, and CH3NCS in density functional theory with asymptotically corrected potentials

In this paper, density functional theory (DFT) with asymptotically corrected potentials is used to investigate CH(3)CN, CH(3)NC, CH(3)SCN, and CH(3)NCS molecules. For the energies of σ* and π* temporary anion states, the stabilized Koopmans' theorem (S-KT) using long-range correction functional and stabilized Koopmans-based (S-KB) approximation using local functional, are adopted. The stabilization procedure is accomplished by varying the exponents of appropriate diffuse functions. Results indicate that the calculations based on asymptotically corrected density functionals can yield better energy results of temporary anion states over conventional DFT methods.     

Effect of the damping function in dispersion corrected density functional theory

It is shown by an extensive benchmark on molecular energy data that the mathematical form of the damping function in DFT-D methods has only a minor impact on the quality of the results. For 12 different functionals, a standard "zero-damping" formula and rational damping to finite values for small interatomic distances according to Becke and Johnson (BJ-damping) has been tested. The same (DFT-D3) scheme for the computation of the dispersion coefficients is used. The BJ-damping requires one fit parameter more for each functional (three instead of two) but has the advantage of avoiding repulsive interatomic forces at shorter distances. With BJ-damping better results for nonbonded distances and more clear effects of intramolecular dispersion in four representative molecular structures are found. For the noncovalently-bonded structures in the S22 set, both schemes lead to very similar intermolecular distances. For noncovalent interaction energies BJ-damping performs slightly better but both variants can be recommended in general. The exception to this is Hartree-Fock that can be recommended only in the BJ-variant and which is then close to the accuracy of corrected GGAs for non-covalent interactions. According to the thermodynamic benchmarks BJ-damping is more accurate especially for medium-range electron correlation problems and only small and practically insignificant double-counting effects are observed. It seems to provide a physically correct short-range behavior of correlation/dispersion even with unmodified standard functionals. In any case, the differences between the two methods are much smaller than the overall dispersion effect and often also smaller than the influence of the underlying density functional.     

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.

     

Failure of time-dependent density functional theory for long-range charge-transfer excited states: the zincbacteriochlorin-bacteriochlorin and bacteriochlorophyll-spheroidene complexes

It is well-known that time-dependent density functional theory (TDDFT) yields substantial errors for the excitation energies of charge-transfer (CT) excited states, when approximate standard exchange-correlation (xc) functionals are used, for example, SVWN, BLYP, or B3LYP. Also, the correct 1/R asymptotic behavior of CT states with respect to a distance coordinate R between the separated charges of the CT state is not reproduced by TDDFT employing these xc-functionals. Here, we demonstrate by analysis of the TDDFT equations that the first failure is due to the self-interaction error in the orbital energies from the ground-state DFT calculation, while the latter is a similar self-interaction error in TDDFT arising through the electron transfer in the CT state. Possible correction schemes, such as inclusion of exact Hartree-Fock or exact Kohn-Sham exchange, as well as aspects of the exact xc-functional are discussed in this context. Furthermore, a practical approach is proposed which combines the benefits of TDDFT and configuration interaction singles (CIS) and which does not suffer from electron-transfer self-interaction. The latter approach is applied to a (1,4)-phenylene-linked zincbacteriochlorin-bacteriochlorin complex and to a bacteriochlorophyll-spheroidene complex, in which CT states may play important roles in energy and electron-transfer processes. The errors of TDDFT alone for the CT states are demonstrated, and reasonable estimates for the true excitation energies of these states are given.     

Modeling temporary anions in density functional theory: calculation of the Fukui function

Two approaches are investigated for modeling electron densities of temporary anions in density functional theory (DFT). Both rely on an artificial binding of the excess electron, in one case by a compact basis set and in the other by a potential wall. The key feature of the calculations is that the degree of binding is controlled in both cases by knowledge of the negative electron affinity of the corresponding neutral, approximated in terms of DFT local functional frontier orbital eigenvalues and vertical ionization potential, A=-(epsilon(LUMO)+epsilon(HOMO))-I. To illustrate the two approaches, Fukui functions for nucleophilic attack are determined in four molecules with increasingly negative electron affinities. They yield very similar results, which are notably different to those determined without artificial electron binding. The use of a potential wall has the attractive feature that large, diffuse basis sets can be used, avoiding the need for a compact basis, tailored to a particular molecule.     

Limitations of density functional theory in application to degenerate states

It is shown that the claims that density functional theory (DFT) can handle orbitally degenerate states are ungrounded. The constraint search formulation of DFT allows one to determine a set of densities and eigenvalues for the degenerate term that, however, are neither observables, nor can they be used to solve the system of coupled equations for the nuclear motions to obtain observables, as in the wave function presentation. A striking example of the failure of the existing versions of DFT to describe degenerate states is provided by the Berry phase problem: the strong dependence of the results on the phase properties of the electronic wave function that are smeared out in the density formulation. The solution of the Jahn-Teller E-e problem illustrates these statements. For nondegenerate states with the full wave function taken in the adiabatic approximation as a product of the electronic and nuclear parts, the formulation of DFT is rigorous if and only if the dependence of the electronic wave function on nuclear coordinates is ignored. This lowers the accuracy of the results, in general, and may lead to erroneous presentation as in the case of molecular systems in strong magnetic fields. © 1997 by John Wiley & Sons, Inc.     

Application of semiempirical long-range dispersion corrections to periodic systems in density functional theory

Ewald summation is used to apply semiempirical long-range dispersion corrections (Grimme, J Comput Chem 2006, 27, 1787; 2004, 25, 1463) to periodic systems in density functional theory. Using the parameters determined before for molecules and the Perdew-Burke-Ernzerhof functional, structure parameters and binding energies for solid methane, graphite, and vanadium pentoxide are determined in close agreement with observed values. For methane, a lattice constant a of 580 pm and a sublimation energy of 11 kJ mol(-1) are calculated. For the layered solids graphite and vanadia, the interlayer distances are 320 pm and 450 pm, respectively, whereas the graphite interlayer energy is -5.5 kJ mol(-1) per carbon atom and layer. Only when adding the semiempirical dispersion corrections, realistic values are obtained for the energies of adsorption of C(4) alkenes in microporous silica (-66 to -73 kJ mol(-1)) and the adsorption and chemisorption (alkoxide formation) of isobutene on acidic sites in the micropores of zeolite ferrierite (-78 to -94 kJ mol(-1)). As expected, errors due to missing self-interaction correction as in the energy for the proton transfer from the acidic site to the alkene forming a carbenium ion are not affected by the dispersion term. The adsorption and reaction energies are compared with the results from M?ller-Plesset second-order perturbation theory with basis set extrapolation.     

Re-examining the properties of the aqueous vapor-liquid interface using dispersion corrected density functional theory

First-principles molecular dynamics simulations, in which the forces are computed from electronic structure calculations, have great potential to provide unique insight into structure, dynamics, electronic properties, and chemistry of interfacial systems that is not available from empirical force fields. The majority of current first-principles simulations are driven by forces derived from density functional theory with generalized gradient approximations to the exchange-correlation energy, which do not capture dispersion interactions. We have carried out first-principles molecular dynamics simulations of air-water interfaces employing a particular generalized gradient approximation to the exchange-correlation functional (BLYP), with and without empirical dispersion corrections. We assess the utility of the dispersion corrections by comparison of a variety of structural, dynamic, and thermodynamic properties of bulk and interfacial water with experimental data, as well as other first-principles and force field-based simulations.     

Double-hybrid density functional theory for excited electronic states of molecules

Double-hybrid density functionals are based on a mixing of standard generalized gradient approximations (GGAs) for exchange and correlation with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional (dubbed B2PLYP) contains only two empirical parameters that describe the mixture of HF and GGA exchange (ax) and of the PT2 and GGA correlation (ac), respectively. Extensive testing has recently demonstrated the outstanding accuracy of this approach for various ground state problems in general chemistry applications. The method is extended here without any further empirical adjustments to electronically excited states in the framework of time-dependent density functional theory (TD-DFT) or the closely related Tamm-Dancoff approximation (TDA-DFT). In complete analogy to the ground state treatment, a scaled second-order perturbation correction to configuration interaction with singles (CIS(D)) wave functions developed some years ago by Head-Gordon et al. [Chem. Phys. Lett. 219, 21 (1994)] is computed on the basis of density functional data and added to the TD(A)-DFTGGA excitation energy. The method is implemented by applying the resolution of the identity approximation and the efficiency of the code is discussed. Extensive tests for a wide variety of molecules and excited states (of singlet, triplet, and doublet multiplicities) including electronic spectra are presented. In general, rather accurate excitation energies (deviations from reference data typically <0.2 eV) are obtained that are mostly better than those from standard functionals. Still, systematic errors are obtained for Rydberg (too low on average by about 0.3 eV) and charge-transfer transitions but due to the relatively large ax parameter (0.53), B2PLYP outperforms most other functionals in this respect. Compared to conventional HF-based CIS(D), the method is more robust in electronically complex situations due to the implicit account of static correlation effects by the GGA parts. The (D) correction often works in the right direction and compensates for the overestimation of the transition energy at the TD level due to the elevated fraction of HF exchange in the hybrid GGA part. Finally, the limitations of the method are discussed for challenging systems such as transition metal complexes, cyanine dyes, and multireference cases.     

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.     

Asymptotic correction of the exchange-correlation kernel of time-dependent density functional theory for long-range charge-transfer excitations

Time-dependent density functional theory (TDDFT) calculations of charge-transfer excitation energies omegaCT are significantly in error when the adiabatic local density approximation (ALDA) is employed for the exchange-correlation kernel fxc. We relate the error to the physical meaning of the orbital energy of the Kohn-Sham lowest unoccupied molecular orbital (LUMO). The LUMO orbital energy in Kohn-Sham DFT--in contrast to the Hartree-Fock model--approximates an excited electron, which is correct for excitations in compact molecules. In CT transitions the energy of the LUMO of the acceptor molecule should instead describe an added electron, i.e., approximate the electron affinity. To obtain a contribution that compensates for the difference, a specific divergence of fxc is required in rigorous TDDFT, and a suitable asymptotically correct form of the kernel fxc(asymp) is proposed. The importance of the asymptotic correction of fxc is demonstrated with the calculation of omegaCT(R) for the prototype diatomic system HeBe at various separations R(He-Be). The TDDFT-ALDA curve omegaCT(R) roughly resembles the benchmark ab initio curve omegaCT CISD(R) of a configuration interaction calculation with single and double excitations in the region R=1-1.5 A, where a sizable He-Be interaction exists, but exhibits the wrong behavior omegaCT(R)相似文献