Next-order correction to the Dirac exchange energy of the free electron gas in the thermodynamic limit and generalized gradient approximations

Abstract

We derive the next order correction to the Dirac exchange energy for the free electron gas in a box with zero boundary conditions in the thermodynamic limit. The correction is of the order of the surface area of the box, and comes from three different contributions: (i) a real-space boundary layer, (ii) a boundary-condition-induced small shift of Fermi momentum and bulk density, and (iii) a long-range electrostatic finite-size correction. Moreover we show that the local density approximation, in addition to capturing the bulk term exactly, also produces a correction of the correct order but not the correct size. Generalized gradient approximation (GGA) corrections are found to be capable of capturing the surface term exactly, provided the gradient enhancement factor satisfies a simple explicit integral constraint. For current GGAs such as B88 and Perdew, Burke and Ernzerhof we find that the new constraint is not satisfied and the size of the surface correction is overestimated by about ten percent. The new constraint might thus be of interest for the design of future exchange functionals.