Hierarchy of model Kohn-Sham potentials for orbital-dependent functionals: a practical alternative to the optimized effective potential method.

We describe a method for constructing a hierarchy of model potentials approximating the functional derivative of a given orbital-dependent exchange-correlation functional with respect to electron density. Each model is derived by assuming a particular relationship between the self-consistent solutions of Kohn-Sham (KS) and generalized Kohn-Sham (GKS) equations for the same functional. In the KS scheme, the functional is differentiated with respect to density, in the GKS scheme--with respect to orbitals. The lowest-level approximation is the orbital-averaged effective potential (OAEP) built with the GKS orbitals. The second-level approximation, termed the orbital-consistent effective potential (OCEP), is based on the assumption that the KS and GKS orbitals are the same. It has the form of the OAEP plus a correction term. The highest-level approximation is the density-consistent effective potential (DCEP), derived under the assumption that the KS and GKS electron densities are equal. The analytic expression for a DCEP is the OCEP formula augmented with kinetic-energy-density-dependent terms. In the case of exact-exchange functional, the OAEP is the Slater potential, the OCEP is roughly equivalent to the localized Hartree-Fock approximation and related models, and the DCEP is practically indistinguishable from the true optimized effective potential for exact exchange. All three levels of the proposed hierarchy require solutions of the GKS equations as input and have the same affordable computational cost.