Ideal glass transitions, shear modulus, activated dynamics, and yielding in fluids of nonspherical objects

Thermoreversible gels of hollow silica nanorod dispersions

Colloidal suspensions of anisotropic particles are ubiquitous in particle-based industries. Consequently, there is a need to quantify the effects of particle shape on equilibrium phases and kinetic state transitions, particularly at lower aspect ratios (L/D ≈ 1-10). We present a new, colloidal system comprised of hollow, octadecyl-coated silica rods with 40 nm diameter with controlled aspect ratio and thermoreversible short-range attractions. Rheology and dynamic light scattering measurements on suspensions of these hollow adhesive hard rods with nominal aspect ratio ≈3 suspended in tetradecane exhibit thermoreversible gelation without complicating effects of gravitational settling. Small angle neutron scattering measurements of the microstructure are analyzed to determine the effective strength of attraction in the form of Baxter sticky parameter. Quantitative agreement is found with simulation predictions of the thermoreversible gel transition as a function of volume fraction, further validating a universal state diagram and providing guidance for the effects of aspect ratio on gelation.

Penetrant shape effects on activated dynamics and selectivity in polymer melts and networks based on self-consistent cooperative hopping theory

We generalize and apply the microscopic self-consistent cooperative hopping theory for activated penetrant dynamics in polymer melts and crosslinked networks to address the role of highly variable non-spherical molecular shape. The focus is on vastly different shaped penetrants that have identical space filling volume in order to isolate how non-spherical shape explicitly modifies dynamics over a wide range of temperature down to the kinetic glass transition temperature. The theory relates intramolecular and intermolecular structure and kinetic constraints, and reveals how different solvation packing of polymer monomers around variable shaped penetrants impact penetrant hopping. A highly shape-dependent penetrant activated relaxation, including alpha time decoupling and trajectory level cooperativity of the hopping process, is predicted in the deeply supercooled regime for relatively larger penetrants which is sensitive to whether the polymer matrix is a melt or heavily crosslinked network. In contrast, for smaller size penetrants or at high/medium temperatures the effect of isochoric penetrant shape is relatively weak. We propose an aspect ratio variable that organizes how penetrant shape influences the activated relaxation times, leading to a (near) collapse or master curve. The relative absolute values of the penetrant relaxation time (inverse hopping rate) in polymer melts versus in crosslinked networks are found to be opposite when compared at a common absolute temperature versus when they are compared at a fixed value of distance from the glass transition based on the variable Tg/T with Tg the glass transition temperature. Quantitative comparison with recent diffusion experiments on chemically complex molecular penetrants of variable shape but fixed volume in crosslinked networks reveals good agreement, and testable new predictions are made. Extension of the theoretical approach to more complex systems of high experimental interest are discussed, including applications to realize selective transport in membrane separation applications.

Dynamic arrest of adhesive hard rod dispersions

The phenomenon of dynamic arrest, more commonly referred to as gel and glass formation, originates as particle motion slows significantly. Current understanding of gels and glasses stems primarily from dispersions of spherical particles, but much less is known about how particle shape affects dynamic arrest transitions. To better understand the effects of particle shape anisotropy on gel and glass formation, we systematically measure the rheology, particle dynamics, and static microstructure of thermoreversible colloidal dispersions of adhesive hard rods (AHR). First, the dynamic arrest transitions are mapped as a function of temperature T, aspect ratio L/D≈ 3 to 7, and volume fraction φ≈ 0.1 to 0.5. The critical gel temperature T_gel and glass volume fraction φ_g are determined from the particle dynamics and rheology. Second, an effective orientation-averaged, short-range attraction between rods is quantified from small-angle scattering measurements and characterized by a reduced temperature τ. Similar τ is found at low rod concentrations, indicating that rod gelation occurs at similar effective attraction strength independent of L/D. Monte Carlo simulations reveal a similar convergence in τ when rods cluster and percolate with an average bond coordination number 〈n_c〉≈ 2.4, supporting the link between physical gelation and rigidity percolation. Lastly, AHR results are mapped onto a dimensionless state diagram to compare with previous predictions of attraction-driven gels, repulsion-driven glasses, and liquid crystal phases.

Short and Soft: Multidomain Organization, Tunable Dynamics, and Jamming in Suspensions of Grafted Colloidal Cylinders with a Small Aspect Ratio

The yet virtually unexplored class of soft colloidal rods with a small aspect ratio is investigated and shown to exhibit a very rich phase and dynamic behavior, spanning from liquid to nearly melt state. Instead of the nematic order, these short and soft nanocylinders alter their organization with increasing concentration from isotropic liquid with random orientation to small domains with preferred local orientation and eventually a multidomain arrangement with a local orientational order. The latter gives rise to a kinetically suppressed state akin to structural glass with detectable terminal relaxation, which, on further increasing concentration, reveals features of hexagonally packed order as in ordered block copolymers. The respective dynamic response comprises four regimes, all above the overlapping concentration of 0.02 g/mL:(I) from 0.03 to 0.1 g/mol, the system undergoes a liquid-to-solidlike transition with a structural relaxation time that grows by 4 orders of magnitude. (II) From 0.1 to 0.2 g/mL, a dramatic slowing-down is observed and is accompanied by an evolution from isotropic to a multidomain structure. (III) Between 0.2 and 0.6 g/mol, the suspensions exhibit signatures of shell interpenetration and jamming, with the colloidal plateau modulus depending linearly on concentration. (IV) At 0.74 g/mL, in the densely jammed state, the viscoelastic signature of hexagonally packed cylinders from microphase-separated block copolymers is detected. These properties set short and soft nanocylinders apart from long colloidal rods (with a large aspect ratio) and provide insights for fundamentally understanding the physics in this intermediate soft colloidal regime and for tailoring the flow properties of nonspherical soft colloids.

Evolutionary drivers of protein shape

Diffusional motion within the crowded environment of the cell is known to be crucial to cellular function as it drives the interactions of proteins. However, the relationships between protein diffusion, shape and interaction, and the evolutionary selection mechanisms that arise as a consequence, have not been investigated. Here, we study the dynamics of triaxial ellipsoids of equivalent steric volume to proteins at different aspect ratios and volume fractions using a combination of Brownian molecular dynamics and geometric packing. In general, proteins are found to have a shape, approximately Golden in aspect ratio, that give rise to the highest critical volume fraction resisting gelation, corresponding to the fastest long-time self-diffusion in the cell. The ellipsoidal shape also directs random collisions between proteins away from sites that would promote aggregation and loss of function to more rapidly evolving nonsticky regions on the surface, and further provides a greater tolerance to mutation.

Deuterium off-magic-angle spinning NMR spectroscopy for pure-quadrupole spectra of paramagnetic solids

We examined a simple two-dimensional ²H nuclear magnetic resonance spectroscopy for paramagnetic solids utilizing off-magic-angle spinning (OMAS). By adding a rotor-synchronized 180° pulse to rotational echo (RE) measurement, the effect of the shift interaction was removed from the indirect dimension. The obtained pure-quadrupole spectrum could be simulated by calculating a quasi-one-dimensional NMR signal without considering the shift interaction. The sensitivity of the proposed method was compared with that of previous static NMR methods.

Glass transition of repulsive charged rods (fd-viruses)

It has recently been shown that suspensions of long and thin charged fibrous viruses (fd) form a glass at low ionic strengths. The corresponding thick electric double layers give rise to long-ranged repulsive electrostatic interactions, which lead to caging and structural arrest at concentrations far above the isotropic-nematic coexistence region. Structural arrest and freezing of the orientational texture are found to occur at the same concentration. In addition, various types of orientational textures are equilibrated below the glass transition concentration, ranging from a chiral-nematic texture with a large pitch (of about 100 μm), an X-pattern, and a tightly packed domain texture, consisting of helical domains with a relatively small pitch (of about 10 μm) and twisted boundaries. The dynamics of both particles as well as the texture are discussed, below and above the glass transition. Dynamic light scattering correlation functions exhibit two dynamical modes, where the slow mode is attributed to the elasticity of helical domains. On approach of the glass-transition concentration, the slow mode increases in amplitude, while as the amplitudes of the fast and slow mode become equal at the glass transition. Finally, interesting features of the "transient" behaviors of charged fd-rod glass are shown as the initial caging due to structural arrest, the propagation of flow originating from stress release, and the transition to the final metastable glass state. In addition to the intensity correlation function, power spectra are presented as a function of the waiting time, at the zero-frequency limit that may access to the thermal anomalities in a charged system.

Glass transition in suspensions of charged rods: structural arrest and texture dynamics

We report on the observation of a glass transition in suspensions of very long and thin, highly charged colloidal rods (fd-virus particles). Structural particle arrest is found to occur at a low ionic strength due to caging of the charged rods in the potential setup by their neighbors through long-ranged electrostatic interactions. The relaxation time of density fluctuations as probed by dynamic light scattering is found to diverge within a small concentration range. The rod concentration where structural particle arrest occurs is well within the full chiral-nematic state, far beyond the two-phase isotropic-nematic coexistence region. The morphology of the suspensions thus consists of nematic domains with various orientations. We quantify the dynamics of the resulting texture with image-time correlation spectroscopy. Interestingly, the decay times of image correlation functions are found to diverge in a discontinuous fashion at the same concentration of charged rods where structural particle arrest is observed. At the glass-transition concentration, we thus find both structural arrest and freezing of the texture dynamics.

Rheology of dense suspensions of shape anisotropic particles designed to show pH-sensitive anisotropic pair potentials

Here we investigate the flow properties of suspensions of dicolloidal particles composed of interpenetrating spheres where one sphere is rich in polystyrene and the second is rich in poly 2-vinyl pyridine. The synthesis method is designed to create both anisotropic shape and anisotropic interaction potentials that should lead to head to tail clustering. These particles are referred to as copolymer dicolloids (CDCs). The viscoelastic properties of stable and gelled suspensions of CDC particles are compared with analogs composed of homopolymer dicolloids (HDCs), having the same shape but not displaying the anisotropic attractions. After coating the particles with a nonionic surfactant to minimize van der Waals attractions, the flow properties of glassy and gelled suspensions of CDCs and HDCs are studied as a function of volume fraction, ionic strength and pH. Suspensions of HDC particles display a high kinetic arrest volume fraction (φ(g) > 0.5) over a wide range of pH and ionic strength up to [I]=0.5 M, demonstrating that the particles experience repulsive or weakly attractive pair potentials. Suspensions of CDC particles behave in a similar manner at high or low pH when [I]=0.001 M, but gel at a volume fraction of φ(g) < 0.3 and display anomalously large elastic moduli at and above the gel transition point for intermediate pH or for pH=9 when [I]=0.5 M. The gelation processes for the CDC particles are reversible by adjusting the solution pH. Interaction potential anisotropy is evident in the processes, during which the CDC particles yield on increasing oscillatory strain.

Percolation, phase separation, and gelation in fluids and mixtures of spheres and rods

The relationship between kinetic arrest, connectivity percolation, structure and phase separation in protein, nanoparticle, and colloidal suspensions is a rich and complex problem. Using a combination of integral equation theory, connectivity percolation methods, naïve mode coupling theory, and the activated dynamics nonlinear Langevin equation approach, we study this problem for isotropic one-component fluids of spheres and variable aspect ratio rigid rods, and also percolation in rod-sphere mixtures. The key control parameters are interparticle attraction strength and its (short) spatial range, total packing fraction, and mixture composition. For spherical particles, formation of a homogeneous one-phase kinetically stable and percolated physical gel is predicted to be possible, but depends on non-universal factors. On the other hand, the dynamic crossover to activated dynamics and physical bond formation, which signals discrete cluster formation below the percolation threshold, almost always occurs in the one phase region. Rods more easily gel in the homogeneous isotropic regime, but whether a percolation or kinetic arrest boundary is reached first upon increasing interparticle attraction depends sensitively on packing fraction, rod aspect ratio and attraction range. Overall, the connectivity percolation threshold is much more sensitive to attraction range than either the kinetic arrest or phase separation boundaries. Our results appear to be qualitatively consistent with recent experiments on polymer-colloid depletion systems and brush mediated attractive nanoparticle suspensions.

Activated dynamics in dense fluids of attractive nonspherical particles. I. Kinetic crossover, dynamic free energies, and the physical nature of glasses and gels

We apply the center-of-mass versions of naïve mode coupling theory and nonlinear Langevin equation theory to study how short-range attractive interactions modify the onset of localization, activated single-particle dynamics, and the physical nature of the transiently arrested state of a variety of dense nonspherical particle fluids (and the spherical analog) as a function of volume fraction and attraction strength. The form of the dynamic crossover boundary depends on particle shape, but the reentrant glass-fluid-gel phenomenon and the repulsive glass-to-attractive glass crossover always occur. Diverse functional forms of the dynamic free energy are found for all shapes including glasslike, gel-like, a glass-gel form defined by the coexistence of two localization minima and two activation barriers, and a "mixed" attractive glass characterized by a single, very short localization length but an activation barrier located at a large displacement as in repulsive-force caged glasses. For the latter state, particle trajectories are expected to be of a two-step activated form and can be accessed at high attraction strength by increasing volume fraction, or by increasing attraction strength at fixed high enough volume fraction. A new classification scheme for slow dynamics of fluids of dense attractive particles is proposed based on specification of both the nature of the localized state and the particle displacements required to restore ergodicity via activated barrier hopping. The proposed physical picture appears to be in qualitative agreement with recent computer simulations and colloid experiments.

Glass formation and shear elasticity in dense suspensions of repulsive anisotropic particles

Kinetic vitrification, shear elasticity, and the approach to jamming are investigated for repulsive nonspherical colloids and contrasted with their spherical analog. Particle anisotropy dramatically increases the volume fraction for kinetic arrest. The shear modulus of all systems increases roughly exponentially with volume fraction, and a universal collapse is achieved based on either the dynamic crossover or random close packing volume fraction as the key nondimensionalizing quantity. Quantitative comparisons with recent microscopic theories are performed and good agreement demonstrated.

Theory of coupled translational-rotational glassy dynamics in dense fluids of uniaxial particles

The naïve mode coupling theory (NMCT) for ideal kinetic arrest and the nonlinear Langevin equation theory of activated single-particle barrier hopping dynamics are generalized to treat the coupled center-of-mass (CM) translational and rotational motions of uniaxial hard objects in the glassy fluid regime. The key dynamical variables are the time-dependent displacements of the particle center-of-mass and orientational angle. The NMCT predicts a kinetic arrest diagram with three dynamical states: ergodic fluid, plastic glass, and fully nonergodic double glass, the boundaries of which meet at a "triple point" corresponding to a most difficult to vitrify diatomic of aspect ratio approximately 1.43. The relative roles of rotation and translation in determining ideal kinetic arrest are explored by examining three limits of the theory corresponding to nonrotating, pure rotation, and rotationally ergodic models. The ideal kinetic arrest boundaries represent a crossover to activated dynamics described by two coupled stochastic nonlinear Langevin equations for translational and rotational motions. The fundamental quantity is a dynamic free-energy surface, which for small aspect ratios in the high-volume fraction regime exhibits two saddle points reflecting a two-step activated dynamics where relatively rapid rotational dynamics coexists with slower CM translational motions. For large-enough aspect ratios, the dynamic free-energy surface has one saddle point which corresponds to a system-specific coordinated translation-rotation motion. The entropic barriers as a function of the relative amount of rotation versus translation are determined.

Glassy dynamics and kinetic vitrification of isotropic suspensions of hard rods

A nonlinear Langevin equation (NLE) theory for the translational center-of-mass dynamics of hard nonspherical objects has been applied to isotropic fluids of rigid rods. The ideal kinetic glass transition volume fraction is predicted to be a monotonically decreasing function beyond an aspect ratio of two. The functional form of the decrease is weaker than the inverse aspect ratio. Vitrification occurs at lower volume fractions for corrugated tangent bead rods compared to their smooth spherocylinder analogs. The ideal glass transition signals a crossover to activated dynamics, which is estimated to be observable before the nematic phase boundary is encountered if the aspect ratio is less than roughly 25. Calculations of the glassy elastic shear modulus and absolute yield stress reveal a roughly exponential growth with volume fraction. The dependence of entropic barriers and mean barrier hopping times on concentration for rods of variable aspect ratios can be collapsed quite well based on a difference volume fraction variable that quantifies the distance from the ideal glass boundary. Full numerical solution of the NLE theory via stochastic trajectory simulation was performed for tangent bead rods, and the results were compared to their hard sphere analogs. With increasing shape anisotropy the characteristic length scales of the nonequilibrium free energy increase and the magnitude of the localization well and entropic barrier curvatures decreases. These changes result in a significant aspect ratio dependence of dynamical properties and time correlation functions including weaker intermediate time subdiffusive transport, stronger two-step decay of the incoherent dynamic structure factor, longer mean alpha relaxation time, and stronger wavevector-dependent decoupling of relaxation times and the self-diffusion constant. The theoretical results are potentially testable via computer simulation, confocal microscopy, and dynamic light scattering.

Theory of gelation, vitrification, and activated barrier hopping in mixtures of hard and sticky spheres

Naive mode coupling theory (NMCT) and the nonlinear stochastic Langevin equation theory of activated dynamics have been generalized to mixtures of spherical particles. Two types of ideal nonergodicity transitions are predicted corresponding to localization of both, or only one, species. The NMCT transition signals a dynamical crossover to activated barrier hopping dynamics. For binary mixtures of equal diameter hard and attractive spheres, a mixture composition sensitive "glass-melting" type of phenomenon is predicted at high total packing fractions and weak attractions. As the total packing fraction decreases, a transition to partial localization occurs corresponding to the coexistence of a tightly localized sticky species in a gel-like state with a fluid of hard spheres. Complex behavior of the localization lengths and shear moduli exist because of the competition between excluded volume caging forces and attraction-induced physical bond formation between sticky particles. Beyond the NMCT transition, a two-dimensional nonequilibrium free energy surface emerges, which quantifies cooperative activated motions. The barrier locations and heights are sensitive to the relative amplitude of the cooperative displacements of the different species.

Collisions, caging, thermodynamics, and jamming in the barrier hopping theory of glassy hard sphere fluids

An ultralocal limit of the microscopic single particle barrier hopping theory of glassy dynamics is proposed which allows explicit analytic expressions for the characteristic length scales, energy scales, and nonequilibrium free energy to be derived. All properties are shown to be controlled by a single coupling constant determined by the fluid density and contact value of the radial distribution function. This parameter quantifies an effective mean square force exerted on a tagged particle due to collisions with its surroundings. The analysis suggests a conceptual basis for previous surprising findings of multiple inter-relationships between characteristics of the transient localized state, the early stages of cage escape, non-Gaussian or dynamic heterogeneity effects, and the barrier hopping process that defines the alpha relaxation event. The underlying physical picture is also relevant to fluids of nonspherical molecules and sticky colloidal suspensions. The possibility of a unified view of liquid dynamics is suggested spanning the range from dense gases to the zero mobility jammed state.

Ideal vitrification, barrier hopping, and jamming in fluids of modestly anisotropic hard objects

Our recent theory for the glassy dynamics of fluids and suspensions of hard nonspherical objects is applied to several modestly anisotropic shapes. The role of bond length and aspect ratio is studied for diatomics, triatomics, and spherocylinders. As spherical symmetry is broken the ideal kinetic glass transition volume fraction of all objects increases linearly with aspect ratio with the same slope, in surprising agreement with the jamming phase diagram of hard granular ellipsoids. The ideal glass boundary of all shapes is a nonmonotonic function of aspect ratio which is also in qualitative accord with the jamming behavior of spherocylinders and ellipsoids. The maximum glass volume fraction shifts to higher values, and larger aspect ratios, as the object becomes smoother. Suggestions for why the nonequilibrium jamming and kinetic ideal glass formation (dynamical crossover) boundaries are similar are advanced. Beyond the ideal glass volume fraction the nonequilibrium free energy acquires a localization well and entropic barrier. Although its form is highly nonuniversal, if different shapes are compared at constant barrier height then a good collapse is found. Collapse of the volume fraction dependence of the barrier height for different shapes is also predicted for modest shape anisotropy, but increasingly fails as the aspect ratio exceeds 2. For a given volume fraction the mean barrier hopping times are nonmonotonic functions of aspect ratio. The functional form of this dependence, and order of magnitude variation with aspect ratio, is distinct for each object.