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.     

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.     

Rotational diffusion of colloidal particles near confining walls

We study the rotational diffusion of a spherical colloid confined in a narrow channel between parallel plane hard walls. The walls damp translational diffusion much more than rotational diffusion so that there is expected to be little translation-rotation coupling. Using a recent calculation of the nonisotropic rotational mobilities arising from the hydrodynamic interactions with the walls, we set up the rotational Smoluchowski equation for either a particle with a permanent dipole moment or a polarizable particle with axisymmetric polarizabilities subject to an external electric field. Using the Smoluchowski equation dynamics we calculate the time-correlation functions of orientation that are measured in depolarized light scattering for the cases of no external field, external field normal to the walls, and external field parallel to the walls. The decay of correlations is shown to be given by a weighted sum of decaying exponentials and can be characterized by an initial and a mean characteristic decay time. The weights and decay rates of each component and the characteristic decay times are studied numerically for a range of field strengths. The nonisotropic rotational mobilities make these decay times highly sensitive to the distance of the particle from the confining walls. This position dependence can be used as a method of measuring the rotational mobilities or, conversely, the rate of decay of correlations can be used as a probe of particle position between the confining walls.     

Diffusion of spherical particles in microcavities

The diffusive motion of a colloidal particle trapped inside a small cavity filled with fluid is reduced by hydrodynamic interactions with the confining walls. In this work, we study these wall effects on a spherical particle entrapped in a closed cylinder. We calculate the diffusion coefficient along the radial, azimuthal, and axial direction for different particle positions. At all locations the diffusion is smaller than in a bulk fluid and it becomes anisotropic near the container's walls. We present a simple model which reasonably well describes the simulation results for the given dimensions of the cylinder, which are taken from the recent experimental work.     

Adsorption and disjoining pressure isotherms of confined polymers using dissipative particle dynamics

The adsorption and disjoining pressure isotherms of polymers confined by planar walls are obtained using Monte Carlo (MC) simulations in the Grand Canonical (GC) ensemble in combination with the mesoscopic technique known as dissipative particle dynamics (DPD). Two models of effective potentials for the confining surfaces are used: one with both an attractive and a repulsive term and one with a purely repulsive term. As for the polymer, seven-bead linear model of polyethylene glycol (PEG) dissolved in water is used. The results indicate remarkably good agreement between the trends shown by our adsorption isotherms and those obtained from experiments of PEG on oxide surfaces. Additionally, the disjoining pressure isotherm of water shows oscillations, while those of PEG display the same trend for both wall models. Moreover, it is found that the disjoining pressure isotherms are in qualitative agreement with those from experiments on confined linear polymers.     

A modeling approach to describe the adhesion of rough, asymmetric particles to surfaces

A combined theoretical and experimental study of the adhesion of alumina particles and polystyrene latex spheres to silicon dioxide surfaces was performed. A boundary element technique was used to model electrostatic interactions between micron-scale particles and planar surfaces when the particles and surfaces were in contact. This method allows quantitative evaluation of the effects of particle geometry and surface roughness on the electrostatic interaction. The electrostatic interactions are combined with a previously developed model for van der Waals forces in particle adhesion. The combined model accounts for the effects of particle and substrate geometry, surface roughness and asperity deformation on the adhesion force. Predictions from the combined model are compared with experimental measurements made with an atomic force microscope. Measurements are made in aqueous solutions of varying ionic strength and solution pH. While van der Waals forces are generally dominant when particles are in contact with surfaces, results obtained here indicate that electrostatic interactions contribute to the overall adhesion force in certain cases. Specifically, alumina particles with complex geometries were found to adhere to surfaces due to both electrostatic and van der Waals interactions, while polystyrene latex spheres were not affected by electrostatic forces when in contact with various surfaces.     

Evaluation of the dynamic properties of polymer latex particles interacting with quartz interface by evanescent wave dynamic light scattering

The diffusion behavior of polymer latex particles in dispersion near the quartz interface has been estimated by evanescent wave dynamic light scattering (EVDLS) technique. The diffusion coefficient of the particles was measured as a function of the distance between the particle and interface. The apparent diffusion coefficient estimated by EVDLS was small for particles near the interface and increased upon increasing the distance from the interface, and then saturated at a certain value which is close to the value expected for free-motion. The range of the distance over which diffusion was affected by interaction with the interface depended on the added salt concentration. This means that the diffusion of the particle is influenced by an electrostatic interaction between the particle and quartz interface in addition to the hydrodynamic effect near the wall. This range was found to be more than 800?nm at 0?M salt condition but about 400?nm at 10^-4 and 10^-3?M salt conditions. Hence it is appropriate to say that the hydrodynamic effect reaches up to 400?nm and the electrostatic effect is longer ranged, more than 800?nm, for the system studied here. The EVDLS technique is a very powerful tool for quantitative estimations of the dynamic behavior of the particle near the interface and for estimation of the range where the wall effect is dominant. EVDLS will give us an answer to the question of “where is the ‘interface’ and where is the ‘bulk’?”.     

Electronic properties of Ag nanoparticle arrays. A Kelvin probe and high resolution XPS study

The electronic properties of citrate stabilised Ag nanoparticles with sizes ranging from 4 to 35 nm were investigated by the Kelvin probe method and high resolution XPS. Two and three dimensional assemblies of the particles were prepared by electrostatic adsorption from aqueous solution onto poly-l-lysine modified surfaces. The work function of the Ag particles increased from 5.29 +/- 0.05 to 5.53 +/- 0.05 eV as the particle size decreased. These values are approximately 0.8 eV higher than for clean polycrystalline Ag surfaces. The origin of these remarkable high work functions cannot be explained in terms of either citrate induced changes in the surface dipole or image forces in the confined metallic domains. High resolution XPS spectra of the Ag 3d(5/2) core level were characterised by broad bands and a 0.4 eV shift towards lower binding energies for the smallest particles. Comparisons with reported studies on extended Ag surfaces indicate that as-grown particles exhibit partially oxidised surfaces. The behaviour of the work function further suggests that the strength of the Ag-O bonding increases with decreasing particle sizes. These findings are highly relevant to the interpretation of the catalytic properties of Ag nanoparticles.     

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.     

Brush-sheathed particles diffusing at brush-coated surfaces in the thermally responsive PNIPAAm system

Phase-contrast microscopy and particle tracking algorithms are used to study the near-surface diffusion of poly(N-isopropylacrylamide) (PNIPAAm) brush functionalized micron-sized silica microspheres after sedimentation from aqueous suspension onto planar substrates coated with a similar polymer brush above and below the lower critical solution temperature (LCST) of PNIPAAm, 32 degrees C. A small negative charge on the wall and the particles (zeta potential = -6 mV) prevents adhesion above and below the LCST. The near-surface translational diffusion coefficient (D(surface)) is compared to the bulk-phase translational diffusion coefficient (D(bulk)), which was measured by dynamic light scattering. We find that D(surface)/D(bulk) is approximately equal to 0.6 at temperatures T < 32 degrees C but rises abruptly to approximately 0.8-0.9 at T > 32 degrees C. Near-surface diffusion is expected to be slower than bulk diffusion owing to hydrodynamic coupling to the wall, implying reduced hydrodynamic coupling at the higher temperatures, perhaps mediated by enhanced electrostatic repulsion above the LCST transition.     

The rheology and microstructure of charged colloidal suspensions

The effects of electric charge interation and particle correlations on suspension rheology are examined. A one-component fluid analysis using a Smoluchowski equation for the equilibrium structure is applied to charged suspensions of spherical colloids under shear. The frequency dependent modulus and viscosity, predicted as functions of particle and added salt concentrations, are compared with published rheological measurements on model suspensions. Recent improvements in the statistical mechanical theories for the equilibrium microstructure, its nonequilibrium deformation, and the bulk shear stresses are included. The direct electrostatic interaction is found to drive the divergence in the shear viscosity near the liquid-solid phase transition. Extensions of the theory predict the elastic modulus of binary mixtures of charged colloids. Estimates of the primary electroviscous effect, hydrodynamic interactions, and errors in the Yukawa limiting form for the potential and applications of asymptotic theories are presented. Predictions for the rheology based on effective hard-sphere models are found to be reasonable when using a parameter fit from the equilibrium phase behavior. Mean-field mode coupling theories predict larger relaxation times than calculated from the Smoluchowski equation (=SE). A study of binary mixing effects on elasticity shows non-ideal behavior. It is noted that equilibrium structural information can be used to resolve discrepancies between the theoretical predictions and the measured rheology.     

Hydrodynamic crossover in dynamic microparticle adhesion on surfaces of controlled nanoscale heterogeneity

This note documents the crossover from a regime where shear flow hinders microparticle adhesion on collecting surfaces to that where increased flow aids particle capture. Flow generally works against adhesion and successfully hinders particle capture when the net physicochemical attractions between the particles and collector are weak compared with hydrodynamic forces on the particle. Conversely, with strong attractions between particles and collector, flow aids particle capture by increasing the mass transport of particles to the interfacial region. Here, local hydrodynamics still generally oppose adhesion but are insufficient to pull particles off of the surface. Thus, flow actually increases the particle capture rate through the increased transport to the surface. These behaviors are demonstrated using 1 mum silica spheres flowing over electrostatically heterogeneous (length scales near 10 nm) collecting surfaces at shear rates from 22 to 795 s(-1). The net surface charge on the collector is varied systematically from strongly negative (pure silica) to strongly positive (a saturated polycationic overlayer), demonstrating the interplay between physicochemical and hydrodynamic contributions. These results clearly apply to situations where heterogeneous particle-surface interactions are electrostatic in nature; however, qualitatively similar behavior was previously reported for the effect receptor density on bacterial adhesion.     

Generalized Langevin dynamics of a nanoparticle using a finite element approach: thermostating with correlated noise

A direct numerical simulation (DNS) procedure is employed to study the thermal motion of a nanoparticle in an incompressible Newtonian stationary fluid medium with the generalized Langevin approach. We consider both the Markovian (white noise) and non-Markovian (Ornstein-Uhlenbeck noise and Mittag-Leffler noise) processes. Initial locations of the particle are at various distances from the bounding wall to delineate wall effects. At thermal equilibrium, the numerical results are validated by comparing the calculated translational and rotational temperatures of the particle with those obtained from the equipartition theorem. The nature of the hydrodynamic interactions is verified by comparing the velocity autocorrelation functions and mean square displacements with analytical results. Numerical predictions of wall interactions with the particle in terms of mean square displacements are compared with analytical results. In the non-Markovian Langevin approach, an appropriate choice of colored noise is required to satisfy the power-law decay in the velocity autocorrelation function at long times. The results obtained by using non-Markovian Mittag-Leffler noise simultaneously satisfy the equipartition theorem and the long-time behavior of the hydrodynamic correlations for a range of memory correlation times. The Ornstein-Uhlenbeck process does not provide the appropriate hydrodynamic correlations. Comparing our DNS results to the solution of an one-dimensional generalized Langevin equation, it is observed that where the thermostat adheres to the equipartition theorem, the characteristic memory time in the noise is consistent with the inherent time scale of the memory kernel. The performance of the thermostat with respect to equilibrium and dynamic properties for various noise schemes is discussed.     

Charge stabilization in nonpolar solvents

While the important role of electrostatic interactions in aqueous colloidal suspensions is widely known and reasonably well-understood, their relevance to nonpolar suspensions remains mysterious. We measure the interaction potentials of colloidal particles in a nonpolar solvent with reverse micelles. We find surprisingly strong electrostatic interactions characterized by surface potentials, |ezeta|, from 2.0 to 4.4 k(B)T and screening lengths, kappa(-1), from 0.2 to 1.4 microm. Interactions depend on the concentration of reverse micelles and the degree of confinement. Furthermore, when the particles are weakly confined, the values of |ezeta| and kappa extracted from interaction measurements are consistent with bulk measurements of conductivity and electrophoretic mobility. A simple thermodynamic model, relating the structure of the micelles to the equilibrium ionic strength, is in good agreement with both conductivity and interaction measurements. Since dissociated ions are solubilized by reverse micelles, the entropic incentive to charge a particle surface is qualitatively changed from aqueous systems, and surface entropy plays an important role.     

Experimental study of the effect of external electric fields on interfacial dynamics of colloidal particles

The interaction of colloidal particles with a planar surface (i.e., wall) in the presence of an electric field applied parallel to the planar surface is of interest in various microfluidic devices. Evanescent wave-based particle-tracking velocimetry was used to investigate the dynamics of a dilute suspension of polystyrene and silica particles (radii a = 110-463 nm) in a monovalent electrolyte solution with a Debye length of 6.8 nm driven through a microchannel by external electric fields E = 15-31 V/cm over the first 300 nm next to the channel wall. The particle velocity parallel to the wall due to electrophoresis and electroosmosis was in good agreement with the Helmholtz-Smoluchowski relation, and the hydrodynamic interactions between the wall and the particles were negligible, for all particle types. Measurements of the distribution of particles along the wall-normal coordinate, however, suggest that an additional force as great as 30 fN that repels the negatively charged particles away from the wall is induced by nonzero E. The results suggest that the magnitude of this force scales as E(2) and a(2) but is independent of the particle ζ-potential, in agreement with previous theoretical studies. However, estimates of the force assuming that the particles have a Boltzmann distribution were up to 40 times greater than the theoretical predictions, which only considered "remote" particle-wall interactions. These results are, to our knowledge, the first to observe a repulsive wall-normal force due to an applied electric field for near-wall colloidal particles.     

Orientational orders of small anisotropic molecules confined in slit pores

Based on a constant-pressure Monte Carlo molecular simulation, we have studied orientationally ordered transitions of small anisotropic molecules confined in two parallel hard walls. These molecules are modeled by the hard Gaussian overlap model. The molecular elongations of the chosen molecules are so small that the molecules cannot form stable liquid-crystal (LC) phases in the bulk. But in the slit pores, we found, while the distance between two walls of the pores decreases to the molecular scale, an orientationally ordered phase can form. It shows that even hard confining surfaces favor the alignment of the small anisotropic molecules. Thus we conclude that the required molecular elongation for forming LC phases will decrease in confinement. Our results indicate that some non-LC small molecules may form stable LC phases due to the inducement of confining surfaces.     

Anisotropic diffusion in confined colloidal dispersions: the evanescent diffusivity

We employ an analogy to traditional dynamic light scattering to describe the inhomogeneous and anisotropic diffusion of colloid particles near a solid boundary measured via evanescent wave dynamic light scattering. Following this approach, we generate new expressions for the short-time self- and collective diffusivities of colloidal dispersions with arbitrary volume fraction. We use these expressions in combination with accelerated Stokesian dynamics simulations to calculate the diffusivities in the limit of large and small scattering wave numbers for evanescent penetration depths ranging from four particle radii to one-fifth of a particle radius and volume fractions from 10% to 40%. We show that at high volume fractions, and larger penetration depths, the boundaries have little effect on the dynamics of the suspension parallel to the wall since, to a first approximation, the boundary acts hydrodynamically much as another nearby particle. However, near and normal to the wall, the diffusivity shows a strong dependence on penetration depth for all volume fractions. This is due to the lubrication interactions between the particles and the boundary as the particle moves relative to the wall. These results are novel and comprehensive with respect to the range of penetration depth and volume fraction and provide a complete determination of the effect of hydrodynamic interactions on colloidal diffusion adjacent to a rigid boundary.     

硅酸盐岩石微孔中流体混合物扩散系数的分子动力学模拟

用计算机模拟研究了受限于硅酸盐岩石微孔中等摩尔氩氪流体混合物的扩散系数.为验证程序的正确性,分别将宏量条件下氩和氪的分子动力学模拟值与文献实验值进行了比较.然后,对受限于硅酸岩平板壁面中氩氪流体混合物的扩散系数进行了模拟计算.根据计算结果提出了五个不同对比壁间距下氩氪混合流体扩散系数的经验模型,这些模型可应用于工业估算.     

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.