CORRELATION-ENERGY FUNCTIONAL AND ITS HIGH-DENSITY LIMIT OBTAINED FROM A COUPLING-CONSTANT PERTURBATION EXPANSION

A perturbation theory is developed for the correlation energy E(c)[n], of a finite-density system, with respect to the coupling constant a w hich multiplies the electron-electron repulsion operator in H(alpha)=T +alphaV(ee)+SIGMA(i)v(alpha)(r(i)). The external potential v(alpha) is constrained to keep the ground-state density n fixed for all a greate r-than-or-equal-to 0. v(alpha) is given completely in terms of functio nal derivatives at full charge (alpha = 1), from which e(c,2)[n]+lambd a-1e(c,3)[n]+lambda-2e(c,4)[n]+..., where each e(c,j)[n] is expressed in terms of integrals involving Kohn-Sham determinants. Here, n(lambda )(x,y,x)=lambda3n(lambdax,lambday,lambdaz) and lambda=alpha-1. The ide ntification of lim(lambda-->infinity)E(c)[n(lambda)], which is a high- density limit, as the second-order energy e(c,2)[n] allows one to comp ute bounds upon lim(lambda-->infinity)E(c)[n(lambda)]; the bounds impl y that lim(lambda-->infinity)E(c)[n(lambda)]congruent-to E(c)[n] for a large class of small atoms and molecules, and suggest that lim(lambda -->infinity)E(c)[n(lambda)] should be of the same order of magnitude a s E(c)[n] infinite insulators and semiconductors. Approximations to E( c)[n] should reflect all this. In contrast, perhaps the well-known ove rbinding of the local-density approximation (LDA) in molecules and sol ids is due, in part, to the fact that the LDA correlation energy is to o sensitive to a coordinate scaling of n. Indeed, the LDA for E(c)[n(l ambda)] diverges when lambda-->infinity because of the presence of the -1n(lambda) term in the Gell-Mann and Brueckner high-density expressi on for the correlation energy, per particle, of a homogeneous density, which is infinite. In a sense, the derived perturbation expansion tra nsforms the Gell-Mann and Brueckner expression into one that applies s pecifically to an inhomogeneous density which integrates to a finite n umber of electrons.