Free energy and phase equilibria for the restricted primitive model of ionic fluids from Monte Carlo simulations

Coexistence and criticality in size-asymmetric hard-core electrolytes

Liquid-vapor coexistence curves and critical parameters for hard-core 1:1 electrolyte models with diameter ratios lambda = sigma(-)/sigma(+) = 1 to 5.7 have been studied by fine-discretization Monte Carlo methods. Normalizing via the length scale sigma(+/-) = 1 / 2(sigma(+)+sigma(-)), relevant for the low densities in question, both T(*)(c) ( = k(B)T(c)sigma(+/-)/q(2)) and rho(*)(c) ( = rho(c)sigma(3)(+/-)) decrease rapidly (from approximately 0.05 to 0.03 and 0.08 to 0.04, respectively) as lambda increases. These trends, which unequivocally contradict current theories, are closely mirrored by results for tightly tethered dipolar dimers (with T(*)(c) lower by approximately 0%-11% and rho(*)(c) greater by 37%-12%).

Order-disorder transition in the solid phase of a charged hard sphere model

We investigate the solid phases of the restricted primitive model (RPM). Monte Carlo simulations show the existence of an order-disorder transition from a substitutionally disordered face centered cubic lattice (fcc) to a new ordered fcc structure which is proposed as the ground state of the RPM at the close packing density. Our results suggest that the new phase might turn out in a new triple point in the RPM phase diagram involving three solid phases: CsCl, fcc ordered and fcc disordered structures. The order-disorder transition is also studied using the cell theory. The theory shows good agreement with the simulation results and suggests that the transition is weakly first order.

Ion association in model ionic fluids

Monte Carlo simulations are used to investigate the temperature and pressure (density) dependence of ion association in the restricted primitive model. It is shown that at temperatures below the critical temperature T(c) the vapor consists almost exclusively of strongly bound ion pairs at or near contact. Significant ion-pair dissociation begins at temperatures very near T(c). This raises the possibility that compositional fluctuations between strongly bound and free ions influence the critical behavior. We note qualitative similarities between the present results and the Kosterlitz-Thouless transition in the two-dimensional Coulomb gas.