Phase transitions and tricriticality in the lattice restricted primitive model supplemented by short-range interactions

Ginzburg criterion for ionic fluids: the effect of Coulomb interactions

The effect of the Coulomb interactions on the crossover between mean-field and Ising critical behavior in ionic fluids is studied using the Ginzburg criterion. We consider the charge-asymmetric primitive model supplemented by short-range attractive interactions in the vicinity of the gas-liquid critical point. The model without Coulomb interactions exhibiting typical Ising critical behavior is used to calibrate the Ginzburg temperature of the systems comprising electrostatic interactions. Using the collective variables method, we derive a microscopic-based effective Hamiltonian for the full model. We obtain explicit expressions for all the relevant Hamiltonian coefficients within the framework of the same approximation, i.e., the one-loop approximation. Then we consistently calculate the reduced Ginzburg temperature t(G) for both the purely Coulombic model (a restricted primitive model) and the purely nonionic model (a hard-sphere square-well model) as well as for the model parameters ranging between these two limiting cases. Contrary to the previous theoretical estimates, we obtain the reduced Ginzburg temperature for the purely Coulombic model to be about 20 times smaller than for the nonionic model. For the full model including both short-range and long-range interactions, we show that t(G) approaches the value found for the purely Coulombic model when the strength of the Coulomb interactions becomes sufficiently large. Our results suggest a key role of Coulomb interactions in the crossover behavior observed experimentally in ionic fluids as well as confirm the Ising-like criticality in the Coulomb-dominated ionic systems.

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.

Lattice model of equilibrium polymerization. V. Scattering properties and the width of the critical regime for phase separation

Dynamic clustering associated with self-assembly in many complex fluids can qualitatively alter the shape of phase boundaries and produce large changes in the scale of critical fluctuations that are difficult to comprehend within the existing framework of theories of critical phenomena for nonassociating fluids. In order to elucidate the scattering and critical properties of associating fluids, we consider several models of equilibrium polymerization that describe widely occurring types of associating fluids at equilibrium and that exhibit the well defined cluster geometry of linear polymer chains. Specifically, a Flory-Huggins-type lattice theory is used, in conjunction with the random phase approximation, to compute the correlation length amplitude xi(o) and the Ginzburg number Gi corresponding, respectively, to the scale of composition fluctuations and to a parameter characterizing the temperature range over which Ising critical behavior is exhibited. Our calculations indicate that upon increasing the interparticle association energy, the polymer chains become increasingly long in the vicinity of the critical point, leading naturally to a more asymmetric phase boundary. This increase in the average degree of polymerization implies, in turn, a larger xi(o) and a drastically reduced width of the critical region (as measured by Gi). We thus obtain insight into the common appearance of asymmetric phase boundaries in a wide range of "complex" fluids and into the observation of apparent mean field critical behavior even rather close to the critical point.

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.

Flory-Huggins model of equilibrium polymerization and phase separation in the Stockmayer fluid

The competition between chain formation and phase separation in the Stockmayer fluid (SF) of dipolar particles is analyzed using a renormalized Flory-Huggins model of equilibrium polymerization. Calculated critical properties (T(c), phi(c), Z(c)) for the SF compare favorably with simulations over a wide range of the dimensionless dipolar (or "sticking") energy mu*. We find that the polymerization transition preempts phase separation for a large mu*, i.e., (mu*)2 >22.

Thermodynamics of electrolytes on anisotropic lattices

The phase behavior of ionic fluids on simple cubic and tetragonal (anisotropic) lattices has been studied by grand canonical Monte Carlo simulations. Systems with both the true lattice Coulombic potential and continuous-space 1/r electrostatic interactions have been investigated. At all degrees of anisotropy, only coexistence between a disordered low-density phase and an ordered high-density phase with the structure similar to ionic crystal was found, in contrast to recent theoretical predictions. Tricritical parameters were determined to be increasing functions of anisotropy parameters which is consistent with theoretical calculations based on the Debye-Hückel approach. At large anisotropies a two-dimensional-like behavior is observed, from which we estimated the dimensionless tricritical temperature and density for the two-dimensional square lattice electrolyte to be T*(tri)=0.14 and rho*(tri)=0.70.

Effect of space discretization on phase diagrams in ionic systems: a field-theoretic approach

Using Landau-Ginzburg-Wilson field-theoretic methods we determine for what kind of space discretization the transition between charge-disordered and charge-ordered phases in the restricted primitive model is fluctuation-induced first order. We identify this transition with ionic-crystal formation. We predict four types of generic phase diagrams in ionic systems for various kinds of space discretization for low and intermediate densities of ions. Our results also shed light on the simulation results obtained for an off-lattice ionic system over a wide range of densities, including the fcc crystal.