Optimized Effective Potential Research Articles

Avoiding spin contamination and spatial symmetry breaking by exact-exchange-only optimized-effective-potential methods within the symmetrized Kohn-Sham framework.

For open-shell atoms and molecules, Kohn-Sham (KS) methods typically resort to spin-polarized approaches that exhibit spin-contamination and often break spatial symmetries. As a result, the KS Hamiltonian operator and the KS orbitals do not exhibit the space and spin symmetry of the physical electron system. The KS formalism can be symmetrized in a rigorous way only in real space, only in spin space, or both in real and spin space. Within such symmetrized KS frameworks, we present exact-exchange-only optimized-effective-potential (OEP) methods that are free of spin contamination and/or spatial symmetry breaking. The effect of symmetrizations on the total energy and its parts and on the exchange potential is analyzed. The presented exact-exchange-only OEP methods may serve as a starting point for high-level symmetrized KS methods based, e.g., on the adiabatic-connection fluctuation-dissipation theorem.

Time-Dependent Density-Functional Theory for Determining the Electron-Capture Cross Section for Protons Impacting on Atoms and Molecules.

The use of the Time-Dependent Density-Functional Theory (TDDFT) has increased in the atomic collision field. Calculating the electron-capture cross section (ECCS) for protons is an important question in hadrontherapy and plasma physics, among other areas. In previous studies, it was shown that the approach based on the Local Density Approximation (LDA) fails in the 1-50 keV region, requiring the use of the Optimized Effective Potential (OEP) method. In this work, the ECCS values for 1-50 keV protons impacting on isolated hydrogen, carbon, nitrogen, oxygen, and nitrogenous atoms were determined using the TDDFT. It is shown that adding the Self Interaction Correction to the LDA (LDA-Sic) allows obtaining results close to those provided by the OEP and experiments, with the advantage that the LDA-Sic consumes less computational time. In addition, it was demonstrated that it is imperative to include the spin correction for the specific helium and oxygen cases, in order to get good results for the ECCS using the TDDFT. Theoretical results obtained in this work show excellent agreement with experimental values.

Analytical Forces for the Optimized Effective Potential Calculations.

The optimized effective potential (OEP) equation is an ill-conditioned linear system when using finite basis sets. Without any special treatment, the obtained exchange-correlation (XC) potential may have unphysical oscillations. One way to alleviate this problem is to regularize the solutions; however, a regularized XC potential is not the exact solution to the OEP equation. As a result, the system's energy is no longer variational against the Kohn-Sham (KS) potential, and the analytical forces cannot be derived from the Hellmann-Feynman theorem. In this work, we develop a robust and nearly black-box OEP method to ensure that the system's energy is variational against the KS potential. The basic idea is to add a penalty function that regularizes the XC potential to the energy functional. Analytical forces can then be derived based on the Hellmann-Feynman theorem. Another key result is that the impact of the regularization can be much reduced by regularizing the difference between the XC potential and an approximate XC potential rather than regularizing the XC potential. Numerical tests show that forces and the energy differences between systems are not sensitive to the regularization coefficient, which indicates that in practice accurate structural and electronic properties can be obtained without extrapolating the regularization coefficient to zero. We expect this new method to be found useful for calculations that employ advanced, orbital-based functionals, especially for applications that require efficient force calculations.

Dissipation and spontaneous emission in quantum electrodynamical density functional theory based on optimized effective potential: A proof of concept study

We generalize the optimized effective potential (OEP) formalism in the quantum electrodynamical density functional theory (QEDFT) to the case of continuous distribution of photon modes, and study its applicability to dissipative dynamics of electron systems interacting with photons of lossy cavities. Specifically, we test whether this technique is capable of capturing the quantum features of electron-photon interaction related to spontaneous emission and the corresponding energy transfer from the electrons to cavity photons. For this purpose, we analyze a discrete three-site system with one electron coupled to photons of the cavity, which, in fact, is a minimal model allowing to eliminate classical radiation and the corresponding energy loss, but still have nontrivial density dynamics. By considering two typical spectral densities of photon modes, modeling (i) lossy cavity with Lorentzian broadening of photon peaks, and (ii) the Ohmic bath, and several representative dynamical regimes, we find that OEP-QEDFT demonstrates a good qualitative and quantitative performance, especially in the case when the disspation is dominated by one-photon processes.

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.

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self-interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-L\"owdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSIC-KLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Finally the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

Methods to generate reference total and Pauli kinetic potentials

We have derived a new method which allows to compute the full and the Pauli reference kinetic potentials for atoms and molecules in a real space representation. This is done by applying the optimized effective potential (OEP) method to {the} Kohn-Sham non-interacting kinetic energy expression. Additionally, we have also derived a simplified OEP variant based on the common energy denominator approximation which has proven to give much more stable and robust results than the original OEP one. Moreover, we have also proved that at the solution point our approach is formally equivalent to the commonly used Bartolotti-Acharya formula.

Analytical energy gradient for the embedded cluster density approximation.

We recently developed the embedded cluster density approximation (ECDA), which is a local correlation method for scaling up Kohn-Sham (KS) density functional theory calculations using high-level exchange-correlation (XC) functionals. In ECDA, a system's XC energy is obtained by patching locally calculated, high-level XC energy densities over the entire system. Our previous formulation of ECDA is not variational, making it difficult to derive the analytical energy gradient. In this work, we present a fully variational formulation of ECDA and derive the analytical energy gradient. The challenge for making ECDA a variational method is that both partitioning the system's density and solving the system's XC potential are the optimized effective potential (OEP) problems. Simply regularizing these two OEP equations makes ECDA a nonvariational method. We show how to regularize these two OEP problems while still keeping ECDA variational. KS linear responses are involved in the calculations of the system's XC potential and the analytical energy gradients, but are not explicitly constructed. The terms involving the KS linear responses are calculated by solving the Sternheimer equation. The analytical energy gradients are validated with a Si2H6 molecule and are used to relax the geometry of Si6H10. In both examples, the exact exchange is used as the high-level XC functional and is patched over the molecules.

Strong Correlation and Charge Localization in Kohn-Sham Theories with Fractional Orbital Occupations.

We study static correlation and delocalization errors and show that even methods with good energies can yield significant delocalization errors that affect the density, leading to large errors in predicting, e.g., dipole moments. We illustrate this point by comparing existing state-of-the-art approaches with an accurate exchange-correlation functional based on a generalized valence-bond ansatz in which orbitals and fractional occupations are treated as variational parameters via an optimized effective potential (OEP). We show that the OEP exhibits step and peak features that, similar to the exact Kohn-Sham potential of density functional theory (DFT), are crucial to prevent charge delocalization. We further show that the step is missing in common approximations within reduced density matrix functional theory, resulting in delocalization errors comparable to those found in DFT approximations. Finally, we explain the delocalization error as coming from artificial mixing of the ground state with a charge-transfer excited state that is avoided if the occupation numbers exhibit discontinuities.

Potential Functional Embedding Theory with an Improved Kohn-Sham Inversion Algorithm.

Potential functional embedding theory (PFET) is a rigorous theory that can yield a unique, self-consistent embedding potential shared by different subsystems treated at different levels of theory. Application of PFET has been limited by the time-consuming and sometimes unstable optimized effective potential (OEP) procedure. Here, we improve the performance of PFET by replacing the OEP algorithm with a new method to reconstruct the effective Kohn-Sham (KS) potential. We propose a direct, efficient KS inversion algorithm to solve for the effective KS potential and then employ the resulting algorithm in PFET. We benchmark our KS inversion algorithm against the recently reported modified Ryabinkin-Kohut-Staroverov (mRKS) procedure. Numerical examples show that, with sufficiently large basis sets, our KS inversion algorithm generates almost as accurate results as the mRKS procedure does, except in the vicinity of atomic nuclei, and that it requires less computational time. Three types of chemical interactions then were tested using the new KS inversion algorithm in PFET; the energetics computed from the updated formalism compare well to benchmarks.

Density-functional theory for systems with noncollinear spin: Orbital-dependent exchange-correlation functionals and their application to the Hubbard dimer

A new class of orbital-dependent exchange-correlation (xc) potentials for applications in noncollinear spin-density-functional theory is developed. Starting from the optimized effective potential (OEP) formalism for the exact exchange potential - generalized to the noncollinear case - correlation effects are added via a self-consistent procedure inspired by the Singwi-Tosi-Land-Sjolander (STLS) method. The orbital-dependent xc potentials are applied to the Hubbard dimer in uniform and noncollinear magnetic fields and compared to exact diagonalization and to the Bethe-ansatz local spin-density approximation. The STLS gives the overall best performance for total energies, densities and magnetizations, particularly in the weakly to moderately correlated regime.

Renormalized Møller–Plesset Correlation Energy Functional Used in the Optimized Effective Potential Method

We present the renormalized Moller–Plesset correlation energy functional that is used in the optimized effective potential (OEP) method. The correlation energy functional is developed with the aid ...

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.

Application of the Active Space Self-Interaction-Correction Method to Molecular Systems

Within the context of the active space of the self-interactioncorrection (SIC) optimized effective potential (OEP) method (J. Chem. Phys. 2001, 114, 639-651), the effect of the inclusion of the SIC at the level of only use the HOMO orbital is analyzed for a set of small molecules and for a model of an interstitial region surrounded by positively charged groups in a polypeptide; the model is representative of a class of regions occurring in proteins. It is shown, for the molecular systems treated in this work, that the inclusion of the HOMO orbital, within the SIC-OEP, induces a remarkable change on the eigenvalue spectrum. For the interstitial state model, the improvement is systematic as one increase the active space from one to ten orbitals; also, the influence on the local behavior of the interstitial virtual state closest to the Fermi level is important and enhances its regional character. As this method reduces the computational effort to introduce the SIC it seems promising to deal with self-interaction-corrections in systems with many atoms/electrons such as biomolecules

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.

Self-consistent double-hybrid density-functional theory using the optimized-effective-potential method.

We introduce an orbital-optimized double-hybrid (DH) scheme using the optimized-effective-potential (OEP) method. The orbitals are optimized using a local potential corresponding to the complete exchange-correlation energy expression including the second-order Møller-Plesset correlation contribution. We have implemented a one-parameter version of this OEP-based self-consistent DH scheme using the BLYP density-functional approximation and compared it to the corresponding non-self-consistent DH scheme for calculations on a few closed-shell atoms and molecules. While the OEP-based self-consistency does not provide any improvement for the calculations of ground-state total energies and ionization potentials, it does improve the accuracy of electron affinities and restores the meaning of the LUMO orbital energy as being connected to a neutral excitation energy. Moreover, the OEP-based self-consistent DH scheme provides reasonably accurate exchange-correlation potentials and correlated densities.

Accurate Kohn-Sham ionization potentials from scaled-opposite-spin second-order optimized effective potential methods.

One important property of Kohn-Sham (KS) density functional theory is the exact equality of the energy of the highest occupied KS orbital (HOMO) with the negative ionization potential of the system. This exact feature is out of reach for standard density-dependent semilocal functionals. Conversely, accurate results can be obtained using orbital-dependent functionals in the optimized effective potential (OEP) approach. In this article, we investigate the performance, in this context, of some advanced OEP methods, with special emphasis on the recently proposed scaled-opposite-spin OEP functional. Moreover, we analyze the impact of the so-called HOMO condition on the final quality of the HOMO energy. Results are compared to reference data obtained at the CCSD(T) level of theory. © 2016 Wiley Periodicals, Inc.

More realistic band gaps from meta-generalized gradient approximations: Only in a generalized Kohn-Sham scheme

Unlike the local density approximation (LDA) and the generalized gradient approximation (GGA), calculations with meta-generalized gradient approximations (meta-GGA) are usually done according to the generalized Kohn-Sham (gKS) formalism. The exchange-correlation potential of the gKS equation is non-multiplicative, which prevents systematic comparison of meta-GGA bandstructures to those of the LDA and the GGA. We implement the optimized effective potential (OEP) of the meta-GGA for periodic systems, which allows us to carry out meta-GGA calculations in the same KS manner as for the LDA and the GGA. We apply the OEP to several meta-GGAs, including the new SCAN functional [Phys. Rev. Lett. 115, 036402 (2015)]. We find that the KS gaps and KS band structures of meta-GGAs are close to those of GGAs. They are smaller than the more realistic gKS gaps of meta-GGAs, but probably close to the less-realistic gaps in the band structure of the exact KS potential, as can be seen by comparing with the gaps of the EXX+RPA OEP potential. The well-known grid sensitivity of meta-GGAs is much more severe in OEP calculations.

Spin-dependent optimized effective potential formalism for open and closed systems

Orbital-based exchange $(x)$ correlation $(c)$ energy functionals, leading to the optimized effective potential (OEP) formalism of density-functional theory (DFT), are gaining increasing importance in ground-state DFT, as applied to the calculation of the electronic structure of closed systems with a fixed number of particles, such as atoms and molecules. These types of functionals prove also to be extremely valuable for dealing with solid-state systems with reduced dimensionality, such as is the case of electrons trapped at the interface between two different semiconductors, or narrow metallic slabs. In both cases, electrons build a quasi-two-dimensional electron gas, or Q2DEG. We provide here a general DFT-OEP formal scheme valid both for Q2DEGs either isolated (closed) or in contact with a particle bath (open), and show that both possible representations are equivalent, being the choice of one or the other essentially a question of convenience. Based on this equivalence, a calculation scheme is proposed which avoids the noninvertibility problem of the density response function for closed systems. We also consider the case of spontaneously spin-polarized Q2DEGs, and find that far from the region where the Q2DEG is localized, the exact $x\text{-only}$ exchange potential approaches two different, spin-dependent asymptotic limits. As an example, aside from these formal results, we also provide numerical results for a spin-polarized jellium slab, using the new OEP formalism for closed systems. The accuracy of the Krieger-Li-Iafrate approximation has been also tested for the same system, and found to be as good as it is for atoms and molecules.