Test of athermal terms of activity coefficient models by Monte Carlo simulation with hard-core models

The infinite dilution activity coefficients of exactly athermal fluids were calculated by Monte Carlo simulation with hard-core models. The hard-core models used in this work were hard-sphere and hard-spherocylinder models. The Widom test particle method was adopted to calculate the residual chemical potentials of solutes in pure solvent and in pure solute solutions. The infinite dilution activity coefficients of solutes were obtained from the residual chemical potentials of solutes. The infinite dilution activity coefficients calculated by Monte Carlo simulation were compared with those of athermal terms in activity coefficient equations. Staverman–Guggenheim equation overestimates the activity coefficients. The deviations of activity coefficients increase with increasing the hard-core volume of solute. Flory–Huggins equation based on molar volume gives good results for the hard-spherocylinder systems. Elbro-FV equation gives good results for both the hard-sphere and hard-spherocylinder systems.