Activated hopping and dynamical fluctuation effects in hard sphere suspensions and fluids

Single particle Brownian dynamics simulation methods are employed to establish the full trajectory level predictions of our nonlinear stochastic Langevin equation theory of activated hopping dynamics in glassy hard sphere suspensions and fluids. The consequences of thermal noise driven mobility fluctuations associated with the barrier hopping process are determined for various ensemble-averaged properties and their distributions. The predicted mean square displacements show classic signatures of transient trapping and anomalous diffusion on intermediate time and length scales. A crossover to a stronger volume fraction dependence of the apparent nondiffusive exponent occurs when the entropic barrier is of order the thermal energy. The volume fraction dependences of various mean relaxation times and rates can be fitted by empirical critical power laws with parameters consistent with ideal mode-coupling theory. However, the results of our divergence-free theory are largely a consequence of activated dynamics. The experimentally measurable alpha relaxation time is found to be very similar to the theoretically defined mean reaction time for escape from the barrier-dominated regime. Various measures of decoupling have been studied. For fluid states with small or nonexistent barriers, relaxation times obey a simple log-normal distribution, while for high volume fractions the relaxation time distributions become Poissonian. The product of the self-diffusion constant and mean alpha relaxation time increases roughly as a logarithmic function of the alpha relaxation time. The cage scale incoherent dynamic structure factor exhibits nonexponential decay with a modest degree of stretching. A nearly universal collapse of the different volume fraction results occurs if time is scaled by the mean alpha relaxation time. Hence, time-volume fraction superposition holds quite well, despite the presence of stretching and volume fraction dependent decoupling associated with the stochastic barrier hopping process. The relevance of other origins of dynamic heterogeneity (e.g., mesoscopic domains), and comparison of our results with experiments, simulations, and alternative theories, is discussed.     

Theory of nonlinear elasticity, stress-induced relaxation, and dynamic yielding in dense fluids of hard nonspherical colloids

We generalize the microscopic na?ve mode coupling and nonlinear Langevin equation theories of the coupled translation-rotation dynamics of dense suspensions of uniaxial colloids to treat the effect of applied stress on shear elasticity, cooperative cage escape, structural relaxation, and dynamic and static yielding. The key concept is a stress-dependent dynamic free energy surface that quantifies the center-of-mass force and torque on a moving colloid. The consequences of variable particle aspect ratio and volume fraction, and the role of plastic versus double glasses, are established in the context of dense, glass-forming suspensions of hard-core dicolloids. For low aspect ratios, the theory provides a microscopic basis for the recently observed phenomenon of double yielding as a consequence of stress-driven sequential unlocking of caging constraints via reduction of the distinct entropic barriers associated with the rotational and translational degrees of freedom. The existence, and breadth in volume fraction, of the double yielding phenomena is predicted to generally depend on both the degree of particle anisotropy and experimental probing frequency, and as a consequence typically occurs only over a window of (high) volume fractions where there is strong decoupling of rotational and translational activated relaxation. At high enough concentrations, a return to single yielding is predicted. For large aspect ratio dicolloids, rotation and translation are always strongly coupled in the activated barrier hopping event, and hence for all stresses only a single yielding process is predicted.     

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.     

Dynamic yielding, shear thinning, and stress rheology of polymer-particle suspensions and gels

The nonlinear rheological version of our barrier hopping theory for particle-polymer suspensions and gels has been employed to study the effect of steady shear and constant stress on the alpha relaxation time, yielding process, viscosity, and non-Newtonian flow curves. The role of particle volume fraction, polymer-particle size asymmetry ratio, and polymer concentration have been systematically explored. The dynamic yield stress decreases in a polymer-concentration- and volume-fraction-dependent manner that can be described as apparent power laws with effective exponents that monotonically increase with observation time. Stress- or shear-induced thinning of the viscosity becomes more abrupt with increasing magnitude of the quiescent viscosity. Flow curves show an intermediate shear rate dependence of an effective power-law form, becoming more solidlike with increasing depletion attraction. The influence of polymer concentration, particle volume fraction, and polymer-particle size asymmetry ratio on all properties is controlled to a first approximation by how far the system is from the gelation boundary of ideal mode-coupling theory (MCT). This emphasizes the importance of the MCT nonergodicity transition despite its ultimate destruction by activated barrier hopping processes. Comparison of the theoretical results with limited experimental studies is encouraging.     

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. I. Activated relaxation, kinetic vitrification, and fragility

The microscopic nonlinear Langevin equation theory of activated glassy dynamics is applied to dense fluids of spherical particles that interact via a finite range Hertzian contact soft repulsion. The activation barrier and mean alpha relaxation time are predicted to be rich functions of volume fraction and particle stiffness, exhibiting a non-monotonic variation with concentration at high volume fractions. The latter is due to a structural "soft jamming" crossover where the real space local cage order weakens when soft particles significantly overlap. The highly variable dependences of the relaxation time on temperature and volume fraction are reasonably well collapsed onto two distinct master curves that are qualitatively consistent with a recent scaling ansatz and computer simulation study. A kinetic vitrification diagram is constructed and compared to its dynamic crossover analog. Intersection of the dynamic crossover and soft jamming threshold boundaries occurs for particles that are sufficiently soft, implying the nonexistence of a clear activated dynamics regime or kinetic arrest transition for such particles. The isothermal dynamic fragility is predicted to vary over a wide range as a function of particle stiffness, and soft particles behave as strong glasses. Qualitative comparisons with simulations and microgel experiments reveal good agreement.     

The microwave and infrared spectroscopy of benzaldehyde: conflict between theory and experimental deductions

Recently, it has been proposed that ab initio calculations cannot accurately treat molecules comprised of a benzene ring with a pi-conjugated substituent, for example, benzaldehyde. Theoretical predictions of the benzaldehyde barrier to internal rotation are typically a factor of 2 too high in comparison to the experimental values of 4.67 (infared) and 4.90 (microwave) kcal mol(-1). However, both experiments use Pitzer's 1946 model to compute the reduced moment of inertia and employ the experimentally observed torsional frequency to deduce benzaldehyde's rotational barrier. When Pitzer's model is applied to a system with a nonconjugated functional group, such as phenol, the model and theoretical values are in close agreement. Therefore, we conclude the model may not account for conjugation between the substituent and the pi-system of benzene. The experimental values of the benzaldehyde rotational barrier are therefore misleading. The true rotational barrier lies closer to the theoretically extrapolated limit of 7.7 kcal mol(-1), based on coupled cluster theory.     

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.     

Structure, surface excess and effective interactions in polymer nanocomposite melts and concentrated solutions

The Polymer Reference Interaction Site Model (PRISM) theory is employed to investigate structure, effective forces, and thermodynamics in dense polymer-particle mixtures in the one and two particle limit. The influence of particle size, degree of polymerization, and polymer reduced density is established. In the athermal limit, the surface excess is negative implying an entropic dewetting interface. Polymer induced depletion interactions are quantified via the particle-particle pair correlation function and potential of mean force. A transition from (nearly) monotonic decaying, attractive depletion interactions to much stronger repulsive-attractive oscillatory depletion forces occurs at roughly the semidilute-concentrated solution boundary. Under melt conditions, the depletion force is extremely large and attractive at contact, but is proceeded by a high repulsive barrier. For particle diameters larger than roughly five monomer diameters, division of the force by the particle radius results in a nearly universal collapse of the depletion force for all interparticle separations. Molecular dynamics simulations have been employed to determine the depletion force for nanoparticles of a diameter five times the monomer size over a wide range of polymer densities spanning the semidilute, concentrated, and melt regimes. PRISM calculations based on the spatially nonlocal hypernetted chain closure for particle-particle direct correlations capture all the rich features found in the simulations, with quantitative errors for the amplitude of the depletion forces at the level of a factor of 2 or less. The consequences of monomer-particle attractions are briefly explored. Modification of the polymer-particle pair correlations is relatively small, but much larger effects are found for the surface excess including an energetic driven transition to a wetting polymer-particle interface. The particle-particle potential of mean force exhibits multiple qualitatively different behaviors (contact aggregation, steric stabilization, local bridging attraction) depending on the strength and spatial range of the polymer-particle attraction.     

Pressure and temperature dependence of hydrophobic hydration: volumetric, compressibility, and thermodynamic signatures

The combined effect of pressure and temperature on hydrophobic hydration of a nonpolar methanelike solute is investigated by extensive simulations in the TIP4P model of water. Using test-particle insertion techniques, free energies of hydration under a range of pressures from 1 to 3000 atm are computed at eight temperatures ranging from 278.15 to 368.15 K. Corresponding enthalpy, entropy, and heat capacity accompanying the hydration process are estimated from the temperature dependence of the free energies. Partial molar and excess volumes calculated using pressure derivatives of the simulated free energies are consistent with those determined by direct volume simulations; but direct volume determination offers more reliable estimates for compressibility. At 298.15 K, partial molar and excess isothermal compressibilities of methane are negative at 1 atm. Partial molar and excess adiabatic (isentropic) compressibilities are estimated to be also negative under the same conditions. But partial molar and excess isothermal compressibilities are positive at high pressures, with a crossover from negative to positive compressibility at approximately 100-1000 atm. This trend is consistent with experiments on aliphatic amino acids and pressure-unfolded states of proteins. For the range of pressures simulated, hydration heat capacity exhibits little pressure dependence, also in apparent agreement with experiment. When pressure is raised at constant room temperature, hydration free energy increases while its entropic component remains essentially constant. Thus, the increasing unfavorability of hydration under raised pressure is seen as largely an enthalpic effect. Ramifications of the findings of the authors for biopolymer conformational transitions are discussed.     

Theory of correlated two-particle activated glassy dynamics: general formulation and heterogeneous structural relaxation in hard sphere fluids

We generalize the nonlinear Langevin equation theory of activated single particle dynamics to describe the correlated motion of two tagged spherical particles in a glass- or gel-forming fluid as a function of their initial separation. The theory is built on the concept of a two-dimensional dynamic free energy surface which quantifies the forces on two particles moving in a cooperative manner. For the hard sphere fluid, above a threshold volume fraction we generically find two relaxation channels corresponding largely, but not exclusively, to a center-of-mass-like displacement and a radial separation of the two tagged particles. The entropic barriers and mean first passage times are computed and found to systematically vary with volume fraction and initial particle separation; both oscillate as a function of the latter in a manner related to the equilibrium pair correlation function. A dynamic correlation length is estimated as the length scale beyond which the two-particle activated dynamics becomes uncorrelated in space and time, and is found to modestly grow with increasing mean relaxation time. The theory is also applied to a simplified model of cage escape, the elementary step of structural relaxation. Predictions for characteristic relaxation times, translation-relaxation decoupling, and stretched-exponential decay of time correlation functions are obtained. A novel mechanism for understanding why strong decoupling emerges in the activated regime, but stretched nonexponential time correlation functions do not change shape as the mean relaxation time grows, is presented and favorably compared with experiment. The theory may serve as a starting point for constructing a predictive model of multiple correlated caging and hopping (forward and backward) events of a pair of tagged particles.     

Nonlinear elasticity and yielding of nanoparticle glasses

We employ experiment and theory to explore the nonlinear elasticity and yielding of concentrated suspensions of nanoparticles which interact via purely repulsive forces. These glassy suspensions are found to exhibit high exponent power law or simple exponential dependences of the shear elastic modulus and perturbative yield stress on nanoparticle volume fraction, as well as a monotonic decrease of the perturbative yield strain with increasing concentration. Our experimental observations are in good agreement with the predictions of a recently developed microscopic statistical mechanical theory, which describes glassy dynamics based on a nonequilibrium free energy that incorporates local cage correlations and activated barrier hopping processes [(1) Schweizer, K. S.; Saltzman, E. J. J. Chem. Phys. 2003, 119, 1181. (2) Saltzman, E. J.; Schweizer, K. S. J. Chem. Phys. 2003, 119, 1197. (3) Kobelev, V.; Schweizer, K. S. Phy. Rev. E 2005, 71, 021401].     

Re-entrant kinetic arrest and elasticity of concentrated suspensions of spherical and nonspherical repulsive and attractive colloids

We have designed and studied a new experimental colloidal system to probe how the weak shape anisotropy of uniaxial particles and variable repulsive (Coulombic) and attractive (van der Waals) forces influence slow dynamics, shear elasticity, and kinetic vitrification in dense suspensions. The introduction of shape anisotropy dramatically delays kinetic vitrification and reduces the shear elastic modulus of colloidal diatomics relative to their chemically identical spherical analogs. Tuning the interparticle interaction from repulsive, to nearly hard, to attractive by increasing suspension ionic strength reveals a nonmonotonic re-entrant dynamical phase behavior (glass-fluid-gel) and a rich variation of the shear modulus. The experimental results are quantitatively confronted with recent predictions of ideal mode coupling and activated barrier hopping theories of kinetic arrest and elasticity, and good agreement is generally found with a couple of exceptions. The systems created may have interesting materials science applications such as flowable ultrahigh volume fraction suspensions, or responsive fluids that can be reversibly switched between a flowing liquid and a solid nonequilibrium state based on in situ modification of suspension ionic strength.     

Glassy dynamics and mechanical response in dense fluids of soft repulsive spheres. II. Shear modulus, relaxation-elasticity connections, and rheology

We apply the quiescent and mechanically driven versions of nonlinear Langevin equation theory to study how particle softness influences the shear modulus, the connection between shear elasticity and activated relaxation, and nonlinear rheology of the repulsive Hertzian contact model of dense soft sphere fluids. Below the soft jamming threshold, the shear modulus follows a power law dependence on volume fraction over a narrow interval with an apparent exponent that grows with particle stiffness. To a first approximation, the elastic modulus and transient localization length are controlled by a single coupling constant determined by local fluid structure. In contrast to the behavior of hard spheres, an approximately linear relation between the shear modulus and activation barrier is predicted. This connection has recently been observed for microgel suspensions and provides a microscopic realization of the elastic shoving model. Yielding, shear and stress thinning of the alpha relaxation time and viscosity, and flow curves are also studied. Yield strains are relatively weakly dependent on volume fraction and particle stiffness. Shear thinning commences at values of the effective Peclet number far less than unity, a signature of stress-assisted activated relaxation when barriers are high. Apparent power law reduction of the viscosity with shear rate is predicted with a thinning exponent less than unity. In the vicinity of the soft jamming threshold, a power law flow curve occurs over an intermediate reduced shear rate range with an apparent exponent that decreases as fluid volume fraction and/or repulsion strength increase.     

Ergodicity breaking and conformational hysteresis in the dynamics of a polymer tethered at a surface stagnation point

We study the dynamics of long chain polymer molecules tethered to a plane wall and subjected to a stagnation point flow. Using a combination of theory and numerical techniques, including Brownian dynamics (BD), we demonstrate that a chain conformation hysteresis exists even for freely draining (FD) chains. Hydrodynamic interactions (HI) between the polymer and the wall are included in the BD simulations. We find qualitative agreement between the FD and HI simulations, with both exhibiting simultaneous coiled and stretched states for a wide range of fixed flow strengths. The range of state coexistence is understood by considering an equivalent projected equilibrium problem of a two state reaction. Using this formalism, we construct Kramers rate theory (from the inverse mean first passage time for a Markov process) for the hopping transition from coil to stretch and stretch to coil. The activation energy for this rate is found to scale proportionally to chain length or Kuhn step number. Thus, in the limit of infinite chain size the hopping rates at a fixed value of the suitably defined Deborah number approach zero and the states are "frozen." We present the results that demonstrate this "ergodicity breaking."