The liquid–liquid phase transition in ionic solutions: Coexistence curves of tetra-n-butylammonium pricrate in alkyl alcohols

Phase Transitions and Electrochemical Properties of Ionic Liquids and Ionic Liquid-Solvent Mixtures

Recent advances in studies of ionic liquids (IL) and ionic liquid-solvent mixtures are reviewed. Selected experimental, simulation, and theoretical results for electrochemical, thermodynamical, and structural properties of IL and IL-solvent mixtures are described. Special attention is paid to phenomena that are not predicted by the classical theories of the electrical double layer or disagree strongly with these theories. We focus on structural properties, especially on distribution of ions near electrodes, on electrical double layer capacitance, on effects of confinement, including decay length of a dissjoining pressure between confinig plates, and on demixing phase transition. In particular, effects of the demixing phase transition on electrochemical properties of ionic liquid-solvent mixtures for different degrees of confinement are presented.

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.

Monte Carlo simulation of the nonadditive restricted primitive model of ionic fluids: phase diagram and clustering

We report an accurate Monte Carlo calculation of the phase diagram and clustering properties of the restricted primitive model with nonadditive hard-sphere diameters. At high density the positively nonadditive fluid shows more clustering than in the additive model and the negatively nonadditive fluid shows less clustering than in the additive model; at low density the reverse scenario appears. A negative nonadditivity tends to favor the formation of neutrally charged clusters starting from the dipole. A positive nonadditivity favors the pairing of like ions at high density. The critical point of the gas-liquid phase transition moves at higher temperatures and higher densities for a negative nonadditivity and at lower temperatures and lower densities for a positive nonadditivity. The law of corresponding states does not seem to hold strictly. Our results can be used to interpret recent experimental works on room temperature ionic liquids.

The critical behavior of the refractive index near liquid-liquid critical points

The nature of the critical behavior in the refractive index n is revisited in the framework of the complete scaling formulation. A comparison is made with the critical behavior of n as derived from the Lorentz-Lorenz equation. Analogue anomalies to those predicted for the dielectric constant ε, namely, a leading |t|(2β) singularity in the coexistence-curve diameter in the two-phase region and a |t|(1-α) along the critical isopleth in the one phase region, are expected in both cases. However, significant differences as regards the amplitudes of both singularities are obtained from the two approaches. Analysis of some literature data along coexistence in the two-phase region and along the critical isopleth in the one-phase region provide evidence of an intrinsic effect, independent of the density, in the critical anomalies of n. This effect is governed by the shift of the critical temperature with an electric field, which is supposed to take smaller values at optical frequencies than at low frequencies in the Hz to MHz range.

Asymmetric criticality in weakly compressible liquid mixtures

The thermodynamics of asymmetric liquid-liquid criticality is updated by incorporating pressure effects into the complete-scaling formulation earlier developed for incompressible liquid mixtures [C. A. Cerdeirina et al., Chem. Phys. Lett. 424, 414 (2006); J. T. Wang et al., Phys. Rev. E 77, 031127 (2008)]. Specifically, we show that pressure mixing enters into weakly compressible liquid mixtures as a consequence of the pressure dependence of the critical parameters. The theory is used to analyze experimental coexistence-curve data in the mole fraction-temperature, density-temperature, and partial density-temperature planes for a large number of binary liquid mixtures. It is shown how the asymmetry coefficients in the scaling fields are related to the difference in molecular volumes of the two liquid components. The work resolves the question of the so-called "best order parameter" discussed in the literature during the past decades.

Liquid-liquid phase separation in solutions of ionic liquids: phase diagrams, corresponding state analysis and comparison with simulations of the primitive model

Phase diagrams of ionic solutions of the ionic liquid C(18)mim(+)NTF(2)(-) (1-n-octadecyl-3-methyl imidazolium bistrifluormethylsulfonylimide) in decalin, cyclohexane and methylcyclohexane are reported and compared with that of solutions of other imidazolium ionic liquids with the anions NTF(2)(-), Cl(-) and BF4(-) in arenes, CCl(4), alcohols and water. The phase diagrams are analysed presuming Ising criticality and taking into account the asymmetry of the phase diagrams. The resulting parameters are compared with simulation results for equal-sized charged hard spheres in a dielectric continuum, the restricted primitive model (RPM) and the primitive model (PM) that allows for ions of different size. In the RPM temperature scale the critical temperatures vary almost linearly with the dielectric permittivity of the solvent. The RPM critical temperatures of the solutions in non-polar solvents are very similar, somewhat below the RPM value. Correlations with the boiling temperatures of the solvents and a dependence on the length of the side chain of the imidazolium cations show that dispersion interactions modify the phase transition, which is mainly determined by Coulomb forces. Critical concentrations, widths of the phase diagrams and the slopes of the diameter are different for the solutions in protic and aprotic solvents. The phase diagrams of the solutions in alcohols and water get a lower critical solution point when represented in RPM variables.

Phase Separation in Solutions of Room Temperature Ionic Liquids in Hydrocarbons

The room temperature ionic liquids (RTIL) trihexyl-tetradecyl phosphonium chloride (P666 14Cl) and the bromide (P666 14Br) are soluble in hydrocarbons. The investigated solutions in heptane, octane, nonane and decane show liquid–liquid phase separation with an upper critical solution point at ambient temperatures at molar fractions near 0.03 of the salt. Phase diagrams are reported and analysed presuming Ising criticality. The critical temperatures and the critical densities increase with the chain length of the hydrocarbons, where the figures corresponding to the bromides are above that of the chlorides. Scaled by the critical data the phase diagrams show corresponding state behaviour. In accordance with the prediction of the restricted primitive model (RPM), which is a model fluid of equal sized, charged hard spheres in a dielectric continuum, the critical points are located at low temperature and low concentration, when the corresponding state variables of this model are used. However, the critical temperature T c * and the critical density ρc * are well below the figures of the RPM prediction. Comparison is made with the phase diagrams of alcohol solutions of imidazolium ionic liquids and with simulation results of the RPM.

Closed Solubility Loops in Liquid Mixtures

Based on the principle of isomorphic critical behavior of fluid mixtures we derive a simple equation for closed solubility loops in a temperature–concentration diagram at constant pressure. The validity of the simple equation is tested by comparisons with experimental solubility data for liquid mixtures at atmospheric pressure and at elevated pressures. We also elucidate the nature of liquid–liquid phase separation in the vicinity of a double critical point. Specifically, we show how from an expansion around the critical double point one can develop a description of solubility loops as a function of pressure.

Critical behavior on approaching a special critical point in a complex fluid

The critical behavior of osmotic susceptibility is investigated in the re-entrant complex mixture 1-propanol (P)+water (W)+potassium chloride (KCl) through light-scattering measurements. The measurements are performed on approaching a special critical point [i.e., the double critical point (DCP)] in this mixture, along the line of upper critical solution temperatures (T(U)'s), by varying t from the high temperature one-phase region. The light-scattering data analysis emphasizes the need for corrections to the asymptotic Ising behavior and yields very large magnitudes for the correction-to-scaling amplitudes A(1) and A(2), with the first-correction amplitude A(1) being negative, signifying a nonmonotonic crossover behavior of the susceptibility exponent in this mixture. For the T(U) closest to the DCP, the effective susceptibility exponent gamma(eff) displays a nonmonotonic crossover from its nearly doubled three dimensional (3D)-Ising value toward its nearly double mean-field value with an increase in t. While for that far away from the DCP, it displays a nonmonotonic crossover from its single-limit Ising value toward a value slightly lower than its mean-field value of 1 with an increase in t. This feature of the effective susceptibility exponent is interpreted in terms of the possibility of a nonmonotonic crossover to the mean-field value from lower values in the nonasymptotic high t region. The renormalized Ising regime extends over a larger t range for the sample (or T(U)) closest to the DCP when compared to that far away from it. The in-between T(U)'s display a trend toward shrinkage in the renormalized Ising regime as T(U) shifts away from the DCP. Nevertheless, the crossover to the mean-field behavior is completed only beyond t>10(-2) for the T(U)'s studied. The observed crossover behavior is attributed to the presence of strong ion-induced clustering in this mixture, as revealed by various structure probing techniques, while the observed unique trend in the crossover behavior is discussed in terms of the varying influence of the DCP on the critical behavior along the T(U) line. The crossover behavior for the T(U)'s is pronounced and more sharp compared to the T(L)'s (lower critical solution temperatures) [U. K. Pradeep, J. Chem. Phys. 129, 134506 (2008)] in this mixture, although there exists no difference in the growth of the mesoscale clusters in the lower and upper one-phase regions in this mixture. Our observations suggest the need to look at the crossover behavior probably from two perspectives, namely, the dielectric effect and the clustering effect. The effective susceptibility exponent as a function of the field variable t(UL), instead of the conventional variable t, displays a sharp nonmonotonic crossover from its asymptotic 3D-Ising value ( approximately 1.24) toward a value slightly lower than its nonasymptotic mean-field value of 1, as that observed in the t analysis for the T(U) far away from the influence of the DCP.

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.

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.

Universality of ionic criticality: size- and charge-asymmetric electrolytes

Grand-canonical simulations designed to resolve critical universality classes are reported for z:1 hard-core electrolyte models with diameter ratios lambda=a+/a- less than or approximately equal 6. For z=1 Ising-type behavior prevails. Unbiased estimates of Tc(lambda) are within 1% of previous (biased) estimates but the critical densities are approximately 5% lower. Ising character is also established for the 2:1 and 3:1 equisized models, along with critical amplitudes and improved Tc estimates. For z=3, however, strong finite-size effects reduce the confidence level although classical and O (n>or=3) criticality are excluded.

Critical point of electrolyte mixtures

The critical behavior of electrolyte mixtures was studied using grand canonical Monte Carlo simulations. Mixtures consist of large multivalent macroions and small monovalent co- and counterions. The system can be viewed as a binary mixture of macroions (with their counterions) and salt (co- and counterion pair). The primitive model description was used, in which the ions are point charges with a hard core and the solvent is treated as a uniform dielectric continuum. The grand canonical simulations are based on insertions and removals of neutral molecules: macroion with its counterions or coions and a counterion. We propose a distance biasing method that enables direct grand canonical simulations up to charge asymmetry of 10:1. We calculated the critical loci that connect the salt-free state, which consists of only macroions and counterions, with the pure salt state using mixed-field finite-size scaling with no pressure mixing. The critical parameters are determined for macroion to counterion charge asymmetries of 2:1, 3:1, and 10:1. Our results suggest that binary electrolyte mixtures are type-I mixtures, where the two components mix continuously.

Supercritical Water as a Solvent

Water is not restricted to moderate temperatures and low pressures, but can exist up to very high temperatures, far above its critical point at 647 K. In this supercritical regime, water can be gradually compressed from gas-like to liquid-like densities. The resulting dense supercritical states have extraordinary properties which can be tuned by temperature and pressure, and form the basis for innovative technologies. This Review covers the current knowledge of the major properties of supercritical water and its solutions with nonpolar, polar, and ionic compounds, and of the underlying molecular processes.

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.

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.

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.

Discretization dependence of criticality in model fluids: a hard-core electrolyte

Grand-canonical simulations at various levels, zeta=5-20, of fine-lattice discretization are reported for the near-critical 1:1 hard-core electrolyte or restricted primitive model (RPM). With the aid of finite-size scaling analyses, it is shown convincingly that, contrary to recent suggestions, the universal critical behavior is independent of zeta (> or approximately 4), thus the continuum (zeta--> infinity ) RPM exhibits Ising-type (as against classical, self-avoiding walk, XY, etc.) criticality. A general consideration of lattice discretization provides effective extrapolation of the intrinsically erratic zeta dependence, yielding (T*(c),rho*(c)) approximately equal to (0.0493(3),0.075) for the zeta=infinity RPM.

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.

Critical concentration fluctuations of the ionic binary mixture ethylammonium nitrate-n-octanol: an ultrasonic spectrometry study

Between 200 kHz and 130 MHz, the ultrasonic attenuation spectrum of the ionic ethylammonium nitrate--n-octanol mixture of critical composition has been measured at various reduced temperatures (1.5 x 10(-4)<or= t <or=5.7 x 10(-2)). Using static and dynamic light scattering data as well as viscosity and heat capacity data from the literature, the experimental spectra have been evaluated to yield the scaling function, with the background contribution to the spectra as the only adjustable parameter. Agreement, within the limits of experimental error, of the measured scaling function with that of the nonionic binary system ethanol--dodecane and with the theoretical predictions of the Bhattacharjee-Ferrell dynamic scaling model is found. The amplitude of the fluctuation correlation length xi(o) (=0.47 nm) and the amount of the coupling constant /g/ (=1.3) are rather high as compared to nonionic binary critical mixtures. The amplitude of the relaxation rate of order parameter fluctuations Gamma(o)(=2.6 x 10(8) s(-1)) exhibits an unusual small value, likely to the most part a reflection of the high viscosity and thus small diffusion coefficient of the ionic liquid.

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

Heat capacity behavior in the critical region of the ionic binary mixture ethylammonium nitrate-n-octanol

At temperatures between 30 degrees C and 58 degrees C we have recorded the heat capacity of the ionic ethylammonium nitrate-n-octanol mixture of critical composition and also of the constituents. Different samples of the binary mixture have been measured with its upper critical demixing temperature T(c) varying between 41.04 degrees C and 46.87 degrees C, depending on small traces of water within the liquid under test. Almost identical heat capacity profiles result if the data are displayed as a function of the temperature distance to the actual T(c) value. In the homogeneous phase the critical contribution exhibits power-law behavior with the critical exponent alpha=0.11 as theoretically predicted for nonionic liquids and in conformity with our understanding of the ionic criticality as being asymptotically Ising like. The noncritical background part of the heat capacity can be related to the heat capacities of the constituents using a simple mixture relation. In the two-phase regime a series of almost perfectly reproducible events is found which may be taken to indicate the existence of nonequilibrium intermediate states in the ethylammonium nitrate-n-octanol system.

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.