Coexistence and criticality in size-asymmetric hard-core electrolytes

Effect of Bjerrum pairs on the electrostatic properties of an electrolyte solution near charged surfaces: a mean-field approach

In this paper, we investigate the consequences of ion association, coupled with the considerations of finite size effects and orientational ordering of Bjerrum pairs as well as ions and water molecules, on the electric double layer near charged surfaces. Based on the lattice statistical mechanics accounting for finite sizes and dipole moments of ions, Bjerrum pairs and solvent molecules, we consider the formation of Bjerrum pairs and derive the mathematical expressions for Bjerrum pair number density as well as cation/anion number density and water molecule number density. We reveal several significant phenomena. Firstly, it is shown that our approach naturally yields the equilibrium constant for dissociation-association equilibrium between Bjerrum pairs and ions. Secondly, at low surface charge densities, an increase in the bulk concentration of Bjerrum pairs enhances the permittivity and decreases the differential capacitance. Next, for the cases where Bjerrum pairs in an alcohol electrolyte solution have a high value of dipole moment, the Bjerrum pair number density increases with decreasing distance from the charged surface, and the differential capacitance and permittivity are high compared to those for the cases with lower values of Bjerrum-pair dipole moments. Finally, we show that the difference in the concentration and dipole moment of Bjerrum pairs can lead to some variation in osmotic pressure between two similarly charged surfaces.

Theory of electrolytes including steric, attractive, and hydration interactions

We present a continuum theory of electrolytes composed of a waterlike solvent and univalent ions. First, we start with a density functional F for the coarse-grained solvent, cation, and anion densities, including the Debye-Hückel free energy, the Coulombic interaction, and the direct interactions among these three components. These densities fluctuate obeying the distribution ∝exp(-F/k_BT). Eliminating the solvent density deviation in F, we obtain the effective non-Coulombic interactions among the ions, which consist of the direct ones and the solvent-mediated ones. We then derive general expressions for the ion correlation, the apparent partial volume, and the activity and osmotic coefficients up to linear order in the average salt density n_s. Second, we perform numerical analysis using the Mansoori-Carnahan-Starling-Leland model [J. Chem. Phys. 54, 1523 (1971)] for three-component hardspheres. The effective interactions sensitively depend on the cation and anion sizes due to competition between the steric and hydration effects, which are repulsive between small-large ion pairs and attractive between symmetric pairs. These agree with previous experiments and Collins' rule [Biophys. J. 72, 65 (1997)]. We also give simple approximate expressions for the ionic interaction coefficients valid for any ion sizes.

Self-assembly of charged colloidal cubes

In this work, we show how and why the interactions between charged cubic colloids range from radially isotropic to strongly directionally anisotropic, depending on tuneable factors. Using molecular dynamics simulations, we illustrate the effects of typical solvents to complement experimental investigations of cube assembly. We find that in low-salinity water solutions, where cube self-assembly is observed, the colloidal shape anisotropy leads to the strongest attraction along the corner-to-corner line, followed by edge-to-edge, with a face-to-face configuration of the cubes only becoming energetically favorable after the colloids have collapsed into the van der Waals attraction minimum. Analysing the potential of mean force between colloids with varied cubicity, we identify the origin of the asymmetric microstructures seen in experiment.

Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids

Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. The theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the "accessibility" of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.

Critical charge and density coupling in ionic spherical models

We investigate ionic criticality on the basis of a specially devised spherical model that accounts both for Coulomb and nonionic forces in binary systems. We show in detail here the consequences of the entanglement of density and charge correlation functions G_{NN} and G_{ZZ} on criticality and screening. We also show on this soluble model how, because of electroneutrality, the long-range Coulomb interactions do not change the universality class of criticality in the model driven primarily by sufficiently attractive nonionic interactions. Near criticality, G_{NN} and G_{ZZ} are fully decoupled in charge symmetric systems. However, in more realistic nonsymmetric models, charge and density fluctuations couple in leading order so that the charge and density correlation lengths diverge asymptotically in a similar way. Similarly, the Stillinger-Lovett sum rule, which characterizes a conducting fluid, is violated at criticality in nonsymmetric models when the critical-point density-decay exponent η vanishes. In addition, if quantum effects are accounted for semiclassically by incorporating algebraically decaying interactions, G_{ZZ} decays only as a power law in the whole phase space, contrary to the usually expected exponential Debye screening. We expect these results on this soluble toy model to be general and to reveal general mechanisms ruling ionic criticality.

Vapor-liquid phase behavior of a size-asymmetric model of ionic fluids confined in a disordered matrix: The collective-variables-based approach

We develop a theory based on the method of collective variables to study the vapor-liquid equilibrium of asymmetric ionic fluids confined in a disordered porous matrix. The approach allows us to formulate the perturbation theory using an extension of the scaled particle theory for a description of a reference system presented as a two-component hard-sphere fluid confined in a hard-sphere matrix. Treating an ionic fluid as a size- and charge-asymmetric primitive model (PM) we derive an explicit expression for the relevant chemical potential of a confined ionic system which takes into account the third-order correlations between ions. Using this expression, the phase diagrams for a size-asymmetric PM are calculated for different matrix porosities as well as for different sizes of matrix and fluid particles. It is observed that general trends of the coexistence curves with the matrix porosity are similar to those of simple fluids under disordered confinement, i.e., the coexistence region gets narrower with a decrease of porosity and, simultaneously, the reduced critical temperature T_{c}^{*} and the critical density ρ_{i,c}^{*} become lower. At the same time, our results suggest that an increase in size asymmetry of oppositely charged ions considerably affects the vapor-liquid diagrams leading to a faster decrease of T_{c}^{*} and ρ_{i,c}^{*} and even to a disappearance of the phase transition, especially for the case of small matrix particles.

Dielectric effects in the self-assembly of binary colloidal aggregates

Electrostatic interactions play an important role in numerous self-assembly phenomena, including colloidal aggregation. Although colloids typically have a dielectric constant that differs from the surrounding solvent, the effective interactions that arise from inhomogeneous polarization charge distributions are generally neglected in theoretical and computational studies. We introduce an efficient technique to resolve polarization charges in dynamical dielectric geometries, and demonstrate that dielectric effects qualitatively alter the predicted self-assembled structures, with surprising colloidal strings arising from many-body effects.

Thermo-molecular orientation effects in fluids of dipolar dumbbells

Plots of first-order (left) and novel second-order (right) thermomolecular orientation effects in fluids of dipolar dumbbells.

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.

Spatial inhomogeneities in ionic liquids, charged proteins, and charge stabilized colloids from collective variables theory

Effects of size and charge asymmetry between oppositely charged ions or particles on spatial inhomogeneities are studied for a large range of charge and size ratios. We perform a stability analysis of the primitive model of ionic systems with respect to periodic ordering using the collective variables-based theory. We extend previous studies [Ciach et al., Phys. Rev. E 75, 051505 (2007)] in several ways. First, we employ a nonlocal approximation for the reference hard-sphere fluid which leads to the Percus-Yevick pair direct correlation functions for the uniform case. Second, we use the Weeks-Chandler-Anderson regularization scheme for the Coulomb potential inside the hard core. We determine the relevant order parameter connected with the periodic ordering and analyze the character of the dominant fluctuations along the λ lines. We show that the above-mentioned modifications produce large quantitative and partly qualitative changes in the phase diagrams obtained previously. We discuss possible scenarios of the periodic ordering for the whole range of size and charge ratios of the two ionic species, covering electrolytes, ionic liquids, charged globular proteins or nanoparticles in aqueous solutions, and charge-stabilized colloids.

Thermodynamics and diffusion in size-symmetric and asymmetric dense electrolytes

MD simulation results for model size-symmetric and asymmetric electrolytes at high densities and temperatures (well outside the liquid-gas coexistence region) are generated and analyzed focusing on thermodynamic and diffusion properties. An extension of the mean spherical approximation for electrolytes originally derived for charged hard sphere fluids is adapted to these systems by exploiting the separation of short range and Coulomb interaction contributions intrinsic to these theoretical models and is found to perform well for predicting equation of state quantities. The diffusion coefficients of these electrolytes can also be reasonably well predicted using entropy scaling ideas suitably adapted to charged systems and mixtures. Thus, this approach may provide an avenue for studying dense electrolytes or complex molecular systems containing charged species at high pressures and temperatures.

Microscopic approach to the kinetics of pattern formation of charged molecules on surfaces

A microscopic formalism based on computing many-particle densities is applied to the analysis of the diffusion-controlled kinetics of pattern formation in oppositely charged molecules on surfaces or adsorbed at interfaces with competing long-range Coulomb and short-range Lennard-Jones interactions. Particular attention is paid to the proper molecular treatment of energetic interactions driving pattern formation in inhomogeneous systems. The reverse Monte Carlo method is used to visualize the spatial molecular distribution based on the calculated radial distribution functions (joint correlation functions). We show the formation of charge domains for certain combinations of temperature and dynamical interaction parameters. The charge segregation evolves into quasicrystalline clusters of charges, due to the competing long- and short-range interactions. The clusters initially co-exist with a gas phase of charges that eventually add to the clusters, generating "fingers" or line of charges of the same sign, very different than the nanopatterns expected by molecular dynamics in systems with competing interactions in two dimensions, such as strain or dipolar versus van der Waals interactions.

Gas-liquid critical parameters of asymmetric models of ionic fluids

The effects of size and charge asymmetry on the gas-liquid critical parameters of a primitive model (PM) of ionic fluids are studied within the framework of the statistical field theory based on the collective variables method. Recently, this approach has enabled us to obtain the correct trends of the both critical parameters of the equisize charge-asymmetric PM without assuming ionic association. In this paper, we focus on the general case of an asymmetric PM characterized by the two parameters: hard-sphere diameter, lambda=sigma+/sigma-, and charge, z=q+/|q-|, ratios of the two ionic species. We derive an explicit expression for the chemical potential conjugate to the order parameter which includes the effects of correlations up to the third order. Based on this expression we consider the three versions of PM: a monovalent size-asymmetric PM (lambda not equal 1, z=1) , an equisize charge-asymmetric PM (lambda=1, z not equal 1) and a size- and charge-asymmetric PM (lambda not equal 1, z=2) . Similar to simulations, our theory predicts that the critical temperature and the critical density decrease with the increase in size asymmetry. Regarding the effects of charge asymmetry, we obtain the correct trend of the critical temperature with z , while the trend of the critical density obtained in this approximation is inconsistent with simulations, as well as with our previous results found in the higher-order approximation. We expect that the consideration of the higher-order correlations will lead to the correct trend of the critical density with charge asymmetry.

Phase diagrams of charged colloids from thermodynamic integration

We present full phase diagrams (including solid phases) of spherical charged colloids, using Monte Carlo sampling and thermodynamic integration of the Helmholtz free energy. Colloids and their co- and counterions are described by the primitive model for ionic systems that consists of hard-spheres with central point charges, while the solvent is taken into account solely through its dielectric constant. Two systems are considered: (i) a size-asymmetric system of oppositely charged spheres with size ratios q = 0.3 and 0.5 and (ii) a charge- and size-asymmetric system with colloid charge Q = 10 and counterions of charge -1 in the presence of monovalent added salt. In system (i), for both size ratios, the stable solid phase is equivalent to the NaCl crystal where the oppositely charged spheres take the lattice positions of Na and Cl ions. In system (ii), the phase diagram consists of gas-liquid and fluid-solid coexistence regions. We show that added salt stabilizes the fluid phase and shrinks the fluid-solid coexistence region, in agreement with experimental and theoretical results.

Inflation of the screening length induced by Bjerrum pairs

Within a modified Poisson-Boltzmann theory we study the effect of Bjerrum pairs on the typical length scale [Formula: see text] over which electric fields are screened in electrolyte solutions, taking into account a simple association-dissociation equilibrium between free ions and Bjerrum pairs. At low densities of Bjerrum pairs, this length scale is well approximated by the Debye length [Formula: see text], with ρ(s) the free-ion density. At high densities of Bjerrum pairs, however, we find [Formula: see text], which is significantly larger than 1/κ due to the enhanced effective permittivity of the electrolyte, caused by the polarization of Bjerrum pairs. We argue that this mechanism may explain the recently observed anomalously large colloid-free zones between an oil-dispersed colloidal crystal and a colloidal monolayer at the oil-water interface.

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.

Liquid-liquid phase transition in solutions of ionic liquids with halide anions: Criticality and corresponding states

Measurements of the liquid-liquid phase diagrams of solutions of the ionic liquids (ILs) 1-dodecyl-3-methylimidazolium chloride (C12mimCl) in arenes (benzene, toluene, o-xylene, tetraline) and 1-tetradecyl-3-methylimidazolium chloride (C14mimCl) in CCl4 are reported and compared with those of solutions of trihexyl-tetradecyl-phosphonium halides (P666 14Cl, P666 14Br) in hydrocarbons and 1-alkyl-3-methylimidazolium tetrafluoroborates (CnmimBF4) in alcohols and water. The phase diagrams of solutions of tetrapentyl-ammonium bromide (N5555Br) in water and KI in SO2 are also discussed. Except for the KI/SO2 system, which features a lower critical solution point (LCSP), all systems have an upper critical solution point (UCSP) and show corresponding-states behavior. The experimental data are compared with results from simulations and theory concerning the model fluid of charged hard spheres in a dielectric continuum, termed restricted primitive model (RPM). The analysis in terms of of RPM variables shows agreement with the location of the critical point (CP) of the model with noticeable systematic deviations. However, for protic solvents, the CP becomes an LCSP, while in aprotic solvents the CP is a UCSP as expected for Coulomb systems. This indicates that in aprotic solvents, the phase transition is essentially determined by the Coulomb interactions, while in the solutions in protic solvents with hydrogen bonds, both Coulomb and solvophobic interactions are important.

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.

Symmetry, equivalence, and molecular self-assembly

Molecular self-assembly at equilibrium is fundamental to the fields of biological self-organization, the development of novel environmentally responsive polymeric materials, and nanofabrication. Our approach to understanding the principles governing this process is inspired by existing models and measurements for the self-assembly of actin, tubulin, and the ubiquitous icosahedral shell structures of viral capsids. We introduce a family of simple potentials that give rise to the self-assembly of linear polymeric, random surface ("membrane"), tubular ("nanotube"), and hollow icosahedral structures that are similar in many respects to their biological counterparts. The potentials involve equivalent particles and an interplay between directional (dipolar, multipolar) and short-range (van der Waals) interactions. Specifically, we find that the dipolar potential, having a continuous rotational symmetry about the dipolar axis, gives rise to chain formation, while particles with multipolar potentials, having discrete rotational symmetries (square quadrupole or triangular ring of dipoles or "hexapole"), lead to the self-assembly of open sheet, nanotube, and hollow icosahedral geometries. These changes in the geometry of self-assembly are accompanied by significant changes in the kinetics of the organization.

Liquid-gas separation in colloidal electrolytes

The liquid-gas transition of an electroneutral mixture of oppositely charged colloids, studied by Monte Carlo simulations, is found in the low-temperature-low-density region. The critical temperature shows a nonmonotonous behavior as a function of the interaction range, kappa(-1), with a maximum at kappasigma approximately 10, implying an island of coexistence in the kappa-rho plane. The system is arranged in such a way that each particle is surrounded by shells of particles with alternating charge. In contrast with the electrolyte primitive model, both neutral and charged clusters are obtained in the vapor phase.

Melting line of charged colloids from primitive model simulations

We develop an efficient simulation method to study suspensions of charged spherical colloids using the primitive model. In this model, the colloids and the co- and counterions are represented by charged hard spheres, whereas the solvent is treated as a dielectric continuum. In order to speed up the simulations, we restrict the positions of the particles to a cubic lattice, which allows precalculation of the Coulombic interactions at the beginning of the simulation. Moreover, we use multiparticle cluster moves that make the Monte Carlo sampling more efficient. The simulations are performed in the semigrand canonical ensemble, where the chemical potential of the salt is fixed. Employing our method, we study a system consisting of colloids carrying a charge of 80 elementary charges and monovalent co- and counterions. At the colloid densities of our interest, we show that lattice effects are negligible for sufficiently fine lattices. We determine the fluid-solid melting line in a packing fraction eta-inverse screening length kappa plane and compare it with the melting line of charged colloids predicted by the Yukawa potential of the Derjaguin-Landau-Verwey-Overbeek theory. We find qualitative agreement with the Yukawa results, and we do not find any effects of many-body interactions. We discuss the difficulties involved in the mapping between the primitive model and the Yukawa model at high colloid packing fractions (eta>0.2).

Size-asymmetric primitive model at low temperature: description of ion pairing and location of the critical point

We argue that Bjerrum's approach to ion pairing is inappropriate for the size-asymmetric primitive model in the neighborhood of its critical point, and propose a new approach based on the Stillinger-Lovett pairing procedure. The new approach recursively scales up the ion size until linear approximations are suitable for analyzing such a model. To locate the critical point, a residual van der Waals interaction between pairs is added, with an energy cutoff adjusted to match the critical temperature of the restricted primitive model. The locations and downward trends of T(c) and rho(c) with asymmetry are found to compare favorably with simulations.

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.

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.

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.

A modified Poisson–Boltzmann analysis of the capacitance behavior of the electric double layer at low temperatures

The modified Poisson-Boltzmann theory is used to analyze the anomalous behavior of the electric double layer capacitance for small surface charge at low temperatures and densities. Good agreement is found with simulation and recent density-functional theory results. Negative adsorption is also found in line with theory and simulation. An unsatisfactory feature is the relatively poor structure in this region due to the inherent approximations in the theory. This feature is unimportant in relation to the capacitance results but has implications when calculating adsorption properties.

Monte Carlo simulations of oppositely charged macroions in solution

The structure and phase behavior of oppositely charged macroions in solution have been studied with Monte Carlo simulations using the primitive model where the macroions and small ions are described as charged hard spheres. Size and charge symmetric, size asymmetric, and charge asymmetric macroions at different electrostatic coupling strengths are considered, and the properties of the solutions have been examined using cluster size distribution functions, structure factors, and radial distribution functions. At increasing electrostatic coupling, the macroions form clusters and eventually the system displays a phase instability, in analogy to that of simple electrolyte solutions. The relation to the similar cluster formation and phase instability occurring in solutions containing oppositely charged polymers is also discussed.

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.

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

Molecular dynamics simulations of the surface tension of ionic liquids

We report molecular dynamics computer simulations of the surface tension and interfacial thickness of ionic liquid-vapor interfaces modeled with a soft core primitive model potential. We find that the surface tension shows an anomalous oscillatory behavior with interfacial area. This observation is discussed in terms of finite size effects introduced by the periodic boundary conditions employed in computer simulations. Otherwise we show that the thickness of the liquid-vapor interface increases with surface area as predicted by the capillary wave theory. Data on the surface tension of size-asymmetric ionic liquids are reported and compared with experimental data of molten salts. Our data suggest that the surface tensions of size-asymmetric ionic liquids do not follow a corresponding states law.