Role of polydispersity in anomalous interactions in electrostatically levitated colloidal systems

In this paper, we investigate the effects of using inverse analyses developed for monodisperse particles to extract particle-particle and particle-surface potentials from simulated interfacial colloidal configurations having finite-size polydispersity. Forward Monte Carlo simulations are used to generate three-dimensional equilibrium configurations of log normal-distributed polydisperse particles confined by gravity near an underlying surface. Particles remain levitated above the substrate and stabilized against aggregation by repulsive electrostatic Derjaguin-Landau-Verwey-Overbeek pair potentials. An inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method are used to obtain interactions from simulated distribution functions as a function of polydispersity (sigma), relative range of repulsion (kappa a), and projected interfacial concentration (rho). Both inverse analyses successfully recover input potentials for all monodisperse cases, but fail for polydispersities often encountered in experiments. For different conditions (sigma, kappa a, and rho), our results indicate softened short-range repulsion, anomalous long-range attraction, and apparent particle overlaps, which are similar to commonly reported observations in optical microscopy measurements of quasi-two-dimensional interfacial colloidal ensembles. By demonstrating signatures of, and limitations due to, polydispersity when extracting pair potentials from measured distribution functions, our specific goal is to provide a basis to objectively interpret and resolve the effects of polydispersity in optical microscopy experiments.     

Anomalous potentials from inverse analyses of interfacial polydisperse attractive colloidal fluids

This paper investigates effects of using monodisperse inverse analyses to extract particle-particle and particle-surface potentials from simulated interfacial colloidal fluids of polydisperse attractive particles. Effects of polydispersity are investigated as functions of particle concentration and attractive well depth and range for van der Waals and depletion potentials. Forward Monte Carlo simulations are used to generate particle distribution functions for polydisperse interfacial colloidal fluids from which inverted potentials are obtained using an inverse Ornstein-Zernike analysis and an inverse Monte Carlo simulation method. Attractive potentials are successfully recovered for monodisperse colloidal fluids, but polydispersity that is unaccounted for in inverse analyses produces (1) apparent softening of strong forces, (2) anomalous repulsive and attractive interactions, and (3) aphysical particle overlaps. This investigation provides insights into the role of polydispersity in altering the equilibrium structure and corresponding inverted potentials of attractive colloidal fluids near surfaces. These findings should assist the design and interpretation of optical microscopy experiments involving interfacial colloidal fluids similar to the simulated experiments reported here.     

Experimental verification of an exact evanescent light scattering model for TIRM

Total internal reflection microscopy (TIRM) is a method for the precise measurement of interaction potentials between a spherical colloidal particle and a wall. The method is based on single-particle evanescent wave light scattering. The well-established model used to interpret TIRM data is based on an exponential relation between scattering intensity and particle wall distance. However, applying this model for a certain range of experimental parameters leads to significant distortions of the measured potentials. Using a TIRM setup based on a two-wavelength illumination technique, we were able to directly measure the intensity distance relation revealing deviations from an exponential decay. The intensity-distance relations could be compared to scattering simulations taking into account exact experimental parameters and multiple reflections between a particle and the wall. Converging simulation results were independently obtained by the T-matrix method and the discrete sources method (DSM) and show excellent agreement with experiments. Using the new scattering model for data evaluation, we could reconstruct the correct potential shape for distorted interaction potentials as we demonstrate. The comparison of simulations to experiment intrinsically yields a new method to determine absolute particle-wall distances, a highly desired quantity in TIRM experiments.     

Direct measurement of single and ensemble average particle-surface potential energy profiles

This work involves the development of a novel technique that integrates total internal reflection and video microscopy methods to simultaneously measure single particle and ensemble average particle-surface interactions. For the 2 mum silica colloids and glass coverslip used in this study, particle size polydispersity is found to be a dominant factor in determining the distribution of single particle profiles about ensemble average profiles. In conjunction with this observation, chemical and physical nonuniformity are not evident in any of our measurements even with sensitivity to interactions on the order of kT. One advantage of using ensemble averaging in conjunction with time averaging is the ability to dramatically decrease the time required to measure average particle-wall interactions which scales inversely with interfacial particle concentration. A number of experimental issues are addressed in the development of this technique including (1) combining single particle distribution functions, (2) statistical sampling of distribution functions using both time and ensemble averaging, and (3) correcting overlapping scattering signals between adjacent particles. The capabilities of the ensemble averaging technique are also demonstrated to provide unique measurements of particle-surface interactions in metastable systems by selecting only height excursions of levitated particles when calculating potentials. Ultimately, this new technique provides several important advantages over single particle measurements, which provides a foundation for measuring interactions in increasingly complex interfacial systems.     

Charged micelle depletion attraction and interfacial colloidal phase behavior

Ensemble total internal reflection microscopy (TIRM) is used to directly measure the evolution of colloid-surface depletion attraction with increasing sodium dodecyl sulfate (SDS) concentration near the critical micelle concentration (CMC). Measured potentials are well described by a modified Asakura-Oosawa (AO) depletion potential in addition to electrostatic and van der Waals contributions. The modified AO potential includes effects of electrostatic interactions between micelles and surfaces via effective depletant dimensions in an excluded volume term and partitioning in an osmotic pressure term. Directly measured colloid-surface depletion potentials are used in Monte Carlo (MC) simulations to capture video microscopy (VM) measurements of micelle-mediated quasi-two-dimensional phase behavior including fluid, crystal, and gel microstructures. Our findings provide information to develop more rigorous and analytically simple models of depletion attraction in charged micellar systems.     

kT-scale colloidal interactions in high-frequency inhomogeneous AC electric fields. II. Concentrated ensembles

We report nonintrusive optical microscopy measurements of ensembles of polystyrene colloids in inhomogeneous AC electric fields as a function of field frequency and particle size. By using an inverse Monte Carlo (MC) simulation analysis of time-averaged particle microstructures, we sensitively measure induced dipole-dipole interactions on the kT energy scale. Measurements are reported for frequencies when the particle polarizability is greater and less than the medium, as well as the crossover between these conditions when dipole-dipole interactions vanish. By using measured single dipole-field interactions and associated parameters from Part I as input in the inverse analysis, the dipole-dipole interactions in this work are accurately modeled with no adjustable parameters for conditions away from the crossover frequency (i.e., |f(CM)| > 0). As dipolar interactions vanish at the crossover, a single frequency-dependent parameter is introduced to account for microstructures that appear to result from weak AC electro-osmotic flow induced interactions. By connecting quantitative measures of equilibrium microstructures and kT-scale dipole-field and dipole-dipole interactions, our findings provide a basis to understand colloidal assembly in inhomogeneous AC electric fields.     

Density profiles and solvation forces for a Yukawa fluid in a slit pore

The effect of varying wall-particle and particle-particle interactions on the density profiles near a single wall and the solvation forces between two walls immersed in a fluid of particles is investigated by grand canonical Monte Carlo simulations. Attractive and repulsive particle-particle and particle-wall interactions are modeled by a versatile hard-core Yukawa form. These simulation results are compared to theoretical calculations using the hypernetted chain integral equation technique, as well as with fundamental measure density functional theory (DFT), where particle-particle interactions are either treated as a first order perturbation using the radial distribution function or else with a DFT based on the direct-correlation function. All three theoretical approaches reproduce the main trends fairly well, but exhibit inconsistent accuracy, particularly for attractive particle-particle interactions. We show that the wall-particle and particle-particle attractions can couple together to induce a nonlinear enhancement of the adsorption and a related "repulsion through attraction" effect for the effective wall-wall forces. We also investigate the phenomenon of bridging, where an attractive wall-particle interaction induces strongly attractive solvation forces.     

Interpretation of conservative forces from Stokesian dynamic simulations of interfacial and confined colloids

This paper presents Stokesian dynamics simulations of experiments involving one or two charged colloids near either a single charged wall or confined between parallel charged walls. Equilibrium particle-particle and particle-wall interactions are interpreted from dynamic particle trajectories in simulations involving (1) a single particle levitated above a wall, (2) two particles below a wall, and (3) two particles confined between two parallel walls. By specifying only repulsive electrostatic Derjaguin-Landau-Verwey-Overbeek (DLVO) potentials and including multibody hydrodynamics, we successfully recover expected potentials in some cases, while anomalous attraction is observed in other cases. Attraction inferred in the latter simulations displays quantitative agreement with literature measurements when particle dynamics are interpreted using reported analyses. Because anomalous attraction is reproduced in simulations using only electrostatic repulsive DLVO potentials, our results reveal the one-dimensional analyses to be invalid for configurations that are inherently multidimensional via multibody hydrodynamics. Parameters related to experimental sampling of particle dynamics are also found to be critical for obtaining accurate potentials. We explain the anomalous attraction in each experiment using effective potentials, which can be employed in an a priori fashion to assist the confident design of future experiments involving interfacial and confined colloids. Ultimately, our findings reveal the importance of dimensionality and multibody hydrodynamics for understanding nonequilibrium dynamics of colloids near surfaces.     

Interfacial colloidal sedimentation equilibrium. II. Closure-based density functional theory

In Part I [R. E. Beckham and M. A. Bevan, J. Chem. Phys. 127, 164708 (2007)], results were presented for the sedimentation equilibrium of concentrated colloidal dispersions using confocal scanning laser microscopy experiments, Monte Carlo (MC) simulations, and a local density approximation perturbation theory. In this paper, we extended the modeling effort on those systems to include nonlocal density functional theory (DFT), which is capable of predicting the microstructure of the sediment at length scales comparable to the colloidal particle dimension. Specifically, we use a closure-based DFT formulation to predict interfacial colloidal sedimentation equilibrium density profiles. The colloid-colloid and colloid-surface interactions were modeled with DLVO screened electrostatic potentials using parameters taken directly from the experimental work. The DFT profiles were compared to the experimental and MC results from Part I. Good agreement was found for relatively dilute interfacial colloidal fluids, but agreement was less satisfactory as interfacial layering became more pronounced for conditions approaching the onset of interfacial crystallization. We also applied DFT in an inverse sense using the measured colloid density profile to extract the underlying colloid-surface potential; this can be thought of as a microscopic analog to the well-known procedure of using the macroscopic (coarse-grained) density profile to extract the osmotic equation of state. For the dilute interfacial fluid, the inverse DFT calculations reproduced the true colloid-surface potential to within 0.5kT at all elevations.     

Influence of magnetic interactions between clusters on particle orientational characteristics and viscosity of a colloidal dispersion composed of ferromagnetic spherocylinder particles: analysis by means of mean field approximation for a simple shear flow

We have theoretically investigated the particle orientational distribution and viscosity of a dense colloidal dispersion composed of ferromagnetic spherocylinder particles under an applied magnetic field. The mean field approximation has been applied to take into account the magnetic interactions of the particle of interest with the other ones that belong to the neighboring clusters, besides those that belong to its own cluster. The basic equation of the orientational distribution function, which is an integrodifferential equation, has approximately been solved by Galerkin's method and the method of successive approximation. Some of the main results obtained here are summarized as follows. Even when the magnetic interaction between particles is of the order of the thermal energy, the effect of particle-particle interactions on the orientational distribution comes to appear more significant with increasing volumetric fraction of particles; the orientational distribution function exhibits a sharper peak in the direction nearer to the magnetic field one as the volumetric fraction increases. Such a significant inclination of the particle in the field direction induces the large increase in viscosity in the range of larger values of the volumetric fraction. The above-mentioned characteristics of the orientational distribution and viscosity come to appear more significantly when the influence of the applied magnetic field is not so strong compared with that of magnetic particle-particle interactions.     

Effect of anionic surfactants on synthesis and self-assembly of silica colloidal nanoparticles

Effects of the anionic surfactants, sodium dodecyl sulfate and sodium oleate, on the formation and properties of silica colloidal nanoparticles were investigated. At a concentration of approximately 1 x 10(-3) M, adsorption of anionic surfactants increased particle size, monodispersity, and negative surface charge density of synthesized silica particles. As uniformity of particle size and particle-particle interactions increase, colloidal photonic crystals readily self-assemble without extensive washing of the synthesized silica nanoparticles. The photonic crystals diffract light in the visible region according to Bragg's law. The assembled colloidal particle arrays exhibit a face-centered cubic structure in dried thin films. This study offers a new approach for producing ordered colloidal silica thin films.     

Self-diffusion in submonolayer colloidal fluids near a wall

Theoretical expressions are developed to describe self-diffusion in submonolayer colloidal fluids that require only equilibrium structural information as input. Submonolayer colloidal fluids are defined for the purpose of this work to occur when gravity confines colloids near a planar wall surface so that they behave thermodynamically as two dimensional fluids. Expressions for self-diffusion are generalized to consider different colloid and surface interaction potentials and interfacial concentrations from infinite dilution to near fluid-solid coexistence. The accuracy of these expressions is demonstrated by comparing self-diffusion coefficients predicted from Monte Carlo simulated equilibrium particle configurations with standard measures of self-diffusion evaluated from Stokesian Dynamics simulated particle trajectories. It is shown that diffusivities predicted for simulated equilibrium fluid structures via multibody hydrodynamic resistance tensors and particle distribution functions display excellent agreement with values computed from mean squared displacements and autocorrelation functions of simulated tracer particles. Results are obtained for short and long time self-diffusion both parallel and normal to underlying planar wall surfaces in fluids composed of particles having either repulsive electrostatic or attractive van der Waals interactions. The demonstrated accuracy of these expressions for self-diffusion should allow their direct application to experiments involving submonolayer colloidal fluids having a range of interaction potentials and interfacial concentrations.     

Rheological properties and particle behaviors of a nondilute colloidal dispersion composed of ferromagnetic spherocylinder particles subjected to a simple shear flow: analysis by means of mean-field approximation for the two typical external magnetic field directions

We have analyzed the orientational distributions and rheological properties of a nondilute colloidal dispersion composed of ferromagnetic spherocylinder particles subjected to a simple shear flow. In order to understand the effects of the magnetic interactions between the particles, we have applied the mean-field theory to a nondilute colloidal dispersion for the two typical external magnetic field directions, that is, the direction parallel to the shear flow and the direction parallel to the angular velocity vector of the shear flow. The main results are summarized as follows. The particle-particle interactions suppress the Brownian motion of the particles and, therefore, make the particles incline toward the same direction. For the magnetic direction parallel to the shear flow, the influence of the particle-particle interactions makes the peak of the orientational distribution sharper and higher. The viscosity generally increases as the interactions between particles become stronger in the case where the effects of the shear flow and magnetic field are relatively small. For the magnetic direction parallel to the angular velocity vector of the shear flow, the influence of the particle-particle interactions on the orientational distribution appears significantly, when the influences of the shear flow and magnetic field are not so strong that the particles can be aligned sufficiently to form stable chainlike clusters in a certain direction.     

Techniques for measuring surface forces

The forces acting between colloidal particles and between surfaces are of utmost importance for determining the behaviour of dispersed systems and adhesion phenomena. Several techniques are now available for direct measurement of these surface forces. In this review we focus on some of these methods. Two techniques for measuring forces between macroscopic solid surfaces; the interferometric surface force apparatus, known as the SFA, and a novel instrument which is based on a bimorph force sensor, the so-called MASIF, are described in some detail. Forces between a macroscopic surface and a particle can be measured with the atomic force microscope (AFM) using a colloidal probe, or by employing total internal reflection microscipy (TIRM) to monitor the position of a colloidal particle trapped by a laser beam. We also describe two different techniques that can be used for measuring forces between “soft” interfaces, the thin film balance (TFB) for single foam, emulsion and solid/fluid/fluid films, and osmotic stress methods, commonly used for studying interactions in liquid crystalline surfactant phases or in concentrated dispersions. The advantages and limitations of each of these techniques are discussed and typical results are presented.     

Interfacial colloidal crystallization via tunable hydrogel depletants

We demonstrate an approach using temperature-dependent hydrogel depletants to thermoreversibly tune colloidal attraction and interfacial colloidal crystallization. Total internal reflection and video microscopy are used to measure temperature-dependent depletion potentials between approximately 2 microm silica colloids and surfaces as mediated by approximately 0.2 microm poly-N-isopropylacrylamide (PNIPAM) hydrogel particles. Measured depletion potentials are modeled using the Asakura-Oosawa theory while treating PNIPAM depletants as swellable hard spheres. Monte Carlo simulations using the measured potentials predict reversible, quasi-2D crystallization and melting at approximately 27 degrees C in quantitative agreement with video microscopy images of measured microstructures (i.e., radial distribution functions) over the temperature range of interest (20-29 degrees C). Additional measurements of short-time self-diffusivities display excellent agreement with predicted diffusivities by considering multibody hydrodynamic interactions and using a swellable hard sphere model for the PNIPAM solution viscosity. Our findings demonstrate the ability to quantitatively measure, model, and manipulate kT-scale depletion attraction and phase behavior as a means of formally engineering interfacial colloidal crystallization.     

Closure-based density functional theory applied to interfacial colloidal fluids

The behavior of dense colloidal fluids near surfaces can now be probed in great detail with experimental techniques like confocal microscopy. In fact, we are approaching a point where quantitative comparisons of experiment with particle-level theory, such as classical density functional theory (DFT), are appropriate. In a forward sense, we may use a known surface potential to predict a particle density distribution function from DFT; in an inverse sense, we may use an experimentally measured particle density distribution function to predict the underlying surface potential from DFT. In this paper, we tested the ability of the closure-based DFT of Zhou and Ruckenstein (J. Chem. Phys. 2000, 112, 8079-8082) to perform forward and inverse calculations on potential models commonly employed for colloidal particles and surfaces. To reduce sources of uncertainty in this initial study, Monte Carlo simulation results played the role of experimental data. The combination of Rogers-Young and modified-Verlet closures consistently performed well across the different potential models. For a reasonable range of choices of the density, temperature, and potential parameters, the inversion procedure yielded particle-surface potentials to an accuracy on the order of 0.1kT.     

Sonofragmentation of molecular crystals

Possible mechanisms for the breakage of molecular crystals under high-intensity ultrasound were investigated using acetylsalicylic acid (aspirin) crystals as a model compound for active pharmaceutical ingredients. Surprisingly, kinetics experiments ruled out particle-particle collisions as a viable mechanism for sonofragmentation. Two other possible mechanisms (particle-horn and particle-wall collisions) were dismissed on the basis of decoupling experiments. Direct particle-shock wave interactions are therefore indicated as the primary mechanism of sonofragmentation of molecular crystals.     

Two-dimensional crystallization of microspheres by a coplanar AC electric field

The particle-field and particle-particle interactions induced by alternating electric fields can be conveniently used for on-chip assembly of colloidal crystals. Two coplanar electrodes with a millimeter-sized gap between them are used here to assemble two-dimensional crystals from suspensions of either latex or silica microspheres. When an AC voltage is applied, the particles accumulate and crystallize on the surface between the electrodes. Light diffraction and microscopic observations demonstrate that the hexagonal crystal is always oriented with one axis along the direction of the field. The particles disassemble when the field is turned off, and the process can be repeated many times. The diffraction patterns from all consecutively formed crystals are identical. This assembly is driven by forces that depend on the electric field gradient, and a model is proposed involving a combination of dielectrophoresis and induced dipole chaining. The organization of large two-dimensional crystals allows characterization of the electrostatic interactions in the particle ensembles. The process can be controlled via the field strength, the frequency, and the viscosity of the liquid media. It could be used to make rudimentary optical switches or to separate mixtures of particles of different sizes.     

Electrostatically confined nanoparticle interactions and dynamics

We report integrated evanescent wave and video microscopy measurements of three-dimensional trajectories of 50, 100, and 250 nm gold nanoparticles electrostatically confined between parallel planar glass surfaces separated by 350 and 600 nm silica colloid spacers. Equilibrium analyses of single and ensemble particle height distributions normal to the confining walls produce net electrostatic potentials in excellent agreement with theoretical predictions. Dynamic analyses indicate lateral particle diffusion coefficients approximately 30-50% smaller than expected from predictions including the effects of the equilibrium particle distribution within the gap and multibody hydrodynamic interactions with the confining walls. Consistent analyses of equilibrium and dynamic information in each measurement do not indicate any roles for particle heating or hydrodynamic slip at the particle or wall surfaces, which would both increase diffusivities. Instead, lower than expected diffusivities are speculated to arise from electroviscous effects enhanced by the relative extent (kappaa approximately 1-3) and overlap (kappah approximately 2-4) of electrostatic double layers on the particle and wall surfaces. These results demonstrate direct, quantitative measurements and a consistent interpretation of metal nanoparticle electrostatic interactions and dynamics in a confined geometry, which provides a basis for future similar measurements involving other colloidal forces and specific biomolecular interactions.     

Anisotropic and hindered diffusion of colloidal particles in a closed cylinder

Video microscopy and particle tracking were used to measure the spatial dependence of the diffusion coefficient (D(α)) of colloidal particles in a closed cylindrical cavity. Both the height and radius of the cylinder were equal to 9.0 particle diameters. The number of trapped particles was varied between 1 and 16, which produced similar results. In the center of the cavity, D(α) turned out to be 0.75D(0) measured in bulk liquid. On approaching the cylindrical wall, a transition region of about 3 particle diameters wide was found in which the radial and azimuthal components of D(α) decrease to respective values of 0.1D(0) and 0.4D(0), indicating asymmetrical diffusion. Hydrodynamic simulations of local drag coefficients for hard spheres produced very good agreement with experimental results. These findings indicate that the hydrodynamic particle-wall interactions are dominant and that the complete 3D geometry of the confinement needs to be taken into account to predict the spatial dependence of diffusion accurately.