Simple iterative construction of the optimized effective potential for orbital functionals, including exact exchange

Comparison of optimized effective potential with inverse Kohn-Sham method for Hartree-Fock exchange energy

The inverse Kohn-Sham (inv-KS) density-functional theory for the electron density of the Hartree-Fock (HF) wave function was revisited within the context of the optimized effective potential (HF-OEP). First, we clarify the relationship between the inv-KS and the HF-OEP within the framework of the potential-functional theory. The similarities and the differences of the approaches are then discussed on the basis of their methodological details, which motivates comparisons of the wave function provided by each method. Next, the real-space grid implementations of the inv-KS and the HF-OEP are addressed for the comparisons. The total HF energies EHF[{φiinv-KS}] for the wave functions φiinv-KS on the effective potentials optimized by the inv-KS are computed for a set of small molecules. It is found that the mean absolute deviation (MAD) of EHF[{φiinv-KS}] from the HF energy is clearly smaller than the MAD of EHF[{φiOEP}], demonstrating that the inv-KS is advantageous in constructing the detailed structure of the exchange potential υx as compared with the HF-OEP. The inv-KS method is also applied to an ortho-benzyne radical known as a strongly correlated polyatomic molecule. It is revealed that the spin populations on the atomic sites computed by the UHF calculation can be faithfully reproduced by the wave functions on the inv-KS potential.

Density functional theory from spherically symmetric densities: Ground and excited states of Coulomb systems

Recently, Theophilou [J. Chem. Phys. 149, 074104 (2018)] proposed a peculiar version of the density functional theory by showing that the set of spherical averages of the density around the nuclei determines uniquely the external potential in atoms, molecules, and solids. Here, this novel theory is extended to individual excited states. The generalization is based on the method developed in the series of papers by Ayers, Levy, and Nagy [Phys. Rev. A 85, 042518 (2012)]. Generalized Hohenberg-Kohn theorems are proved to the set of spherically symmetric densities using constrained search. A universal variational functional for the sum of the kinetic and electron-electron repulsion energies is constructed. The functional is appropriate for the ground state and all bound excited states. Euler equations and Kohn-Sham equations for the set are derived. The Euler equations can be rewritten as Schrödinger-like equations for the square root of the radial densities, and the effective potentials in them can be expressed in terms of wave function expectation values. The Hartree plus exchange-correlation potentials can be given by the difference of the interacting and the non-interacting effective potentials.

Performance of Density Functionals and Semiempirical 3c Methods for Small Gold-Thiolate Clusters

In light of the recent surge in computational studies of gold thiolate clusters, we present a comparison of popular density functionals (DFAs) and three-part corrected methods (3c-methods) on their performance by taking a data set named as AuSR18 consisting of 18 isomers of Au_n(SCH₃)_m (m ≤ n = 1-3). We have compared the efficiency and accuracy of the DFAs and 3c-methods in geometry optimization with RI-SCS-MP2 as the reference method. Similarly, the performance for accurate and efficient energy evaluation was compared with DLPNO-CCSD(T) as the reference method. The lowest energy structure among the isomers of the largest stoichiometry from our data set, AuSR18, i.e., Au₃(SCH₃)₃, is considered to evaluate the computational time for SCF and gradient evaluations. Alongside this, the numbers of optimization steps to locate the most stable minima of Au₃(SCH₃)₃ are compared to assess the efficiency of the methods. A comparison of relevant bond lengths with the reference geometries was made to estimate the accuracy in geometry optimization. Some methods, such as LC-BLYP, ωB97M-D3BJ, M06-2X, and PBEh-3c, could not locate many of the minima found by most of the other methods; thus, the versatility in locating various minima is also an important criterion in choosing a method for the given project. To determine the accuracy of the methods, we compared the relative energies of the isomers in each stoichiometry and the interaction energy of the gold core with the ligands. The dependence of basis set size and relativistic effects on energies are also compared. The following are some of the highlights. TPSS has shown accuracy, while mPWPW shows comparable speed and accuracy. For the relative energies of the clusters, the hybrid range-separated DFAs are the best option. CAM-B3LYP excels, whereas B3LYP performs poorly. Overall, LC-BLYP is a balanced performer considering both the geometry and relative stability of the structures, but it lacks diversity. The 3c-methods, although fast, are less impressive in relative stability.

Scaled σ-functionals for the Kohn-Sham correlation energy with scaling functions from the homogeneous electron gas

The recently introduced σ-functionals constitute a new type of functionals for the Kohn-Sham (KS) correlation energy. σ-Functionals are based on the adiabatic-connection fluctuation-dissipation theorem, are computationally closely related to the well-known direct random phase approximation (dRPA), and are formally rooted in many-body perturbation theory along the adiabatic connection. In σ-functionals, the function of the eigenvalues σ of the Kohn-Sham response matrix that enters the coupling constant and frequency integration in the dRPA is replaced by another function optimized with the help of reference sets of atomization, reaction, transition state, and non-covalent interaction energies. σ-Functionals are highly accurate and yield chemical accuracy of 1 kcal/mol in reaction or transition state energies, in main group chemistry. A shortcoming of σ-functionals is their inability to accurately describe processes involving a change of the electron number, such as ionizations or electron attachments. This problem is attributed to unphysical self-interactions caused by the neglect of the exchange kernel in the dRPA and σ-functionals. Here, we tackle this problem by introducing a frequency- and σ-dependent scaling of the eigenvalues σ of the KS response function that models the effect of the exchange kernel. The scaling factors are determined with the help of the homogeneous electron gas. The resulting scaled σ-functionals retain the accuracy of their unscaled parent functionals but in addition yield very accurate ionization potentials and electron affinities. Moreover, atomization and total energies are found to be exceptionally accurate. Scaled σ-functionals are computationally highly efficient like their unscaled counterparts.

Numerically stable inversion approach to construct Kohn-Sham potentials for given electron densities within a Gaussian basis set framework

We present a Kohn-Sham (KS) inversion approach to construct KS exchange-correlation potentials corresponding to given electron densities. This method is based on an iterative procedure using linear response to update potentials. All involved quantities, i.e., orbitals, potentials, and response functions, are represented by Gaussian basis functions. In contrast to previous KS inversion methods relying on Gaussian basis sets, the method presented here is numerically stable even for standard basis sets from basis set libraries due to a preprocessing of the auxiliary basis used to represent an exchange-correlation charge density that generates the exchange-correlation potential. The new KS inversion method is applied to reference densities of various atoms and molecules obtained by full configuration interaction or CCSD(T) (coupled cluster singles doubles perturbative triples). The considered examples encompass cases known to be difficult, such as stretched hydrogen or lithium hydride molecules or the beryllium isoelectronic series. For the stretched hydrogen molecule, potentials of benchmark quality are obtained by employing large basis sets. For the carbon monoxide molecule, we show that the correlation potential from the random phase approximation (RPA) is in excellent qualitative and quantitative agreement with the correlation potential from the KS inversion of a CCSD(T) reference density. This indicates that RPA correlation potentials, in contrast to those from semi-local density-functionals, resemble the exact correlation potential. Besides providing exchange-correlation potentials for benchmark purposes, the proposed KS inversion method may be used in density-partition-based quantum embedding and in subsystem density-functional methods because it combines numerical stability with computational efficiency.

Numerically stable optimized effective potential method with standard Gaussian basis sets

We present a numerically stable optimized effective potential (OEP) method based on Gaussian basis sets. The key point of the approach is a sequence of preprocessing steps of the auxiliary basis set used to represent exchange or correlation potentials, the Kohn-Sham (KS) response function, and the right-hand side of the OEP equation in conjunction with a representation of exchange or correlation potentials via exchange or correlation charge densities whose electrostatic potentials generate the potentials. Due to the preprocessing, standard Gaussian basis sets from basis set libraries can be used in OEP calculations. As examples, we present numerical stable computational setups based on aux-cc-pwCVXZ basis sets with X = T, Q, 5 for the orbitals and aux-cc-pVDZ/mp2fit and aux-cc-pVTZ/mp2fit auxiliary basis sets and use them to calculate KS exchange potentials with the exact exchange-only KS method for various atoms and molecules. The resulting exchange potentials not only are numerically stable and physically reasonable but also show convergence with increasing quality of the orbital basis sets. The effect of incorporating exact conditions that the KS exchange potential has to obey is discussed. Moreover, it is briefly demonstrated that the presented approach not only works for KS exchange potentials but equally well for correlation potentials within the direct random phase approximation. Besides for OEP methods, the introduced preprocessing of auxiliary basis sets should also be beneficial in procedures to calculate back effective KS potentials from given electron densities.

Optimized effective potentials from the random-phase approximation: Accuracy of the quasiparticle approximation

The optimized effective potential (OEP) method presents an unambiguous way to construct the Kohn-Sham potential corresponding to a given diagrammatic approximation for the exchange-correlation functional. The OEP from the random-phase approximation (RPA) has played an important role ever since the conception of the OEP formalism. However, the solution of the OEP equation is computationally fairly expensive and has to be done in a self-consistent way. So far, large scale solid state applications have, therefore, been performed only using the quasiparticle approximation (QPA), neglecting certain dynamical screening effects. We obtain the exact RPA-OEP for 15 semiconductors and insulators by direct solution of the linearized Sham-Schlüter equation. We investigate the accuracy of the QPA on Kohn-Sham bandgaps and dielectric constants, and comment on the issue of self-consistency.

A range-separated generalized Kohn-Sham method including a long-range nonlocal random phase approximation correlation potential

Based on our recently published range-separated random phase approximation (RPA) functional [Kreppel et al., "Range-separated density-functional theory in combination with the random phase approximation: An accuracy benchmark," J. Chem. Theory Comput. 16, 2985-2994 (2020)], we introduce self-consistent minimization with respect to the one-particle density matrix. In contrast to the range-separated RPA methods presented so far, the new method includes a long-range nonlocal RPA correlation potential in the orbital optimization process, making it a full-featured variational generalized Kohn-Sham (GKS) method. The new method not only improves upon all other tested RPA schemes including the standard post-GKS range-separated RPA for the investigated test cases covering general main group thermochemistry, kinetics, and noncovalent interactions but also significantly outperforms the popular G₀W₀ method in estimating the ionization potentials and fundamental gaps considered in this work using the eigenvalue spectra obtained from the GKS Hamiltonian.

Introductory lecture: when the density of the noninteracting reference system is not the density of the physical system in density functional theory

A major challenge in density functional theory (DFT) is the development of density functional approximations (DFAs) to overcome errors in existing DFAs, leading to more complex functionals. For such functionals, we consider roles of the noninteracting reference systems. The electron density of the Kohn-Sham (KS) reference with a local potential has been traditionally defined as being equal to the electron density of the physical system. This key idea has been applied in two ways: the inverse calculation of such a local KS potential for the reference from a given density and the direct calculation of density and energy based on given DFAs. By construction, the inverse calculation can yield a KS reference with the density equal to the input density of the physical system. In application of DFT, however, it is the direct calculation of density and energy from a DFA that plays a central role. For direct calculations, we find that the self-consistent density of the KS reference defined by the optimized effective potential (OEP), is not the density of the physical system, when the DFA is dependent on the external potential. This inequality holds also for the density of generalized KS (GKS) or generalized OEP reference, which allows a nonlocal potential, when the DFA is dependent on the external potential. Instead, the density of the physical system, consistent with a given DFA, is given by the linear response of the total energy with respect to the variation of the external potential. This is a paradigm shift in DFT on the use of noninteracting references: the noninteracting KS or GKS references represent the explicit computational variables for energy minimization, but not the density of the physical system for external potential-dependent DFAs. We develop the expressions for the electron density so defined through the linear response for general DFAs, demonstrate the results for orbital functionals and for many-body perturbation theory within the second-order and the random-phase approximation, and explore the connections to developments in DFT.

Lieb-Oxford bound and pair correlation functions for density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem

Compliance with the Lieb-Oxford bound for the indirect Coulomb energy and for the exchange-correlation energy is investigated for a number of density-functional methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem to treat correlation. Furthermore, the correlation contribution to the pair density resulting from these methods is compared with highly accurate reference values for the helium atom and for the hydrogen molecule at several bond distances. For molecules, the Lieb-Oxford bound is obeyed by all considered methods. For the homogeneous electron gas, it is violated by all methods for low electron densities. The simplest considered ACFD method, the direct random phase approximation (dRPA), violates the Lieb-Oxford bound much earlier than more advanced ACFD methods that, in addition to the simple Hartree kernel, take into account the exchange kernel and an approximate correlation kernel in the calculation of the correlation energy. While the dRPA yields quite poor correlation contributions to the pair density, those from more advanced ACFD methods are physically reasonable but still leave room for improvements, particularly in the case of the stretched hydrogen molecule.

First Ionization Energy as the Asymptotic Limit of the Average Local Electron Energy

The first vertical ionization energy of an atom or molecule is encoded in the rate of exponential decay of the exact natural orbitals. For natural orbitals represented in terms of Gaussian basis functions, this property does not hold even approximately. We show that it is nevertheless possible to deduce the first ionization energy from the long-range behavior of Gaussian-basis-set wave functions by evaluating the asymptotic limit of a quantity called the average local electron energy (ALEE), provided that the most diffuse functions of the basis set have a suitable shape and location. The ALEE method exposes subtle qualitative differences between seemingly analogous Gaussian basis sets and complements the extended Koopmans theorem by being robust in situations where the one-electron reduced density matrix is ill-conditioned.

Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems

Over the last few years, extraordinary advances in experimental and theoretical tools have allowed us to monitor and control matter at short time and atomic scales with a high degree of precision. An appealing and challenging route toward engineering materials with tailored properties is to find ways to design or selectively manipulate materials, especially at the quantum level. To this end, having a state-of-the-art ab initio computer simulation tool that enables a reliable and accurate simulation of light-induced changes in the physical and chemical properties of complex systems is of utmost importance. The first principles real-space-based Octopus project was born with that idea in mind, i.e., to provide a unique framework that allows us to describe non-equilibrium phenomena in molecular complexes, low dimensional materials, and extended systems by accounting for electronic, ionic, and photon quantum mechanical effects within a generalized time-dependent density functional theory. This article aims to present the new features that have been implemented over the last few years, including technical developments related to performance and massive parallelism. We also describe the major theoretical developments to address ultrafast light-driven processes, such as the new theoretical framework of quantum electrodynamics density-functional formalism for the description of novel light-matter hybrid states. Those advances, and others being released soon as part of the Octopus package, will allow the scientific community to simulate and characterize spatial and time-resolved spectroscopies, ultrafast phenomena in molecules and materials, and new emergent states of matter (quantum electrodynamical-materials).

Ab Initio Optimized Effective Potentials for Real Molecules in Optical Cavities: Photon Contributions to the Molecular Ground State

We introduce a simple scheme to efficiently compute photon exchange-correlation contributions due to the coupling to transversal photons as formulated in the newly developed quantum-electrodynamical density-functional theory (QEDFT).1-5 Our construction employs the optimized-effective potential (OEP) approach by means of the Sternheimer equation to avoid the explicit calculation of unoccupied states. We demonstrate the efficiency of the scheme by applying it to an exactly solvable GaAs quantum ring model system, a single azulene molecule, and chains of sodium dimers, all located in optical cavities and described in full real space. While the first example is a two-dimensional system and allows to benchmark the employed approximations, the latter two examples demonstrate that the correlated electron-photon interaction appreciably distorts the ground-state electronic structure of a real molecule. By using this scheme, we not only construct typical electronic observables, such as the electronic ground-state density, but also illustrate how photon observables, such as the photon number, and mixed electron-photon observables, for example, electron-photon correlation functions, become accessible in a density-functional theory (DFT) framework. This work constitutes the first three-dimensional ab initio calculation within the new QEDFT formalism and thus opens up a new computational route for the ab initio study of correlated electron-photon systems in quantum cavities.

Time-Local Equation for the Exact Optimized Effective Potential in Time-Dependent Density Functional Theory

A long-standing challenge in the time-dependent density functional theory is to efficiently solve the exact time-dependent optimized effective potential (TDOEP) integral equation derived from orbital-dependent functionals, especially for the study of nonadiabatic dynamics in time-dependent external fields. In this Letter, we formulate a completely equivalent time-local TDOEP equation that admits a unique real-time solution in terms of time-dependent Kohn-Sham and effective memory orbitals. The time-local formulation is numerically implemented, with the incorporation of exponential memory loss to address the unaccounted for correlation component in the exact-exchange-only functional, to enable the study of the many-electron dynamics of a one-dimensional hydrogen chain. It is shown that the long time behavior of the electric dipole converges correctly and the zero-force theorem is fulfilled in the current implementation.

Atoms and molecules in cavities, from weak to strong coupling in quantum-electrodynamics (QED) chemistry

In this work, we provide an overview of how well-established concepts in the fields of quantum chemistry and material sciences have to be adapted when the quantum nature of light becomes important in correlated matter-photon problems. We analyze model systems in optical cavities, where the matter-photon interaction is considered from the weak- to the strong-coupling limit and for individual photon modes as well as for the multimode case. We identify fundamental changes in Born-Oppenheimer surfaces, spectroscopic quantities, conical intersections, and efficiency for quantum control. We conclude by applying our recently developed quantum-electrodynamical density-functional theory to spontaneous emission and show how a straightforward approximation accurately describes the correlated electron-photon dynamics. This work paves the way to describe matter-photon interactions from first principles and addresses the emergence of new states of matter in chemistry and material science.

A local Fock-exchange potential in Kohn-Sham equations

We derive and employ a local potential to represent the Fock exchange operator in electronic single-particle equations. This local Fock-exchange (LFX) potential is very similar to the exact exchange (EXX) potential in density functional theory (DFT). The practical software implementation of the two potentials (LFX and EXX) yields robust and accurate results for a variety of systems (semiconductors, transition metal oxides) where Hartree-Fock and popular approximations of DFT typically fail. This includes examples traditionally considered qualitatively inaccessible to calculations that omit correlation.

Real-space grids and the Octopus code as tools for the development of new simulation approaches for electronic systems

Real-space grids are a powerful alternative for the simulation of electronic systems. One of the main advantages of the approach is the flexibility and simplicity of working directly in real space where the different fields are discretized on a grid, combined with competitive numerical performance and great potential for parallelization. These properties constitute a great advantage at the time of implementing and testing new physical models. Based on our experience with the Octopus code, in this article we discuss how the real-space approach has allowed for the recent development of new ideas for the simulation of electronic systems. Among these applications are approaches to calculate response properties, modeling of photoemission, optimal control of quantum systems, simulation of plasmonic systems, and the exact solution of the Schrödinger equation for low-dimensionality systems.

Optimized Effective Potential for Quantum Electrodynamical Time-Dependent Density Functional Theory

We propose an orbital exchange-correlation functional for applying time-dependent density functional theory to many-electron systems coupled to cavity photons. The time nonlocal equation for the electron-photon optimized effective potential (OEP) is derived. In the static limit our OEP energy functional reduces to the Lamb shift of the ground state energy. We test the new approximation in the Rabi model. It is shown that the OEP (i) reproduces quantitatively the exact ground-state energy from the weak to the deep strong coupling regime and (ii) accurately captures the dynamics entering the ultrastrong coupling regime. The present formalism opens the path to a first-principles description of correlated electron-photon systems, bridging the gap between electronic structure methods and quantum optics for real material applications.

Optimized effective potential method and application to static RPA correlation

The optimized effective potential (OEP) method is a promising technique for calculating the ground state properties of a system within the density functional theory. However, it is not widely used as its computational cost is rather high and, also, some ambiguity remains in the theoretical framework. In order to overcome these problems, we first introduced a method that accelerates the OEP scheme in a static RPA-level correlation functional. Second, the Krieger-Li-Iafrate (KLI) approximation is exploited to solve the OEP equation. Although seemingly too crude, this approximation did not reduce the accuracy of the description of the magnetic transition metals (Fe, Co, and Ni) examined here, the magnetic properties of which are rather sensitive to correlation effects. Finally, we reformulated the OEP method to render it applicable to the direct RPA correlation functional and other, more precise, functionals. Emphasis is placed on the following three points of the discussion: (i) level-crossing at the Fermi surface is taken into account; (ii) eigenvalue variations in a Kohn-Sham functional are correctly treated; and (iii) the resultant OEP equation is different from those reported to date.

Fock-exchange for periodic structures in the real-space formalism and the KLI approximation

The calculation of Fock-exchange interaction is an important task in the computation of molecule and solid properties.

Accurate and efficient approximation to the optimized effective potential for exchange

We devise an efficient practical method for computing the Kohn-Sham exchange-correlation potential corresponding to a Hartree-Fock electron density. This potential is almost indistinguishable from the exact-exchange optimized effective potential (OEP) and, when used as an approximation to the OEP, is vastly better than all existing models. Using our method one can obtain unambiguous, nearly exact OEPs for any reasonable finite one-electron basis set at the same low cost as the Krieger-Li-Iafrate and Becke-Johnson potentials. For all practical purposes, this solves the long-standing problem of black-box construction of OEPs in exact-exchange calculations.

Piecewise linearity of approximate density functionals revisited: implications for frontier orbital energies

In the exact Kohn-Sham density-functional theory, the total energy versus the number of electrons is a series of linear segments between integer points. However, commonly used approximate density functionals produce total energies that do not exhibit this piecewise-linear behavior. As a result, the ionization potential theorem, equating the highest occupied eigenvalue with the ionization potential, is grossly disobeyed. Here, we show that, contrary to conventional wisdom, most of the required piecewise linearity of an arbitrary approximate density functional can be restored by careful consideration of the ensemble generalization of density-functional theory. Furthermore, the resulting formulation introduces the desired derivative discontinuity to any approximate exchange-correlation functional, even one that is explicitly density dependent. This opens the door to calculations of the ionization potential and electron affinity, even without explicit electron removal or addition. All these advances are achieved while neither introducing empiricism nor changing the underlying functional form. The power of the approach is demonstrated on benchmark systems using the local density approximation as an illustrative example.

Time-dependent density-functional theory in massively parallel computer architectures: the OCTOPUS project

Octopus is a general-purpose density-functional theory (DFT) code, with a particular emphasis on the time-dependent version of DFT (TDDFT). In this paper we present the ongoing efforts to achieve the parallelization of octopus. We focus on the real-time variant of TDDFT, where the time-dependent Kohn-Sham equations are directly propagated in time. This approach has great potential for execution in massively parallel systems such as modern supercomputers with thousands of processors and graphics processing units (GPUs). For harvesting the potential of conventional supercomputers, the main strategy is a multi-level parallelization scheme that combines the inherent scalability of real-time TDDFT with a real-space grid domain-partitioning approach. A scalable Poisson solver is critical for the efficiency of this scheme. For GPUs, we show how using blocks of Kohn-Sham states provides the required level of data parallelism and that this strategy is also applicable for code optimization on standard processors. Our results show that real-time TDDFT, as implemented in octopus, can be the method of choice for studying the excited states of large molecular systems in modern parallel architectures.

ALGEBRAIC EQUATION AND ITERATIVE OPTIMIZATION FOR THE OPTIMIZED EFFECTIVE POTENTIAL IN DENSITY FUNCTIONAL THEORY

We further develop our recent direct method for the optimized effective potential (OEP) in density functional theory (DFT) [Yang and Wu, Phys. Rev. Lett.89, 143002 (2002)]. First, we show that the stationary condition in our optimization approach leads to a proper nonlinear algebraic equation for the OEP in a finite basis set, which differs from other finite basis set approaches. Then by constructing an approximate second derivative matrix of the energy functional in conjunction with the use of the Newton method, we significantly accelerate the convergence of the iterative optimization for OEP. Enhancement of the method is made in using the Tikhonov regularization method for the inversion of the second derivative matrix when it is singular or nearly singular and the direct inversion in the iterative space. It is shown that under a fixed stepsize condition, the optimization approach is equivalent to the self-consistent solution to the nonlinear algebraic equation for OEP. Because the approximate second derivatives are easy to compute and the iteration numbers are small now, the computation costs of OEP become comparable to that of regular DFT calculations as shown by calculations of some molecules, small and larger ones. We show how to find balanced results between energies and potentials when choosing a basis set for potentials.

Density-Functional Theory with Orbital-Dependent Functionals: Exact-exchange Kohn-Sham and Density-Functional Response Methods

Kohn-Sham methods with orbital-dependent exchange-correlation functionals, in particular exact-exchange methods, are discussed. The numerical stability of optimized-effective-potential methods to construct the exchange-correlation potential in the case of orbital-dependent energy functionals is considered. A numerically stable exact-exchange Kohn-Sham method is briefly presented. A new ansatz for time-dependent density-functional methods in the response regime is considered and an exact-exchange time-dependent density-functional method based on it is discussed and it is demonstrated that this methods is able to treat qualitatively correctly charge-transfer excitations.

Self-interaction Free Relativistic Spin-density Functional Theory

We review the progress in the recent development of the relativistic optimized effective potential (ROEP) method within spin-density functional theory. The ROEP-equations for spin-polarized systems are derived and their application to open-shell atoms using the exact exchange approximation for the exchange correlation-functional is presented. Further, we expand the ROEP framework to treat extended systems within the KKR-multiple scattering formalism. We illustrate the application of the theory to open-shell free atoms and solids.

Orbital-dependent Representation of Correlation Energy Functional

Exchange-correlation energy functionals depending on the Kohn-Sham orbitals and eigenvalues resolve some of the well-known deficiencies of the local density and generalized gradient approximations. Such functionals can be derived in first-principles fashion by use of standard many-body techniques, if the Kohn-Sham Hamiltonian is utilized as non-interacting reference Hamiltonian. In this way one can establish an exact relation for the exchange-correlation functional, which opens several routes for the derivation of approximate functionals. Straightforward expansion in powers of the electron-electron coupling constant gives, to first order, the exact exchange of density functional theory and, to second order, a correlation functional which has the same structure as the second order Møller-Plesset term. This simplest first-principles correlation functional reproduces both the shell structure in the exact correlation potential and dispersion forces. On the other hand, it overestimates all correlation effects and is variationally instable for systems with a very small HOMO-LUMO gap. Both deficiences can be resolved by partial resummation of the Kohn-Sham perturbation expansion. In this contribution two such resummations are discussed, a minimum form designed to address the variational instability at no computational cost (based on the Epstein-Nesbet diagrams) and a more systematic variant (based on the ring diagrams).

Direct mapping between exchange potentials of Hartree-Fock and Kohn-Sham schemes as origin of orbital proximity

It is found that for closed-l-shell atoms, the exact local exchange potential v(x)(r) calculated in the exchange-only Kohn-Sham (KS) scheme of the density functional theory (DFT) is very well represented within the region of every atomic shell by each of the suitably shifted potentials obtained with the nonlocal Fock exchange operator for the individual Hartree-Fock (HF) orbitals belonging to this shell. This newly revealed property is not related to the well-known steplike shell structure in the response part of v(x)(r), but it results from specific relations satisfied by the HF orbital exchange potentials. These relations explain the outstanding proximity of the occupied HF and exchange-only KS orbitals as well as the high quality of the Krieger-Li-Iafrate and localized HF (or, equivalently, common-energy-denominator) approximations to the DFT exchange potential v(x)(r). Another highly accurate representation of v(x)(r) is given by the continuous piecewise function built of shell-specific exchange potentials, each defined as the weighted average of the shifted orbital exchange potentials corresponding to a given shell. The constant shifts added to the HF orbital exchange potentials, to map them onto v(x)(r), are nearly equal to the differences between the energies of the corresponding KS and HF orbitals. It is discussed why these differences are positive and grow when the respective orbital energies become lower for inner orbitals.