Absence of crystals in the phase behavior of hollow microgels.

Solutions of microgels have been widely used as model systems to gain insight into atomic condensed matter and complex fluids. We explore the thermodynamic phase behavior of hollow microgels, which are distinguished from conventional colloids by a central cavity. Small-angle neutron and x-ray scattering are used to probe hollow microgels in crowded environments. These measurements reveal an interplay among deswelling, interpenetration, and faceting and an unusual absence of crystals. Monte Carlo simulations of model systems confirm that, due to the cavity, solutions of hollow microgels more readily form a supercooled liquid than for microgels with a cross-linked core.

Poisson-Boltzmann theory of charged colloids: limits of the cell model for salty suspensions.

Thermodynamic properties of charge-stabilized colloidal suspensions and polyelectrolyte solutions are commonly modelled by implementing the mean-field Poisson-Boltzmann (PB) theory within a cell model. This approach models a bulk system by a single macroion, together with counterions and salt ions, confined to a symmetrically shaped, electroneutral cell. While easing numerical solution of the nonlinear PB equation, the cell model neglects microion-induced interactions and correlations between macroions, precluding modelling of macroion ordering phenomena. An alternative approach, which avoids the artificial constraints of cell geometry, exploits the mapping of a macroion-microion mixture onto a one-component model of pseudo-macroions governed by effective interparticle interactions. In practice, effective-interaction models are usually based on linear-screening approximations, which can accurately describe strong nonlinear screening only by incorporating an effective (renormalized) macroion charge. Combining charge renormalization and linearized PB theories, in both the cell model and an effective-interaction (cell-free) model, we compute osmotic pressures of highly charged colloids and monovalent microions, in Donnan equilibrium with a salt reservoir, over a range of concentrations. By comparing predictions with primitive model simulation data for salt-free suspensions, and with predictions from nonlinear PB theory for salty suspensions, we chart the limits of both the cell model and linear-screening approximations in modelling bulk thermodynamic properties. Up to moderately strong electrostatic couplings, the cell model proves accurate for predicting osmotic pressures of deionized (counterion-dominated) suspensions. With increasing salt concentration, however, the relative contribution of macroion interactions to the osmotic pressure grows, leading predictions from the cell and effective-interaction models to deviate. No evidence is found for a liquid-vapour phase instability driven by monovalent microions. These results may guide applications of PB theory to colloidal suspensions and other soft materials.

Electroneutrality and phase behavior of colloidal suspensions.

Several statistical mechanical theories predict that colloidal suspensions of highly charged macroions and monovalent microions can exhibit unusual thermodynamic phase behavior when strongly deionized. Density-functional, extended Debye-Hückel, and response theories, within mean-field and linearization approximations, predict a spinodal phase instability of charged colloids below a critical salt concentration. Poisson-Boltzmann cell model studies of suspensions in Donnan equilibrium with a salt reservoir demonstrate that effective interactions and osmotic pressures predicted by such theories can be sensitive to the choice of reference system, e.g., whether the microion density profiles are expanded about the average potential of the suspension or about the reservoir potential. By unifying Poisson-Boltzmann and response theories within a common perturbative framework, it is shown here that the choice of reference system is dictated by the constraint of global electroneutrality. On this basis, bulk suspensions are best modeled by density-dependent effective interactions derived from a closed reference system in which the counterions are confined to the same volume as the macroions. Lower-dimensional systems (e.g., monolayers, clusters), depending on the strength of macroion-counterion correlations, may be governed instead by density-independent effective interactions tied to an open reference system with counterions dispersed throughout the reservoir, possibly explaining the observed structural crossover in colloidal monolayers and anomalous metastability of colloidal crystallites.

Phase separation in charge-stabilized colloidal suspensions: influence of nonlinear screening.

The phase behavior of charge-stabilized colloidal suspensions is modeled by a combination of response theory for electrostatic interparticle interactions and variational theory for free energies. Integrating out degrees of freedom of the microions (counterions, salt ions), the macroion-microion mixture is mapped onto a one-component system governed by effective macroion interactions. Linear response of microions to the electrostatic potential of the macroions results in a screened-Coulomb (Yukawa) effective pair potential and a one-body volume energy, while nonlinear response modifies the effective interactions [A. R. Denton, Phys. Rev. E 70, 031404 (2004)]. The volume energy and effective pair potential are taken as input to a variational free energy, based on thermodynamic perturbation theory. For both linear and first-order nonlinear effective interactions, a coexistence analysis applied to aqueous suspensions of highly charged macroions and monovalent microions yields bulk separation of macroion-rich and macroion-poor phases below a critical salt concentration, in qualitative agreement with predictions of related linearized theories [R. van Roij, M. Dijkstra, and J.-P. Hansen, Phys. Rev. E 59, 2010 (1999); P. B. Warren, J. Chem. Phys. 112, 4683 (2000)]. It is concluded that nonlinear screening can modify phase behavior but does not necessarily suppress bulk phase separation of de-ionized suspensions.

Effective electrostatic interactions in suspensions of polyelectrolyte brush-coated colloids.

Effective electrostatic interactions between colloidal particles, coated with polyelectrolyte brushes and suspended in an electrolyte solvent, are described via linear-response theory. The inner cores of the macroions are modeled as hard spheres, the outer brushes as spherical shells of continuously distributed charge, the microions (counterions and salt ions) as point charges, and the solvent as a dielectric continuum. The multicomponent mixture of macroions and microions is formally mapped onto an equivalent one-component suspension by integrating out from the partition function the microion degrees of freedom. Applying second-order perturbation theory and a random-phase approximation, analytical expressions are derived for the effective pair interaction and a one-body volume energy, which is a natural by-product of the one-component reduction. The combination of an inner core and an outer shell, respectively impenetrable and penetrable to microions, allows the interactions between macroions to be tuned by varying the core diameter and brush thickness. In the limiting cases of vanishing core diameter and vanishing shell thickness, the interactions reduce to those derived previously for star polyelectrolytes and charged colloids, respectively.

Nonlinear screening and effective electrostatic interactions in charge-stabilized colloidal suspensions.

A nonlinear response theory is developed and applied to electrostatic interactions between spherical macroions, screened by surrounding microions, in charge-stabilized colloidal suspensions. The theory describes leading-order nonlinear response of the microions (counterions, salt ions) to the electrostatic potential of the macroions and predicts microion-induced effective many-body interactions between macroions. A linear response approximation [A.R. Denton, Phys. Rev. E 62, 3855 (2000)] yields an effective pair potential of screened-Coulomb (Yukawa) form, as well as a one-body volume energy, which contributes to the free energy. Nonlinear response generates effective many-body interactions and essential corrections to both the effective pair potential and the volume energy. By adopting a random-phase approximation (RPA) for the response functions, and thus neglecting microion correlations, practical expressions are derived for the effective pair and triplet potentials and for the volume energy. Nonlinear screening is found to weaken repulsive pair interactions, induce attractive triplet interactions, and modify the volume energy. Numerical results for monovalent microions are in good agreement with available ab initio simulation data and demonstrate that nonlinear effects grow with increasing macroion charge and concentration and with decreasing salt concentration. In the dilute limit of zero macroion concentration, leading-order nonlinear corrections vanish. Finally, it is shown that nonlinear response theory, when combined with the RPA, is formally equivalent to the mean-field Poisson-Boltzmann theory and that the linear response approximation corresponds, within integral-equation theory, to a linearized hypernetted-chain closure.

Counterion penetration and effective electrostatic interactions in solutions of polyelectrolyte stars and microgels.

Counterion distributions and effective electrostatic interactions between spherical macroions in polyelectrolyte solutions are calculated via second-order perturbation (linear response) theory. By modeling the macroions as continuous charge distributions that are permeable to counterions, analytical expressions are obtained for counterion profiles and effective pair interactions in solutions of star-branched and microgel macroions. The counterions are found to penetrate stars more easily than microgels, with important implications for screening of bare macroion interactions. The effective pair interactions are Yukawa in form for separated macroions, but are softly repulsive and bounded for overlapping macroions. A one-body volume energy, which depends on the average macroion concentration, emerges naturally in the theory and contributes to the total free energy.

Solid-phase structures of the dzugutov pair potential.

In recent computer simulations of a simple monatomic system interacting via the Dzugutov pair potential, freezing of the fluid into an equilibrium dodecagonal quasicrystal has been reported [M. Dzugutov, Phys. Rev. Lett. 70, 2924 (1993)]. Here, using a combination of molecular dynamics simulation and thermodynamic perturbation theory, we conduct a detailed analysis of the relative stabilities of solid-phase structures of the Dzugutov-potential system. At low pressures, the most stable structure is found to be a bcc crystal, which gives way at higher pressures to a fcc crystal. Although a dodecagonal quasicrystal and a sigma-phase crystal compete with the bcc crystal for stability, they always remain metastable.