Diffusion and sedimentation in colloidal suspensions using multiparticle collision dynamics with a discrete particle model

We study self-diffusion and sedimentation in colloidal suspensions of nearly hard spheres using the multiparticle collision dynamics simulation method for the solvent with a discrete mesh model for the colloidal particles (MD+MPCD). We cover colloid volume fractions from 0.01 to 0.40 and compare the MD+MPCD simulations to experimental data and Brownian dynamics simulations with free-draining hydrodynamics (BD) as well as pairwise far-field hydrodynamics described using the Rotne-Prager-Yamakawa mobility tensor (BD+RPY). The dynamics in MD+MPCD suggest that the colloidal particles are only partially coupled to the solvent at short times. However, the long-time self-diffusion coefficient in MD+MPCD is comparable to that in experiments, and the sedimentation coefficient in MD+MPCD is in good agreement with that in experiments and BD+RPY, suggesting that MD+MPCD gives a reasonable description of hydrodynamic interactions in colloidal suspensions. The discrete-particle MD+MPCD approach is convenient and readily extended to more complex shapes, and we determine the long-time self-diffusion coefficient in suspensions of nearly hard cubes to demonstrate its generality.

Theory of dynamic arrest in colloidal mixtures

We present a first-principles theory of dynamic arrest in colloidal mixtures based on the multicomponent self-consistent generalized Langevin equation theory of colloid dynamics [M. A. Chávez-Rojo and M. Medina-Noyola, Phys. Rev. E 72, 031107 (2005); M. A. Chávez-Rojo and M. Medina-Noyola, Phys. Rev. E76, 039902 (2007)]. We illustrate its application with a description of dynamic arrest in two simple model colloidal mixtures: namely, hard-sphere and repulsive Yukawa binary mixtures. Our results include observation of the two patterns of dynamic arrest, one in which both species become simultaneously arrested and the other involving the sequential arrest of the two species. The latter case gives rise to mixed states in which one species is arrested while the other species remains mobile. We also derive the ("bifurcation" or fixed-point") equations for the nonergodic parameters of the system, which takes the surprisingly simple form of a system of coupled equations for the localization length of the particles of each species. The solution of this system of equations indicates unambiguously which species is arrested (finite localization length) and which species remains ergodic (infinite localization length). As a result, we are able to draw the entire ergodic-nonergodic phase diagram of the binary hard-sphere mixture.

Dynamic equivalence between soft- and hard-core Brownian fluids

In this work, we demonstrate the dynamic equivalence between the members of the family of Brownian fluids whose particles interact through strongly repulsive radially symmetric soft-core potentials. We specifically consider pair potentials proportional to inverse powers of (r/sigma). This equivalence is the dynamic extension of the static equivalence between all these pair potentials and the hard-sphere fluid, assumed in the treatment of soft-core reference potentials in the classical (Weeks-Chandler-Andersen or Barker-Henderson) perturbation theories of simple liquids. In contrast with the strict hard-sphere Brownian system, in the case of soft-sphere potentials the conventional Brownian dynamics algorithm is indeed well defined. We find that, except for small values of nu, and/or very short times, the dynamic properties of all these systems collapse into a single universal curve, upon a well-defined rescaling of the time and distance variables. This family of systems includes the hard-sphere limit. This observation permits a conceptually simple, new, and accurate Brownian dynamics algorithm to simulate the dynamic properties of the hard-sphere model dispersion without hydrodynamic interactions. Such an algorithm consists of the straightforward rescaling of the Brownian-dynamics simulated properties of any of the dynamically equivalent soft-sphere systems.

Diffusion of hard disks and rodlike molecules on surfaces

We study the submonolayer diffusion of hard disks and rodlike molecules on smooth surfaces through numerical simulations and theoretical arguments. We concentrate on the behavior of the various diffusion coefficients as a function of the two-dimensional (2D) number density rho in the case where there are no explicit surface-particle interactions. For the hard disk case, we find that while the tracer diffusion coefficient D(T)(rho) decreases monotonically up to the freezing transition, the collective diffusion coefficient D(C)(rho) is wholly determined by the inverse compressibility which increases rapidly on approaching freezing. We also study memory effects associated with tracer diffusion, and present theoretical estimates of D(T)(rho) from the mode-mode coupling approximation. In the case of rigid rods with short-range repulsion and no orientational ordering, we find behavior very similar to the case of disks with the same repulsive interaction. Both D(T)(rho) and the angular diffusion coefficient D(R)(rho) decrease with rho. Also in this case D(C)(rho) is determined by inverse compressibility and increases rapidly close to freezing. This is in contrast to the case of flexible chainlike molecules in the lattice-gas limit, where D(C)(rho) first increases and then decreases as a function of the density due to the interplay between compressibility and mobility.