Kinetic arrest, dynamical transitions, and activated relaxation in dense fluids of attractive nonspherical colloids

Medium-Range Structural Order as the Driver of Activated Dynamics and Complexity Reduction in Glass-Forming Liquids

We analyze in depth the Elastically Collective Nonlinear Langevin Equation theory of activated dynamics in metastable liquids to establish that the predicted inter-relationships between the alpha relaxation time, local cage and collective elastic barriers, dynamic localization length, and shear modulus are causally related within the theory to the medium range order (MRO) static correlation length. The latter grows exponentially with density for metastable hard sphere fluids and as a nonuniversal inverse power law with temperature for supercooled liquids under isobaric conditions. The physical origin of predicted connections between the alpha time and other metrics of cage order and the thermodynamic inverse dimensionless compressibility is fully established. It is discovered that although kinetic constraints from the real space first coordination shell are important for the alpha time, they are of secondary importance compared to the consequences of the more universal MRO correlations in both the modestly and deeply metastable regimes. This understanding sheds new light on the theoretical basis for, and prior successes of, the predictive mapping of chemically complex thermal liquids to effective hard sphere fluids based on matching their dimensionless compressibilities, a scheme we call "complexity reduction". In essence, the latter is equivalent to the physical requirement that the thermal liquid MRO correlation equals that of its effective hard sphere analog. The mapping alone is shown to provide a remarkable level of quantitative predictive power for the glass transition temperature Tg of 21 molecular and polymer liquids. Predictions for the chemically specific absolute magnitude and growth with cooling of the MRO correlation length are obtained and lie in the window of 2-6 nm at Tg. Dynamic heterogeneity, elastic facilitation, and beyond pair structure issues are briefly discussed. Future opportunities to theoretically analyze the equilibrated deep glass regime are outlined.

Monte Carlo simulation of cylinders with short-range attractions

Cylindrical or rod-like particles are promising materials for the applications of fillers in nanocomposite materials and additives to control rheological properties of colloidal suspensions. Recent advances in particle synthesis allows for cylinders to be manufactured with short-ranged attractions to study the gelation as a function of packing fraction, aspect ratio and attraction strength. In order to aid in the analysis of small-angle scattering experiments of rod-like particles, computer simulation methods were used to model these particles with specialized Monte Carlo algorithms and tabular superquadric potentials. The attractive interaction between neighboring rods increases with the amount of locally-accessible surface area, thus leading to patchy-like interactions. We characterize the clustering and percolation of cylinders as the attractive interaction increases from the homogenous fluid at relatively low attraction strength, for a variety of aspect ratios and packing fractions. Comparisons with the experimental scattering results are also presented, which are in agreement.

Microscopic Theory for the Role of Attractive Forces in the Dynamics of Supercooled Liquids

We formulate a microscopic, no adjustable parameter, theory of activated relaxation in supercooled liquids directly in terms of the repulsive and attractive forces within the framework of pair correlations. Under isochoric conditions, attractive forces can nonperturbatively modify slow dynamics, but at high enough density their influence vanishes. Under isobaric conditions, attractive forces play a minor role. High temperature apparent Arrhenius behavior and density-temperature scaling are predicted. Our results are consistent with recent isochoric simulations and isobaric experiments on a deeply supercooled molecular liquid. The approach can be generalized to treat colloidal gelation and glass melting, and other soft matter slow dynamics problems.

Theory of kinetic arrest, elasticity, and yielding in dense binary mixtures of rods and spheres

We extend the quiescent and stressed versions of naïve mode coupling theory to treat the dynamical arrest, shear modulus, and absolute yielding of particle mixtures where one or more species is a nonrotating nonspherical object. The theory is applied in detail to dense isotropic "chemically matched" mixtures of variable aspect ratio rods and spheres that interact via repulsive and short range attractive site-site pair potentials. A remarkably rich ideal kinetic arrest behavior is predicted with up to eight "dynamical phases" emerging: an ergodic fluid, partially localized states where the spheres remain fluid but the rods can be a gel, repulsive glass or attractive glass, doubly localized glasses and gels, a porous rod gel plus sphere glass, and a narrow window where a type of rod glass and gel localization coexist. Dynamical complexity increases with rod length and the introduction of attractive forces between all species which both enhance gel network formation. Multiple dynamic reentrant features and triple points are predicted, and each dynamic phase has unique particle localization characteristics and mechanical properties. Orders of magnitude variation of the linear shear modulus and absolute yield stress are found as rod length, mixture composition and the detailed nature of interparticle attractions are varied. The interplay of total (high) mixture packing fraction and composition at fixed temperature is also briefly studied. The present work provides a foundation to study more complex rod-sphere mixtures of both biological and synthetic interest that include physical features such as interaction site size asymmetry, rod-sphere specific attractions, and/or Coulomb repulsion.