We construct analytic formulas that represent the coupling-constant-averaged pair distribution function $g_{\mathrm{xc}}{(r}_s,\zeta {,k}_{F}u)$ of a three-dimensional nonrelativistic ground-state electron gas constrained to a uniform density with density parameter $r_s=(9\mathrm{}\pi {/4)}^{1/3}{/k}_F$ and relative spin polarization $\zeta$ over the whole range $0<\mathrm{}{r}_s<\infty$ and $-1<\mathrm{}\zeta <1,$ with energetically unimportant long range $\left(\overrightarrow{u}\infty \right)$ oscillations averaged out. The pair distribution function $g_{\mathrm{xc}}$ at the physical coupling constant is then given by differentiation with respect to $r_s.$ Our formulas are constructed using only known theoretical constraints plus the correlation energy ${\epsilon }_c{(r}_s,\zeta ),$ and accurately reproduce the $g_{\mathrm{xc}}$ of the quantum Monte Carlo method and of the fluctuation-dissipation theorem with the Richardson-Ashcroft dynamical local-field factor. Our $g_{\mathrm{xc}}$ is correct even in the high-density ${(r}_s\to 0)$ and low-density ${(r}_s\to \infty )$ limits. When the spin resolution of ${\epsilon }_c$ into $\uparrow \uparrow ,$ $\downarrow \downarrow ,$ and $\uparrow \downarrow$ contributions is known, as it is in the high- and low-density limits, our formulas also yield the spin resolution of $g_{\mathrm{xc}}.$ Because of these features, our formulas may be useful for the construction of density functionals for nonuniform systems. We also analyze the kinetic energy of correlation into contributions from density fluctuations of various wave vectors. The exchange and long-range correlation parts of our $g_{\mathrm{xc}}{(r}_s,\zeta {,k}_{F}u)-1$ are analytically Fourier transformable, so that the static structure factor $S_{\mathrm{xc}}{(r}_s,\zeta {,k/k}_F)$ is easily evaluated.

• Received 10 June 2002

DOI:https://doi.org/10.1103/PhysRevB.66.165118