Construction of integrable model Kohn-Sham potentials by analysis of the structure of functional derivatives

A directly approximated exchange-correlation potential should, by construction, be a functional derivative of some density functional in order to avoid unphysical results. Using generalized gradient approximations (GGAs) as an example, we show that functional derivatives of explicit density functionals have a very rigid inner structure, the knowledge of which allows one to build the entire functional derivative from a small part. Based on this analysis, we develop a method for direct construction of integrable Kohn-Sham potentials. As an illustration, we transform the model potential of van Leeuwen and Baerends (which is not a functional derivative) into a semilocal exchange potential that has a parent GGA, yields accurate energies, and is free from the artifacts inherent in existing semilocal potential approximations.