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

Stability of the homogeneous steady state for a model of a confined quasi-two-dimensional granular fluid

A linear stability analysis of the hydrodynamic equations of a model for confined quasi-two-dimensional granular gases is carried out. The stability analysis is performed around the homogeneous steady state (HSS) reached eventually by the system after a transient regime. In contrast to previous studies (which considered dilute or quasielastic systems), our analysis is based on the results obtained from the inelastic Enskog kinetic equation, which takes into account the (nonlinear) dependence of the transport coefficients and the cooling rate on dissipation and applies to moderate densities. As in earlier studies, the analysis shows that the HSS is linearly stable with respect to long enough wavelength excitations.

Assessment of the Wolf method using the Stillinger-Lovett sum rules: From strong electrolytes to weakly charged colloidal dispersions

The Ewald method has been the cornerstone in molecular simulations for modeling electrostatic interactions of charge-stabilized many-body systems. In the late 1990s, Wolf and collaborators developed an alternative route to describe the long-range nature of electrostatic interactions; from a computational perspective, this method provides a more efficient and straightforward way to implement long-range electrostatic interactions than the Ewald method. Despite these advantages, the validity of the Wolf potential to account for the electrostatic contribution in charged fluids remains controversial. To alleviate this situation, in this contribution, we implement the Wolf summation method to both electrolyte solutions and charged colloids with moderate size and charge asymmetries in order to assess the accuracy and validity of the method. To this end, we verify that the proper selection of parameters within the Wolf method leads to results that are in good agreement with those obtained through the standard Ewald method and the theory of integral equations of simple liquids within the so-called hypernetted chain approximation. Furthermore, we show that the results obtained with the original Wolf method do satisfy the moment conditions described by the Stillinger-Lovett sum rules, which are directly related to the local electroneutrality condition and the electrostatic screening in the Debye-Hückel regime. Hence, the fact that the solution provided by the Wolf method satisfies the first and second moments of Stillinger-Lovett proves, for the first time, the reliability of the method to correctly incorporate the electrostatic contribution in charge-stabilized fluids. This makes the Wolf method a powerful alternative compared to more demanding computational approaches.

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.

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.

Control of Selective Ion Transfer across Liquid-Liquid Interfaces: A Rectifying Heterojunction Based on Immiscible Electrolytes

The current rectification displayed by solid-state p-n semiconductor diodes relies on the abundance of electrons and holes near the interface between the p-n junction. In analogy to this electronic device, we propose here the construction of a purely ionic liquid-state electric rectifying heterojunction displaying an excess of monovalent cations and anions near the interface between two immiscible solvents with different dielectric properties. This system does not need any physical membrane or material barrier to show preferential ion transfer but relies on the ionic solvation energy between the two immiscible solvents. We construct a simple device, based on an oil/water interface, displaying an asymmetric behavior of the electric current as a function of the polarity of an applied electric field. This device also exhibits a region of negative differential conductivity, analogous to that observed in brain and heart cells via voltage clamp techniques. Computer simulations and mean field theory calculations for a model of this system show that the application of an external electric field is able to control the bulk concentrations of the ionic species in the immiscible liquids in a manner that is asymmetric with respect to the polarity or direction of the applied electric field. These properties make possible to enhance or suppress selective ion transport at liquid-liquid interfaces with the application of an external electric field or electrostatic potential, mimicking the function of biological ion channels, thus creating opportunities for varied applications.

A simulation assessment of the thermodynamics of dense ion-dipole mixtures with polarization

Molecular dynamics (MD) simulations are employed to ascertain the relative importance of various electrostatic interaction contributions, including induction interactions, to the thermodynamics of dense, hot ion-dipole mixtures. In the absence of polarization, we find that an MD-constrained free energy term accounting for the ion-dipole interactions, combined with well tested ionic and dipolar contributions, yields a simple, fairly accurate free energy form that may be a better option for describing the thermodynamics of such mixtures than the mean spherical approximation (MSA). Polarization contributions induced by the presence of permanent dipoles and ions are found to be additive to a good approximation, simplifying the thermodynamic modeling. We suggest simple free energy corrections that account for these two effects, based in part on standard perturbative treatments and partly on comparisons with MD simulation. Even though the proposed approximations likely need further study, they provide a first quantitative assessment of polarization contributions at high densities and temperatures and may serve as a guide for future modeling efforts.

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.

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.

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.

Attractions between like-charged surfaces with dumbbell-shaped counterions

We study the effect of dumbbell-like counterions on the interactions between similarly charged surfaces. Via a systematic study using Monte Carlo simulations and field theory, we fully consider electrostatic correlations and ion structure and find that their intricate coupling determines the equilibrium phase behaviors. In particular, an energetic bridging mechanism is revealed to cause surface attractions for a finite range of surface separations, even in the Poisson-Boltzmann limit.

Molecular simulation of the swelling of polyelectrolyte gels by monovalent and divalent counterions

Permanently crosslinked polyelectrolyte gels are known to undergo discontinuous first-order volume phase transitions, the onset of which may be caused by a number of factors. In this study we examine the volumetric properties of such polyelectrolyte gels in relation to the progressive substitution of monovalent counterions by divalent counterions as the gels are equilibrated in solvents of different dielectric qualities. We compare the results of coarse-grained molecular dynamics simulations of polyelectrolyte gels with previous experimental measurements by others on polyacrylate gels. The simulations show that under equilibrium conditions there is an approximate cancellation between the electrostatic contribution and the counterion excluded-volume contribution to the osmotic pressure in the gel-solvent system; these two contributions to the osmotic pressure have, respectively, energetic and entropic origins. The finding of such a cancellation between the two contributions to the osmotic pressure of the gel-solvent system is consistent with experimental observations that the swelling behavior of polyelectrolyte gels can be described by equations of state for neutral gels. Based on these results, we show and explain that a modified form of the Flory-Huggins model for nonionic polymer solutions, which accounts for neither electrostatic effects nor counterion excluded-volume effects, fits both experimental and simulated data for polyelectrolyte gels. The Flory-Huggins interaction parameters obtained from regression to the simulation data are characteristic of ideal polymer solutions, whereas the experimentally obtained interaction parameters, particularly that associated with the third virial coefficient, exhibit a significant departure from ideality, leading us to conclude that further enhancements to the simulation model, such as the inclusion of excess salt, the allowance for size asymmetric electrolytes, or the use of a distance-dependent solvent dielectricity model, may be required. Molecular simulations also reveal that the condensation of divalent counterions onto the polyelectrolyte network backbone occurs preferentially over that of monovalent counterions.

Faceting ionic shells into icosahedra via electrostatics

Shells of various viruses and other closed packed structures with spherical topology exhibit icosahedral symmetry because the surface of a sphere cannot be tiled without defects, and icosahedral symmetry yields the most symmetric configuration with the minimum number of defects. Icosahedral symmetry is different from icosahedral-shaped structures, which include some large viruses, cationic-anionic vesicles, and fullerenes. We present a faceting mechanism of ionic shells into icosahedral shapes that breaks icosahedral symmetry resulting from different arrangements of the charged components among the facets. These self-organized ionic structures may favor the formation of flat domains on curved surfaces. We show that icosahedral shapes without rotational symmetry can have lower energy than spheres with icosahedral symmetry caused by preferred bending directions in the planar ionic lattice. The ability to create icosahedral shapes without icosahedral symmetry may lead to the design of new functional materials. The electrostatically driven faceting mechanism we present here suggests that we can design faceted polyhedra with diverse symmetries by coassembling oppositely charged molecules of different stoichiometric ratios.

Disappearance of the gas-liquid phase transition for highly charged colloids

We calculate the full phase diagram of spherical charged colloidal particles using Monte Carlo free energy calculations. The system is described using the primitive model, consisting of explicit hard-sphere colloids and point counterions in a uniform dielectric continuum. We show that the gas-liquid critical point becomes metastable with respect to a gas-solid phase separation at colloid charges Q > or =20 times the counterion charge. Approximate free energy calculations with only one and four particles in the fluid and solid phases, respectively, are used to determine the critical line for highly charged colloids up to Q=2000. We propose the scaling law T*(c) approximately Q(1/2) for this critical temperature.

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

Lattice model of equilibrium polymerization. V. Scattering properties and the width of the critical regime for phase separation

Dynamic clustering associated with self-assembly in many complex fluids can qualitatively alter the shape of phase boundaries and produce large changes in the scale of critical fluctuations that are difficult to comprehend within the existing framework of theories of critical phenomena for nonassociating fluids. In order to elucidate the scattering and critical properties of associating fluids, we consider several models of equilibrium polymerization that describe widely occurring types of associating fluids at equilibrium and that exhibit the well defined cluster geometry of linear polymer chains. Specifically, a Flory-Huggins-type lattice theory is used, in conjunction with the random phase approximation, to compute the correlation length amplitude xi(o) and the Ginzburg number Gi corresponding, respectively, to the scale of composition fluctuations and to a parameter characterizing the temperature range over which Ising critical behavior is exhibited. Our calculations indicate that upon increasing the interparticle association energy, the polymer chains become increasingly long in the vicinity of the critical point, leading naturally to a more asymmetric phase boundary. This increase in the average degree of polymerization implies, in turn, a larger xi(o) and a drastically reduced width of the critical region (as measured by Gi). We thus obtain insight into the common appearance of asymmetric phase boundaries in a wide range of "complex" fluids and into the observation of apparent mean field critical behavior even rather close to the critical point.

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.

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.

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.

Equilibrium polymerization in the Stockmayer fluid as a model of supermolecular self-organization

A diverse range of molecular self-organization processes arises from a competition between directional and isotropic van der Waals intermolecular interactions. We conduct Monte Carlo simulations of the Stockmayer fluid (SF) with a large dipolar interaction as a minimal self-organization model and focus on basic thermodynamic properties that are needed to characterize the polymerization transition that occurs in this fluid. In particular, we determine the polymerization transition lines from the maximum in the specific heat, C(v), and the inflection point in the extent of polymerization, Phi. We also characterize the geometry (radius of gyration R(g), chain length L, chain topology) of the clusters that form in this associating fluid as a function of temperature, T, and concentration, rho . The pressure, P, and the second virial coefficient, B2, were determined, since these properties contain essential information about the strength of the isotropic (van der Waals) interactions. Our simulations indicate that the locations of the polymerization lines are quantitatively consistent with a model of equilibrium polymerization with the enthalpy of polymerization ("sticking energy") fixed by the minimum in the intermolecular potential. The polymerization transition in the SF is accompanied by a topological transition from predominantly linear to ring polymers upon cooling that is driven by the minimization of the dipolar energy of the clusters. We also find that the basic interaction parameters describing polymerization and phase separation in the SF can be estimated based on the existing theory of equilibrium polymerization, but the theory must be refined to account for ring formation in order to accurately describe the configurational properties of this model self-organizing fluid.

Molecular multivalent electrolytes: microstructure and screening lengths

We study small rod-like molecular electrolytes solutions with their corresponding atomic counterions. The asymptotic length scales (decay length and wavelength) of the structural correlations are analyzed using the formalism of the dressed interaction site theory (DIST). The correlation functions are determined using the reference interaction site model equation complemented with a mixed approach in which the hypernetted-chain closure is used for the repulsive interactions, and the mean spherical approximation is used for the attractive interactions. The results from this scheme are in good agreement with the Monte Carlo computer simulations reported here. The asymptotic properties of the correlation functions of this molecular system are compared against those corresponding to two related simple (atomic) electrolyte models. The main conclusion is that the molecular structure of the ions lowers by two orders of magnitude the concentration at which the transition from monotonic to oscillatory decay occurs.

Phase diagram of charged dumbbells: a random phase approximation approach

The phase diagram of the charged hard dumbbell system (hard spheres of opposite unit charge fixed at contact) is obtained with the use of the random phase approximation (RPA). The effect of the impenetrability of charged spheres on charge-charge fluctuations is described by introduction of a modified electrostatic potential. The correlations of ions in a pair are included via a correlation function in the RPA. The coexistence curve is in good agreement with Monte Carlo simulations. The relevance of the theory to the restricted primitive model is discussed.

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.

Flory-Huggins model of equilibrium polymerization and phase separation in the Stockmayer fluid

The competition between chain formation and phase separation in the Stockmayer fluid (SF) of dipolar particles is analyzed using a renormalized Flory-Huggins model of equilibrium polymerization. Calculated critical properties (T(c), phi(c), Z(c)) for the SF compare favorably with simulations over a wide range of the dimensionless dipolar (or "sticking") energy mu*. We find that the polymerization transition preempts phase separation for a large mu*, i.e., (mu*)2 >22.

Monte Carlo simulation of a coarse-grained model of polyelectrolyte networks

The structure and properties of a coarse-grained model of a polyelectrolyte network is studied by means of Monte Carlo simulations. Counterions are treated explicitly, and permanent tetrafunctional cross-linking sites are annealed. The resulting pressure-density relationships exhibit a strong dependence on the strength of electrostatic interactions. A discontinuous volume change is observed when electrostatic interactions are strong. The structure of the model networks is examined at various conditions, and it is found to be considerably different from that of noncross-linked polyelectrolytes.