A Monte Carlo finite size scaling study of charged hard-sphere criticality

Self-Consistent Description of Vapor-Liquid Interface in Ionic Fluids

Inhomogeneity of ion correlation widely exists in many physicochemical, soft matter, and biological systems. Here, we apply the modified Gaussian renormalized fluctuation theory to study the classic example of the vapor-liquid interface of ionic fluids. The ion correlation is decomposed into a short-range contribution associated with the local electrostatic environment and a long-range contribution accounting for the spatially varying ionic strength and dielectric permittivity. For symmetric salt, both the coexistence curve and the interfacial tension predicted by our theory are in quantitative agreement with simulation data reported in the literature. Furthermore, we provide the first theoretical prediction of interfacial structure for asymmetric salt, highlighting the importance of capturing local charge separation.

Improving the efficiency of Monte Carlo simulations of ions using expanded grand canonical ensembles

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions. We introduce an expanded ensemble approach which may be used in conjunction with these existing methods to further improve efficiency. Individual ions are inserted or deleted in one Monte Carlo trial rather than simultaneous insertion/deletion of neutral sets. We show how charge neutrality is maintained and show rigorous quantitative agreement between the conventional and the proposed expanded ensemble approaches, but with up to an order of magnitude increase in efficiency at high densities. The expanded ensemble approach is also more straightforward to implement than simultaneous insertion/deletion of neutral sets, and its implementation is demonstrated within open source software.

Smectic phases in ionic liquid crystals

Ionic liquid crystals (ILCs) are anisotropic mesogenic molecules which carry charges and therefore combine properties of liquid crystals, e.g. the formation of mesophases, and of ionic liquids, such as low melting temperatures and tiny triple-point pressures. Previous density functional calculations have revealed that the phase behavior of ILCs is strongly affected by their molecular properties, i.e. their aspect ratio, the loci of the charges, and their interaction strengths. Here, we report new findings concerning the phase behavior of ILCs as obtained by density functional theory and Monte Carlo simulations. The most important result is the occurrence of a novel, wide smectic-A phase [Formula: see text], at low temperature, the layer spacing of which is larger than that of the ordinary high-temperature smectic-A phase [Formula: see text]. Unlike the ordinary smectic S _A phase, the structure of the [Formula: see text] phase consists of alternating layers of particles oriented parallel to the layer normal and oriented perpendicular to it.

Critical point of gas-liquid type phase transition and phase equilibrium functions in developed two-component plasma model

A two-component plasma model, which we called a "shelf Coulomb" model has been developed in this work. A Monte Carlo study has been undertaken to calculate equations of state, pair distribution functions, internal energies, and other thermodynamics properties. A canonical NVT ensemble with periodic boundary conditions was used. The motivation behind the model is also discussed in this work. The "shelf Coulomb" model can be compared to classical two-component (electron-proton) model where charges with zero size interact via a classical Coulomb law. With important difference for interaction of opposite charges: electrons and protons interact via the Coulomb law for large distances between particles, while interaction potential is cut off on small distances. The cut off distance is defined by an arbitrary ɛ parameter, which depends on system temperature. All the thermodynamics properties of the model depend on dimensionless parameters ɛ and γ = βe(2)n(1/3) (where β = 1/kBT, n is the particle's density, kB is the Boltzmann constant, and T is the temperature) only. In addition, it has been shown that the virial theorem works in this model. All the calculations were carried over a wide range of dimensionless ɛ and γ parameters in order to find the phase transition region, critical point, spinodal, and binodal lines of a model system. The system is observed to undergo a first order gas-liquid type phase transition with the critical point being in the vicinity of ɛ(crit) ≈ 13(T(*)(crit) ≈ 0.076), γ(crit) ≈ 1.8(v(*)(crit) ≈ 0.17), P(*)(crit) ≈ 0.39, where specific volume v* = 1/γ(3) and reduced temperature T(*) = ɛ(-1).

Liquid-vapor transition of systems with mean field universality class

We have considered a system where the interaction, v(r)=v(IS)(r)+xi(2)v(MF)(r), is given as a linear combination of two potentials, each of which being characterized with a well-defined critical behavior: for v(IS)(r) we have chosen the potential of the restricted primitive model which is known to belong to the three-dimensional Ising universality class, while for v(MF)(r) we have considered a long-range interaction in the Kac [J. Math. Phys. 4, 216 (1963)] limit, displaying mean field (MF) behavior. We study the performance of two theoretical approaches and of computer simulations in the critical region for this particular system and give a detailed comparison between theories and simulation of the critical region and the location of the critical point. Having shown by theoretical arguments that the system belongs to the MF universality class for any positive value of xi and shows nonclassical behavior only for xi=0, we examine to which extent theoretical approximations and simulation can reproduce this behavior. While in this limiting case theoretical approaches are known to fail, we find good agreement for the critical properties between the theoretical approaches and the simulations for xi(2) larger than 0.05.

Volume terms for charged colloids: a grand-canonical treatment

We present a study of thermodynamic properties of suspensions of charged colloids on the basis of linear Poisson-Boltzmann theory. We calculate the effective Hamiltonian of the colloids by integrating out the ionic degrees of freedom grand canonically. This procedure not only yields the well-known pairwise screened-Coulomb interaction between the colloids, but also additional volume terms that affect the phase behavior and the thermodynamic properties, such as the osmotic pressure. These calculations are greatly facilitated by the grand-canonical character of our treatment of the ions and allow for relatively fast computations compared to earlier studies in the canonical ensemble. Moreover, the present derivation of the volume terms are relatively simple, make a direct connection with Donnan equilibrium, yield an explicit expression for the effective screening constant, and allow for extensions to include, for instance, nonlinear effects.

Criticality in charge-asymmetric hard-sphere ionic fluids

Phase separation and criticality are analyzed in z:1 charge-asymmetric ionic fluids of equisized hard spheres by generalizing the Debye-Hückel approach combined with ionic association, cluster solvation by charged ions, and hard-core interactions, following lines developed by Fisher and Levin for the 1:1 case (i.e., the restricted primitive model). Explicit analytical calculations for 2:1 and 3:1 systems account for ionic association into dimers, trimers, and tetramers and subsequent multipolar cluster solvation. The reduced critical temperatures, Tc* (normalized by z), decrease with charge asymmetry, while the critical densities increase rapidly with . The results compare favorably with simulations and represent a distinct improvement over all current theories such as the mean spherical approximation, symmetric Poisson-Boltzmann theory, etc. For z not equal to 1, the interphase Galvani (or absolute electrostatic) potential difference, Deltaphi(T), between coexisting liquid and vapor phases is calculated and found to vanish as absolute value (T-Tc) beta when T-->Tc-with, since our approximations are classical, beta = (1/2). Above Tc, the compressibility maxima and so-called k-inflection loci (which aid the fast and accurate determination of the critical parameters) are found to exhibit a strong z dependence.

Phase coexistence in a polydisperse charged hard-sphere fluid: Polymer mean spherical approximation

We have reconsidered the phase behavior of a polydisperse mixture of charged hard spheres (CHSs) introducing the concept of minimal size neutral clusters. We thus take into account ionic association effects observed in charged systems close to the phase boundary where the properties of the system are dominated by the presence of neutral clusters while the amount of free ions or charged clusters is negligible. With this concept we clearly pass beyond the simple level of the mean spherical approximation (MSA) that we have presented in our recent study of a polydisperse mixture of CHS [Yu. V. Kalyuzhnyi, G. Kahl, and P. T. Cummings, J. Chem. Phys. 120, 10133 (2004)]. Restricting ourselves to a 1:1 and possibly size-asymmetric model we treat the resulting polydisperse mixture of neutral, polar dimers within the framework of the polymer MSA, i.e., a concept that--similar as the MSA--readily can be generalized from the case of a mixture with a finite number of components to the polydisperse case: again, the model belongs to the class of truncatable free-energy models so that we can map the formally infinitely many coexistence equations onto a finite set of coupled, nonlinear equations in the generalized moments of the distribution function that characterizes the system. This allows us to determine the full phase diagram (in terms of binodals as well as cloud and shadow curves), we can study fractionation effects on the level of the distribution functions of the coexisting daughter phases, and we propose estimates on how the location of the critical point might vary in a polydisperse mixture with an increasing size asymmetry and polydispersity.

Phase behavior of the restricted primitive model of ionic fluids with association in slitlike pores. Density-functional approach

We present results of investigations of condensation of restricted primitive model of electrolyte solutions with association between oppositely charged ions confined to slitlike pores. The associative interaction leads to the formation of ionic pairs. It is accounted for by incorporating the first-order thermodynamic perturbation theory into the free energy functional. In order to elucidate the role of association, the phase diagrams are compared with those obtained by us recently [O. Pizio et al., J. Chem. Phys. 121, 11957 (2004)] for the restricted primitive model. The inclusion of the association into the theory leads to lowering the critical temperature for the fluid confined to pores with uncharged and with charged walls. We have observed that the average fraction of bonded ions is high along the coexistence envelope.

Influence of Confinement on the Electrostatic Interaction between Charged Colloids: a (N,V,T) Monte Carlo Study within Hyperspherical Geometry

Monte Carlo simulations within closed hyperspherical geometry are used to analyze the ionic distribution around two confined charged colloids to determine the origin of the net attraction recently reported in the literature. A scaling procedure is used to compare our numerical results obtained with small ideal colloids with the conclusion of the measurements performed with large silica colloids. Although no electrostatic attraction is detected under confinement, our simulations exhibit a significant reduction of the electrostatic repulsion between charged colloids confined between two weakly charged walls. After rescaling to reproduce the electrostatic repulsion between large confined colloids, our numerical results are qualitatively consistent with the reported attraction because we reasonably expect a reduction of the electrostatic force between such confined colloids below the order of magnitude of their van der Waals attraction.

Oppositely charged colloidal binary mixtures: A colloidal analog of the restricted primitive model

The equilibrium phase diagram of a colloidal system composed of 1:1 mixture of positive and negative particles with equal charge is studied by means of Monte Carlo simulations. The system is the colloidal analog of the restricted primitive model (RPM) for ionic fluids. A liquid-gas transition is found in the low-temperature-low-density region, similar to the liquid-gas transition in the RPM. The fluid-crystal transition is also studied, and the liquid phase is shown to be stable in a narrow range of temperatures. In the liquid, the pair distribution function shows alternating layers of particles with opposite sign of charge surrounding every particle. In the vapor phase, clusters of particles are observed, again in agreement with the RPM. However, a decreasing distribution of clusters is obtained, instead of the discrimination between charged and neutral clusters found in the RPM.

Phase coexistence in polydisperse charged hard-sphere fluids: Mean spherical approximation

Taking advantage of the availability of the analytic solution of the mean spherical approximation for a mixture of charged hard spheres with an arbitrary number of components we show that the polydisperse fluid mixture of charged hard spheres belongs to the class of truncatable free energy models, i.e., to those systems where the thermodynamic properties can be represented by a finite number of (generalized) moments of the distribution function that characterizes the mixture. Thus, the formally infinitely many equations that determine the parameters of the two coexisting phases can be mapped onto a system of coupled nonlinear equations in these moments. We present the formalism and demonstrate the power of this approach for two systems; we calculate the full phase diagram in terms of cloud and shadow curves as well as binodals and discuss the distribution functions of the coexisting daughter phases and their charge distributions.

Dipole fluid as a basic model for the equation of state of ionic liquids in the vicinity of their critical point

A model of dipole fluid for ionic liquids similar to molten NaCl is proposed. The estimates for the critical parameters are obtained with the help of the van der Waals equation of state. The influence of the rotation on the characteristics of a dipole pair and the location of the critical point is discussed. The dissociation of such a fluid near the critical point is considered.

Dipolar origin of the gas-liquid coexistence of the hard-core 1:1 electrolyte model

We present a systematic study of the effect of the ion pairing on the gas-liquid phase transition of hard-core 1:1 electrolyte models. We study a class of dipolar dimer models that depend on a parameter R(c), the maximum separation between the ions that compose the dimer. This parameter can vary from sigma(+/-) that corresponds to the tightly tethered dipolar dimer model to R(c)--> infinity that corresponds to the Stillinger-Lovett description of the free ion system. The coexistence curve and critical point parameters are obtained as a function of R(c) by grand-canonical Monte Carlo techniques. Our results show that this dependence is smooth but nonmonotonic and converges asymptotically towards the free ion case for relatively small values of R(c). This fact allows us to describe the gas-liquid transition in the free ion model as a transition between two dimerized fluid phases. The role of the unpaired ions can be considered as a perturbation of this picture.

Critical behavior of ionic liquids

The renormalization of the Landau-Ginzburg Hamiltonian for a system with Coulombic interactions caused by spatially inhomogeneous polarizational effects is discussed. It is shown that for ionic liquids with a strong dependence of the degree of dissociation on density, the nonclassical fluctuation region is significantly narrowed. The essential role of the association of ions is noted.

Universality class of criticality in the restricted primitive model electrolyte

The 1:1 equisized hard-sphere electrolyte or restricted primitive model has been simulated via grand-canonical fine-discretization Monte Carlo. Newly devised unbiased finite-size extrapolation methods using loci in the temperature-density or (T,rho) plane of isothermal rho(2-k) vs pressure inflections, of Q identical with<m(2)>(2)/<m(4)> maxima, and of canonical and C(V) criticality, yield estimates of (T(c),rho(c)) to +/-(0.04,3)%. Extrapolated exponents and Q ratio are (gamma,nu,Q(c)) = [1.24(3), 0.63(3); 0.624(2)], which support Ising (n = 1) behavior with (1.23(9), 0.630(3); 0.623(6)), but exclude classical, XY (n = 2), self-avoiding walk (n = 0), and n = 1 criticality with potentials varphi(r)>Phi/r(4.9) when r-->infinity.

Phase equilibria of charge-, size-, and shape-asymmetric model electrolytes

The low-temperature phase behavior of two 2:1 hard-core electrolyte models has been investigated by Monte Carlo simulations. In the first model, both bivalent cations and monovalent anions are spherical, and the charges are located at the ion's centers; in the second model, bivalent cations are modeled as rigid dimers composed of two tangent hard spheres, each carrying a positive charge at the center. It is found that the critical temperature and the critical density are strongly affected by the size asymmetry and the shape of the ions. The results presented in this work provide insights into the behavior of charged colloidal suspensions and polyelectrolytes, where large, symmetric or asymmetric ionic species carrying like charges can attract each other and give rise to thermodynamically unstable conditions.

Phase transitions in 2:1 and 3:1 hard-core model electrolytes

Critical temperatures, T(c), densities, rho(c), and coexistence curves for 2:1 and 3:1 hard-core model electrolytes have been found by fine-discretization Monte Carlo simulation. The size ratio of + and - ions strongly affects T(c) and rho(c); the trends contradict most current theories. Large multivalent ions screened by small monovalent counterions exhibit normal gas-liquid transitions of direct relevance to phase separation in charge-stabilized colloids. Conversely, extrapolation suggests the absence of such transitions for sufficiently small multivalent ions with large monovalent counterions.

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.

Criticality in confined ionic fluids

A theory of a confined two-dimensional electrolyte is presented. The positive and negative ions, interacting by a 1/r potential, are constrained to move on an interface separating two solvents with dielectric constants epsilon(1) and epsilon(2). It is shown that the Debye-Hückel type of theory predicts that this two-dimensional Coulomb fluid should undergo a phase separation into a coexisting liquid (high-density) and gas (low-density) phases. We argue, however, that the formation of polymerlike chains of alternating positive and negative ions can prevent this phase transition from taking place.

Precise simulation of criticality in asymmetric fluids

Extensive grand canonical Monte Carlo simulations have been performed for the hard-core square-well fluid with interaction range b=1.5 sigma. The critical exponent for the correlation length has been estimated in an unbiased fashion as nu=0.63+/-0.03 via finite-size extrapolations of the extrema of properties measured along specially constructed, asymptotically critical loci that represent pseudosymmetry axes. The subsequent location of the critical point achieves a precision of five parts in 10(4) for Tc and about 0.3% for the critical density rhoc. The effective exponents gamma+(eff) and beta(eff) indicate Ising-type critical-point values to within 2% and 5.6%, respectively, convincingly distinguishing the universality class from the "nearby" XY and n=0 (self-avoiding walk) classes. Simulations of the heat capacity CV(T,rho) and d2psigma/dT2, where psigma is the vapor pressure below Tc, suggest a negative but small Yang-Yang anomaly, i.e., a specific-heat-like divergence in the corresponding chemical potential derivative (d2 musigma/dT2) that requires a revision of the standard asymptotic scaling description of asymmetric fluids.

Phase equilibria of size-asymmetric primitive model electrolytes

The low-temperature phase coexistence of size-asymmetric primitive model electrolyte solutions has been investigated by means of Monte Carlo simulations. Binodal curves and critical parameters are reported as a function of size ratio lambda = sigma(+)/sigma(-) in the range 0.05 to 1. Both the critical temperature and the critical density are found to decrease as lambda decreases. These trends are in conflict with available theoretical predictions. For highly asymmetric systems, the depression of the critical parameters is accompanied by the formation of large chainlike and ringlike structures, thereby giving rise to considerable finite-size effects.

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

Criticality and crossover in accessible regimes

The near-critical behavior of (d = 3)-dimensional Ising-model ferromagnets or simple lattice gases with equivalent first, second, and third nearest-neighbor interactions is studied through Monte Carlo simulations using histogram reweighting techniques and comparisons with series expansions. By carefully analyzing numerical data from relatively small finite systems using scaling and extrapolation methods, it is demonstrated that one can reliably estimate critical exponents, critical temperatures, and universal amplitude ratios, thereby distinguishing convincingly between different "nearby" universality classes and revealing systematic crossover effects. This study is preparatory to extending similar techniques to study criticality in continuum fluid models lacking symmetries, with Coulomb interactions, etc.