Many-body interactions and the melting of colloidal crystals

Electrostatic interactions between charge regulated spherical macroions

We study the interaction between two charge regulating spherical macroions with dielectric interior and dissociable surface groups immersed in a monovalent electrolyte solution. The charge dissociation is modelled via the Frumkin-Fowler-Guggenheim isotherm, which allows for multiple adsorption equilibrium states. The interactions are derived from the solutions of the mean-field Poisson-Boltzmann type theory with charge regulation boundary conditions. For a range of conditions we find symmetry breaking transitions from symmetric to asymmetric charge distribution exhibiting annealed charge patchiness, which results in like-charge attraction even in a univalent electrolyte-thus fundamentally modifying the nature of electrostatic interactions in charge-stabilized colloidal suspensions.

Jellium and cell model for titratable colloids with continuous size distribution

A good understanding and determination of colloidal interactions is paramount to comprehend and model the thermodynamic and structural properties of colloidal suspensions. In concentrated aqueous suspensions of colloids with a titratable surface charge, this determination is, however, complicated by the density dependence of the effective pair potential due to both the many-body interactions and the charge regulation of the colloids. In addition, colloids generally present a size distribution which results in a virtually infinite combination of colloid pairs. In this paper, we develop two methods and describe the corresponding algorithms to solve this problem for arbitrary size distributions. An implementation in Nim is also provided. The methods, inspired by the seminal work of Torres et al., [J. Chem. Phys. 128, 154906 (2008)] are based on a generalization of the cell and renormalized jellium models to polydisperse suspensions of spherical colloids with a charge regulating boundary condition. The latter is described by the one-pK-Stern model. The predictions of the models are confronted to the equations of state of various commercially available silica dispersions. The renormalized Yukawa parameters (effective charges and screening lengths) are also calculated. The importance of size and charge polydispersity as well as the validity of these two models is discussed in light of the results.

Decoupling conservative forces and hydrodynamic interactions between optically trapped spheres

Characterizing the interactions between colloidal particles is important, both from a fundamental perspective as well as due to its technological importance. However, current methods to measure the interaction forces between two colloids have significant limitations. Here we describe a method that exploits the fluctuation spectra of two optically trapped microspheres in order to extract, and decouple, the conservative forces acting between them and their hydrodynamic coupling. We demonstrate the proposed method with two silica microspheres, and find good agreement between our results and previous predictions for the hydrodynamic and electrostatic interactions between the spheres.

Renormalized jellium model for colloidal mixtures

In an attempt to quantify the role of polydispersity in colloidal suspensions, we present an efficient implementation of the renormalized jellium model for a mixture of spherical charged colloids. The different species may have different size, charge, and density. Advantage is taken from the fact that the electric potential pertaining to a given species obeys a Poisson's equation that is species independent; only boundary conditions do change from one species to the next. All species are coupled through the renormalized background (jellium) density, that is determined self-consistently. The corresponding predictions are compared to the results of Monte Carlo simulations of binary mixtures, where Coulombic interactions are accounted for exactly, at the primitive model level (structureless solvent with fixed dielectric permittivity). An excellent agreement is found.

Ion density deviations in semipermeable ionic microcapsules

By implementing the nonlinear Poisson-Boltzmann theory in a cell model, we theoretically investigate the influence of polyelectrolye gel permeability on ion densities and pH deviations inside the cavities of ionic microcapsules. Our calculations show that variations in permeability of a charged capsule shell cause a redistribution of ion densities within the capsule, which ultimately affects the pH deviation and Donnan potential induced by the electric field of the shell. We find that semipermeable capsules can induce larger pH deviations inside their cavities that can permeable capsules. Furthermore, with increasing capsule charge, the influence of permeability on pH deviations progressively increases. Our theory, while providing a self-consistent method for modeling the influence of permeability on fundamental properties of ionic microgels, makes predictions of practical significance for the design of microcapsules loaded with fluorescent dyes, which can serve as biosensors for diagnostic purposes.

Two-Dimensional Nanoparticle Supracrystals: A Model System for Two-Dimensional Melting

In a Langmuir trough, successive compression cycles can drive a two-dimensional (2D) nanoparticle supracrystal (NPSC) closer to its equilibrium structure. Here, we show a series of equilibrated 2D NPSCs consisting of gold NPs of uniform size, varying solely in the length of their alkanethiol ligands. The ordering of the NPSC is governed by the ligand length, thus providing a model system to investigate the nature of 2D melting in a system of NPs. As the ligand length increases the supracrystal transitions from a crystalline to a liquid-like phase with evidence of a hexatic phase at an intermediate ligand length. The phase change is interpreted as an entropy-driven phenomenon associated with steric constraints between ligand shells. The density of topological defects scales with ligand length, suggesting an equivalence between ligand length and temperature in terms of melting behavior. On the basis of this equivalence, the experimental evidence indicates a two-stage 2D melting of NPSCs.

Effective charges and virial pressure of concentrated macroion solutions

The stability of colloidal suspensions is crucial in a wide variety of processes, including the fabrication of photonic materials and scaffolds for biological assemblies. The ionic strength of the electrolyte that suspends charged colloids is widely used to control the physical properties of colloidal suspensions. The extensively used two-body Derjaguin-Landau-Verwey-Overbeek (DLVO) approach allows for a quantitative analysis of the effective electrostatic forces between colloidal particles. DLVO relates the ionic double layers, which enclose the particles, to their effective electrostatic repulsion. Nevertheless, the double layer is distorted at high macroion volume fractions. Therefore, DLVO cannot describe the many-body effects that arise in concentrated suspensions. We show that this problem can be largely resolved by identifying effective point charges for the macroions using cell theory. This extrapolated point charge (EPC) method assigns effective point charges in a consistent way, taking into account the excluded volume of highly charged macroions at any concentration, and thereby naturally accounting for high volume fractions in both salt-free and added-salt conditions. We provide an analytical expression for the effective pair potential and validate the EPC method by comparing molecular dynamics simulations of macroions and monovalent microions that interact via Coulombic potentials to simulations of macroions interacting via the derived EPC effective potential. The simulations reproduce the macroion-macroion spatial correlation and the virial pressure obtained with the EPC model. Our findings provide a route to relate the physical properties such as pressure in systems of screened Coulomb particles to experimental measurements.

Quantitative assessment of the accuracy of the Poisson-Boltzmann cell model for salty suspensions

The cell model is a ubiquitous, fast, and relatively easily implemented model used to estimate the osmotic pressure of a colloidal dispersion. It has been shown to yield accurate approximations of the pressure in dispersions with a low salt content. It is generally accepted that it performs well when long-ranged interactions are involved and the structure of the dispersion is solidlike. The aim of the present work is to determine quantitatively the error committed by assuming the pressure computed with the cell model is the real osmotic pressure of a dispersion. To this end, cell model pressures are compared to a correct estimation of the actual pressures obtained from Poisson-Boltzmann Brownian dynamics simulations including many-body electrostatics and the thermal motion of the colloids. The comparison is performed for various colloidal sizes and charges, salt contents, and volume fractions. It is demonstrated that the accuracy of the cell model predictions is a function of only the average intercolloid distance scaled by Debye's length κd̅ and the normalized colloidal charge. The cell model is accurate for κd̅ < 1 and not reliable for κd̅ > 5 independently of the colloidal charge. In the 1 < κd̅ < 5 range, covering a wide set of experimental conditions, the colloidal surface charge has a large influence on the error associated with the cell approximation. The results presented in this article should provide a useful reference to determine a priori if the cell model can be expected to predict accurately an equation of state for a given set of physicochemical parameters.

Pattern formation and coarse-graining in two-dimensional colloids driven by multiaxial magnetic fields

We study the pattern formation in a two-dimensional system of superparamagnetic colloids interacting via spatially coherent induced interactions driven by an external precessing magnetic field. On the pair level, upon changing the opening angle of the external field, the interactions smoothly vary from purely repulsive (opening angle equal to zero) to purely attractive (time-averaged pair interactions at an opening angle of 90°). In the experiments, we observed ordered hexagonal crystals at the repulsive end and coarsening frothlike structures for purely attractive interactions. In both of these limiting cases, the dense colloidal systems can be sufficiently accurately described by assuming pairwise additivity of the interaction potentials. However, for a range of intermediate angles, pronounced many-body depolarization effects compete with the direct induced interactions, resulting in inherently anisotropic effective interactions. Under such conditions, we observed the decay of hexagonal order with the concomitant formation of short chains and percolated networks of chains coexisting with free colloids. In order to describe and investigate these systems theoretically, we developed a coarse-grained model of a binary mixture of patchy and nonpatchy particles with the ratio of patchy and nonpatchy colloids as the order parameter. Combining genetic algorithms with Monte Carlo simulations, we optimized the model parameters and quantitatively reproduced the experimentally observed sequence of colloidal structures. The results offer new insight into the anisotropy induced by the many-body effects. At the same time, they allow for a very efficient description of the system by means of a pairwise-additive Hamiltonian, whereupon the original, one-component system features a two-component mixture of isotropic and patchy colloids.

Poisson-Boltzmann description of interaction forces and aggregation rates involving charged colloidal particles in asymmetric electrolytes

Forces and aggregation rates involving spherical particles are studied numerically within the theory of Derjaguin, Landau, Verwey, and Overbeek (DLVO) for asymmetric and mixed electrolytes. Thereby, the double layer interactions are treated at the Debye-Hückel (DH) and Poisson-Boltzmann (PB) levels. The DH model is applicable for weakly charged systems, and effects of ion valence enter only implicitly through the ionic strength. The PB model is necessary for more highly charged systems, and depends on the actual ionic composition. One finds that forces in asymmetric electrolytes at fixed ionic strength weaken when the valence of the counterions is increased or when the valence of the coions is decreased. In symmetric electrolytes, the effect of counterions is more important than the one of the coions. For weakly charged systems, the critical coagulation concentration (CCC) decreases with the square of the valence in symmetric electrolytes, while this decrease is weaker in asymmetric ones. With increasing charge density, the dependence of the CCC on the valence becomes stronger, but the classical Schulze-Hardy decrease with the sixths power of the valence is only recovered for unrealistically high charge densities. Mixtures of electrolytes are treated within the same framework, and one observes that already small amounts of multivalent ions affect the system considerably. An empirical mixing rule is proposed to describe the calculated CCCs.

Colloidal ionic complexes on periodic substrates: ground-state configurations and pattern switching

We theoretically and numerically studied ordering of "colloidal ionic clusters" on periodic substrate potentials such as those generated by optical trapping. Each cluster consists of three charged spherical colloids: two negatively and one positively charged. The substrate is a square or rectangular array of traps, each confining one such cluster. By varying the lattice constant from large to small, the observed clusters are first rodlike and form ferro- and antiferrolike phases, then they bend into a bananalike shape, and finally they condense into a percolated structure. Remarkably, in a broad parameter range between single-cluster and percolated structures, we have found stable supercomplexes composed of six colloids forming grapelike or rocketlike structures. We investigated the possibility of macroscopic pattern switching by applying external electrical fields.

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.

Stable centred tetragonal phases in the hard core Yukawa system

The hard core Yukawa system, consisting of particles with an impenetrable core and interacting via a repulsive Yukawa pair potential, is known to solidify either in a bcc or fcc crystal. Using optimization strategies based on genetic algorithms we give evidence that at zero temperature an additional centred tetragonal phase structure in the range of high packing fractions emerges.

Many-body electrostatic forces between colloidal particles at vanishing ionic strength

Electrostatic forces between small groups of colloidal particles are measured using blinking optical tweezers. When the electrostatic screening length is longer than the interparticle separation, forces are found to be non-pairwise-additive. Both pair and multiparticle forces are well described by the linearized Poisson-Boltzmann equation with constant potential boundary conditions. These findings may play an important role in understanding the structure and stability of a wide variety of systems, from micron-sized particles in oil to aqueous nanocolloids.

Ion-mediated interactions in suspensions of oppositely charged nanoparticles

The structure of oppositely charged spherical nanoparticles (polyions), dispersed in ionic solutions with continuous solvent (primitive model), is investigated by Monte Carlo (MC) simulations, within explicit and implicit microion representations, over a range of polyion valences and densities, and microion concentrations. Systems with explicit microions are explored by semigrand canonical MC simulations, and allow density-dependent effective polyion pair potentials v(alphabeta) (eff)(r) to be extracted from measured partial pair distribution functions. Implicit microion MC simulations are based on pair potentials of mean force v(alphabeta) ((2))(r) computed by explicit microion simulations of two charged polyions, in the low density limit. In the vicinity of the liquid-gas separation expected for oppositely charged polyions, the implicit microion representation leads to an instability against density fluctuations for polyion valences mid R:Zmid R: significantly below those at which the instability sets in within the exact explicit microion representation. Far from this instability region, the v(alphabeta) ((2))(r) are found to be fairly close to but consistently more repulsive than the effective pair potentials v(alphabeta) (eff)(r). This is corroborated by additional calculations of three-body forces between polyion triplets, which are repulsive when one polyion is of opposite charge to the other two. The explicit microion MC data were exploited to determine the ratio of salt concentrations c and c(o) within the dispersion and the reservoir (Donnan effect). c/c(o) is found to first increase before finally decreasing as a function of the polyion packing fraction.

The renormalized jellium model for spherical and cylindrical colloids

Starting from a mean-field description for a dispersion of highly charged spherical or (parallel) rodlike colloids, we introduce the simplification of a homogeneous background to include the contribution of other polyions to the static field created by a tagged polyion. The charge of this background is self-consistently renormalized to coincide with the polyion effective charge, the latter quantity thereby exhibiting a nontrivial density dependence, which directly enters into the equation of state through a simple analytical expression. The good agreement observed between the pressures calculated using the renormalized jellium and Monte Carlo simulations confirms the relevance of the renormalized jellium model for theoretical and experimental purposes and provides an alternative to the Poisson-Boltzmann cell model since it is free of some of the intrinsic limitations of this approach.

Nonlinear screening and gas-liquid separation in suspensions of charged colloids

We calculate phase diagrams of charged colloidal spheres (valency Z and radius a) in a 1:1 electrolyte from multicentered nonlinear Poisson-Boltzmann theory. Our theory takes into account charge renormalization of the colloidal interactions and volume terms due to many-body effects. For valencies as small as Z = 1 and as large as 10(4) we find a gas-liquid spinodal instability in the colloid-salt phase diagram provided Z lambdaB/a > or similar 24+/-1, where lambdaB is the Bjerrum length.

Effective charges along the melting line of colloidal crystals

The shear modulus G of charged colloidal crystals was measured at several constant particle densities n and varying salt concentrations c up to the melting salt concentration cM using torsional resonance spectroscopy. Far from the phase boundary the samples are polycrystalline and the shear modulus stays roughly constant as a function of c. Upon approaching the melting transition an increasing amount of wall based crystal material is formed surrounding a shrinking polycrystalline core and G drops nearly linearly. When the transition is complete G again stays constant. The morphologic transitions may be scaled upon a single master curve. For the polycrystalline morphology, the elastic data are evaluated in terms of a pairwise additive screened Coulomb interaction yielding a particle effective charge Z(G)*. Under de-ionized conditions Z(0,G)* is independent of n and significantly lower than expected from charge renormalization theory. With increasing salt concentration Z(G)* increases. The increase becomes more pronounced at larger n. By extrapolation we further obtain the melting line effective elasticity charge Z(M,G)*. Z(M,G)* shows a steplike increase with increasing nM and cM to values consistent with charge renormalization theory. Interestingly, the increase coincides semi-quantitatively with the one expected from the universal melting line for charged spheres, thus facilitating a consistent description of phase behavior and elasticity over an extended range of the phase diagram.

Effective triplet interactions in nematic colloids

Three-body effective interactions emerging between parallel cylindrical rods immersed in a nematic liquid crystals are calculated within the Landau-de Gennes free-energy description. Collinear, equilateral and midplane configurations of the three colloidal particles are considered. In the last two cases the effective triplet interaction is of the same magnitude and range as the pair one.

Direct numerical simulations of electrophoresis of charged colloids

We propose a numerical method to simulate electrohydrodynamic phenomena in charged colloidal dispersions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier-Stokes equations so that the electrohydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters; however, our simulation predicts notable deviations in the case of dense dispersions.

Volume terms for charged colloids: a grand-canonical treatment

We present a study of thermodynamic properties of suspensions of charged colloids on the basis of linear Poisson-Boltzmann theory. We calculate the effective Hamiltonian of the colloids by integrating out the ionic degrees of freedom grand canonically. This procedure not only yields the well-known pairwise screened-Coulomb interaction between the colloids, but also additional volume terms that affect the phase behavior and the thermodynamic properties, such as the osmotic pressure. These calculations are greatly facilitated by the grand-canonical character of our treatment of the ions and allow for relatively fast computations compared to earlier studies in the canonical ensemble. Moreover, the present derivation of the volume terms are relatively simple, make a direct connection with Donnan equilibrium, yield an explicit expression for the effective screening constant, and allow for extensions to include, for instance, nonlinear effects.

Self-diffusion in a liquid complex plasma

Self-diffusion has been experimentally studied in a two-dimensional underdamped liquid complex (dusty) plasma. It was found that the self-diffusion coefficient D increases linearly with the temperature T: D/omega(E)a2 = (0.019 +/- 0.007)(T/T(m) - 1), where T(m), omega(E), and a are the melting temperature, the Einstein frequency, and the mean particle separation, respectively. No superdiffusion was observed, whereas a subdiffusion occurred at temperatures close to melting.

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

Effect of counterion species on colloidal crystal

The effect of counterion species on the colloidal crystal structure in a dispersion was carefully investigated as a function of the degree of neutralization (alpha) by the ultra-small-angle X-ray scattering technique. The nearest neighbor interparticle distance (2D(exp)) first increased with decreasing alpha, and then decreased after passing through the maximum. This behavior was confirmed for K(+), Li(+), Ca(2+), TMA(+) (tetramethylammonium) as a counterion, and Na(+) in our previous report (Harada, T.; Matsuoka, H.; Ikeda, T.; Yamaoka, H. Langmuir 2000, 16, 1612). However, the alpha value of the maximum position (alpha(max)) largely depended on the counterion species, and it was in the order K(+) < Na(+) < TMA(+) approximately Li(+). This behavior was well characterized by the specific features of each ion: the alpha(max) map could be well superimposed in the Stokes radius-crystal ion radius relationship of counterions. The alpha(max) dependence on Stokes radius was very similar to that of the B coefficient by Jones and Dole except in the case of Ca(2+). In principle, the smaller the value for B, the smaller alpha(max), indicating that a water structure breaker such as K(+) can more easily destroy the colloidal crystal structure. In other words, the effect of the counterion species on colloidal crystal stability follows the Hofmeister series. Including Ca(2+), the relationship was linear for the alpha(max) values plotted as a function of the limiting equivalent conductivity of small ions. A counterion with larger conductivity would be a stronger breaker for the colloidal crystal structure.

Electrostatic double layer force between a sphere and a planar substrate in the presence of previously deposited spherical particles

A finite element model of the electrostatic double layer interaction between an approaching colloidal particle and a small region of a charged planar surface containing four previously deposited particles is presented. The electrostatic interaction force experienced by the approaching particle is obtained by solving the Poisson-Boltzmann equation with appropriate boundary conditions representing this complex geometry. The interaction forces obtained from the detailed three-dimensional finite element simulations suggest that for the many-body scenario addressed here, the electrostatic double layer repulsion experienced by the approaching particle is less than the corresponding sphere-plate interaction due to the presence of the previously deposited particles. The reduction in force is quite significant when the screening length of the electric double layer becomes comparable to the particle radius (kappaa approximately 1). The results also suggest that the commonly used technique of pairwise addition of binary interactions can grossly overestimate the net electrostatic double layer interaction forces in such situations. The simulation methodology presented here can form a basis for investigating the influence of several previously deposited particles on the electrostatic repulsion experienced by a particle during deposition onto a substrate.

Fcc-bcc transition for Yukawa interactions determined by applied strain deformation

Calculations of the work required to transform between bcc and fcc phases yield a high-precision bcc-fcc transition line for monodisperse point Yukawa (screened-Coulomb) systems. Our results agree qualitatively but not quantitatively with recently published simulations and phenomenological criteria for the bcc-fcc transition. In particular, the bcc-fcc-fluid triple point lies at a higher inverse screening length than previously reported.

Three-body interactions in colloidal systems

We present a direct measurement of three-body interactions in a colloidal system comprised of three charged colloidal particles. Two of the particles have been confined by means of a scanned laser tweezers to a line-shaped optical trap where they diffused due to thermal fluctuations. Upon the approach of a third particle, attractive three-body interactions have been observed. The results are in qualitative agreement with additionally performed nonlinear Poisson-Boltzmann calculations, which also allow us to investigate the microionic density distributions in the neighborhood of the interacting colloidal particles.

Direct measurement of three-body interactions amongst charged colloids

Three-body interactions amongst three charged colloidal particles are measured in a deionized aqueous solution. Two of the particles are confined to an optical line trap while the third one is approached by means of a focused laser beam. From the observed particle configurations we extract the three-body potential which is found to be attractive and roughly of the same magnitude and range as the pair interactions. In addition, numerical calculations are performed which show qualitative agreement with the experimental results.