Comparison of integral equation theories of the liquid state

The Ornstein-Zernike equation is a powerful tool in liquid state theory for predicting structural and thermodynamic properties of fluids. Combined with a suitable closure, it has been shown to reproduce, e.g., the static structure factor, pressure, and compressibility of liquids to a great degree of accuracy. However, out of the multitude of closures that exist for the Ornstein-Zernike equation, it is hard to predict a priori which closure will give the most accurate predictions for the system at hand. To alleviate this problem, we compare the predictive power of many closures on a curated set of representative benchmark systems, including those with hard-sphere, inverse power-law, Gaussian core, and Lennard-Jones particles, in three and two dimensions. For example, we find that the well-known and highly used Percus-Yevick closure gives significantly worse predictions than lesser-known closures of equal complexity in all cases studied. We anticipate that the trends observed in our results will aid in making more informed decisions regarding closure choices. To facilitate the adoption of more modern closure theories, we also have packaged, documented, and distributed the code necessary to numerically solve the equations for a given closure and pair interaction potential.