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 alsobeen computed.