Accurate and efficient methods for modeling colloidal mixtures in an explicit solvent using molecular dynamics

Most simulations of colloidal suspensions treat the solvent implicitly or as a continuum. However as particle size decreases to the nanometer scale, this approximation fails and one needs to treat the solvent explicitly. Due to the large number of smaller solvent particles, such simulations are computationally challenging. Additionally, as the ratio of nanoparticle size to solvent size increases, commonly-used molecular dynamics algorithms for neighbor finding and parallel communication become inefficient. Here we present modified algorithms that enable fast single processor performance and reasonable parallel scalability for mixtures with a wide range of particle size ratios. The methods developed are applicable for any system with widely varying force distance cutoffs, independent of particle sizes and independent of the interaction potential. As a demonstration of the new algorithm's effectiveness, we present results for the pair correlation function and diffusion constant for mixtures where colloidal particles interact via integrated potentials. In these systems, with nanoparticles 20 times larger than the surrounding solvent particles, our parallel molecular dynamics code runs more than 100 times faster using the new algorithms.

[1]  R. Vaia,et al.  Influence of polymer matrix composition and architecture on polymer nanocomposite formation: Coarse-grained molecular dynamics simulation , 2003 .

[2]  Ron O Dror,et al.  The midpoint method for parallelization of particle simulations. , 2006, The Journal of chemical physics.

[3]  D. Heyes,et al.  The Newtonian viscosity of concentrated stabilized dispersions: Comparisons with the hard sphere fluid , 2004 .

[4]  P. Strating,et al.  Brownian Dynamics Simulation of a Hard-Sphere Suspension , 1999 .

[5]  L. Verlet.  Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[6]  R. Glowinski,et al.  A fictitious domain approach to the direct numerical simulation of incompressible viscous flow past moving rigid bodies: application to particulate flow , 2001 .

[7]  Sharon C. Glotzer,et al.  Origin of Particle Clustering in a Simulated Polymer Nanocomposite and its Impact on Rheology , 2003 .

[8]  Witold Dzwinel,et al.  A discrete-particle model of blood dynamics in capillary vessels. , 2003, Journal of colloid and interface science.

[9]  A. Malevanets,et al.  Solute molecular dynamics in a mesoscale solvent , 2000 .

[10]  R. C. Ball,et al.  Lubrication breakdown in hydrodynamic simulations of concentrated colloids , 1995 .

[11]  D. Bedrov,et al.  A molecular dynamics simulation study of the viscoelastic properties of polymer nanocomposites , 2002 .

[12]  The Effect of Nanoparticle Shape on Polymer-Nanocomposite Rheology and Tensile Strength , 2007, cond-mat/0701607.

[13]  D. Ermak,et al.  Brownian dynamics with hydrodynamic interactions , 1978 .

[14]  Venkat Ganesan,et al.  A coarse-grained explicit solvent simulation of rheology of colloidal suspensions. , 2005, The Journal of chemical physics.

[15]  P. Español,et al.  FLUID PARTICLE MODEL , 1998 .

[16]  Sharon C. Glotzer,et al.  Molecular dynamics simulation of a polymer melt with a nanoscopic particle , 2002 .

[17]  N. Quirke,et al.  Wetting of nanoparticles and nanoparticle arrays , 2002 .

[18]  R. Gee,et al.  Atomistic simulations of spinodal phase separation preceding polymer crystallization , 2006, Nature materials.

[19]  Pawel Keblinski,et al.  Molecular dynamics simulations of polymer transport in nanocomposites. , 2005, The Journal of chemical physics.

[20]  R Everaers,et al.  Interaction potentials for soft and hard ellipsoids. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  M. Dijkstra,et al.  Phase diagram of highly asymmetric binary hard-sphere mixtures. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[22]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[23]  J. Koelman,et al.  Simulating microscopic hydrodynamic phenomena with dissipative particle dynamics , 1992 .

[24]  Juan J de Pablo,et al.  Fast computation of many-particle hydrodynamic and electrostatic interactions in a confined geometry. , 2007, Physical review letters.

[25]  G. Grest,et al.  Dynamics of entangled linear polymer melts: A molecular‐dynamics simulation , 1990 .

[26]  H. Herrmann,et al.  Simulation of claylike colloids. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  R. Hockney,et al.  Quiet high resolution computer models of a plasma , 1974 .

[28]  Courtenay T. Vaughan,et al.  Parallel Transient Dynamics Simulations: Algorithms for Contact Detection and Smoothed Particle Hydrodynamics , 1998, J. Parallel Distributed Comput..

[29]  J. Brady,et al.  Structure, diffusion and rheology of Brownian suspensions by Stokesian Dynamics simulation , 2000, Journal of Fluid Mechanics.

[30]  Kremer,et al.  Molecular dynamics simulation for polymers in the presence of a heat bath. , 1986, Physical review. A, General physics.

[31]  A. Ladd,et al.  Lattice-Boltzmann Simulations of Particle-Fluid Suspensions , 2001 .

[32]  John F. Brady,et al.  Accelerated Stokesian Dynamics simulations , 2001, Journal of Fluid Mechanics.

[33]  Sivakumar R. Challa,et al.  Molecular simulations of lubrication and solvation forces. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[34]  Robin Ball,et al.  Shear thickening and order–disorder effects in concentrated colloids at high shear rates , 2000 .

[35]  D. Gersappe.  Molecular mechanisms of failure in polymer nanocomposites. , 2002, Physical review letters.

[36]  Peter S. Lomdahl,et al.  MOLECULAR DYNAMICS COMES OF AGE: 320 BILLION ATOM SIMULATION ON BlueGene/L , 2006 .

[37]  J. Padding,et al.  Hydrodynamic interactions and Brownian forces in colloidal suspensions: coarse-graining over time and length scales. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  A. Malevanets,et al.  Mesoscopic model for solvent dynamics , 1999 .

[39]  H. C. Hamaker.  The London—van der Waals attraction between spherical particles , 1937 .

[40]  M. Eleftheriou,et al.  Blue Matter: Approaching the Limits of Concurrency for Classical Molecular Dynamics , 2006, ACM/IEEE SC 2006 Conference (SC'06).