Long-time self-diffusion coefficients of suspensions

Diffusion, density, and defects on spheres

We simulate and model diffusion of spherical colloids of radius, a, on spherical surfaces of radius, R, as a function of relative size and surface concentration. Using Brownian dynamics simulations, we quantify diffusion and microstructure at different concentrations ranging from single particles to dense crystalline states. Self-diffusion and structural metrics (pair distribution, local density, and topological charge) are indistinguishable between spheres and planes for all concentrations up to dense liquid states. For concentrations approaching and greater than the freezing transition, smaller spheres with higher curvature show increased diffusivities and nonuniform density/topological defect distributions, which differ qualitatively from planar surfaces. The total topological charge varies quadratically with sphere radius for dense liquid states and linearly with sphere radius for dense crystals with icosahedrally organized grain scars. Between the dense liquid and dense crystal states on spherical surfaces is a regime of fluctuating and interacting defect clusters. We show local density governs self-diffusion in dense liquids on flat and spherical surfaces via the pair distribution. In contrast, dynamic topological defects couple to finite diffusivities through freezing and in low density crystal states on spherical surfaces, where neither exist on flat surfaces.

Dynamics of wire frame glasses in two dimensions

The dynamics of wire frame particles in concentrated suspension are studied by means of a 2D model and compared to those of rod-like particles. The wire frames have bent or branched structures constructed from infinitely thin, rigid rods. In the model, a particle is surrounded by diffusing points that it cannot cross. We derive a formal expression for the mean squared displacement (MSD) and, by using a self-consistent approximation, we find markedly different dynamics for wire frames and rods. For wire frames, there exists a critical concentration of points above which they become frozen with the long time MSD reaching a plateau. Rods, on the other hand, always diffuse by reptation. We also study the rheology through the elastic stress, and more striking differences are found: the initial magnitude of the stress for wire frames is much larger than for rods, scaling such as the square of the point concentration, and above the critical concentration, the stress for wire frames appears to persist indefinitely while for rods it always decays.

Geometric percolation of hard-sphere dispersions in shear flow

We combine a heuristic theory of geometric percolation and the Smoluchowski theory of colloid dynamics to predict the impact of shear flow on the percolation threshold of hard spherical colloidal particles, and verify our findings by means of molecular dynamics simulations. It appears that the impact of shear flow is subtle and highly non-trivial, even in the absence of hydrodynamic interactions between the particles. The presence of shear flow can both increase and decrease the percolation threshold, depending on the criterion used for determining whether or not two particles are connected and on the Péclet number. Our approach opens up a route to quantitatively predict the percolation threshold in nanocomposite materials that, as a rule, are produced under non-equilibrium conditions, making comparison with equilibrium percolation theory tenuous. Our theory can be adapted straightforwardly for application in other types of flow field, and particles of different shape or interacting via other than hard-core potentials.

An alternative, dynamic density functional-like theory for time-dependent density fluctuations in glass-forming fluids

We propose an alternative theory for the relaxation of density fluctuations in glass-forming fluids. We derive an equation of motion for the density correlation function that is local in time and is similar in spirit to the equation of motion for the average non-uniform density profile derived within the dynamic density functional theory. We identify the Franz-Parisi free energy functional as the non-equilibrium free energy for the evolution of the density correlation function. An appearance of a local minimum of this functional leads to a dynamic arrest. Thus, the ergodicity breaking transition predicted by our theory coincides with the dynamic transition of the static approach based on the same non-equilibrium free energy functional.

Dynamics of interacting Brownian particles: A diagrammatic formulation

We present a diagrammatic formulation of a theory for the time dependence of density fluctuations in equilibrium systems of interacting Brownian particles. To facilitate derivation of the diagrammatic expansion, we introduce a basis that consists of orthogonalized many-particle density fluctuations. We obtain an exact hierarchy of equations of motion for time-dependent correlations of orthogonalized density fluctuations. To simplify this hierarchy we neglect contributions to the vertices from higher-order cluster expansion terms. An iterative solution of the resulting equations can be represented by diagrams with three- and four-leg vertices. We analyze the structure of the diagrammatic series for the time-dependent density correlation function and obtain a diagrammatic interpretation of reducible and irreducible memory functions. The one-loop self-consistent approximation for the latter function coincides with mode-coupling approximation for Brownian systems that was derived previously using a projection operator approach.

Self-diffusion in a system of interacting Langevin particles

The behavior of the self-diffusion constant of Langevin particles interacting via a pairwise interaction is considered. The diffusion constant is calculated approximately within a perturbation theory in the potential strength about the bare diffusion constant. It is shown how this expansion leads to a systematic double expansion in the inverse temperature beta and the particle density rho. The one-loop diagrams in this expansion can be summed exactly and we show that this result is exact in the limit of small beta and rhobeta constants. The one-loop result can also be resummed using a semiphenomenological renormalization group method which has proved useful in the study of diffusion in random media. In certain cases the renormalization group calculation predicts the existence of a diverging relaxation time signaled by the vanishing of the diffusion constant, possible forms of divergence coming from this approximation are discussed. Finally, at a more quantitative level, the results are compared with numerical simulations, in two dimensions, of particles interacting via a soft potential recently used to model the interaction between coiled polymers.

Self-consistent generalized Langevin equation for colloid dynamics

We present a general self-consistent theory of colloid dynamics which, for a system without hydrodynamic interactions, allows us to calculate F(k,t), and its self-diffusion counterpart F(S)(k,t), given the effective interaction pair potential u(r) between colloidal particles, and the corresponding equilibrium static structural properties. This theory is build upon the exact results for F(k,t) and F(S)(k,t) in terms of a hierarchy of memory functions, derived from the application of the generalized Langevin equation formalism, plus the proposal of Vineyard-like connections between F(k,t) and F(S)(k,t) through their respective memory functions, and a closure relation between these memory functions and the time-dependent friction function Delta zeta(t). As an illustrative application, we present and analyze a selection of numerical results of this theory in the short- and intermediate-time regimes, as applied to a two-dimensional repulsive Yukawa Brownian fluid. For this system, we find that our theory accurately describes the dynamic properties contained in F(k,t) in a wide range of conditions, including strongly correlated systems, at the longest times available from our computer simulations.

Overdamped van hove function of colloidal suspensions

The generalized-hydrodynamic theory for collective diffusion of a monodisperse colloidal suspension is developed in the framework of the Onsager-Machlup theory of time-dependent fluctuations. The time evolution of the intermediate scattering function F(k,t) is derived as a contraction of the description involving the instantaneous particle number concentration, the particle current, and the stress tensor of the Brownian fluid as state variables. We show that the proper overdamped limit of this equation requires the explicit separation of the stress tensor in its mutually orthogonal kinetic and configurational contributions. Analogous results also follow for the self-intermediate scattering function F(s)(k,t). We show that neglecting the non-Markovian part of the configurational stress tensor memory, one recovers the single exponential memory approximation (based on sum rules derived from the Smoluchowski equation) for both F(s)(k,t) and F(k,t). We suggest simple approximate manners to relate the collective and the self-memory functions, leading to Vineyard-like approximate relations between F(s)(k,t) and F(k,t).