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.     

Electrophoretic Motion of a Spherical Particle with a Thick Double Layer in Bounded Flows

The electrophoretic mobility of a spherical colloidal particle with low zeta potential near a solid charged boundary is calculated numerically for arbitrary values of the double layer thickness by a generalization of Teubner's method to the case of bounded flow. Three examples are considered: a sphere near a nonconducting planar wall with electric field parallel to the wall, near a perfectly conducting planar wall with electric field perpendicular to the wall, and on the axis of a cylindrical pore with electric field parallel to the axis. The results are compared with recent analytical calculations using the method of reflections. For the case of a charged sphere near a neutral surface, the reflection results are quite good, provided there is no double layer overlap, in which case there can be extra effects for constant potential particles that are entirely missed by the analytical expressions. For a neutral sphere near a charged surface, the reflection results are less successful. The main reason is that the particle feels the profile of the electroosmotic flow, an effect ignored by construction in the method of reflections. The general case is a combination of these, so that the reflections are more reliable when the electrophoretic motion dominates the electroosmotic flow. The effect on particle mobility of particle-wall interactions follows the trend expected on geometric grounds in that sphere-plane interactions are stronger than sphere-sphere interactions and the effect on a sphere in a cylindrical pore is stronger still. In the latter case, particle mobility can fall by more than 50% for thick double layers and a sphere half the diameter of the pore. The agreement between numerical results and analytical results follows the same trend, being worst for the sphere in a pore. Nevertheless, the reflections can be reliable for some geometries if there is no double layer overlap. This is demonstrated for a specific example where reflection results have previously been compared with experiments on protein mobility through a membrane (J. Ennis et al., 1996, J. Membrane Sci. 119, 47). Copyright 1999 Academic Press.     

On the effect of induced electro-osmosis on a cylindrical particle next to a surface

The effect of induced electro-osmosis on a cylindrical particle positioned next to a planar surface (wall) is studied theoretically both under the thin double layer approximation utilizing the Smoluchowski slip velocity approximation and under thick electric double layer conditions by solving the Poisson-Nernst-Planck (PNP) equations. The imposed, undisturbed electric field is parallel to the planar surface. The induced hydrodynamic and electrostatic forces are calculated as functions of the particle's and the medium's dielectric constants and the distance between the particle and the surface. The resultant force acting on the particle is directed normal to and away from the wall. The presence of such a repulsive force may adversely affect the interactions between macromolecules suspended in solution and wall-immobilized molecules and may be significant to near-wall particle imaging velocimetry (PIV) in electrokinetic flows.     

DC-electrokinetic motion of colloidal cylinder(s) in the vicinity of a conducting wall

The boundary effects on DC-electrokinetic behavior of colloidal cylinder(s) in the vicinity of a conducting wall is investigated through a computational model. The contribution of the hydrodynamic drag, gravity, electrokinetic (i.e., electrophoretic and dielectrophoretic), and colloidal forces (i.e., forces due to the electrical double layer and van der Waals interactions) are incorporated in the model. The contribution of electrokinetic and colloidal forces are included by introducing the resulting forces as an external force acting on the particle(s). The colloidal forces are implemented with the prescribed expressions from the literature, and the electrokinetic force is obtained by integrating the corresponding Maxwell stress tensor over the particles' surfaces. The electrokinetic slip-velocity together with the thin electrical double layer assumption is applied on the surfaces. The position and velocity of the particles and the resulting electric and flow fields are obtained and the physical insight for the behavior of the colloidal cylinders are discussed in conjunction with the experimental observations in the literature.     

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.     

Electrodipping force acting on solid particles at a fluid interface

We report experimental results which show that the interfacial deformation around glass particles (radius, 200-300 microm) at an oil-water (or air-water) interface is dominated by an electric force, rather than by gravity. It turns out that this force, called for brevity "electrodipping," is independent of the electrolyte concentration in the water phase. The force is greater for oil-water than for air-water interfaces. Under our experimental conditions, it is due to charges at the particle-oil (instead of particle-water) boundary. The derived theoretical expressions, and the experiment, indicate that this electric force pushes the particles into water. To compute exactly the electric stresses, we solved numerically the electrostatic boundary problem, which reduces to a set of differential equations. Convenient analytical expressions are also derived. Both the experimental and the calculated meniscus profile, which are in excellent agreement, exhibit a logarithmic dependence at long distances. This gives rise to a long-range electric-field-induced capillary attraction between the particles, detected by other authors. Deviation from the logarithmic dependence is observed at short distances from the particle surface due to the electric pressure difference across the meniscus. The latter effect gives rise to an additional short-range contribution to the capillary interaction between two floating particles. The above conclusions are valid for either planar or spherical fluid interfaces, including emulsion drops. The electrodipping force, and the related long-range capillary attraction, can engender two-dimensional aggregation and self-assembly of colloidal particles. These effects could have implications for colloid science and the development of new materials.     

Collective mechanical behavior of multilayer colloidal arrays of hollow nanoparticles

The collective mechanical behavior of multilayer colloidal arrays of hollow silica nanoparticles (HSNP) is explored under spherical nanoindentation through a combination of experimental, numerical, and theoretical approaches. The effective indentation modulus E(ind) is found to decrease with an increasing number of layers in a nonlinear manner. The indentation force versus penetration depth behavior for multilayer hollow particle arrays is predicted by an approximate analytical model based on the spring stiffness of the individual particles and the multipoint, multiparticle interactions as well as force transmission between the layers. The model is in good agreement with experiments and with detailed finite element simulations. The ability to tune the effective indentation modulus, E(ind), of the multilayer arrays by manipulating particle geometry and layering is revealed through the model, where E(ind) = (0.725m(-3/2) + 0.275)E(mon) and E(mon) is the monolayer modulus and m is number of layers. E(ind) is seen to plateau with increasing m to E(ind_plateau) = 0.275E(mon) and E(mon) scales with (t/R)(2), t being the particle shell thickness and R being the particle radius. The scaling law governing the nonlinear decrease in indentation modulus with an increase in layer number (E(ind) scaling with m(-3/2)) is found to be similar to that governing the indentation modulus of thin solid films E(ind_solid) on a stiff substrate (where E(ind_solid) scales with h(-1.4) and also decreases until reaching a plateau value) which also decreases with an increase in film thickness h. However, the mechanisms underlying this trend for the colloidal array are clearly different, where discrete particle-to-particle interactions govern the colloidal array behavior in contrast to the substrate constraint on deformation, which governs the thickness dependence of the continuous thin film indentation modulus.     

Colloidal dynamics near a wall studied by evanescent wave light scattering: experimental and theoretical improvements and methodological limitations

The dynamics of colloidal spheres near to a wall is studied with an evanescent wave scattering setup that allows for an independent variation of the components of the scattering wave vector normal and parallel to the wall. The correlation functions obtained with this novel instrumentation are interpreted on the basis of an expression for their short time behavior that includes hydrodynamic interactions between the colloidal spheres and the wall. The combination of the evanescent wave scattering setup and the exact expression for the short time behavior of correlation functions allows for an unambiguous measurement of the particle mobility parallel and normal to the wall by means of light scattering. It is possible to measure the viscous wall drag effect on the dynamics of particles with radii as small as 27 nm, where, however, the method reaches its limits due to the low scattering intensities of such small particles.     

Dynamic adhesion behavior of micrometer-scale particles flowing over patchy surfaces with nanoscale electrostatic heterogeneity

The dynamic adhesion behavior of micrometer-scale silica particles is investigated numerically for a low Reynolds number shear flow over a planar collecting wall with randomly distributed electrostatic heterogeneity at the 10-nanometer scale. The hydrodynamic forces and torques on a particle are coupled to spatially varying colloidal interactions between the particle and wall. Contact and frictional forces are included in the force and torque balances to capture particle skipping, rolling, and arrest. These dynamic adhesion signatures are consistent with experimental results and are reminiscent of motion signatures observed in cell adhesion under flowing conditions, although for the synthetic system the particle–wall interactions are controlled by colloidal forces rather than physical bonds between cells and a functionalized surface. As the fraction of the surface (Θ) covered by the cationic patches is increased from zero, particle behavior sequentially transitions from no contact with the surface to skipping, rolling, and arrest, with the threshold patch density for adhesion (Θ_crit) always greater than zero and in quantitative agreement with experimental results. The ionic strength of the flowing solution determines the extent of the electrostatic interactions and can be used to tune selectively the dynamic adhesion behavior by modulating two competing effects. The extent of electrostatic interactions in the plane of the wall, or electrostatic zone of influence, governs the importance of spatial fluctuations in the cationic patch density and thus determines if flowing particles contact the wall. The distance these interactions extend into solution normal to the wall determines the strength of the particle–wall attraction, which governs the transition from skipping and rolling to arrest. The influence of Θ, particle size, Debye length, and shear rate is quantified through the construction of adhesion regime diagrams, which delineate the regions in parameter space that give rise to different dynamic adhesion signatures and illustrate selective adhesion based on particle size or curvature. The results of this study are suggestive of novel ways to control particle–wall interactions using randomly distributed surface heterogeneity.     

Vertical motion of a charged colloidal particle near an AC polarized electrode with a nonuniform potential distribution: theory and experimental evidence

Electroosmotic flow in the vicinity of a colloidal particle suspended over an electrode accounts for observed changes in the average height of the particle when the electrode passes alternating current at 100 Hz. The main findings are (1) electroosmotic flow provides sufficient force to move the particle and (2) a phase shift between the purely electrical force on the particle and the particle's motion provides evidence of an E2 force acting on the particle. The electroosmotic force in this case arises from the boundary condition applied when faradaic reactions occur on the electrode. The presence of a potential-dependent electrode reaction moves the likely distribution of electrical current at the electrode surface toward uniform current density around the particle. In the presence of a particle the uniform current density is associated with a nonuniform potential; thus, the electric field around the particle has a nonzero radial component along the electrode surface, which interacts with unbalanced charge in the diffuse double layer on the electrode to create a flow pattern and impose an electroosmotic-flow-based force on the particle. Numerical solutions are presented for these additional height-dependent forces on the particle as a function of the current distribution on the electrode and for the time-dependent probability density of a charged colloidal particle near a planar electrode with a nonuniform electrical potential boundary condition. The electrical potential distribution on the electrode, combined with a phase difference between the electric field in solution and the electrode potential, can account for the experimentally observed motion of particles in ac electric fields in the frequency range from approximately 10 to 200 Hz.     

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.     

Dielectrophoretic force on a sphere near a planar boundary

The small gap distance separating a spherical colloidal particle in electrophoretic motion from a planar nonconducting surface is a required parameter for calculating its electrophoretic mobility. In the presence of an externally applied electric field, this gap distance is determined by balancing the van der Waals, electrical double layer interaction, and gravitational forces with a dielectrophoretic (DEP) force. Here, the DEP force was determined analytically by integration of the Maxwell stress over the surface of the particle. The account of this force showed that its previous omission from the analysis always resulted in underpredicted gap distances. Furthermore, the DEP force dominated under conditions of low particle density or high electric field strength and led to much higher gap distances on the order of a few microns. In one particular case, a combination of low particle density and small particle size produced two possible equilibrium gap distances for the particle. However, the particle was unstable in the second equilibrium position when subjected to small perturbations. In general, larger particles had smaller gap sizes. The effects of four other parameters on gap distance were studied, and gap distances were found to increase with lower particle density, higher electric field strength, higher particle and wall zeta potentials, and lower Hamaker constants. Retardation effects on van der Waals attraction were considered.     

Depletion force between anisometric colloidal particles

A simple mathematical model for the depletion force between two arbitrarily shaped large convex colloidal particles immersed in a suspension of small spherical particles is proposed. Using differential geometry, the interaction potential is expressed in terms of the mean and Gaussian curvature of the particle surfaces. The accuracy of theoretical results is tested by Monte Carlo simulations for parallel and nonparallel circular cylinders. The agreement between theoretical results and simulated data is very good if the density of the depletion agent is not too high.     

Polarization and interactions of colloidal particles in ac electric fields

Micrometer-sized polystyrene particles form two-dimensional crystals in alternating current (ac) electric fields. The induced dipole-dipole interaction is the dominant force that drives this assembly. We report measurements of forces between colloidal particles in ac electric fields using optical tweezers and find good agreement with the point dipole model. The magnitude of the pair interaction forces depends strongly on the bulk solution conductivity and decreases as the ionic strength increases. The forces also decrease with increasing field frequency. The salt and frequency dependences are consistent with double layer polarization with a characteristic relaxation frequency omega(CD) approximately a(2)/D, where a is the particle radius and D is the ion diffusivity. This enables us to reinterpret the order-disorder transition reported for micrometer-sized polystyrene particles [Lumsdon et al., Langmuir 20, 2108 (2004)], including the dependence on particle size, frequency, and ionic strength. These results provide a rational framework for identifying assembly conditions of colloidal particles in ac fields over a wide range of parameters.     

Electrostatic double-layer interaction between spherical particles inside a rough capillary

Predictions of electrostatic double-layer interaction forces between two similarly charged spherical colloidal particles inside an infinitely long "rough" capillary are presented. A simple model of a rough cylindrical surface is proposed, which assumes the capillary wall to be a periodic function of axial position. The periodic roughness of the wall is characterized by the wavelength and amplitude of the undulations. The electrostatic double-layer interaction force between two spherical particles located axially inside this rough capillary is determined by solving the nonlinear Poisson-Boltzmann equation employing finite element analysis. The effect of surface roughness of the cylindrical enclosure on the interaction force between two particles is extensively studied on the basis of this model. The simulations are carried out for dimensionless amplitudes (amplitude/particle radii) ranging from 0.05 to 0.15 and scaled wavelengths (wavelength/particle radii) ranging from 0.4 to 4.0. The interaction force between the particles is significantly modified by the proximity of the rough capillary wall. Generally, the interaction force for rough capillaries oscillates around the corresponding interaction force in a smooth capillary depending on the magnitudes of the scaled amplitude and wavelength of the roughness. The influence of roughness on the electrostatic interactions becomes more pronounced when the surface potential of the cylinder wall is different from the sphere surface potentials. When the cylinder and the particle surfaces have large potential differences, the axial force experienced by a particle is dominated by the capillary roughness. There are dramatic oscillations of the force, which alternately becomes repulsive and attractive as the particle moves from the crest to the trough of the rough capillary wall. These results suggest that manipulation of colloidal particles in narrow microchannels may be subject to significant force variations owing to the roughness inherent in microfabricated channels etched on metal films.     

Observations of the near-wall accumulation of suspended particles due to shear and electroosmotic flow in opposite directions

On the basis of previous studies, the particles in a dilute (volume fractions φ_∞ < 4 × 10^–3) suspension in combined Poiseuille and electroosmotic “counterflow” at flow Reynolds numbers Re ≤ 1 accumulate, then assemble into structures called “bands,” within ∼6 μm of the channel wall. The experimental studies presented here use a small fraction of tracer particles labeled with a different fluorophore from the majority “bulk” particles to visualize the dynamics of individual particles in a φ_∞ = 1.7 × 10^–3 suspension. The results at two different near-wall shear rates and three electric field magnitudes E show that the near-wall particles are concentrated about 150-fold when the bands start to form, and are then concentrated about 200-fold to a maximum near-wall volume fraction of ∼0.34. The growth in the near-wall particles during this accumulation stage appears to be exponential. This near-wall particle accumulation is presumably driven by a wall-normal “lift” force. The observations of how the particles accumulate near the wall are compared with recent analyses that predict that suspended particles subject to shear flow and a dc electric field at small particle Reynolds numbers experience such a lift force. A simple model that assumes that the particles are subject to this lift force and Stokes drag suggests that the force driving particles toward the wall, of O(10^–17 N), is consistent with the time scales for particle accumulation observed in the experiments.     

Electrokinetic lift force for a charged particle moving near a charged wall — a modified theory and experiment

We present the refined theory of the electrokinetic lift force for a charged particle moving at a charged wall at the distance much larger than the double layer thickness. The theory is based on the lubrication approximation for the solution of the Stokes equation for the flow around a long cylinder moving near a solid wall. The “thin double layer” approximation is used to solve the ionic balance and electro-osmotic flow equations. The electrokinetic lift force is then obtained by integration of the viscous stress tensor as well as the Maxwell stress tensor over the particle surface. The resulting lift force for the cylinder translating, rotating at the wall as well as for the stationary cylinder in the wall shear flow, is considered. Following this, we apply the Derjaguin approximation to transform the obtained results to the sphere–wall geometry and we compare our theoretical predictions with the measurements of the electrokinetic lift force performed in the “colloidal particle collider” apparatus for the latex particles suspended in the glycerol–water solutions. Our theoretical results for the electrokinetic lift force exceeds by several orders of magnitude one obtained from the previously developed theory and are in a good agreement with experimental findings.     

Electromagnetic interactions in plasmonic nanoparticle arrays

Single two-dimensional planar silver arrays and one-dimensional linear gold chains of nanoparticles were investigated by dark-field surface plasmon spectroscopy and studied as a function of interparticle distance, particle size, and number of particles. In agreement with recent theoretical predictions, a red shift of the surface plasmon resonance occurring in two-dimensional arrays was found for lattice spacings below 200 nm. This red shift is associated with a significant broadening of the resonance and is attributed to the onset of near-field interactions. We found that the relative contributions of the long-range and short-range interactions in two-dimensional arrays of particles are fundamentally different to those occurring in individual linear chains.     

Adsorption and order formation of colloidal nanoparticles on a substrate: a Brownian dynamics study

Adsorption process and order formation of electrostatically stabilized colloidal particles with a radius of 50 nm onto a planar surface with countercharge are examined. We perform Brownian dynamics simulations with a new three-dimensional cell model, in which the particle-particle and particle-substrate interactions are modeled based on the DLVO theory. The simulations yield the following results: (1) a larger bulk concentration would be required for larger kappaa to reach order formation to compensate for the decrease in the bulk potential; (2) the phase transition from a disordered to an ordered structure of the adsorbed particles on the substrate is considered to be of the Kirkwood-Alder type of transition through the examination of the two-dimensional pressure of the adsorbed particles; (3) the adsorbed particles are found to form a hexagonally ordered array, only if what we call "one-directional average force" acting on an adsorbed particle exceeds a critical value, which is independent of the ionic strength, or the interaction potentials. The critical value of the one-directional average force is interpreted as the force needed to keep an ordered structure by localizing adsorbed particles at fixed positions. In addition, the critical force is used to develop a new model to estimate the surface coverage at the order-disorder transition and it is demonstrated that the new model gives better estimation than other models previously reported.