Why the Ising and continuous-space models of ionic systems exhibit essentially different critical behavior

The line tension of two-dimensional ionic fluids

Pressure tensor components are very useful in the calculation of the tension associated with a liquid-vapor interface. In this work, we present expressions for the pressure tensor components of two-dimensional ionic fluids, modeled at the level of the primitive model. As an application, we carried out molecular dynamics simulations of liquid-vapor interfaces to calculate the line tension of the 1:1 two-dimensional ionic fluid, whose liquid-vapor coexistence curve had already been obtained in a previous work. The pressure tensor components were validated by simulating states of one phase and reproducing the scalar pressure, previously obtained from bulk simulations and reported in the literature. The effects on the line tension and the coexisting densities, originated by the choice of the Ewald parameters, the cutoff radius, and the interfacial length were also evaluated.

Universal sequence of ordered structures obtained from mesoscopic description of self-assembly

Mesoscopic theory for soft-matter systems that combines density functional and statistical field theory is derived by a systematic coarse-graining procedure. For particles interacting with spherically symmetric potentials of arbitrary form, the grand-thermodynamic potential consists of two terms. The first term is associated with microscopic length-scale fluctuations, and has the form of the standard density functional. The second term is associated with mesoscopic length-scale fluctuations, and has the form known from the statistical field theory. For the correlation function between density fluctuations in mesoscopic regions, a pair of equations similar to the Ornstein-Zernicke equation with a new closure is obtained. In the special case of weak ordering on the mesoscopic length scale, the theory takes a form similar to either the Landau-Ginzburg-Wilson (LGW) or the Landau-Brazovskii (LB) field theory, depending on the form of the interaction potential. Microscopic expressions for the phenomenological parameters that appear in the Landau-type theories are obtained. Within the framework of this theory, we obtain a universal sequence of phases: disordered, bcc, hexagonal, lamellar, inverted hexagonal, inverted bcc, disordered, for increasing density well below the close-packing density. The sequence of phases agrees with experimental observations and with simulations of many self-assembling systems. In addition to the above phases, more complex phases may appear depending on the interaction potentials. For a particular form of the short-range attraction long-range repulsion potential, we find the bicontinuous gyroid phase ( Ia3d symmetry) that may be related to a network-forming cluster of colloids in a mixture of colloids and nonadsorbing polymers.

Field theory for size- and charge-asymmetric primitive model of ionic systems: mean-field stability analysis and pretransitional effects

The primitive model of ionic systems is investigated within a field-theoretic description for the whole range of diameter-, lambda , and charge, Z ratios of the two ionic species. Two order parameters (OP) are identified. The relation of the OP's to physically relevant quantities is nontrivial. Each OP is a linear combination of the charge density and the number-density waves. Instabilities of the disordered phase associated with the two OP's are determined in the mean-field approximation (MF). In MF a gas-liquid separation occurs for any Z and lambda is not equal to 1 . In addition, an instability with respect to various types of periodic ordering of the two kinds of ions is found. Depending on lambda and Z , one or the other transition is metastable in different thermodynamic states. For sufficiently large size disparity we find a sequence of fluid-crystal-fluid transitions for the increasing volume fraction of ions, in agreement with experimental observations. The instabilities found in MF represent weak ordering of the most probable instantaneous states, and are identified with structural loci associated with pretransitional effects.

Nonasymptotic Critical Behavior of a Ternary Ionic System

Refractive indices n and salt concentrations ms of coexisting phases of the ternary system 1,4-dioxane + water + potassium chloride were measured along the liquid-liquid-solid coexistence curve near the liquid-liquid critical end point. Refractive index measurements were carried out in the range 0.689 x 10-3 < t = (T - Tc)/Tc < 0.118 while salt concentrations were determined for the temperature range 1.84 x 10-3 < t < 8.07 x 10-2. From these experimental results, compositions fD (mass fraction of dioxane on a salt-free basis) and densities rho of coexisting phases were obtained. The shape of the coexistence curve was analyzed using alternatively n, ms, fD, and rho as order parameters. In all cases, the obtained coexistence curve displays, asymptotically, Ising behavior. Outside the asymptotic critical domain, n, ms, and rho show significant deviations of the effective critical exponent from its Ising value, while the concentration variable fD requires no corrections to simple scaling. On the basis of the present results, we conclude that this system shows no indication of multicritical behavior.

Simulation of symmetric tricritical behavior in electrolytes

Despite extensive experimental, theoretical, and simulation efforts, a unified description of ionic phase transitions and criticality has not yet emerged. In this work, we investigate the phase behavior of the restricted primitive model of electrolyte solutions on the simple cubic lattice using grand canonical Monte Carlo simulations and finite-size scaling techniques. The phase diagram of the system is distinctly different from its continuum-space analog. We find order-disorder transitions for reduced temperatures T* < or approximately 0.51, where the ordered structures resemble those of the NaCl crystal. The order-disorder transition is continuous for 0.15 < or approximately T* < or approximately 0.51 and becomes first order at lower temperatures. The line of first-order transitions is a line of three-phase coexistence between a disordered and two ordered phases. The line of continuous, second-order transitions meets this line of triple points at a tricritical point at T* approximately equal to 0.1475. We locate the line of continuous transitions, and the line of triple points using finite-size scaling techniques. The tricritical temperature is estimated by extrapolation of the size-dependent tricritical temperatures obtained from a sixth-order Landau expansion of the free energy. Our calculated phase diagram is in qualitative agreement with mean-field theories.

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.

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)].

Criticality, tricriticality, and crystallization in discretized models of electrolytes

The Restricted Primitive Model of ionic systems is studied within a field-theoretic approach in order to provide a theoretic basis for the qualitative difference in the phase diagrams obtained in simulations for sigma/a=1,sigma/a=2, and sigma/a>/=3 ( a is the lattice constant and sigma is the ion diameter). The evolution of the phase diagrams from the case sigma/a=1 to the case sigma/a=square root[2] [nearest-neighbor ( NN ) occupancy excluded] is studied in the model with NN repulsion, 0</=J</= infinity, supplementing Coulomb forces. The boundary of stability of the charge-disordered phase with respect to short-wavelength charge fluctuations and the tricritical point are found in a mean-field ( MF ) approximation. Next, the effect of fluctuations is studied and we find that for J exceeding a particular value J0 a fluctuation-induced first-order phase transition should be expected instead of the continuous transition found in MF. At J= J0 the line of continuous transitions splits into two lines enclosing the two-phase region, whose thickness increases from zero when J>/= J0 increases. We argue that this transition corresponds to formation of a bcc ionic crystal. For high densities the ions form an fcc crystal, for which we find a fluctuation-induced first-order charge-ordered-charge-disordered transition, in agreement with recent simulation studies. Our results also shed light on the simulation results obtained for an off-lattice ionic system, for which a schematic phase diagram is constructed.

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.

Hierarchical reference theory study of the lattice restricted primitive model

A three dimensional model of point charges, named lattice restricted primitive model (LRPM), is investigated by using the hierachical reference theory of fluids. This approach, which generalizes the momentum renormalization group technique, is shown to capture the physics of the model and provides a quantitative description of the phase diagram. The comparison with recent numerical simulations and with other theoretical approaches is discussed both for the LRPM and for the Blume-Capel model, which can be seen as a screened version of LRPM. The nonuniversal crossover region close to the tricritical point is also discussed.

Critical behavior of ionic solids

Phase transitions of lattice models of ionic crystals are studied by computer simulation. The nature of order-disorder transitions on different crystal structures is established and compared with the behavior of related Ising models. It is found that for both, continuous and first order transitions the basic features seem to be similar to those of Ising systems.