Accuracy of the mean activity coefficient expression (Hansen-Vieillefosse-Belloni equation), valid within the hypernetted chain (HNC) approximation, was tested in a wide concentration range against new Monte Carlo (MC) data for +1:−1 and +2:−2 primitive model electrolytes. The expression has an advantage that the excess chemical potential can be obtained directly, without invoking the time consuming Gibbs-Duhem calculation. We found the HNC results for the mean activity coefficient to be in good agreement with the machine calculations performed for the same model. In addition, the thermodynamic consistency of the HNC approximation was tested. The mean activity coefficients, calculated via the Gibbs-Duhem equation, seem to follow the MC data slightly better than the Hansen-Vieillefosse-Belloni expression. For completeness of the calculation, the HNC excess internal energies and osmotic coefficients are also presented. These results are compared with the calculations based on other theories commonly used to describe electrolyte solutions, such as the mean spherical approximation, Pitzer's extension of the Debye-Hückel theory, and the Debye-Hückel limiting law.

The Replica Ornstein-Zernike (ROZ) equations were used to study the adsorption of ions from electrolyte mixtures. The adsorbent was represented as a quenched primitive model +1:−1 size symmetric electrolyte, while the mobile particles were ions differing in charge and/or size. The ROZ equations in hypernetted-chain (HNC) approximation were tested against new Monte Carlo results in the grand canonical ensemble; good agreement between the two methods was obtained. The ROZ/HNC theory was then used to study the exclusion coefficients as a function of size and/or charge asymmetry of the annealed ions.

Very recently the effect of equisized charged hard sphere solutes in a mixture with core-softened fluid model on the structural and thermodynamic anomalies of the system has been explored in detail by using Monte Carlo simulations and integral equations theory [J. Chem. Phys., 2012, 137, 244502]. Our objective of the present short work is to complement this study by considering univalent ions of unequal diameters in a mixture with the same soft-core fluid model. Specifically, we are interested in the analysis of changes of the temperature of maximum density (TMD) lines with ion concentration for three model salt solutes, namely sodium chloride, potassium chloride and rubidium chloride models. We resort to Monte Carlo simulations for this purpose. Our discussion also involves the dependences of the pair contribution to excess entropy and of constant volume heat capacity on the temperature of maximum density line. Some examples of the microscopic structure of mixtures in question in terms of pair distribution functions are given in addition.