Collective diffusion in a two-dimensional Brownian fluid

Coupled dynamics in binary mixtures of model colloidal Yukawa systems

The dynamical behavior of binary mixtures consisting of highly charged colloidal particles is studied by means of Brownian dynamics simulations. We investigate differently sized, but identically charged particles with nearly identical interactions between all species in highly dilute suspensions. Different short-time self-diffusion coefficients induce, mediated by electrostatic interactions, a coupling of both self and collective dynamics of differently sized particles: the long-time self-diffusion coefficients of a larger species are increased by the presence of a more mobile, smaller species and vice versa. Similar coupling effects are observed in collective dynamics where both time constant and functional form of intermediate scattering functions' initial decay are influenced by the presence of a differently sized species. We provide a systematic analysis of coupling effects in dependence on the ratio of sizes, number densities, and the strength of electrostatic interactions.

Dynamic arrest within the self-consistent generalized Langevin equation of colloid dynamics

This paper presents a recently developed theory of colloid dynamics as an alternative approach to the description of phenomena of dynamic arrest in monodisperse colloidal systems. Such theory, referred to as the self-consistent generalized Langevin equation (SCGLE) theory, was devised to describe the tracer and collective diffusion properties of colloidal dispersions in the short- and intermediate-time regimes. Its self-consistent character, however, introduces a nonlinear dynamic feedback, leading to the prediction of dynamic arrest in these systems, similar to that exhibited by the well-established mode coupling theory of the ideal glass transition. The full numerical solution of this self-consistent theory provides in principle a route to the location of the fluid-glass transition in the space of macroscopic parameters of the system, given the interparticle forces (i.e., a nonequilibrium analog of the statistical-thermodynamic prediction of an equilibrium phase diagram). In this paper we focus on the derivation from the same self-consistent theory of the more straightforward route to the location of the fluid-glass transition boundary, consisting of the equation for the nonergodic parameters, whose nonzero values are the signature of the glass state. This allows us to decide if a system, at given macroscopic conditions, is in an ergodic or in a dynamically arrested state, given the microscopic interactions, which enter only through the static structure factor. We present a selection of results that illustrate the concrete application of our theory to model colloidal systems. This involves the comparison of the predictions of our theory with available experimental data for the nonergodic parameters of model dispersions with hard-sphere and with screened Coulomb interactions.

Self-consistent generalized Langevin equation for colloidal mixtures

A self-consistent theory of collective and tracer diffusion in colloidal mixtures is presented. This theory is based on exact results for the partial intermediate scattering functions derived within the framework of the generalized Langevin equation formalism, plus a number of conceptually simple and sensible approximations. The first of these consists of a Vineyard-like approximation between collective and tracer diffusion, which writes the collective dynamics in terms of the memory function related to tracer diffusion. The second consists of interpolating this only unknown memory function between its two exact limits at small and large wave vectors; for this, a phenomenologically determined, but not arbitrary, interpolating function is introduced: a Lorentzian with its inflection point located at the first minimum of the partial static structure factor. The small wave-vector exact limit involves a time-dependent friction function, for which we take a general approximate result, previously derived within the generalized Langevin equation formalism. This general result expresses the time-dependent friction function in terms of the partial intermediate scattering functions, thus closing the system of equations into a fully self-consistent scheme. This extends to mixtures a recently proposed self-consistent theory developed for monodisperse suspensions [Yeomans-Reyna and Medina-Noyola, Phys. Rev. E 64, 066114 (2001)]. As an illustration of its quantitative accuracy, its application to a simple model of a binary dispersion in the absence of hydrodynamic interactions is reported.

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.

Influence of hydrodynamics on many-particle diffusion in 2D colloidal suspensions

We study many-particle diffusion in 2D colloidal suspensions with full hydrodynamic interactions through a novel mesoscopic simulation technique. We focus on the behaviour of the effective scaled tracer and collective-diffusion coefficients DT(rho)/D0 and DC(rho)/D0 respectively, where D0 is the single-particle diffusion coefficient, as a function of the density of the colloids rho. At low Schmidt numbers Sc - 1, we find that hydrodynamics has essentially no effect on the behaviour of DT (rho)/D0. At larger Sc, DT (rho)/D0 seems to be enhanced at all densities, although the differences compared to the case without hydrodynamics are rather minor. The collective-diffusion coefficient, on the other hand, is much more strongly coupled to hydrodynamical conservation laws and is distinctly different from the purely dissipative case without hydrodynamic interactions.

Self-consistent theory of collective Brownian dynamics: theory versus simulation

A recently developed theory of collective diffusion in colloidal suspensions is tested regarding the quantitative accuracy of its description of the dynamics of monodisperse model colloidal systems without hydrodynamic interactions. The idea is to exhibit the isolated effects of the direct interactions, which constitute the main microscopic relaxation mechanism, in the absence of other effects, such as hydrodynamic interactions. Here we compare the numerical solution of the fully self-consistent theory with the results of Brownian dynamics simulation of the van Hove function G(r,t) and/or the intermediate scattering function F(k,t) of four simple model systems. Two of them are representative of short-ranged soft-core repulsive interactions [(sigma/r)(mu), with mu>>1], in two and in three dimensions. The other two involve long-ranged repulsive forces in two (dipolar, r(-3) potential) and in three (screened Coulomb, or repulsive Yukawa interactions) dimensions. We find that the theory, without any sort of adjustable parameters or rescaling prescriptions, provides an excellent approximate description of the collective dynamics of these model systems, particularly in the short- and intermediate-time regimes. We also compare our results with those of the single exponential approximation and with the competing mode-mode coupling theory.

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).

Vineyard-like approximations for colloid dynamics

In this paper we propose a hierarchy of higher-order Vineyard-like approximations for colloidal systems. These consist of approximate expressions for the intermediate scattering function F(k,t) in terms of the self-intermediate scattering function F(s)(k,t) (or some memory function associated with it), and of other static structural properties of the suspension. In order to assess the accuracy of the proposed approximations, we perform Brownian dynamics simulations in a simple model system (a two-dimensional Yukawa Brownian fluid), in which we determine F(k,t), F(s)(k,t), and the required static structural properties. We study proposals for "second-order" and "third-order" Vineyard-like approximations. We find that the detailed structure of the relationship between the corresponding collective and self-memory functions turns out to be most important, as quantified by our simulation results.