Exact-exchange Kohn-Sham potential, surface energy, and work function of jellium slabs

Exact-exchange self-consistent calculations of the Kohn-Sham potential, surface energy, and work function of jellium slabs are reported in the framework of the optimized effective potential (OEP) scheme of density functional theory. In the vacuum side of the jellium surface and at a distance z that is larger than the slab thickness, the exchange-only Kohn-Sham potential is found to be imagelike (∼−e2/z) but with a coefficient that differs from that of the classical image potential Vim(z)=−e2/4z. The three OEP contributions to the surface energy (kinetic, electrostatic, and exchange) are found to oscillate as a function of the slab thickness, as occurs in the case of the corresponding calculations based on the use of single-particle orbitals and energies obtained in the local-density approximation (LDA). The OEP work function presents large quantum size effects that are absent in the LDA and which reflect the intrinsic derivative discontinuity of the exact Kohn-Sham potential.