Combining molecular dynamics with Lattice Boltzmann: A hybrid method for the simulation of (charged) colloidal systems

Bridging the gap between molecular dynamics and hydrodynamics in nanoscale Brownian motions

Through molecular dynamics simulations, we examined the hydrodynamic behavior of the Brownian motion of fullerene particles based on molecular interactions. The solvation free energy and velocity autocorrelation function (VACF) were calculated by using the Lennard-Jones (LJ) and Weeks-Chandler-Andersen (WCA) potentials for the solute-solvent and solvent-solvent interactions and by changing the size of the fullerene particles. We also measured the diffusion constant of the fullerene particles and the shear viscosity of the host fluid, and then the hydrodynamic radius aHD was quantified from the Stokes-Einstein relation. The aHD value exceeds that of the gyration radius of the fullerene when the solvation free energy exhibits largely negative values using the LJ potential. In contrast, aHD is similar to the size of bare fullerene when the solvation free energy is positive using the WCA potential. Furthermore, the VACF of the fullerene particles is directly comparable with the analytical expressions utilizing the Navier-Stokes equations both in incompressible and compressible forms. A hydrodynamic long-time tail t-3/2 is demonstrated for timescales longer than the kinematic time of the momentum diffusion over the particle size. However, the VACF at shorter timescales deviates from the hydrodynamic description, particularly for smaller fullerene particles and for the LJ potential. This occurs even though the compressible effect is considered when characterizing the decay of the VACF around the sound-propagation timescale over the particle size. These results indicate that the nanoscale Brownian motion is influenced by the solvation structure around the solute particles originating from the molecular interaction.

Physical foundation of the fluid particle dynamics method for colloid dynamics simulation

Colloid dynamics is significantly influenced by many-body hydrodynamic interactions mediated by a suspending fluid. However, theoretical and numerical treatments of such interactions are extremely difficult. To overcome this situation, we developed a fluid particle dynamics (FPD) method [H. Tanaka and T. Araki, Phys. Rev. Lett., 2000, 35, 3523], which is based on two key approximations: (i) a colloidal particle is treated as a highly viscous particle and (ii) the viscosity profile is described by a smooth interfacial profile function. Approximation (i) makes our method free from the solid-fluid boundary condition, significantly simplifying the treatment of many-body hydrodynamic interactions while satisfying the incompressible condition without the Stokes approximation. Approximation (ii) allows us to incorporate an extra degree of freedom in a fluid, e.g., orientational order and concentration, as an additional field variable. Here, we consider two fundamental problems associated with these approximations. One is the introduction of thermal noise and the other is the incorporation of coupling of the colloid surface with an order parameter introduced into a fluid component, which is crucial when considering colloidal particles suspended in a complex fluid. Here, we show that our FPD method makes it possible to simulate colloid dynamics properly while including full hydrodynamic interactions, inertia effects, incompressibility, thermal noise, and additional degrees of freedom of a fluid, which may be relevant for wide applications in colloidal and soft matter science.

Simulation of phase separation with temperature-dependent viscosity using lattice Boltzmann method

This paper presents an exploration of the phase separation behavior and pattern formation in a binary fluid with temperature-dependent viscosity via a coupled lattice Boltzmann method (LBM). By introducing a viscosity-temperature relation into the LBM, the coupling effects of the viscosity-temperature coefficient [Formula: see text] , initial viscosity [Formula: see text] and thermal diffusion coefficient [Formula: see text] , on the phase separation were successfully described. The calculated results indicated that an increase in initial viscosity and viscosity-temperature coefficient, or a decrease in the thermal diffusion coefficient, can lead to the orientation of isotropic growth fronts over a wide range of viscosity. The results showed that droplet-type phase structures and lamellar phase structures with domain orientation parallel or perpendicular to the walls can be obtained in equilibrium by controlling the initial viscosity, thermal diffusivity, and the viscosity-temperature coefficient. Furthermore, the dataset was rearranged for growth kinetics of domain growth and thermal diffusion fronts in a plot by the spherically averaged structure factor and the ratio of separated and continuous phases. The analysis revealed two different temporal regimes: spinodal decomposition and domain growth stages, which further quantified the coupled effects of temperature and viscosity on the evolution of temperature-dependent phase separation. These numerical results provide guidance for setting optimum temperature ranges to obtain expected phase separation structures for systems with temperature-dependent viscosity.

The raspberry model for hydrodynamic interactions revisited. I. Periodic arrays of spheres and dumbbells

The so-called "raspberry" model refers to the hybrid lattice-Boltzmann and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally developed by Lobaskin and Dünweg [New J. Phys. 6, 54 (2004)], wherein discrete surface points are used to achieve fluid-particle coupling. This technique has been used in many simulation studies on the behavior of colloids. However, there are fundamental questions with regards to the use of this model. In this paper, we examine the accuracy with which the raspberry method is able to reproduce Stokes-level hydrodynamic interactions when compared to analytic expressions for solid spheres in simple-cubic crystals. To this end, we consider the quality of numerical experiments that are traditionally used to establish these properties and we discuss their shortcomings. We show that there is a discrepancy between the translational and rotational mobility reproduced by the simple raspberry model and present a way to numerically remedy this problem by adding internal coupling points. Finally, we examine a non-convex shape, namely, a colloidal dumbbell, and show that the filled raspberry model replicates the desired hydrodynamic behavior in bulk for this more complicated shape. Our investigation is continued in de Graaf et al. [J. Chem. Phys. 143, 084108 (2015)], wherein we consider the raspberry model in the confining geometry of two parallel plates.

Effects of structural properties of the Stern layer on the electrophoretic migration of a highly charged spherical macroion

The electrophoretic migration of a highly charged spherical macroion suspended in an aqueous solution of NaCl is studied using the molecular dynamic method. The objective is to examine the effects of the colloidal surface charge density on the electrophoretic mobility (μ) of the spherical macroion. The bare charge and the size of the macroion are varied separately to induce changes in the colloidal surface charge density. Our results indicate that μ depends on colloidal surface charge density in a nonmonotonic manner, but that this relationship is independent of the way the surface charge density is varied. It is found that an increase in colloidal surface charge density may lead to the formation of new sublayers in the Stern layer. The μ profile is also found to have a local maximum for a bare charge at which a new sublayer is formed in the Stern layer, and a local minimum for a bare charge at which the outer sublayer becomes relatively dense. Finally, the electrophoretic flow caused by the migration of the spherical macroion is studied to find that one decisive factor causing the electrophoretic flow is the ability of the macroion to carry anions in the electrolyte solution.

Computer simulations of single particles in external electric fields

Applying electric fields is an attractive way to control and manipulate single particles or molecules, e.g., in lab-on-a-chip devices. However, the response of nanosize objects in electrolyte solution to external fields is far from trivial. It is the result of a variety of dynamical processes taking place in the ion cloud surrounding charged particles and in the bulk electrolyte, and it is governed by an intricate interplay of electrostatic and hydrodynamic interactions. Already systems composed of one single particle in electrolyte solution exhibit a complex dynamical behaviour. In this review, we discuss recent coarse-grained simulations that have been performed to obtain a molecular-level understanding of the dynamic and dielectric response of single particles and single macromolecules to external electric fields. We address both the response of charged particles to constant fields (DC fields), which can be characterized by an electrophoretic mobility, and the dielectric response of both uncharged and charged particles to alternating fields (AC fields), which is described by a complex polarizability. Furthermore, we give a brief survey of simulation algorithms and highlight some recent developments.

AC-field-induced polarization for uncharged colloids in salt solution: a dissipative particle dynamics simulation

We study the response of a spherical colloid under alternating electric fields (AC fields) by mesoscopic simulation method, accounting in full for hydrodynamic and electrostatic interactions. We focus on a special case of uncharged colloids. The main polarization mechanism is the "volume polarization", where ionic currents are deflected by the core of the uncharged colloid. Specifically, we compute the polarizability of a single colloid and systematically investigate the effect of AC field frequency and salt concentration. The simulation results are compared with predictions from classical Maxwell-Wagner theory and electrokinetic theory.

Dielectric response of nanoscopic spherical colloids in alternating electric fields: a dissipative particle dynamics simulation

We study the response of single nanosized spherical colloids in electrolyte solution to an alternating electric field (AC field) by computer simulations. We use a coarse-grained mesoscopic simulation approach that accounts in full for hydrodynamic and electrostatic interactions as well as for thermal fluctuations. The solvent is modeled as a fluid of single dissipative particle dynamics (DPD) beads, and the colloidal particle is modeled as a rigid body made of DPD beads. We compute the mobility and the polarizability of a single colloid and investigate systematically the effect of amplitude and frequency of the AC fields. Even though the thickness of the Debye layer is not 'thin' compared to the radius of the colloid, and the thermal fluctuations are significant, the results are in good agreement with the theoretical prediction of the Maxwell-Wagner-O'Konski theory, especially for uncharged colloids.

The role of effective charges in the electrophoresis of highly charged colloids

We study the variation of electrophoretic mobility μ of highly charged spherical colloidal macroions for varying surface charge density σ on the colloid using computer simulations of the primitive model for charged colloids. Hydrodynamic interactions between ions are incorporated by coupling the primitive model of charged colloids to the lattice Boltzmann model (LB) of the fluid. In the highly charged regime, the mobility μ of the colloid is known to decrease with the increase of bare charge Q of the colloid; the aim of this paper is to investigate the cause of this. We have identified that the two main factors contributing to the decrease of μ are counterion charge condensation on the highly charged colloid and an increase in effective friction of the macroion-counterion complex due to the condensed counterions. Thus the established O'Brien and White theory, which identified the dipolar force originating from distortion of the electric double layer as the cause of decreasing μ, seems to break down for the case of highly charged colloids with σ in the range of 30-400 µC cm (- 2). To arrive at our conclusions, we counted the number of counterions q0 moving along with the spherical macroion. We observe in our simulations that q0 increases with the increase of bare charge Q, such that the effective charge Qeff = Q - q0 remains approximately constant. Interestingly for our nanometer-sized charged colloid, we observe that, if surface charge density σ of the colloid is increased by decreasing the radius RM of the colloid but fixed bare charge Q, the effective charge Q - q0 decreases with the increase of σ. This behavior is qualitatively different when σ is increased by increasing Q keeping RM fixed. Our observations address a controversy about the effective charge of a strongly charged macroion: some studies claim that effective charge is independent of the bare charge (Alexander et al 1984 J. Chem. Phys. 80 5776; Trizac et al 2003 Langmuir 19 4027) whereas others claim that Qeff decreases with increasing bare charge Q (dos Santos 2009 J. Chem. Phys. 130 124110; Diehl et al 2004 J. Chem. Phys. 121 12100; Groot et al 1991 J. Chem. Phys. 95 9191) at relatively high values of σ.

Hydrodynamics of confined colloidal fluids in two dimensions

We apply a hybrid molecular dynamics and mesoscopic simulation technique to study the dynamics of two-dimensional colloidal disks in confined geometries. We calculate the velocity autocorrelation functions and observe the predicted t;{-1} long-time hydrodynamic tail that characterizes unconfined fluids, as well as more complex oscillating behavior and negative tails for strongly confined geometries. Because the t;{-1} tail of the velocity autocorrelation function is cut off for longer times in finite systems, the related diffusion coefficient does not diverge but instead depends logarithmically on the overall size of the system. The Langevin equation gives a poor approximation to the velocity autocorrelation function at both short and long times.

Simulating (electro)hydrodynamic effects in colloidal dispersions: smoothed profile method

Previously, we have proposed a direct simulation scheme for colloidal dispersions in a Newtonian solvent (Phys. Rev. E 71, 036707 (2005)). An improved formulation called the "Smoothed Profile (SP) method" is presented here in which simultaneous time-marching is used for the host fluid and colloids. The SP method is a direct numerical simulation of particulate flows and provides a coupling scheme between the continuum fluid dynamics and rigid-body dynamics through utilization of a smoothed profile for the colloidal particles. Moreover, the improved formulation includes an extension to incorporate multicomponent fluids, allowing systems such as charged colloids in electrolyte solutions to be studied. The dynamics of the colloidal dispersions are solved with the same computational cost as required for solving non-particulate flows. Numerical results which assess the hydrodynamic interactions of colloidal dispersions are presented to validate the SP method. The SP method is not restricted to particular constitutive models of the host fluids and can hence be applied to colloidal dispersions in complex fluids.

Influence of hydrodynamic interactions on lane formation in oppositely charged driven colloids

The influence of hydrodynamic interactions on lane formation of oppositely charged driven colloidal suspensions is investigated using Brownian dynamics computer simulations performed on the Rotne-Prager level of the mobility tensor. Two cases are considered, namely sedimentation and electrophoresis. In the latter case the Oseen contribution to the mobility tensor is screened due to the opposite motion of counterions. The simulation results are compared to that resulting from simple Brownian dynamics where hydrodynamic interactions are neglected. For sedimentation, we find that hydrodynamic interactions strongly disfavor laning. In the steady state of lanes, a macroscopic phase separation of lanes is observed. This is in marked contrast to the simple Brownian case where a finite size of lanes was obtained in the steady state. For strong Coulomb interactions between the colloidal particles a lateral square lattice of oppositely driven lanes is stable similar to the simple Brownian dynamics. In an electric field, on the other hand, the behavior is found in qualitative and quantitative accordance with the case of neglected hydrodynamics.

Lattice-Boltzmann simulation of the sedimentation of charged disks

We report a series of lattice-Boltzmann simulations of the sedimentation velocity of charged disks. In these simulations, we explicitly account for the hydrodynamic and electrostatic forces on disks and on their electrical double layer. By comparing our results with those for spheres with equal surface and charge, we can clarify the effect of the particle shape on the sedimentation process. We find that disks and spheres exhibit a different dependence of the sedimentation velocity on the Debye screening length. An analysis of the behavior of highly charged disks (beyond the scope of the linearized Poisson-Boltzmann equation) shows that, in that regime, the charge dependence of the sedimentation velocity of disks and spheres is similar. This suggests that, at high charge, the effective hydrodynamic shape of the disks becomes more spherical.

Electrophoretic properties of highly charged colloids: A hybrid molecular dynamics∕lattice Boltzmann simulation study

Using computer simulations, the electrophoretic motion of a positively charged colloid (macroion) in an electrolyte solution is studied in the framework of the primitive model. In this model, the electrolyte is considered as a system of negatively and positively charged microions (counterions and coions, respectively) that are immersed into a structureless medium. Hydrodynamic interactions are fully taken into account by applying a hybrid simulation scheme, where the charged ions (i.e., macroion and electrolyte), propagated via molecular dynamics, are coupled to a lattice Boltzmann (LB) fluid. In a recent electrophoretic experiment by Martin-Molina et al. [J. Phys. Chem. B 106, 6881 (2002)], it was shown that, for multivalent salt ions, the mobility mu initially increases with charge density sigma, reaches a maximum, and then decreases with further increase of sigma. The aim of the present work is to elucidate the behavior of mu at high values of sigma. Even for the case of monovalent microions, a decrease of mu with sigma is found. A dynamic Stern layer is defined that includes all the counterions that move with the macroion while subjected to an external electrical field. The number of counterions in the Stern layer, q(0), is a crucial parameter for the behavior of mu at high values of sigma. In this case, the mobility mu depends primarily on the ratio q(0)/Q (with Q the valency of the macroion). The previous contention that the increase in the distortion of the electric double layer (EDL) with increasing sigma leads to the lowering of mu does not hold for high sigma. In fact, it is shown that the deformation of the EDL decreases with the increase of sigma. The role of hydrodynamic interactions is inferred from direct comparisons to Langevin simulations where the coupling to the LB fluid is switched off. Moreover, systems with divalent counterions are considered. In this case, at high values of sigma the phenomenon of charge inversion is found.

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

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 coarse-graining 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.

Direct numerical simulations of electrophoresis of charged colloids

We propose a numerical method to simulate electrohydrodynamic phenomena in charged colloidal dispersions. This method enables us to compute the time evolutions of colloidal particles, ions, and host fluids simultaneously by solving Newton, advection-diffusion, and Navier-Stokes equations so that the electrohydrodynamic couplings can be fully taken into account. The electrophoretic mobilities of charged spherical particles are calculated in several situations. The comparisons with approximation theories show quantitative agreements for dilute dispersions without any empirical parameters; however, our simulation predicts notable deviations in the case of dense dispersions.