Reconstruction of Density Functionals from Kohn-Sham Potentials by Integration along Density Scaling Paths

J Chem Theory Comput. 2009 Apr 14;5(4):699-707. doi: 10.1021/ct800514z. Epub 2009 Mar 10.

Abstract

We demonstrate by specific examples that if a Kohn-Sham exchange-correlation potential is given explicitly in terms of the electron density and its derivatives, then one can easily reconstruct the parent density functional by evaluating analytically (or numerically with one-dimensional quadratures) the van Leeuwen-Baerends line integral (Phys. Rev. A 1995, 51, 170-178) along a path of (coordinate)-scaled densities. The choice of a density scaling path amounts to defining the gauge of the resultant exchange-correlation energy density. The well-known Levy-Perdew virial relation for exchange potentials can be viewed as an analytical line integral along the electron-number-conserving uniform density scaling path. Energies obtained from model exchange-correlation potentials should be interpreted with caution because the reconstructed density functional is unique (up to a gauge transformation) only if the model Kohn-Sham potential is a functional derivative.