A still open issue in many-body theory is the asymptotic behavior of the exchange-correlation energy and potential in the vacuum region of a metal surface. Here we report a numerical study of the position-dependent exchange-correlation energy for jel
lium slabs, as obtained by combining the formally exact adiabatic-connection-fluctuation-dissipation theorem with either time-dependent density-functional theory or an inhomogeneous Singwi-Tosi-Land-Sjolander approach. We find that the inclusion of correlation allows to obtain well-converged semi-infinite-jellium results (independent of the slab thickness) that exhibit an image-like asymptotic behavior close to the classical image potential $V_{im}(z)=-e^2/4z$.
We construct and study several semilocal density functional approximations for the positive Kohn-Sham kinetic energy density. These functionals fit the kinetic energy density of the Airy gas and they can be accurate for integrated kinetic energies of
atoms, molecules, jellium clusters and jellium surfaces. We find that these functionals are the most accurate ones for atomization kinetic energies of molecules and for fragmentation of jellium clusters. We also report that local and semilocal kinetic energy functionals can show binding when the density of a spin unrestricted Kohn-Sham calculation is used.
We present a detailed study of the coupling-constant-averaged exchange-correlation hole density at a jellium surface, which we obtain in the random-phase approximation (RPA) of many-body theory. We report contour plots of the exchange-only and exchan
ge-correlation hole densities, the integration of the exchange-correlation hole density over the surface plane, the on-top correlation hole, and the energy density. We find that the on-top correlation hole is accurately described by local and semilocal density-functional approximations. We also find that for electrons that are localized far outside the surface the main part of the corresponding exchange-correlation hole is localized at the image plane.
We propose a generalized gradient approximation (GGA) for the angle- and system-averaged exchange-correlation hole of a many-electron system. This hole, which satisfies known exact constraints, recovers the PBEsol (Perdew-Burke-Ernzerhof for solids)
exchange-correlation energy functional, a GGA that accurately describes the equilibrium properties of densely packed solids and their surfaces. We find that our PBEsol exchange-correlation hole describes the wavevector analysis of the jellium exchange-correlation surface energy in agreement with a sophisticated time-dependent density-functional calculation (whose three-dimensional wavevector analysis we report here).
Local and semilocal density-functional approximations for the exchange-correlation energy fail badly in the zero-thickness limit of a quasi-two-dimensional electron gas, where the density variation is rapid almost everywhere. Here we show that a full
y nonlocal fifth-rung functional, the inhomogeneous Singwi-Tosi-Land-Sjolander (STLS) approach, which employs both occupied and unoccupied Kohn-Sham orbitals, recovers the true two-dimensional STLS limit and appears to be remarkably accurate for any thickness of the slab (and thus for the dimensional crossover). We also show that this good behavior is only partly due to the use of the full exact exchange energy.
We propose a simple dynamic exchange-correlation kernel of the uniform electron gas. We model the reduction of the electron-electron interaction due to short-range exchange-correlation effects by introducing a frequency-dependent error-function effec
tive interaction. By imposing the fulfillment of the compresibility and the third-frequency-moment sum rules, as well as the correct asymptotic behavior at large wave vectors, we find an accurate and simple dynamic exchange-correlation kernel that accurately reproduces the wave-vector analysis and the imaginary-frequency analysis of the correlation energy of the uniform electron gas.
