Transport properties of suspensions–critical assessment of Beenakker-Mazur method

Short- and long-time diffusion and dynamic scaling in suspensions of charged colloidal particles

We report on a comprehensive theory-simulation-experimental study of collective and self-diffusion in concentrated suspensions of charge-stabilized colloidal spheres. In theory and simulation, the spheres are assumed to interact directly by a hard-core plus screened Coulomb effective pair potential. The intermediate scattering function, f_c(q, t), is calculated by elaborate accelerated Stokesian dynamics (ASD) simulations for Brownian systems where many-particle hydrodynamic interactions (HIs) are fully accounted for, using a novel extrapolation scheme to a macroscopically large system size valid for all correlation times. The study spans the correlation time range from the colloidal short-time to the long-time regime. Additionally, Brownian Dynamics (BD) simulation and mode-coupling theory (MCT) results of f_c(q, t) are generated where HIs are neglected. Using these results, the influence of HIs on collective and self-diffusion and the accuracy of the MCT method are quantified. It is shown that HIs enhance collective and self-diffusion at intermediate and long times. At short times self-diffusion, and for wavenumbers outside the structure factor peak region also collective diffusion, are slowed down by HIs. MCT significantly overestimates the slowing influence of dynamic particle caging. The dynamic scattering functions obtained in the ASD simulations are in overall good agreement with our dynamic light scattering (DLS) results for a concentration series of charged silica spheres in an organic solvent mixture, in the experimental time window and wavenumber range. From the simulation data for the time derivative of the width function associated with f_c(q, t), there is indication of long-time exponential decay of f_c(q, t), for wavenumbers around the location of the static structure factor principal peak. The experimental scattering functions in the probed time range are consistent with a time-wavenumber factorization scaling behavior of f_c(q, t) that was first reported by Segrè and Pusey [Phys. Rev. Lett. 77, 771 (1996)] for suspensions of hard spheres. Our BD simulation and MCT results predict a significant violation of exact factorization scaling which, however, is approximately restored according to the ASD results when HIs are accounted for, consistent with the experimental findings for f_c(q, t). Our study of collective diffusion is amended by simulation and theoretical results for the self-intermediate scattering function, f_s(q, t), and its non-Gaussian parameter α₂(t) and for the particle mean squared displacement W(t) and its time derivative. Since self-diffusion properties are not assessed in standard DLS measurements, a method to deduce W(t) approximately from f_c(q, t) is theoretically validated.

Generalization of Clausius-Mossotti approximation in application to short-time transport properties of suspensions

In 1983, Felderhof, Ford, and Cohen gave microscopic explanation of the famous Clausius-Mossotti formula for the dielectric constant of nonpolar dielectric. They based their considerations on the cluster expansion of the dielectric constant, which relates this macroscopic property with the microscopic characteristics of the system. In this article, we analyze the cluster expansion of Felderhof, Ford, and Cohen by performing its resummation (renormalization). Our analysis leads to the ring expansion for the macroscopic characteristic of the system, which is an expression alternative to the cluster expansion. Using similarity of structures of the cluster expansion and the ring expansion, we generalize (renormalize) the Clausius-Mossotti approximation. We apply our renormalized Clausius-Mossotti approximation to the case of the short-time transport properties of suspensions, calculating the effective viscosity and the hydrodynamic function with the translational self-diffusion and the collective diffusion coefficient. We perform calculations for monodisperse hard-sphere suspensions in equilibrium with volume fraction up to 45%. To assess the renormalized Clausius-Mossotti approximation, it is compared with numerical simulations and the Beenakker-Mazur method. The results of our renormalized Clausius-Mossotti approximation lead to comparable or much less error (with respect to the numerical simulations) than the Beenakker-Mazur method for the volume fractions below ϕ≈30% (apart from a small range of wave vectors in hydrodynamic function). For volume fractions above ϕ≈30%, the Beenakker-Mazur method gives in most cases lower error than the renormalized Clausius-Mossotti approximation.

Short-time diffusion in concentrated bidisperse hard-sphere suspensions

Diffusion in bidisperse Brownian hard-sphere suspensions is studied by Stokesian Dynamics (SD) computer simulations and a semi-analytical theoretical scheme for colloidal short-time dynamics, based on Beenakker and Mazur's method [Physica A 120, 388-410 (1983); 126, 349-370 (1984)]. Two species of hard spheres are suspended in an overdamped viscous solvent that mediates the salient hydrodynamic interactions among all particles. In a comprehensive parameter scan that covers various packing fractions and suspension compositions, we employ numerically accurate SD simulations to compute the initial diffusive relaxation of density modulations at the Brownian time scale, quantified by the partial hydrodynamic functions. A revised version of Beenakker and Mazur's δγ-scheme for monodisperse suspensions is found to exhibit surprisingly good accuracy, when simple rescaling laws are invoked in its application to mixtures. The so-modified δγ scheme predicts hydrodynamic functions in very good agreement with our SD simulation results, for all densities from the very dilute limit up to packing fractions as high as 40%.

Rotational self-diffusion in suspensions of charged particles: simulations and revised Beenakker-Mazur and pairwise additivity methods

We present a comprehensive joint theory-simulation study of rotational self-diffusion in suspensions of charged particles whose interactions are modeled by the generic hard-sphere plus repulsive Yukawa (HSY) pair potential. Elaborate, high-precision simulation results for the short-time rotational self-diffusion coefficient, D(r), are discussed covering a broad range of fluid-phase state points in the HSY model phase diagram. The salient trends in the behavior of D(r) as a function of reduced potential strength and range, and particle concentration, are systematically explored and physically explained. The simulation results are further used to assess the performance of two semi-analytic theoretical methods for calculating D(r). The first theoretical method is a revised version of the classical Beenakker-Mazur method (BM) adapted to rotational diffusion which includes a highly improved treatment of the salient many-particle hydrodynamic interactions. The second method is an easy-to-implement pairwise additivity (PA) method in which the hydrodynamic interactions are treated on a full two-body level with lubrication corrections included. The static pair correlation functions required as the only input to both theoretical methods are calculated using the accurate Rogers-Young integral equation scheme. While the revised BM method reproduces the general trends of the simulation results, it significantly underestimates D(r). In contrast, the PA method agrees well with the simulation results for D(r) even for intermediately concentrated systems. A simple improvement of the PA method is presented which is applicable for large concentrations.

Dynamics of suspensions of hydrodynamically structured particles: analytic theory and applications to experiments

We present an easy-to-use analytic toolbox for the calculation of short-time transport properties of concentrated suspensions of spherical colloidal particles with internal hydrodynamic structure, and direct interactions described by a hard-core or soft Hertz pair potential. The considered dynamic properties include self-diffusion and sedimentation coefficients, the wavenumber-dependent diffusion function determined in dynamic scattering experiments, and the high-frequency shear viscosity. The toolbox is based on the hydrodynamic radius model (HRM) wherein the internal particle structure is mapped on a hydrodynamic radius parameter for unchanged direct interactions, and on an existing simulation data base for solvent-permeable and spherical annulus particles. Useful scaling relations for the diffusion function and self-diffusion coefficient, known to be valid for hard-core interaction, are shown to apply also for soft pair potentials. We further discuss extensions of the toolbox to long-time transport properties including the low-shear zero-frequency viscosity and the long-time self-diffusion coefficient. The versatility of the toolbox is demonstrated by the analysis of a previous light scattering study of suspensions of non-ionic PNiPAM microgels [Eckert et al., J. Chem. Phys., 2008, 129, 124902] in which a detailed theoretical analysis of the dynamic data was left as an open task. By the comparison with Hertz potential based calculations, we show that the experimental data are consistently and accurately described using the Verlet-Weis corrected Percus-Yevick structure factor as input, and for a solvent penetration length equal to three percent of the excluded volume radius. This small amount of solvent permeability of the microgel particles has a significant dynamic effect at larger concentrations.

Adsorption, organization, and rheology of catanionic layers at the air/water interface

We have investigated the adsorption and organization at the air/water interface of catanionic molecules released from a dispersion of solid-like catanionic vesicles composed of myristic acid and cetyl trimethylammonium chloride at the 2:1 ratio. These vesicles were shown recently to be promising foam stabilizers. Using Brewster angle microscopy, we observed the formation of a catanionic monolayer at the air/water interface composed of liquid-condensed domains in a liquid-expanded matrix. Further adsorption of catanionic molecules forced them to pack, thereby forming a very dense monolayer that prevented further vesicle rupture by avoiding contact of the vesicles with air. Moreover, confocal fluorescence microscopy revealed the presence of layers of intact vesicles that were progressively creaming toward this catanionic monolayer; the surface tension of the vesicle dispersion remained constant upon creaming. The catanionic monolayer behaved as a soft glassy material, an amorphous solid with time- and temperature-dependent properties. Using interfacial oscillatory rheology, we found that the monolayer relaxed mechanical stresses in seconds and melted at a temperature very close to the melting transition temperature of the vesicle bilayers. These results have potential application in the design of smart foams that have temperature-tunable stability.