Kohn-Sham Orbital Energies Research Articles

Koopmans: An Open-Source Package for Accurately and Efficiently Predicting Spectral Properties with Koopmans Functionals.

Over the past decade we have developed Koopmans functionals, a computationally efficient approach for predicting spectral properties with an orbital-density-dependent functional framework. These functionals impose a generalized piecewise linearity condition to the entire electronic manifold, ensuring that orbital energies match the corresponding electron removal/addition energy differences (in contrast to semilocal DFT, where a mismatch between the two lies at the heart of the band gap problem and, more generally, the unreliability of Kohn-Sham orbital energies). This strategy has proven to be very powerful, yielding molecular orbital energies and solid-state band structures with comparable accuracy to many-body perturbation theory but at greatly reduced computational cost while preserving a functional formulation. This paper reviews the theory of Koopmans functionals, discusses the algorithms necessary for their implementation, and introduces koopmans, an open-source package that contains all of the code and workflows needed to perform Koopmans functional calculations and obtain reliable spectral properties of molecules and materials.

Fractional-Electron and Transition-Potential Methods for Core-to-Valence Excitation Energies Using Density Functional Theory.

Methods for computing X-ray absorption spectra based on a constrained core hole (possibly containing a fractional electron) are examined. These methods are based on Slater's transition concept and its generalizations, wherein core-to-valence excitation energies are determined using Kohn-Sham orbital energies. Methods examined here avoid promoting electrons beyond the lowest unoccupied molecular orbital, facilitating robust convergence. Variants of these ideas are systematically tested, revealing a best-case accuracy of 0.3-0.4 eV (with respect to experiment) for K-edge transition energies. Absolute errors are much larger for higher-lying near-edge transitions but can be reduced below 1 eV by introducing an empirical shift based on a charge-neutral transition-potential method, in conjunction with functionals such as SCAN, SCAN0, or B3LYP. This procedure affords an entire excitation spectrum from a single fractional-electron calculation, at the cost of ground-state density functional theory and without the need for state-by-state calculations. This shifted transition-potential approach may be especially useful for simulating transient spectroscopies or in complex systems where excited-state Kohn-Sham calculations are challenging.

Ionization potentials and fundamental gaps in atomic systems from the Ensemble-DFT approach.

Calculations in Kohn-Sham density functional theory crucially rely on high-quality approximations for the exchange-correlation (xc) functional. Standard local and semi-local approximations fail to predict the ionization potential (IP) and the fundamental gap, departing from the Kohn-Sham orbital energies, due to the deviation of the total energy from piecewise-linearity and the absence of the derivative discontinuity. The ensemble generalization procedure introduced in Phys. Rev. Lett. 110, 126403 (2013) restores, to a large extent, these features in any approximate xc functional and improves its ability to predict the IP and the fundamental gap with negligible additional computational effort. In this work we perform an extensive study of atoms and first ions across the Periodic Table, generalizing the local spin-density and the Perdew-Burke-Ernzerhof approximations. By applying the ensemble generalization to a variety of systems, with s-, p-, and d-character, we assess the accuracy of the method and identify important trends. In particular, we find that the accuracy of our approach heavily depends on the character of the frontier orbitals: when d-orbitals are involved, the performance is far less accurate. Possible sources of error are discussed and ways for further improvement are outlined.

Localized Resolution of Identity Approach to the Analytical Gradients of Random-Phase Approximation Ground-State Energy: Algorithm and Benchmarks.

We develop and implement a formalism which enables calculating the analytical gradients of particle-hole random-phase approximation (RPA) ground-state energy with respect to the atomic positions within the atomic orbital basis set framework. Our approach is based on a localized resolution of identity (LRI) approximation for evaluating the two-electron Coulomb integrals and their derivatives, and the density functional perturbation theory for computing the first-order derivatives of the Kohn-Sham (KS) orbitals and orbital energies. Our implementation allows one to relax molecular structures at the RPA level using both Gaussian-type orbitals (GTOs) and numerical atomic orbitals (NAOs). Benchmark calculations against previous implementations show that our approach delivers adequate numerical precision, highlighting the usefulness of LRI in the context of RPA gradient evaluations. A careful assessment of the quality of RPA geometries for small molecules reveals that post-KS RPA systematically overestimates the bond lengths. We furthermore optimized the geometries of the four low-lying water hexamers-cage, prism, cyclic, and book isomers, and determined the energy hierarchy of these four isomers using RPA. The obtained RPA energy ordering is in good agreement with that yielded by the coupled cluster method with single, double and perturbative triple excitations, despite that the dissociation energies themselves are appreciably underestimated. The underestimation of the dissociation energies by RPA is well corrected by the renormalized single excitation correction.

Machine Learning of Quasiparticle Energies in Molecules and Clusters.

We present a Δ-machine learning approach for the prediction of GW quasiparticle energies (ΔMLQP) and photoelectron spectra of molecules and clusters, using orbital-sensitive representations (OSRs) based on molecular Cartesian coordinates in kernel ridge regression-based supervised learning. Coulomb matrix, bag-of-bond, and bond-angle-torsion representations are made orbital-sensitive by augmenting them with atom-centered orbital charges and Kohn–Sham orbital energies, both of which are readily available from baseline calculations at the level of density functional theory (DFT). We first illustrate the effects of different constructions of the OSRs on the prediction of frontier orbital energies of 22k molecules of the QM8 data set and show that it is possible to predict the full photoelectron spectrum of molecules within the data set using a single model with a mean absolute error below 0.1 eV. We further demonstrate that the OSR-based ΔMLQP captures the effects of intra- and intermolecular conformations in application to water monomers and dimers. Finally, we show that the approach can be embedded in multiscale simulation workflows, by studying the solvatochromic shifts of quasiparticle and electron–hole excitation energies of solvated acetone in a setup combining molecular dynamics, DFT, the GW approximation, and the Bethe–Salpeter equation. Our findings suggest that the ΔMLQP model allows us to predict quasiparticle energies and photoelectron spectra of molecules and clusters with GW accuracy at DFT cost.

ICMMS-2: Metal−organic Frameworks for Hydrogen Storage: Theoretical Prospective

The adsorption of H2 in dehydrated and hydrated Materials of Institute Lavoisier (MIL-101) was investigated theoretically. The effect of terminal water molecules on adsorption as one of the more vital MIL-n trivalent chromium-based porous carboxylates in metal-organic frameworks application in the renewable energy field was also theoretically studied. The MIL-101 structures were optimized for geometry and energy minimization was performed. The calculations were carried out using density functional theory approach with B3LYP functional and mixed basis set of Lanl2DZ and 6-31G(d, p) for Cr and light atoms (C, H, O, F), respectively, as implemented in the Gaussian 09 program package. The spin and atomic charges distribution on the Cr metal atoms, adsorbate, and water molecules are calculated using natural bond orbital (NBO). The density of states (DOS) for the clusters was obtained using Gaussian smearing of Kohn–Sham orbital energies. The natural bond orbital (NBO) for molecular orbital analysis and atomic charge calculations were utilized. For the dehydrated MIL-101, more adsorbate molecules were found near the exposed Cr2 sites than the fluorine saturated Cr1 sites. Furthermore, terminal water molecules in the hydrated MIL-101 made more interaction sites and enhanced adsorption.

Density Functional Prediction of Quasiparticle, Excitation, and Resonance Energies of Molecules With a Global Scaling Correction Approach.

Molecular quasiparticle and excitation energies determine essentially the spectral characteristics measured in various spectroscopic experiments. Accurate prediction of these energies has been rather challenging for ground-state density functional methods, because the commonly adopted density function approximations suffer from delocalization error. In this work, by presuming a quantitative correspondence between the quasiparticle energies and the generalized Kohn–Sham orbital energies, and employing a previously developed global scaling correction approach, we achieve substantially improved prediction of molecular quasiparticle and excitation energies. In addition, we also extend our previous study on temporary anions in resonant states, which are associated with negative molecular electron affinities. The proposed approach does not require any explicit self-consistent field calculation on the excited-state species, and is thus highly efficient and convenient for practical purposes.

Charge-Transfer Excitation Energies Expressed as Orbital Energies of Kohn-Sham Density Functional Theory with Long-Range Corrected Functionals.

Previously proposed theoretical schemes for estimating one-electron excitation energies using Kohn-Sham (KS) solutions with long-range corrected (LC) functionals are applied to the charge-transfer (CT) excitations of the ethylene···tetrafluoroethylene (C2H4-C2F4) system, and the CT complex between an aromatic donor (Ar = benzene, toluene, o-xylene, naphthalene, anthracene, and various meso-substituted anthracenes) and the tetracyanoethylene (TCNE) acceptor. The CT excited state is described well as a single-electron excitation between specific orbitals of donor and acceptor. Thus, CT excitation energies are well approximated by the orbital energies because of the satisfaction of the Koopmans-type theorem and the asymptotic behavior of the LC functional. We have examined three computational schemes: scheme 1 employs the orbital energies for the neutral and cationic systems, scheme 2 utilizes orbital energies of just the cation, and in scheme 3, because the electron affinity of a molecule is the ionization energy of its anion, a scale factor is applied to enforce this identity. The present schemes reproduce the correct asymptotic behavior of CT excitation energy of C2H4···C2F4 for the long intermolecular distances and give good agreement with accurate ab initio results. Calculated CT excitation energies for Ar-TCNE are compared with those of TD-DFT and ΔSCF methods. Scheme 1 with the optimal range-separation parameter μ accurately reproduces CT excitation energies for all Ar-TCNE systems and gives good agreement with the best TD-DFT calculations and experiment. Scheme 1, scheme 3, and TD-DFT show similar tendencies with respect to the variation in μ. Scheme 2 and ΔSCF approaches are rather insensitive to changes in μ, but both considerably underestimate the CT excitation energies for these systems. KS orbital energies are physically meaningful and they are practically useful; if the range-separation parameter is tuned, then good results can be obtained.

Rapid Prediction of Ultraviolet-Visible Spectra from Conventional (Non-Time-Dependent) Density Functional Theory Calculations.

We present an approximate approach for the simulation of UV/vis spectra using conventional [non-time-dependent (non-TD)] DFT computations. It uses Kohn-Sham orbitals and orbital energies to estimate both the excitation energies and the associated oscillator strengths. For a wide range of systems from small molecules to large molecular dyes used in electrochromic and solar-cell applications, reasonable UV/vis spectra are generated, each with just two conventional DFT computations. The accuracy is generally comparable to what one would expect from TD-DFT calculations. In comparison to TD-DFT, the protocol of the present study provides an intuitive and notably more rapid means for simulating electronic absorption properties. It enables efficient screening of materials for a wide range of relevant applications.

Density Functional Extension to Excited-State Mean-Field Theory.

We investigate an extension of excited-state mean-field theory in which the energy expression is augmented with density functional components in an effort to include the effects of weak electron correlations. The approach remains variational and entirely time independent, allowing it to avoid some of the difficulties associated with linear response and the adiabatic approximation. In particular, all of the electrons' orbitals are relaxed state specifically, and there is no reliance on Kohn-Sham orbital energy differences, both of which are important features in the context of charge transfer. Preliminary testing shows clear advantages for single-component charge transfer states, but the method, at least in its current form, is less reliable for states in which multiple particle-hole transitions contribute significantly.

Approximating Quasiparticle and Excitation Energies from Ground State Generalized Kohn-Sham Calculations.

Quasiparticle energies and fundamental band gaps in particular are critical properties of molecules and materials. It was rigorously established that the generalized Kohn-Sham HOMO and LUMO orbital energies are the chemical potentials of electron removal and addition and thus good approximations to band edges and fundamental gaps from a density functional approximation (DFA) with minimal delocalization error. For other quasiparticle energies, their connection to the generalized Kohn-Sham orbital energies has not been established but remains highly interesting. We provide the comparison of experimental quasiparticle energies for many finite systems with calculations from the GW Green function and localized orbitals scaling correction (LOSC), a recently developed correction to semilocal DFAs, which has minimal delocalization error. Extensive results with over 40 systems clearly show that LOSC orbital energies achieve slightly better accuracy than the GW calculations with little dependence on the semilocal DFA, supporting the use of LOSC DFA orbital energies to predict quasiparticle energies. This also leads to the calculations of excitation energies of the N-electron systems from the ground state DFA calculations of the ( N - 1)-electron systems. Results show good performance with accuracy similar to TDDFT and the delta SCF approach for valence excitations with commonly used DFAs with or without LOSC. For Rydberg states, good accuracy was obtained only with the use of LOSC DFA. This work highlights the pathway to quasiparticle and excitation energies from ground density functional calculations.

Density functional approximations for orbital energies and total energies of molecules and solids.

The relation of Kohn-Sham (KS) orbital energies to ionization energies and electron affinities is different in molecules and solids. In molecules, the local density approximation (LDA) and generalized gradient approximations (GGA) approximate the exact ionization energy (I) and affinity (A) rather well with self-consistently calculated (total energy based) ILDFA and ALDFA, respectively. The highest occupied molecular orbital (HOMO) energy and lowest unoccupied molecular orbital (LUMO) energy, however, differ significantly (by typically 4-6 eV) from these quantities, ϵHLDFA(mol)>-I(mol)≈-ILDFA(mol), ϵLLDFA(mol)<-A(mol)≈-ALDFA(mol). In solids, these relations are very different, due to two effects. The (almost) infinite extent of a solid makes the difference of orbital energies and (L)DFA calculated ionization energy and affinity disappear: in the solid state limit, ϵH(L)DFA(solid)=-I(L)DFA(solid) and ϵL(L)DFA(solid)=-A(L)DFA(solid). Slater's relation ∂E/∂ni = ϵi for local density functional approximations (LDFAs) [and Hartree-Fock (HF) and hybrids] is useful to prove these relations. The equality of LDFA orbital energies and LDFA calculated -ILDFA and -ALDFA in solids does not mean that they are good approximations to the exact quantities. The LDFA total energies of the ions with a delocalized charge are too low, hence ILDFA(solid) < I and ALDFA(solid) > A, due to the local-approximation error, also denoted delocalization error, of LDFAs in extended systems. These errors combine to make the LDFA orbital energy band gap considerably smaller than the exact fundamental gap, ϵLLDFA(solid)-ϵHLDFA(solid)=ILDFA(solid)-ALDFA(solid)<I-A (the LDFA band gap problem). These results for density functional approximations are compared to exact KS and to HF and hybrids. For the exact KS HOMO energy, one has ϵHKS=-I. The exact KS LUMO energy does not approximate the experimental -A (neither in molecules nor in solids), but is considerably below, which is the main reason for the exact KS HOMO-LUMO energy gap being considerably below the fundamental gap I - A (the exact KS band gap problem).

Range-Separated Double-Hybrid Functional from Nonempirical Constraints.

On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.

Assessing GW Approaches for Predicting Core Level Binding Energies

Here we present a systematic study on the performance of different GW approaches: G0W0, G0W0 with linearized quasiparticle equation (lin-G0W0), and quasiparticle self-consistent GW (qsGW), in predicting core level binding energies (CLBEs) on a series of representative molecules comparing to Kohn-Sham (KS) orbital energy-based results. KS orbital energies obtained using the PBE functional are 20-30 eV lower in energy than experimental values obtained from X-ray photoemission spectroscopy (XPS), showing that any Koopmans-like interpretation of KS core level orbitals fails dramatically. Results from qsGW lead to CLBEs that are closer to experimental values from XPS, yet too large. For the qsGW method, the mean absolute error is about 2 eV, an order of magnitude better than plain KS PBE orbital energies and quite close to predictions from ΔSCF calculations with the same functional, which are accurate within ∼1 eV. Smaller errors of ∼0.6 eV are found for qsGW CLBE shifts, again similar to those obtained using ΔSCF PBE. The computationally more affordable G0W0 approximation leads to results less accurate than qsGW, with an error of ∼9 eV for CLBEs and ∼0.9 eV for their shifts. Interestingly, starting G0W0 from PBE0 reduces this error to ∼4 eV with a slight improvement on the shifts as well (∼0.4 eV). The validity of the G0W0 results is however questionable since only linearized quasiparticle equation results can be obtained. The present results pave the way to estimate CLBEs in periodic systems where ΔSCF calculations are not straightforward although further improvement is clearly needed.

Electronic structure and optical properties of boron difluoride naphthaloyl- and anthracenoylacetonates

Electronic structure and optical properties of boron difluoride naphthaloyl- (I) and anthracenoylacetonates (II) have been studied by the methods of photoelectron, absorption, luminescence spectroscopy as well as by quantum chemistry (DFT and TDDFT). In contrast with compound (I), the correlation of energies and types of transitions in excited states between compound (II) and anthracene molecule are ascertained. In addition, the bathochromic shift of the luminescence band maximum and decrease of luminescence quantum yield are detected. The decrease of luminescence quantum yield correlates with the reduction of C-C bond order between chelate ring and substituting group and with the decrease of the excited state lifetime. The good agreement between UV/X-ray photoelectron spectra data and Kohn-Sham orbital energies and structure confirms the reliability of the results of quantum-chemical modeling.

Excitation energies from particle-particle random phase approximation with accurate optimized effective potentials.

The optimized effective potential (OEP) that gives accurate Kohn-Sham (KS) orbitals and orbital energies can be obtained from a given reference electron density. These OEP-KS orbitals and orbital energies are used here for calculating electronic excited states with the particle-particle random phase approximation (pp-RPA). Our calculations allow the examination of pp-RPA excitation energies with the exact KS density functional theory (DFT). Various input densities are investigated. Specifically, the excitation energies using the OEP with the electron densities from the coupled-cluster singles and doubles method display the lowest mean absolute error from the reference data for the low-lying excited states. This study probes into the theoretical limit of the pp-RPA excitation energies with the exact KS-DFT orbitals and orbital energies. We believe that higher-order correlation contributions beyond the pp-RPA bare Coulomb kernel are needed in order to achieve even higher accuracy in excitation energy calculations.

Optimal Tuning of Range-Separated Hybrids for Solvated Molecules with Time-Dependent Density Functional Theory.

The applicability range of density functional theory (DFT) can be improved with no additional parametrization by imposing some exact conditions. Enforcing equality between the orbital energy of the highest occupied Kohn-Sham orbital and ionization energy, determined from the total energy difference between neutral and ionized states (ΔKS), leads to the concept of optimally tuned range-separated hybrid functionals. Here, we present an alternative tuning scheme for range-separated hybrid functionals based on enforcing the equality between the ΔKS ionization energy and the ionization energy calculated by means of the time-dependent DFT using the concept of ionization as an excitation to the distant center (OT-IEDC scheme). The scheme can be naturally applied to solvated systems described either within the explicit solvation or dielectric continuum models. We test the performance of the scheme on a benchmark set of molecules. We further show that the scheme allows for reliably modeling liquid phase photoemission spectra.

Performance of the G0W0 Method in Predicting the Electronic Gap of TiO2 Nanoparticles.

Using a relativistic all-electron description and numerical atomic-centered orbital basis set, the performance of the G0W0 method on the electronic band gap of (TiO2)n nanoparticles (n = 1-20) is investigated. Results are presented for G0W0 on top of hybrid (PBE0 and a modified version with 12.5% of Fock exchange) functionals. The underestimation of the electronic band gap from Kohn-Sham orbital energies is corrected by the quasiparticle energies from the G0W0 method, which are consistent with the variational ΔSCF approach. A clear correlation between both methods exists regardless of the hybrid functional employed. In addition, the vertical ionization potential and electron affinity from quasiparticle energies show a systematic correlation with the ΔSCF calculated values. On the other hand, the shape of the nanoparticles promotes some deviations on the electronic band gap. In conclusion, this study shows the following: (i) A systematic correlation exists between band gaps, ionization potentials, and electron affinities of TiO2 nanoparticles as predicted from variational ΔSCF and G0W0 methods. (ii) The G0W0 approach can be successfully used to study the electronic band gap of realistic size nanoparticles at an affordable computational cost with a ΔSCF accuracy giving results that are directly related with those from photoemission spectroscopy. (iii) The quasiparticle energies are explicitly required to shed light on the photocatalytic properties of TiO2. (iv) The G0W0 approach emerges as an accurate method to investigate the photocatalytic properties of both nanoparticles and extended semiconductors.

Density Functional Approach and Random Matrix Theory in Proteogenesis

We study the energy-level statistics of amino acids by random matrix theory. The molecular orbital and the Kohn–Sham orbital energies are calculated using ab initio and density-functional formalisms for 20 different amino acids. To generate statistical data, we performed a multipoint calculation on 10000 molecular structures produced via a molecular dynamics simulation. For the valence orbitals, the energy-level statistics exhibit repulsion, but the universality in the random matrix cannot be determined. For the unoccupied orbitals, the energy-level statistics indicate an intermediate distribution between the Gaussian orthogonal ensemble and the semi-Poisson statistics for all 20 different amino acids. These amino acids are considered to be in a type of critical state.

Systematic treatment of spin-reactivity indicators in conceptual density functional theory

The spin-resolved conceptual density functional theory is presented, emphasizing the finite-difference approximations for spin-reactivity indicators and especially the spin-hardness. We derive simple approximations for these spin-reactivity indicators in terms of Kohn–Sham orbital energies and discuss their accuracy and validity. We also show that the second derivative of the energy with respect to the spin-abundance, N S = ½M S , is nonnegative, contradicting the generally negative values obtained for one of the most popular formulas used to approximate this quantity. A new approximation for the spin-hardness is proposed and assessed based on information for spin-philicities and singlet–triplet gaps. One conclusion of our work is that the tendency in the literature to mix linear models (e.g., one-sided derivatives for chemical potentials and Fukui functions) with quadratic models (e.g., two-sided derivatives for hardnesses and dual descriptors) not only is theoretically unsound, but often leads to chemically absurd results.