Criticality in charge-asymmetric hard-sphere ionic fluids

Closure to the PRISM equation derived from nonlinear response theory

Nonlinear response theory is employed to derive a closure to the polymer reference interaction site model equation. The closure applies to a liquid of neutral polymers at melt densities. It can be considered a molecular generalization of the mean spherical approximation (MSA) closure of Lebowitz and Percus to the atomic Ornstein-Zernike (OZ) equation and is similar in some aspects to the reference "molecular" MSA (R-MMSA) closure of Schweizer and Yethiraj to PRISM. For a model binary blend of freely-jointed chains, the new closure predicts an unmixing critical temperature, Tc, via the susceptibility route that scales linearly with molecular weight, N, in agreement with Flory theory. Predictions for Tc of the new closure differ greatest from those of the R-MMSA at intermediate N, the latter being about 40% higher than the former there, but at large N, both theories give about the same values. For an isotopic blend of polyethylene, the new and R-MMSA closures predict a Tc about 25% higher than the experimental value, which is only moderately less accurate than the prediction of atomic OZ-MSA theory for Tc of methane. In this way, the derivation and its consequences help to identify the ingredients in a theory needed to properly model the equilibrium properties of a polymeric liquid at both short and long lengthscales.

Ion partitioning and permeation in charged low-T* membranes

Understanding ion transport in membrane materials is key to engineering and development of desalination and water purification technologies as well as electro-membrane applications. To date, modeling of ion transport has mainly relied on mean-field approaches, originally intended for weak inter-ionic interactions, i.e., high reduced temperature T*. This condition is violated in many membranes, which could explain disagreement between predicted trends and experiments. The paper highlights observed discrepancies and develops a new approach based on the concept of ion association, more adequate in the low-T^⁎ limit. The new model addresses ion binding and mobility consistently within the same physical picture, applied to different types of single and mixed salts. The resulting relations show a significantly weaker connection between ion partitioning and permeability than the standard ones. Estimates using primitive model (PM) of ions in a homogeneous dielectric suggest that non-PM mechanisms, originating from the molecular structure of the ion-solvating environment, might enhance ion association in membranes. PM analysis also predicts that ion solvation and association must be rigidly related, yet non-PM effects may decouple these phenomena and allow a crossover to non-trivial regimes consistent with experiments and simulations. Despite the crude nature of the presented approach and some questions remaining open, it appears to explain most available experimental data and presents a step towards predictive modeling of ion-selective membrane separations in water-, environment- and energy-related applications.

Selectivity and polarization in water channel membranes: lessons learned from polymeric membranes and CNTs

The aspects of ion exclusion and concentration polarization are highlighted as critical for achieving high selectivity in an artificial water channel.

Development of the modern theory of polymeric complex coacervation

Oppositely charged polymers can undergo the process of complex coacervation, which refers to a liquid-liquid phase separation driven by electrostatic attraction. These materials have demonstrated considerable promise as the basis for complex, self-assembled materials. In this review, we provide a broad overview of the theoretical tools used to understand the physical properties of polymeric coacervates. In particular, we discuss historic theories (Voorn-Overbeek, Random Phase Approximation), and then describe recent developments in the field (Field Theoretic, Counterion Release, Molecular Simulation, and Polymer Reference Interaction Site Model methods). We provide context for these methods, and map out the patchwork of theoretical models that are used to describe a diverse array of coacervate systems. We use this review of the literature to clarify a number of important theoretical challenges remaining in our physical understanding of complex coacervation.

Multiply associating electrolytes in the binding mean spherical approximation: thermodynamic properties and speciation

Ionic solutions exhibiting multiple association are described within the binding mean spherical approximation (BiMSA). This model is based on the Wertheim formalism, in the framework of the primitive model at the McMillan-Mayer level. The cation and the anion form the various complexes according to stepwise complexation-equilibria. Analytic expressions for the Helmholtz energy, the internal energy, the speciation, and for the osmotic and activity coefficients are given considering a binary solution with an arbitrary number of association sites on one type of ion (polyion) and one site on the ions of opposite sign (counterions). As an alternative, mean field expressions, as developed in SAFT-type theories, are also presented. The result obtained from the latter approximate method exhibits a reasonable agreement with those from BiMSA for the speciation, and a remarkable one for the osmotic coefficient.

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.

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.

Charge fluctuations and correlation lengths in finite electrolytes

Fluctuations of the charge Q_{Lambda} inside a subdomain Lambda embedded in an electrolyte contained in a finite cubical box of dimensions LxLxL with periodic boundary conditions are investigated. When Lambda is an LxL "slab" of width W , asymptotically exact expressions for the mean-square fluctuation Q_{Lambda};{2} are obtained in terms of the Lebowitz length xi_{L}(T,rho) and of the "true" or asymptotic screening/decay length xi_{Z,infinity}(T,rho) together, when the charge correlation decay is oscillatory, with the characteristic wavelength lambda_{Z}(T,rho) . In finite systems, the normalized charge fluctuations exhibit threefold scaling behavior in the ratios xi_{Z,infinity}W , Wlambda_{Z} , and WL . This enables one to estimate all the correlation lengths away from criticality quite precisely from finite-size grand canonical Monte Carlo simulations. The results for xi_{Z,infinity},lambda_{Z} , and xi_{L} are presented for the restricted primitive model or hard-sphere 1:1 electrolyte for densities rho less, similar1.3rho_{c} and T greater, similar4T_{c} . The fitted values compare favorably with the expectations of generalized Debye-Hückel theory [Lee and Fisher, Europhys. Lett. 39, 611 (1997)]; specifically, if xi_{D} proportional, variant(Trho);{12} is the Debye length, we find xi_{Z,infinity}<xi_{D}<xi_{L} , although, since ion-pairing is neglected, the charge oscillations set in only for densities approximately 1.9 times larger than predicted.

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.

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.