Choice of reference, influence of non-additivity, and present challenges in thermodynamic perturbation theory for mixtures

This work revisits the fundamentals of thermodynamic perturbation theory for fluid mixtures. The choice of reference and governing assumptions can profoundly influence the accuracy of the perturbation theory. The statistical associating fluid theory for variable range interactions of the generic Mie form equation of state is used as a basis to evaluate three choices of hard-sphere reference fluids: single component, additive mixture, and non-additive mixture. Binary mixtures of Lennard-Jones fluids are investigated, where the ratios of σ (the distance where the potential is zero) and the ratios of ϵ (the well depth) are varied. By comparing with Monte Carlo simulations and results from the literature, we gauge the accuracy of different theories. A perturbation theory with a single-component reference gives inaccurate predictions when the σ-ratio differs significantly from unity but is otherwise applicable. Non-additivity becomes relevant in phase-equilibrium calculations for fluids with high ϵ-ratios or when the mixing rule of σ incorporates non-additivity through an adjustable parameter. This can be handled in three ways: by using a non-additive hard-sphere reference, by incorporating an extra term in the additive hard-sphere reference, or with a single-component reference when the σ-ratio is close to unity. For σ- and ϵ-ratios that differ significantly from unity, the perturbation theories overpredict the phase-equilibrium pressures regardless of reference. This is particularly pronounced in the vicinity of the critical region for mixtures with high ϵ-ratios. By comparing with Monte Carlo simulations where we compute the terms in the perturbation theory directly, we find that the shortcomings of the perturbation theory stem from an inaccurate representation of the second- and third-order perturbation terms, a₂ and a₃. As mixtures with molecules that differ significantly in size and depths of their interaction potentials are often encountered in industrial and natural applications, further development of the perturbation theory based on these results is an important future work.