Search

A dissipative particle dynamics (DPD) model is developed and demonstrated for studying dynamics in colloidal rod suspensions. The solvent is modeled as conventional DPD particles, while individual rods are represented by a rigid linear chain consisting of overlapping solid spheres, which interact with solvent particles through a hard repulsive potential. The boundary condition on the rod surface is controlled using a surface friction between the solid spheres and the solvent particles. In this work, this model is employed to study the diffusion of a single colloid in the DPD solvent and compared with theoretical predictions. Both the translational and rotational diffusion coefficients obtained at a proper surface friction show good agreement with calculations based on the rod size defined by the hard repulsive potential. In addition, the system-size dependence of the diffusion coefficients shows that the Navier–Stokes hydrodynamic interactions are correctly included in this DPD model. Comparing our results with experimental measurements of the diffusion coefficients of gold nanorods, we discuss the ability of the model to correctly describe dynamics in real nanorod suspensions. Our results provide a clear reference point from which the model could be extended to enable the study of colloid dynamics in more complex situations or for other types of particles.

Here, we present an implicit-solvent model for dynamic simulations of hard-rod and rod-polymer suspensions. Individual rods are represented by a rigid linear chain consisting of overlapping spheres which interact through a pseudohard-core potential based on the cut-and-shifted Mie (generalized Lennard-Jones) potential with exponents (50, 49). In the rod-polymer suspensions, the polymers are modeled as freely interpenetrable spheres with respect to each other, while there is the pseudohard-core repulsion between the polymer and rod spheres. Dynamic simulations with this model are carried out with a dissipative particle dynamics (DPD) thermostat-each sphere is put in a larger DPD sphere and thus interacts with others via additional pairwise frictional and random forces-which captures the effects of Brownian forces due to the solvent while conserving local momentum. The phase behavior of these models, obtained from continuous compression and expansion simulations, reproduces previous predictions based on theoretical calculations and Monte Carlo simulations. Our method is suited to study dynamic processes in these suspensions, including nucleation and self-assembly, and can be readily extended to colloidal particles of different shapes and chemistry.

In order to better control the assembly of nanorods, knowledge of the pathways by which they form ordered structures is desirable. In this paper, we characterize crystal nucleation in suspensions of spherocylindrical rods with aspect ratio L/D = 2.3 in the presence of both small and large polymer depletants. Using a combination of Langevin dynamics and Monte Carlo simulations, together with biased sampling techniques, we show that the preferred pathway always involves the formation of monolayer assemblies irrespective of the volume fraction of the initial isotropic phase and the diameter of the depletants. This includes the previously neglected case of nucleation from the colloidal liquid phase and shows that the presence of depletion attraction can alter nucleation pathways even when the initial phase is dense.

While growth and dissolution of surface nanobubbles has been widely studied in recent years, their stability under pressure changes or a temperature increase has not received the same level of scrutiny. Here, we present theoretical predictions based on classical theory for pressure and temperature thresholds ($p_c$ and $T_c$) at which unstable growth occurs for the case of air nanobubbles on a solid surface in water. We show that bubbles subjected to pinning have much lower $p_c$ and higher $T_c$ compared to both unpinned and bulk bubbles of similar size, indicating that pinned bubbles can withstand a larger tensile stress (negative pressure) and higher temperatures. The values of $p_c$ and $T_c$ obtained from many-body dissipative particle dynamics (MDPD) simulations of quasi-two-dimensional (quasi-2D) surface nanobubbles are consistent with the theoretical predictions, provided that the lateral expansion during growth is taken into account. This suggests that the modified classical thermodynamic description is valid for pinned bubbles as small as several nanometers. While some discrepancies still exist between our theoretical results and previous experiments, further experimental data is needed before a comprehensive understanding of the stability of surface nanobubbles can be achieved., Comment: 12 pages, 9 figures

This paper describes the mechanism of defect-mediated relaxation in a dodecagonal square-triangle random tiling phase exhibited by a simulated binary mixture of soft discs in 2D. We examine the internal transitions within the elementary mobile defect (christened the 'zipper') that allow it to move, as well as the mechanisms by which the zipper is created and annihilated. The structural relaxation of the random tiling phase is quantified and we show that this relaxation is well described by a model based on the distribution of waiting times for each atom to be visited by the diffusing zipper. This system, representing one of the few instances where a well defined mobile defect is capable of structural relaxation, can provide a valuable test case for general theories of relaxation in complex and disordered materials., Comment: 30 pages, 16 figures, to appear in J. Chem. Phys

The use of the isoconfigurational ensemble to explore structure-dynamic correlations in supercooled liquids is examined. The statistical error of the dynamic propensity and its spatial distribution are determined. The authors present the spatial distribution of the particle non-Gaussian parameter as a measure of the intermittency with which particles exhibit their propensity for motion. The ensemble average of the direction of particle motion is introduced to establish the anisotropy of the dynamic propensity., Comment: Published - see below or http://link.aip.org/link/?JCPSA6/126/154503/1