Three-body interactions in colloidal systems

Compact Expansion of a Repulsive Suspension

Short-range repulsion governs the dynamics of matter from atoms to animals. Using theory, simulations, and experiments, we find that an ensemble of repulsive particles spreads compactly with a sharp boundary, in contrast to the diffusive spreading of Brownian particles. Starting from the pair interactions, at high densities, the many-body dynamics follow nonlinear diffusion with a self-similar expansion, growing as t^{1/4}; At longer times, thermal motion dominates with the classic t^{1/2} expansion. A logarithmic growth controlled by nearest-neighbor interactions connects the two self-similar regimes.

Thermally active nanoparticle clusters enslaved by engineered domain wall traps

The stable assembly of fluctuating nanoparticle clusters on a surface represents a technological challenge of widespread interest for both fundamental and applied research. Here we demonstrate a technique to stably confine in two dimensions clusters of interacting nanoparticles via size-tunable, virtual magnetic traps. We use cylindrical Bloch walls arranged to form a triangular lattice of ferromagnetic domains within an epitaxially grown ferrite garnet film. At each domain, the magnetic stray field generates an effective harmonic potential with a field tunable stiffness. The experiments are combined with theory to show that the magnetic confinement is effectively harmonic and pairwise interactions are of dipolar nature, leading to central, strictly repulsive forces. For clusters of magnetic nanoparticles, the stationary collective states arise from the competition between repulsion, confinement and the tendency to fill the central potential well. Using a numerical simulation model as a quantitative map between the experiments and theory we explore the field-induced crystallization process for larger clusters and unveil the existence of three different dynamical regimes. The present method provides a model platform for investigations of the collective phenomena emerging when strongly confined nanoparticle clusters are forced to move in an idealized, harmonic-like potential.

Forces between solid surfaces in aqueous electrolyte solutions

This review addresses experimental findings obtained with direct force measurements between two similar or dissimilar solid surfaces in aqueous electrolyte solutions. Interpretation of these measurements is mainly put forward in terms of the classical theory of Derjaguin, Landau, Verwey, and Overbeek (DLVO). This theory invokes a superposition of attractive van der Waals forces and repulsive double layer forces. DLVO theory is shown to be extremely reliable, even in the case of multivalent ions. However, such a description is only successful, when appropriate surface charge densities, charge regulation characteristics, and ion pairing or complexation equilibria in solution are considered. Deviations from DLVO theory only manifest themselves at distances of typically below few nm. More long-ranged non-DLVO forces can be observed in some situations, particularly, in concentrated electrolyte solutions, in the presence of strongly adsorbed layers, or for hydrophobic surfaces. The latter forces probably originate from patch-charge surface heterogeneities, which can be induced by ion-ion correlation effects, charge fluctuations, or other types of surface heterogeneities.

Modeling the Electrostatics of Hollow Shell Suspensions: Ion Distribution, Pair Interactions, and Many-Body Effects

Electrostatic interactions play a key role in hollow shell suspensions as they determine their structure, stability, thermodynamics, and rheology and also the loading capacity of small charged species for nanoreservoir applications. In this work, fast, reliable modeling strategies aimed at predicting the electrostatics of hollow shells for one, two, and many colloids are proposed and validated. The electrostatic potential inside and outside a hollow shell with a finite thickness and a specific permittivity is determined analytically in the Debye-Hückel (DH) limit. An expression for the interaction potential between two such hollow shells is then derived and validated numerically. It follows a classical Yukawa form with an effective charge depending on the shell geometry, permittivity, and inner and outer surface charge densities. The predictions of the Ornstein-Zernike (OZ) equation with this pair potential to determine equations of state are then evaluated by comparison to results obtained with a Brownian dynamics algorithm coupled to the resolution of the linearized Poisson-Boltzmann and Laplace equations (PB-BD simulations). The OZ equation based on the DLVO-like potential performs very well in the dilute regime as expected, but also quite well, and more surprisingly, in the concentrated regime in which full spheres exhibit significant many-body effects. These effects are shown to vanish for shells with small thickness and high permittivity. For highly charged hollow shells, we propose and validate a charge renormalization procedure. Finally, using PB-BD simulations, we show that the cell model predicts the ion distribution inside and outside hollow shells accurately in both electrostatically dilute and concentrated suspensions. We then determine the shell loading capacity as a function of salt concentration, volume fraction, and surface charge density for nanoreservoir applications such as drug delivery, sensing, or smart coatings.

Potential of mean force between like-charged nanoparticles: Many-body effect

Ion-mediated interaction is important for the properties of polyelectrolytes such as colloids and nucleic acids. The effective pair interactions between two polyelectrolytes have been investigated extensively, but the many-body effect for multiple polyelectrolytes still remains elusive. In this work, the many-body effect in potential of mean force (PMF) between like-charged nanoparticles in various salt solutions has been comprehensively examined by Monte Carlo simulation and the nonlinear Poisson-Boltzmann theory. Our calculations show that, at high 1:1 salt, the PMF is weakly repulsive and appears additive, while at low 1:1 salt, the additive assumption overestimates the repulsive many-body PMF. At low 2:2 salt, the pair PMF appears weakly repulsive while the many-body PMF can become attractive. In contrast, at high 2:2 salt, the pair PMF is apparently attractive while the many-body effect can cause a weaker attractive PMF than that from the additive assumption. Our microscopic analyses suggest that the elusive many-body effect is attributed to ion-binding which is sensitive to ion concentration, ion valence, number of nanoparticles and charges on nanoparticles.

Three-body critical Casimir forces

Within mean-field theory we calculate universal scaling functions associated with critical Casimir forces for a system consisting of three parallel cylindrical colloids immersed in a near-critical binary liquid mixture. For several geometrical arrangements and boundary conditions at the surfaces of the colloids we study the force between two colloidal particles in the direction normal to their axes, analyzing the influence of the presence of a third particle on that force. Upon changing temperature or the relative positions of the particles we observe interesting features such as a change of sign of this force caused by the presence of the third particle. We determine the three-body component of the forces acting on one of the colloids by subtracting the pairwise forces from the total force. The three-body contribution to the total critical Casimir force turns out to be more pronounced for small surface-to-surface distances between the colloids as well as for temperatures close to criticality. Moreover, we compare our results with similar ones for other physical systems such as three atoms interacting via van der Waals forces.

Pairwise interactions of colloids in two-dimensional geometric confinement

We present the pairwise interaction behaviour of colloids confined to two-dimensional (2D) colloidal cages using optical laser tweezers. A single probe particle inside hexagonal cage particles at a planar oil-water interface is allowed to diffuse freely and the spring constant is extracted from its trajectories. To evaluate the effect of multibody interactions, the pair interactions between the probe particle and each cage particle are directly measured by using optical tweezers. Based on pairwise additivity, Monte Carlo simulations are used to compare the values of the spring constant obtained from experiments and simulations. We find that the multibody interactions negligibly occur and thus the particle interactions confined to such colloidal cages are highly pairwise. This work demonstrates that the use of the pairwise assumption in numerical simulations is rational when interparticle repulsive interactions are sufficiently strong, such as the particle interactions at fluid-fluid interfaces.

Micro-mutual-dipolar model for rapid calculation of forces between paramagnetic colloids

Typically, the force between paramagnetic particles in a uniform magnetic field is calculated using either dipole-based models or the Maxwell stress tensor combined with Laplace's equation for magnetostatics. Dipole-based models are fast but involve many assumptions, leading to inaccuracies in determining forces for clusters of particles. The Maxwell stress tensor yields an exact force calculation, but solving Laplace's equation is very time consuming. Here, we present a more elaborate dipole-based model: the micro-mutual-dipolar model. Our model has a time complexity that is similar to that of other dipole-based models but is much more accurate especially when used to calculate the force of small aggregates. Using this model, we calculate the force between two paramagnetic spheres in a uniform magnetic field and a circular rotational magnetic field and compare our results with those of other models. The forces for three-particle and ten-particle systems dispersed in two-dimensional (2D) space are examined using the same model. We also apply this model to calculate the force between two paramagnetic disks dispersed in 2D space. The micro-mutual-dipolar model is demonstrated to be useful for force calculations in dynamic simulations of small clusters of particles for which both accuracy and efficiency are desirable.

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.

Numerical calculation of interaction forces between paramagnetic colloids in two-dimensional systems

Typically the force between paramagnetic particles in a uniform magnetic field is described using the dipolar model, which is inaccurate when particles are in close proximity to each other. Instead, the exact force between paramagnetic particles can be determined by solving a three-dimensional Laplace's equation for magnetostatics under specified boundary conditions and calculating the Maxwell stress tensor. The analytical solution to this multi-boundary-condition Laplace's equation can be obtained by using a solid harmonics expansion in conjunction with the Hobson formula. However, for a multibody system, finite truncation of the Hobson formula does not lead to convergence of the expansion at all points, which makes the approximation physically unrealistic. Here we present a numerical method for solving this Laplace's equation for magnetostatics. This method uses a smoothed representation to replace all the boundary conditions. A two-step propagation is used to dramatically accelerate the calculation without losing accuracy. Using this method, we calculate the force between two paramagnetic particles in a uniform and a rotational external field and compare our results with other models. Furthermore, the many-body effects for three-particle, ten-particle, and 24-particle systems are examined using the same method. We also calculate the interaction between particles with different magnetic susceptibilities and particle diameters. The Laplace's equation solver method described in this article that is used to determine the force between paramagnetic particles is shown to be very useful for dynamic simulations for both two-particle systems and a large cluster of particles.

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.

Acousto-optically generated potential energy landscapes: potential mapping using colloids under flow

Optical potential energy landscapes created using acousto-optical deflectors are characterized via solvent-driven colloidal particles. The full potential energy of both single optical traps and complex landscapes composed of multiple overlapping traps are determined using a simple force balance argument. The potential of a single trap is shown to be well described by a Gaussian trap with stiffness found to be consistent with those obtained by a thermal equilibrium method. We also obtain directly the depth of the well, which (as with stiffness) varies with laser power. Finally, various complex systems ranging from double-well potentials to random landscapes are generated from individually controlled optical traps. Predictions of these landscapes as a sum of single Gaussian wells are shown to be a good description of experimental results, offering the potential for fully controlled design of optical landscapes, constructed from single optical traps.

Stochastic transport of particles across single barriers

Transport phenomena of interacting particles are of high interest for many applications in biology and mesoscopic systems. Here we present measurements on colloidal particles, which are confined in narrow channels on a substrate and interact with a barrier, which impedes the motion along the channel. The substrate of the particle is tilted in order for the particles to be driven towards the barrier and, if the energy gained by the tilt is large enough, surpass the barrier by thermal activation. We therefore study the influence of this barrier as well as the influence of particle interaction on the particle transport through such systems. All experiments are supported with Brownian dynamics simulations in order to complement the experiments with tests of a large range of parameter space which cannot be accessed in experiments.

Fast parallel interferometric 3D tracking of numerous optically trapped particles and their hydrodynamic interaction

Multi-dimensional, correlated particle tracking is a key technology to reveal dynamic processes in living and synthetic soft matter systems. In this paper we present a new method for tracking micron-sized beads in parallel and in all three dimensions - faster and more precise than existing techniques. Using an acousto-optic deflector and two quadrant-photo-diodes, we can track numerous optically trapped beads at up to tens of kHz with a precision of a few nanometers by back-focal plane interferometry. By time-multiplexing the laser focus, we can calibrate individually all traps and all tracking signals in a few seconds and in 3D. We show 3D histograms and calibration constants for nine beads in a quadratic arrangement, although trapping and tracking is easily possible for more beads also in arbitrary 2D arrangements. As an application, we investigate the hydrodynamic coupling and diffusion anomalies of spheres trapped in a 3 × 3 arrangement.

Electrostatics at the nanoscale

Electrostatic forces are amongst the most versatile interactions to mediate the assembly of nanostructured materials. Depending on experimental conditions, these forces can be long- or short-ranged, can be either attractive or repulsive, and their directionality can be controlled by the shapes of the charged nano-objects. This Review is intended to serve as a primer for experimentalists curious about the fundamentals of nanoscale electrostatics and for theorists wishing to learn about recent experimental advances in the field. Accordingly, the first portion introduces the theoretical models of electrostatic double layers and derives electrostatic interaction potentials applicable to particles of different sizes and/or shapes and under different experimental conditions. This discussion is followed by the review of the key experimental systems in which electrostatic interactions are operative. Examples include electroactive and "switchable" nanoparticles, mixtures of charged nanoparticles, nanoparticle chains, sheets, coatings, crystals, and crystals-within-crystals. Applications of these and other structures in chemical sensing and amplification are also illustrated.

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.

Interferometric 3D tracking of several particles in a scanning laser focus

High-Speed tracking of several particles allows measuring dynamic long-range interactions relevant to biotechnology and colloidal physics. In this paper we extend the successful technique of 3D back-focal plane interferometry to oscillating laser beams and show that two or more particles can be trapped and tracked with a precision of a few nanometers in all three dimensions. The tracking rate of several kHz is only limited by the scan speed of the beam steering device. Several tests proof the linearity and orthogonality of our detection scheme, which is of interest to optical tweezing applications and various metrologies. As an example we show the position cross-correlations of three diffusing particles in a scanning line optical trap.

Adsorption of poly(amido amine) (PAMAM) dendrimers on silica: importance of electrostatic three-body attraction

Adsorption of poly(amido amine) (PAMAM) dendrimers to silicon oxide surfaces was studied as a function of pH, ionic strength, and dendrimer generation. By combining optical reflectometry and atomic force microscopy (AFM), the adsorbed layers can be fully characterized and an unequivocal determination of the adsorbed mass becomes possible. For early stages, the adsorption process is transport limited and of first order with respect to the dendrimer solution concentration. For later stages, the surface saturates and the adsorbed dendrimers form loose but correlated liquidlike surface structures. This correlation is evidenced by a peak in the pair correlation function determined by AFM. The maximum adsorbed amount increases with increasing ionic strength and pH. The increase with the ionic strength is explained by the random sequential adsorption (RSA) model and electrostatic repulsion between the dendrimers. The adsorbing dendrimers interact by the repulsive screened Coulomb potential, whose range decreases with increasing ionic strength and thus leads to increasing adsorbed densities. The pH increase is interpreted as an effect of the substrate and is quantitatively explained by the extended three-body RSA model. This model stipulates the importance of a three-body interaction acting between two adsorbing dendrimers and the charged substrate. The presence of the charged substrate weakens the repulsion between the adsorbing dendrimers and thus leads to higher surface densities. This effect can be interpreted as an additional attractive three-body interaction, which acts in addition to the usual two-body repulsion and originates from the additional screening of the Coulomb repulsion by the counterions accumulating in the diffuse layer.

Nonlinear Poisson-Boltzmann equation in spherical symmetry

The Poisson-Boltzmann problem in spherical symmetry has been considered using the distribution of a self-consistent potential around a charged grain in a thermal collisional plasma as an example. The qualitative patterns of all possible solutions have been presented and a study of their asymptotics has been carried out. It has been demonstrated that for large potentials it is possible to neglect the curvature of the grain surface and to use the solution of the plane problem. It has also been demonstrated that the electrical interaction of the grains is possible only at distances smaller than eight screening lengths.

Nonlinear effects in charge stabilized colloidal suspensions

Molecular-dynamics simulations are used to study the effective interactions in charged stabilized colloidal suspensions. We focus on highly charged macroions in the limit of low salt concentrations. Within this regime, nonlinear corrections to Debye-Hückel (DH) theory have to be considered. For non-bulk-like systems, such as isolated pairs or triples of macroions, we show that nonlinear effects can become relevant, which cannot be described by the charge renormalization concept [S. Alexander, J. Chem. Phys. 80, 5776 (1984)]. For an isolated pair of macroions, we find an almost perfect qualitative agreement between our simulation data and DH theory. However, on a quantitative level, neither DH theory nor the charge renormalization concept can be confirmed in detail. This seems mainly to be related to the fact that for small ion concentrations, microionic layers can strongly overlap, whereas, simultaneously, excluded volume effects are less important. In the case of isolated triples, where we compare between coaxial and triangular geometries, we find attractive corrections to pairwise additivity in the limit of small macroion separations and salt concentrations. These triplet interactions arise if all three microionic layers around the macroions exhibit a significant overlap. In contrast to the case of two isolated colloids, the charge distribution around a macroion in a triple is found to be anisotropic.

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.

Structure of highly asymmetric hard-sphere mixtures: An efficient closure of the Ornstein-Zernike equations

A simple modification of the reference hypernetted chain (RHNC) closure of the multicomponent Ornstein-Zernike equations with bridge functions taken from Rosenfeld's hard-sphere bridge functional is proposed. Its main effect is to remedy the major limitation of the RHNC closure in the case of highly asymmetric mixtures--the wide domain of packing fractions in which it has no solution. The modified closure is also much faster, while being of similar complexity. This is achieved with a limited loss of accuracy, mainly for the contact value of the big sphere correlation functions. Comparison with simulation shows that inside the RHNC no-solution domain, it provides a good description of the structure, while being clearly superior to all the other closures used so far to study highly asymmetric mixtures. The generic nature of this closure and its good accuracy combined with a reduced no-solution domain open up the possibility to study the phase diagram of complex fluids beyond the hard-sphere model.

Influence of Confinement on the Electrostatic Interaction between Charged Colloids: a (N,V,T) Monte Carlo Study within Hyperspherical Geometry

Monte Carlo simulations within closed hyperspherical geometry are used to analyze the ionic distribution around two confined charged colloids to determine the origin of the net attraction recently reported in the literature. A scaling procedure is used to compare our numerical results obtained with small ideal colloids with the conclusion of the measurements performed with large silica colloids. Although no electrostatic attraction is detected under confinement, our simulations exhibit a significant reduction of the electrostatic repulsion between charged colloids confined between two weakly charged walls. After rescaling to reproduce the electrostatic repulsion between large confined colloids, our numerical results are qualitatively consistent with the reported attraction because we reasonably expect a reduction of the electrostatic force between such confined colloids below the order of magnitude of their van der Waals attraction.