Solutions to the Percus-Yevick Equation

The radial distribution function for a classical fluid of particles interacting with the Lennard-Jones potential has been computed by solving the Percus-Yevick integral equation numerically. The solutions and the quantities, p/nkT and E/NkT, are compared with those obtained by Wood and Parker using Monte Carlo techniques. The radial distribution functions agree to better than 15% beyond the first points where they are unity while the thermodynamic quantities differ by, at most, 3% from the range of Monte Carlo values for cases where the system is believed to be in a liquid state. The quantity K=−(1/V)(∂V/∂p)N,T,has also_been computed.