Hydrodynamic interactions and Brownian forces in colloidal suspensions: Coarse-graining over time and length scales

Abstract

We describe in detail how to implement a coarse-grained hybrid molecular dynamics and stochastic rotation dynamics simulation technique that captures the combined effects of Brownian and hydrodynamic forces in colloidal suspensions. The importance of carefully tuning the simulation parameters to correctly resolve the multiple time and length scales of this problem is emphasized. We systematically analyze how our coarsegraining scheme resolves dimensionless hydrodynamic numbers such as the Reynolds number Re, which indicates the importance of inertial effects, the Schmidt number Sc, which indicates whether momentum transport is liquidlike or gaslike, the Mach number, which measures compressibility effects, the Knudsen number, which describes the importance of noncontinuum molecular effects, and the Peclet number, which describes the relative effects of convective and diffusive transport. With these dimensionless numbers in the correct regime the many Brownian and hydrodynamic time scales can be telescoped together to maximize computational efficiency while still correctly resolving the physically relevant processes. We also show how to control a number of numerical artifacts, such as finite-size effects and solvent-induced attractive depletion interactions. When all these considerations are properly taken into account, the measured colloidal velocity autocorrelation functions and related self-diffusion and friction coefficients compare quantitatively with theoretical calculations. By contrast, these calculations demonstrate that, notwithstanding its seductive simplicity, the basic Langevin equation does a remarkably poor job of capturing the decay rate of the velocity autocorrelation function in the colloidal regime, strongly underestimating it at short times and strongly overestimating it at long times. Finally, we discuss in detail how to map the parameters of our method onto physical systems and from this extract more general lessons-keeping in mind that there is no such thing as a free lunch-that may be relevant for other coarse-graining schemes such as lattice Boltzmann or dissipative particle dynamics.

Translational and rotational friction on a colloidal rod near a wall

“…In our units these choices mean that the fluid viscosity takes the value η = 2.5m/a 0 t 0 and the kinematic viscosity is ν = 0.5a 2 0 /t 0 . The Schmidt number Sc, which measures the rate of momentum (vorticity) diffusion relative to the rate of mass transfer, is given by Sc = ν/D f ≈ 5, where D f is the fluid self-diffusion constant [26][27][28]. In a gas Sc ∼ 1, momentum is mainly transported by moving particles, whereas in a liquid Sc is much larger and momentum is primarily transported by interparticle collisions.…”

Toward Hydrodynamics with Solvent Free Lipid Models: STRD Martini

“…(More extensive discussions of SRD (25,27) and MPC (28) can be found elsewhere.) These methods represent the fluid with a number of tracer particles whose dynamics conserve mesoscopic hydrodynamic flows.…”

“…At the scales relevant to molecular transport in membranes, water is a low-Reynolds-number, nearly incompressible solvent (34,35). With an appropriate choice of parameters, SRD dynamics capture this hydrodynamic regime (27). Following the approach of Huang et al (36), analysis of the velocity correlations (Fig.…”

Swimming at Low Reynolds Number: From Sheets to the African Trypanosome

“…To couple our swimmers to the MPCD fluid, we include their vertices in the collision step. In addition, for the model trypanosome the fluid particles are not allowed to go through the cell body by applying a stochastic bounce back rule [23,25], which implements the no-slip boundary condition at the surface. However, for our sheet swimmers, we let the fluid particles move through the sheet during the streaming step.…”