Generation of functional derivatives in Kohn-Sham density-functional theory Original Research Article

A general procedure is presented for obtaining the form of exchange-correlation potentials arising from semilocal electron density functionals. The procedure is implemented as a Mathematica program, which can be used to obtain functional derivatives both in terms of the density and its gradients, and also as an explicit function of position for the case that the density is specified.