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

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.

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

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

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.

Computer simulations of the restricted primitive model at very low temperature and density

The problem of successfully simulating ionic fluids at low temperature and low density states is well known in the simulation literature: using conventional methods, the system is not able to equilibrate rapidly due to the presence of strongly associated cation-anion pairs. In this paper we present a numerical method for speeding up computer simulations of the restricted primitive model (RPM) at low temperatures (around the critical temperature) and at very low densities (down to 10(-10)σ(-3), where σ is the ion diameter). Experimentally, this regime corresponds to typical concentrations of electrolytes in nonaqueous solvents. As far as we are aware, this is the first time that the RPM has been equilibrated at such extremely low concentrations. More generally, this method could be used to equilibrate other systems that form aggregates at low concentrations.

Avoiding unphysical kinetic traps in Monte Carlo simulations of strongly attractive particles

We introduce a "virtual-move" Monte Carlo algorithm for systems of pairwise-interacting particles. This algorithm facilitates the simulation of particles possessing attractions of short range and arbitrary strength and geometry, an important realization being self-assembling particles endowed with strong, short-ranged, and angularly specific ("patchy") attractions. Standard Monte Carlo techniques employ sequential updates of particles and can suffer from low acceptance rates when attractions are strong. In this event, collective motion can be strongly suppressed. Our algorithm avoids this problem by proposing simultaneous moves of collections (clusters) of particles according to gradients of interaction energies. One particle first executes a "virtual" trial move. We determine which of its neighbors move in a similar fashion by calculating individual bond energies before and after the proposed move. We iterate this procedure and update simultaneously the positions of all affected particles. Particles move according to an approximation of realistic dynamics without requiring the explicit computation of forces and without the step size restrictions required when integrating equations of motion. We employ a size- and shape-dependent damping of cluster movements, motivated by collective hydrodynamic effects neglected in simple implementations of Brownian dynamics. We discuss the virtual-move algorithm in the context of other Monte Carlo cluster-move schemes and demonstrate its utility by applying it to a model of biological self-assembly.

Vapor-liquid coexistence in fluids of charged hard dumbbells

Vapor-liquid coexistence in fluids of charged hard dumbbells, each made up of two oppositely charged hard spheres with diameters sigma and separation d, has been studied using grand-canonical Monte Carlo simulations. In the limit d/sigma-->0, and with the temperature scaled accordingly, the system corresponds to dipolar hard spheres. For separations in the range 0.3<d/sigma<or=1 the coexisting vapor phase contains compact clusters. For separations in the range 0.1<or=d/sigma<0.3 the coexistence is between a chainlike vapor and a networklike liquid. Finite-size effects preclude the simulation of the coexistence in systems with d/sigma<0.1, but extrapolations of the results to d/sigma-->0 yield estimates of the apparent critical parameters for dipolar hard spheres.

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.

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.

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.

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.

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.

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.

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.