Criticality, tricriticality, and crystallization in discretized models of electrolytes

Distribution of ions near a charged selective surface in critical binary solvents

Near-critical binary mixtures containing ionic solutes near a charged wall preferentially adsorbing one component of the solvent are studied. Within the Landau-Ginzburg approach extended to include electrostatic interactions and the chemical preference of ions for one component of the solvent, we obtain a simple form for the leading-order correction to the Debye-Hückel theory result for the charge density profile. Our result shows that critical adsorption influences significantly distribution of ions near the wall. This effect may have important implications for the screening of electrostatic interactions between charged surfaces immersed in binary near-critical solvents.

Field-theoretic description of ionic crystallization in the restricted primitive model

Effects of charge-density fluctuations on a phase behavior of the restricted primitive model are studied within a field-theoretic formalism. We focus on a lambda line of continuous transitions between charge-ordered and charge-disordered phases that is observed in several mean-field theories, but is absent in simulation results. In our study the RPM is reduced to a phi(6) theory, and a fluctuation contribution to a grand thermodynamic potential is obtained by generalizing the Brazovskii approach. We find that in a presence of fluctuations the lambda -line disappears. Instead, a fluctuation-induced first-order transition to a charge-ordered phase appears in the same region of a phase diagram, where the liquid-ionic-crystal transition is obtained in simulations. Our results indicate that the charge-ordered phase should be identified with an ionic crystal.

Universality class of the critical point in the restricted primitive model of ionic systems

A coarse-grained description of the restricted primitive model is considered in terms of the local charge- and number-density fields. Exact reduction to a one-field theory is derived, and exact expressions for the number-density correlation functions in terms of higher-order correlation functions for the charge-density are given. It is shown that in continuum space the singularity of the charge-density correlation function associated with short-wavelength charge-ordering disappears when charge-density fluctuations are included by following the Brazovskii approach. The related singularity of the individual Feynman diagrams contributing to the number-density correlation functions is cured when all the diagrams are segregated into disjoint sets according to their topological structure. By performing a resummation of all diagrams belonging to each set a regular expression represented by a secondary diagram is obtained. The secondary diagrams are again segregated into disjoint sets, and the series of all the secondary diagrams belonging to a given set is represented by a hyperdiagram. A one-to-one correspondence between the hyperdiagrams contributing to the number-density vertex functions, and diagrams contributing to the order-parameter vertex functions in a certain model system belonging to the Ising universality class is demonstrated. Corrections to scaling associated with irrelevant operators that are present in the model-system Hamiltonian, and other corrections specific to the RPM are also discussed.

Phase behavior of the lattice restricted primitive model with nearest neighbor exclusion

The global phase behavior of the lattice restricted primitive model with nearest neighbor exclusion has been studied by grand canonical Monte Carlo simulations. The phase diagram is dominated by a fluid (or charge-disordered solid) to charge-ordered solid transition that terminates at the maximum density rho*(max)= sqrt 2 and reduced temperature T* approximately equal to 0.29. At that point, there is a first-order phase transition between two phases of the same density, one charge-ordered, and the other charge-disordered. The liquid-vapor transition for the model is metastable, lying entirely within the fluid-solid phase envelope.

Thermodynamic properties of lattice hard-sphere models

Thermodynamic properties of several lattice hard-sphere models were obtained from grand canonical histogram- reweighting Monte Carlo simulations. Sphere centers occupy positions on a simple cubic lattice of unit spacing and exclude neighboring sites up to a distance sigma. The nearestneighbor exclusion model, sigma = radical2, was previously found to have a second-order transition. Models with integer values of sigma = 1 or 2 do not have any transitions. Models with sigma = radical3 and sigma = 3 have weak first-order fluid-solid transitions while those with sigma = 2 radical2, 2 radical3, and 3 radical2 have strong fluid-solid transitions. Pressure, chemical potential, and density are reported for all models and compared to the results for the continuum, theoretical predictions, and prior simulations when available.

Phase diagrams in the lattice restricted primitive model: from order-disorder to gas-liquid phase transition

The phase behavior of the lattice restricted primitive model (RPM) for ionic systems with additional short-range nearest neighbor (NN) repulsive interactions has been studied by grand canonical Monte Carlo simulations. We obtain a rich phase behavior as the NN strength is varied. In particular, the phase diagram is very similar to that of the continuum RPM for high NN strength. Specifically, we have found both gas-liquid phase separation, with associated Ising critical point, and a first-order liquid-solid transition. We discuss how the line of continuous order-disorder transitions present for the low NN strength changes into continuum-space behavior as one increases the NN strength and compare our findings with recent theoretical results by Ciach and Stell [Phys. Rev. Lett. 91, 060601 (2003)].

Effect of charge-density fluctuations on the order-disorder transition in the lattice restricted primitive models of electrolytes

Lattice restricted primitive models, where equally charged ions of a diameter sigma are located at sites of a simple cubic lattice with a lattice constant a , are studied within the field-theoretic approach. We focus on the transition between charge-disordered and charge-ordered phases for sigma/a=1 and sigma/a=2 . By using renormalization-group methods we show that at high concentrations of ions the transition is continuous for sigma/a=1 , while for sigma/a=2 it is only first order, as found previously in simulations. For sigma/a=1 the effect of charge-density fluctuations on the positions of the continuous transition and the tricritical point (TCP) is determined within a formalism developed in this work. The temperature and the concentration of ions at the TCP agree very well with simulation results.