Optimized and transferable nonlocal separable ab initio pseudopotentials

The purpose of the present work is to obtain a better atomic pseudopotential with respect to convergence and computational efficiency while retaining reasonable transferability in the context of electronic-structure calculations for solids using a plane-wave basis set. We introduce a systematic procedure for generating optimized nonlocal pseudopotentials by minimizing the high Fourier components of the pseudo-wave-functions with the constraints of normalization and continuity of first and second derivatives of the wave function at the core radius. This is based on the recent ideas of Rappe et al. (RRKJ) [Phys. Rev. B 41, 1227 (1990)], but overcomes certain difficulties which we have found with the RRKJ scheme. For computational efficiency this optimized nonlocal pseudopotential is transformed into a Kleinman-Bylander (KB) form. To ensure the transferability we first compare the logarithmic derivative of the all-electron wave function with that of the final KB form of the optimized nonlocal pseudopotential over a wide range of energies. We then test the KB form of the potential in a number of atomic environments. The structural properties of ZnS are calculated to demonstrate the reliability of our optimized nonlocal separable ab initio pseudopotential and its total-energy convergence.