Computations of Radial Distribution Functions for a Classical Electron Gas

The radial distribution functions g for a classical electron gas computed using the Percus-Yevick (PY) equation, convolution hypernetted chain (CHNC) equation, and the Broyles-Sahlin (BS) method, have been compared with the Debye-H\"uckel (DH) theory. The quantities E≡ŪNkt and P≡pn̄kT have been computed from these g's. Computations have been made for values of θ of 20, 10, 5, 3, and 1; θ=kTaq², where a is the ion sphere radius. The PY and BS results show the best agreement, particularly at θ<3. The BS method has been of particular value in this study of a long-range potential. In the range of θ studied, g never exceeds one, that is, there is no oscillatory behavior of g.