Tracer-diffusion in binary colloidal hard-sphere suspensions

Binary colloidal glasses: linear viscoelasticity and its link to the microscopic structure and dynamics

We study the relation between the microscopic structure and dynamics and the macroscopic rheological response of glass-forming colloidal suspensions, namely binary colloidal hard-sphere mixtures with large size asymmetry (1 : 5) that span a large range of mixture compositions close to the glass transition. The dynamical shear moduli are measured by oscillatory rheology and the structure and dynamics on the single-particle level by confocal microscopy. The data are compared with Brownian Dynamics simulations and predictions from mode-coupling theory based on the Percus-Yevick approximation. Experiments, simulations and theory consistently observe a strong decrease of the intermediate-frequency mechanical moduli combined with faster dynamics at intermediate mixing ratios and hence a non-monotonic dependence of these parameters but a localization of the large particles which decreases monotonically as the fraction of small particles is increased. We find that the Generalized-Stokes Einstein relation applied to the mean square displacements of the two components leads to a reasonable estimate of the shear moduli of the mixtures and hence links the rheological response to the particle dynamics which in turn reflects the microscopic structure.

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.

Viscosity scaling in concentrated dispersions and its impact on colloidal aggregation

Viscosity scaling in concentrated dispersions is identified using Brownian dynamics simulations, and its impact on colloidal aggregation is quantified.

Diffusion of a sphere in a dilute solution of polymer coils

We calculate the short time and the long time diffusion coefficients of a spherical tracer particle in a polymer solution in the low density limit by solving the Smoluchowski equation for a two-particle system and applying a generalized Einstein relation (fluctuation dissipation theorem). The tracer particle as well as the polymer coils are idealized as hard spheres with a no-slip boundary condition for the solvent but the hydrodynamic radius of the polymer coils is allowed to be smaller than the direct-interaction radius. We take hydrodynamic interactions up to 11th order in the particle distance into account. For the limit of small polymers, the expected generalized Stokes-Einstein relation is found. The long time diffusion coefficient also roughly obeys the generalized Stokes-Einstein relation for larger polymers whereas the short time coefficient does not. We find good qualitative and quantitative agreement to experiments.