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

Structure of colloids and macroionic solutions with soft,long-range interactions

The structure of strongly interacting, charged dispersions is examined using the Ornstein-Zernike integral-equation formalism with thermodynamic consistency conditions (i.e., with the so-called HMSA closure conditions). It is shown that for highly charged dispersions interacting through screened Coulombic (Yukawa) interactions this integral-equation approach predicts the static structure factorS(q) in excellent agreement with Monte Carlo results and that it is better than the rescaled mean-spherical approximation. The Sogaimi potential (which predicts Coulombic attraction under identical physicochemical conditions) is also considered here as a model potential for representing soft, long-ranged interactions with an attractive component. None of the integralequation formalisms, including the HMSA theory and the reference hypernetted-chain (RHNC) theory, leads to a sufficient level of thermodynamic consistency for the latter potential, and the results deviate noticeably from the Monte Carlo results. We further demonstrate that fitting the experimentally observed structure factors in the neighborhood and beyond the primary peak inS(q) could lead to inaccurate conclusions concerning the nature of interparticle forces, particularly in the case of soft, long-range interactions.     

Interaction of similarly charged surfaces mediated by nanoparticles

We study the interaction between two like charged surfaces embedded in a solution of oppositely charged multivalent rod-like counterions.The counterions consist of two rigidly bonded point charges,each of valency Z.The strength of the electrostatic coupling increases with increasing surface charge density or valency of the charges.The system is analyzed by employing a self-consistent field theory,which treats the short and long range interactions of the counterions within different approximations.We find that in the weak coupling limit,the interactions are only repulsive.In the intermediate coupling regime,the multivalent rod-like counterions can mediate attractive interactions between the surfaces. For sufficiently long rods,bridging contributes to the attractive interaction.In the strong coupling limit,the charge correlations can contribute to the attractive interactions at short separations between the charged surfaces.Two minima can then appear in the force curve between surfaces.     

Effective electrostatic interactions in solutions of polyelectrolyte stars with rigid rodlike arms

In solutions of star-branched polyelectrolytes, electrostatic interactions between charged arms on neighboring stars can compete with intrastar interactions and rotational entropy to induce anisotropy in the orientational distribution of arms. We explore the influence of arm orientational anisotropy on effective star-star interactions for model stars comprising rigid rodlike arms with evenly spaced charged monomers interacting via an effective screened-Coulomb (Yukawa) potential. Monte Carlo simulation and density-functional theory are used to compute the arm orientational distributions and effective pair potentials between weakly charged stars. For comparison, a torque balance analysis is performed to obtain the configuration and energy of the ground state, in which the torque vanishes on each arm of the two-star system. The degree of anisotropy is found to increase with the strength of electrostatic interactions and proximity of the stars. As two stars begin to overlap, the forward arms are pushed back by interstar arm-arm repulsion, but partially interdigitate due to rotational entropy. At center-center separations approaching complete overlap, the arms relax to an isotropic distribution. For nonoverlapping stars, anisotropy-induced changes in the intra- and interstar arm-arm interactions largely cancel and the effective pair interactions are then well approximated by a simple Yukawa potential, as predicted by linear-response theory for a continuum model of isotropic stars [A. R. Denton, Phys. Rev. E 67, 11804 (2003)]. For overlapping stars, the effective pair interactions in the simple rigid-arm-Yukawa model agree closely with simulations of a molecular model that includes flexible arms and explicit counterions [A. Jusufi et al., Phys. Rev. Lett. 88, 018301 (2002); J. Chem. Phys. 116, 11011 (2002)].     

Macroion solutions in the cell model studied by field theory and Monte Carlo simulations

Aqueous solutions of charged spherical macroions with variable dielectric permittivity and their associated counterions are examined within the cell model using a field theory and Monte Carlo simulations. The field theory is based on separation of fields into short- and long-wavelength terms, which are subjected to different statistical-mechanical treatments. The simulations were performed by using a new, accurate, and fast algorithm for numerical evaluation of the electrostatic polarization interaction. The field theory provides counterion distributions outside a macroion in good agreement with the simulation results over the full range from weak to strong electrostatic coupling. A low-dielectric macroion leads to a displacement of the counterions away from the macroion.     

Phase separation in mixtures of Yukawa and charged Yukawa particles from Gibbs ensemble Monte Carlo simulations and the mean spherical approximation

The phase equilibrium of mixtures of Yukawa and charged Yukawa particles is studied by means of Gibbs ensemble Monte Carlo (GEMC) simulation method and the mean spherical approximation (MSA). The strength of the Coulomb energy compared to that of the Yukawa attraction is characterized by a coupling constant. For low coupling constants a classical vapor--liquid phase separation appears with a good agreement between GEMC and the MSA. For high coupling constant, a phase separation between a salt poor and a salt rich phase occurs that resembles the phase equilibrium behavior of the solvent primitive model.     

Charge renormalization of charged spheres based on thermodynamic properties

At strong electrostatic coupling, counterions are accumulated in the vicinity of the surface of the charged particle with intrinsic charge Z. In order to explain the behavior of highly charged particles, effective charge Z(*) is therefore invoked in the models based on Debye-Huckel approximation, such as the Derjaguin-Landau-Verwey-Overbeek potential. For a salt-free colloidal suspension, we perform Monte Carlo simulations to obtain various thermodynamic properties omega in a spherical Wigner-Seitz cell. The effect of dielectric discontinuity is examined. We show that at the same particle volume fraction, counterions around a highly charged sphere with Z may display the same value of omega as those around a weakly charged sphere with Z(*), i.e., omega(Z)=omega(Z(*)). There exists a maximally attainable value of omega at which Z=Z(*). Defining Z(*) as the effective charge, we find that the effective charge passes through a maximum and declines again due to ion-ion correlation as the number of counterions is increased. The effective charge is even smaller if one adopts the Debye-Huckel expression omega(DH). Our results suggest that charge renormalization can be performed by chemical potential, which may be observed in osmotic pressure measurements.     

Use of the Yukawa fluid as the reference system in perturbation theory

The Weeks-Chandler-Anderson (WCA) perturbation theory is studied utilising recent results for the Yukawa model fluid. Replacing the attractive tail of the Lennard-Jones potential with a Yukawa tail, where the Yukawa parameters are chosen using a least squares fit, it is shown that accurates field dstribution functions can be generated via the EXP approximation of the WCA optimized cluster theory. The comparative case and accuracy with which the correlation functions for the Yukawa fluid can be compared render this a very useful method for studying the equilibrium properties of simple liquids.     

Thermodynamic and structural properties of repulsive hard-core Yukawa fluid: integral equation theory, perturbation theory and Monte Carlo simulations

The thermodynamic and structural properties of purely repulsive hard-core Yukawa particles in the fluid state are determined through Monte Carlo simulation and modeled using perturbation theory and integral equation theory in the mean spherical approximation (MSA). Systems of particles with Yukawa screening lengths of 1.8, 3.0, and 5.0 are examined with results compared to variations of MSA and perturbation theory. Thermodynamic properties were predicted well by both theories in the fluid region up to the fluid-solid phase boundary. Further, we found that a simplified exponential version of the MSA is the most accurate at predicting radial distribution function at contact. Radial distribution function of repulsive hard-core Yukawa particles are also reported. The results show that methods based on MSA and perturbation theory that are typically applied to the attractive hard-core Yukawa potential can also be extended to the purely repulsive hard-core Yukawa potential.     

Synergistic effects in mixtures of oppositely charged surfactants as calculated from the Poisson-Boltzmann theory: a comparison between theoretical predictions and experiments

Critical micelle concentrations in mixtures of an anionic surfactant and a cationic amphiphilic drug have been investigated using a model-independent procedure to quantify observed synergistic effects. Experimental results were compared with a theory based on the Poisson-Boltzmann mean field approximation of a charged interface with a diffuse layer of counterions. Explicit expressions for the activity coefficients from which the critical micelle concentration can be calculated and quantitatively predicted have been derived and excellent agreement between experimental data and theory was obtained. As a result, we demonstrate that it is possible to rationalize and predict the magnitude of synergism in mixtures of oppositely charged surfactants in the presence of added salt.     

Effective multipoles and Yukawa electrostatics in dressed molecule theory

In this paper we derive the multipolar expansion of the screened Coulomb potential in electrolyte solutions with molecular solvent. The solute and solvent molecules can have arbitrary sizes, shapes, and internal charge distributions. We use the exact statistical mechanical definition of renormalized charge distributions coming from "dressed molecule theory" to determine the effective multipoles of a molecule immersed in an electrolyte. The effects of many-body correlations are fully included in our formally exact theory. We restrict ourselves to sufficiently dilute solutions so the screened Coulomb potential decays for large distances like a Yukawa function, exp(-kappa r)/r, where r is the distance and 1/kappa is the decay length (it is normally different from the Debye length). The resulting "Yukawa electrostatics" differ in many respects from ordinary, unscreened electrostatics. The "Yukawa charge" of a molecule (the lowest order moment in the multipolar expansion) is in general not equal to its Coulombic charge and it is not the integral of the renormalized charge distribution of the molecule. Moreover, as shown in this paper, the multipolar expansion of the Yukawa potential does not correspond, contrary to the case of the Coulomb potential, to its asymptotic expansion for large r. As a consequence, the charge term in the multipolar expansion is not the leading term in the asymptotic expansion. Instead, for large r values, multipoles of all orders contribute to the leading asymptotic term. Thus, the electrostatic potential from, for example, an electroneutral solvent molecule in an electrolyte solution has generally the same range as that from an ion. The proper asymptotic expansion for electrostatic interactions in electrolytes is derived. It is briefly shown how the multipole expansion formalism can also be applied in the Poisson-Boltzmann approximation for primitive model electrolytes.     

Counterion condensation on charged spheres, cylinders, and planes

We use the framework of counterion condensation theory, in which deviations from linear electrostatics are ascribed to charge renormalization caused by collapse of counterions from the ion atmosphere, to explore the possibility of condensation on charged spheres, cylinders, and planes immersed in dilute solutions of simple salt. In the limit of zero concentration of salt, we obtain Zimm-Le Bret behavior: a sphere condenses none of its counterions regardless of surface charge density, a cylinder with charge density above a threshold value condenses a fraction of its counterions, and a plane of any charge density condenses all of its counterions. The response in dilute but nonzero salt concentrations is different. Spheres, cylinders, and planes all exhibit critical surface charge densities separating a regime of counterion condensation from states with no condensed counterions. The critical charge densities depend on salt concentration, except for the case of a thin cylinder, which exhibits the invariant criticality familiar from polyelectrolyte theory.     

Recent study on counterion-mediated attraction between colloidal particles

Recent experimental results were reviewed. The 1D- and 2D-USAXS studies gave higher orders of Bragg diffraction for single crystals of colloidal silica particles, allowing more accurate determinations of the lattice constant, lattice symmetry, and direction. The closest interparticle spacing thus determined was confirmed to be smaller than the average spacing. The most closely packed planes ((110) planes for bcc) of negatively charged particles were found to be parallel to the likewise negatively charged capillary surface, inconsistently with the accepted double layer interaction theory but consistently with a recent experimental finding of positive adsorption. Shaking caused disruption of the single crystals but newly formed microcrystals retained the lattice constant and the preference of the (110) planes. The liquid-solid-liquid transition, a re-entrant phase transition, was found for silica particles and latex particles at given particle volume fraction and salt concentration, when the charge density of particles was varied. It was demonstrated that the purely repulsive Yukawa potential and the concept of renormalized charge cannot account for the re-entrant behavior. The Monte-Carlo simulation using the Sogami potential, which contains short-range repulsion and long-range attraction, was found to account for the fcc–bcc transition, which was earlier claimed to be explainable only by the Yukawa potential. Furthermore, the homogeneous-inhomogeneous phase transition and void formation could be accounted for by the simulation using the Sogami potential; the Yukawa potential could not reproduce the experiments. Attention was drawn to the experimental conditions in direct measurements of interparticle forces; only short interparticle distance and low charge density particles were covered, which make it practically impossible to detect the long-range counterion-mediated attraction. It is hoped that, by technical improvements, these shortcomings may be made up and quantitative argument become possible on the attraction.     

Theory of inhomogeneous weakly charged polyelectrolytes

A theory of inhomogeneous multicomponent systems containing weakly charged polyelectrolytes is developed. The theory treats the polymer conformation and the electrostatics simultaneously using a functional integral representation of the partition function. A mean‐field approximation to the theory leads to two sets of coupled mean‐field equations: a Poisson‐Boltzmann type equation describing the electrostatic potential, and a set of self‐consistent field equations describing the equilibrium densities. Asymptotic forms of the theory at weak and strong segregation limits are derived. The theory can be used to study the interfacial properties, microphase structures, and adsorptions of a variety of weakly charged polyelectrolyte systems. As a simple example, the interface between the polymer‐rich and polymer‐poor phases of a polyelectrolyte solution is studied.     

Molecular dynamics study of charged dendrimers in salt-free solution: effect of counterions

Polyamidoamine dendrimers, being protonated under physiological conditions, represent a promising class of nonviral, nanosized vectors for drug and gene delivery. We performed extensive molecular dynamics simulations of a generic model dendrimer in a salt-free solution with dendrimer's terminal beads positively charged. Solvent molecules as well as counterions were explicitly included as interacting beads. We find that the size of the charged dendrimer depends nonmonotonically on the strength of electrostatic interactions demonstrating a maximum when the Bjerrum length equals the diameter of a bead. Many other structural and dynamic characteristics of charged dendrimers are also found to follow this pattern. We address such a behavior to the interplay between repulsive interactions of the charged terminal beads and their attractive interactions with oppositely charged counterions. The former favors swelling at small Bjerrum lengths and the latter promotes counterion condensation. Thus, counterions can have a dramatic effect on the structure and dynamics of charged dendrimers and, under certain conditions, cannot be treated implicitly.     

The density distributions of the counterions and the coions confined in two similarly charged plates

By using the field-theoretic method, we established a unified systematic formulation of a model of counterions and coions confined in two similarly charged plates, and calculated the density distributions of counterions and coions with various coupling parameters by the two methods: Poisson-Boltzmann (PB) approach and the strong coupling (SC) theory, respectively. We also performed Monte Carlo simulations, and obtained the density distributions of counterions and coions with several different coupling parameters. Comparing our theoretical results with those from Monte Carlo simulation, we find that the PB approach is valid when the coupling parameter Xi is smaller than 1, but, as Xi > or = 1, the results by the PB approach deviate from the corresponding Monte Carlo simulation data, and the deviation gets larger with the coupling parameter increasing. This shows that the PB approach is completely invalid when the coupling parameter is equal to 1 or larger than 1. For the latter case, the development trend of the distribution curve calculated by SC theory agrees with that from Monte Carlo simulation as the coupling parameter increases. This demonstrates that the SC theory can give a qualitative available explanation on the density distribution of the counterions in the system in which the coupling parameters are strictly confined.     

Improved radial distribution functions for Coulomb charged fluid based on first-order mean spherical approximation

In this work, we solve the Ornstein-Zernike equation in a simple, analytical, and consistent manner to obtain the like and unlike radial distribution functions (RDFs) for charged fluids. To improve mean spherical approximation (MSA) solutions, the direct correlation functions both for the density and charge contributions are modified with the Yukawa potential, respectively. On the basis of the contact values of RDFs and excess internal energy of the system, we construct correlated equations to cope with the potential parameters. Thus obtained equations are solved with the first-order MSA method. The resulting like and unlike RDFs are in good agreement with molecular simulation data within a wide range of densities and temperatures.     

Entropic effects in the electric double layer of model colloids with size-asymmetric monovalent ions

The structure of the electric double layer of charged nanoparticles and colloids in monovalent salts is crucial to determine their thermodynamics, solubility, and polyion adsorption. In this work, we explore the double layer structure and the possibility of charge reversal in relation to the size of both counterions and coions. We examine systems with various size-ratios between counterions and coions (ion size asymmetries) as well as different total ion volume fractions. Using Monte Carlo simulations and integral equations of a primitive-model electric double layer, we determine the highest charge neutralization and electrostatic screening near the electrified surface. Specifically, for two binary monovalent electrolytes with the same counterion properties but differing only in the coion's size surrounding a charged nanoparticle, the one with largest coion size is found to have the largest charge neutralization and screening. That is, in size-asymmetric double layers with a given counterion's size the excluded volume of the coions dictates the adsorption of the ionic charge close to the colloidal surface for monovalent salts. Furthermore, we demonstrate that charge reversal can occur at low surface charge densities, given a large enough total ion concentration, for systems of monovalent salts in a wide range of ion size asymmetries. In addition, we find a non-monotonic behavior for the corresponding maximum charge reversal, as a function of the colloidal bare charge. We also find that the reversal effect disappears for binary salts with large-size counterions and small-size coions at high surface charge densities. Lastly, we observe a good agreement between results from both Monte Carlo simulations and the integral equation theory across different colloidal charge densities and 1:1-electrolytes with different ion sizes.