CONSTRUCTION AND APPLICATION OF AN ACCURATE LOCAL SPIN-POLARIZED KOHN-SHAM POTENTIAL WITH INTEGER DISCONTINUITY - EXCHANGE-ONLY THEORY

An accurate spin-polarized exchange-only Kohn-Sham (KS) potential is constructed from a consideration of the optimized-effective-potential (OEP) method. A detailed analysis of the OEP integral equation for the exchange-only case results in a set of conditions which are manifestly satisfied by the exact OEP; these conditions are employed to construct an approximate OEP, V(x-sigma) and therefore an approximate KS exchange-only potential as a functional of KS orbitals. Further, it is shown that this V(x-sigma) can be derived analytically based on a simple approximation of the Green's functions in the OEP integral equation. The constructed potential, although approximate, contains many of the key analytic features of the exact KS potential: it reduces to the exact KS result in the homogeneous-electron-gas limit, approaches -1/r as r --> infinity, yields highest occupied-orbital energy eigenvalues epsilon(m-sigma) that satisfy Koopmans's theorem, and exhibits an integer discontinuity when considered as a function of fractional occupancy of the highest-energy occupied single-particle state of a given spin projection-sigma. In addition epsilon(m-sigma) nearly exactly satisfies Janak's theorem. The approximate OEP is a simple but remarkably accurate representation of the exact, numerically derived exchange-only OEP. Detailed numerical results obtained by employing V(x-sigma) as the exchange-only potential for ten atoms with closed subshells yield total energies, Hartree potentials, single-particle expectation values, and epsilon(m) which are in excellent agreement with both exact OEP and Hartree-Fock (HF) results and represent a significant improvement over the results obtained by employing other exchange-only potentials. Similarly, the properties of alkali-metal atoms are calculated including the separate spin-up and spin-down densities to obtain results in excellent agreement with those of spin-unrestricted OEP and HF methods. Finally, we demonstrate the accuracy of V(x-sigma) by calculating the total energy, epsilon(m up), and epsilon(m down) as a function of fractional filling f, of the highest occupied single-particle orbital for the magnesium atom (Z = 12) from N = 9-12 electrons and find excellent agreement with both spin-unrestricted OEP and HF results even when epsilon(m-sigma) is strongly dependent on f. In addition we display the integer discontinuity in V(x-sigma) when the highest-energy spin subshell begins to be filled.