Numerical Solutions of the Percus–Yevick Equation for the Lennard-Jones (6–12) and Hard-Sphere Potentials

Numerical solutions of the Percus–Yevick equation for the Lennard-Jones (6–12) potential are presented largely for high densities and the liquid region in a form particularly useful to workers who have access only to small computers which cannot be used for solving the PY equation. The solutions, given in a form suitable for calculating either the radial distribution functions or direct correlation functions, are deposited with the National Auxiliary Publications Service (Document No. NAPS-00720). A table of densities and temperatures at which the PY solutions are available is given. Energies and pressures, including the compressibility pressure for temperatures below the critical, are also tabulated and comparisons are made with molecular dynamics results. Solutions for the hard-sphere potential are also presented completing the tables of Throop and Bearman J. Chem. Phys. 43, 2408 (1965) for large intermolecular separations.