Gibbs ensemble Monte Carlo of nonadditive hard-sphere mixtures.

In this article, we perform Gibbs ensemble Monte Carlo (GEMC) simulations of liquid-liquid phase coexistence in nonadditive hard-sphere mixtures (NAHSMs) for different size ratios and non-additivity parameters. The simulation data are used to provide a benchmark to a number of theoretical and mixed theoretical/computer simulation approaches which have been adopted in the past to study phase equilibria in NAHSMs, including the method of the zero of the Residual Multi-Particle Entropy, Integral Equation Theories (IETs), and classical Density Functional Theory (DFT). We show that while the entropic criterium is quite accurate in predicting the location of phase equilibrium curves, IETs and DFT provide at best a semi-quantitative reproduction of GEMC demixing curves.