Large-amplitude jumps and non-Gaussian dynamics in highly concentrated hard sphere fluids

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.

Fluctuating diffusivity emerges even in binary gas mixtures

Diffusivity in some soft matter and biological systems changes with time, called the fluctuating diffusivity. In this work, we propose a novel origin for fluctuating diffusivity based on stochastic simulations of binary gas mixtures. In this system, the fraction of one component is significantly small, and the mass of the minor component molecule is different from that of the major component. The minor component exhibits fluctuating diffusivity when its mass is sufficiently smaller than that of the major component. We elucidate that this fluctuating diffusivity is caused by the time scale separation between the relaxation of the velocity direction and the speed of the minor component molecule.

Efficient d-dimensional molecular dynamics simulations for studies of the glass-jamming transition

We develop an algorithm suitable for parallel molecular dynamics simulations in d spatial dimensions and describe its implementation in C++. All routines work in arbitrary d; the maximum simulated d is limited only by available computing resources. These routines include several that are particularly useful for studies of the glass-jamming transition, such as SWAP Monte Carlo and FIRE energy minimization. The scalings of simulation runtimes with the number of particles N and number of simulation threads n_{threads} are comparable to popular molecular dynamics codes such as LAMMPS. The efficient parallel implementation allows simulation of systems that are much larger than those employed in previous high-dimensional glass-transition studies. As a demonstration of the code's capabilities, we show that supercooled d=6 liquids can possess dynamics that are substantially more heterogeneous and experience a breakdown of the Stokes-Einstein relation that is substantially stronger than previously reported, owing at least in part to the much smaller system sizes employed in earlier simulations.

Theory of Transient Localization, Activated Dynamics, and a Macromolecular Glass Transition in Ring Polymer Liquids

We construct a segmental scale force level theory for the center-of-mass diffusion constant and corresponding relaxation time for globally compact unconcatenated ring polymer solutions and melts (degree of polymerization N). The approach is based on slowly decaying macromolecular scale intermolecular force dynamic correlations as the origin of their unusual dynamics. Unentangled Rouse, weakly caged, and activated regimes are predicted. The barrier of the activated regime scales linearly with N and as a power law of concentration, which drives a kinetic glass transition on the radius-of-gyration scale. The values of N at the two dynamic crossovers (Rouse to weakly caged, weakly caged to activated) are proportional, with nonuniversality entering mainly via macromolecular volume fraction and dimensionless compressibility. Quantitative comparisons with simulation data reveal good agreement. Aspects of intermediate time dynamics are analyzed, and predictions are made for the conditions required to observe a macromolecular glass transition in the laboratory and on the computer.

Supercooled melt structure and dynamics of single-chain nanoparticles: A computer simulation study

By using coarse-grained molecular dynamics simulations, we have investigated the structure and dynamics of supercooled single-chain cross-linked nanoparticle (SCNP) melts having a range of cross-linking degrees ϕ. We find a nearly linear increase in glass-transition temperature (T_g) with increasing ϕ. Correspondingly, we have also experimentally synthesized a series of polystyrene-based SCNPs and have found that the measured T_g estimated from differential scanning calorimetry is qualitatively consistent with the trend predicted by our simulation estimates. Experimentally, an increase in T_g as large as ΔT_g = 61 K for ϕ = 0.36 is found compared with their linear chain counterparts, indicating that the changes in dynamics with cross-links are quite appreciable. We attribute the increase in T_g to the enlarged effective hard-core volume and the corresponding reduction in the free volume of the polymer segments. Topological constraints evidently frustrate the local packing. In addition, the introduction of intra-molecular cross-linking bonds slows down the structural relaxation and simultaneously enhances the local coupling motion on the length scales within SCNPs. Consequently, a more pronounced dynamical heterogeneity (DH) is observed for larger ϕ, as quantified by measuring the dynamical correlation length through the four-point susceptibility parameter, χ₄. The increase in DH is directly related to the enhanced local cooperative motion derived from intra-molecular cross-linking bonds and structural heterogeneity derived from the cross-linking process. These results shed new light on the influence of intra-molecular topological constraints on the segmental dynamics of polymer melts.

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

We employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant-matrix pair correlation function and the penetrant-matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

Machine Learning to Reveal Nanoparticle Dynamics from Liquid-Phase TEM Videos

Liquid-phase transmission electron microscopy (TEM) has been recently applied to materials chemistry to gain fundamental understanding of various reaction and phase transition dynamics at nanometer resolution. However, quantitative extraction of physical and chemical parameters from the liquid-phase TEM videos remains bottlenecked by the lack of automated analysis methods compatible with the videos' high noisiness and spatial heterogeneity. Here, we integrate, for the first time, liquid-phase TEM imaging with our customized analysis framework based on a machine learning model called U-Net neural network. This combination is made possible by our workflow to generate simulated TEM images as the training data with well-defined ground truth. We apply this framework to three typical systems of colloidal nanoparticles, concerning their diffusion and interaction, reaction kinetics, and assembly dynamics, all resolved in real-time and real-space by liquid-phase TEM. A diversity of properties for differently shaped anisotropic nanoparticles are mapped, including the anisotropic interaction landscape of nanoprisms, curvature-dependent and staged etching profiles of nanorods, and an unexpected kinetic law of first-order chaining assembly of concave nanocubes. These systems representing properties at the nanoscale are otherwise experimentally inaccessible. Compared to the prevalent image segmentation methods, U-Net shows a superior capability to predict the position and shape boundary of nanoparticles from highly noisy and fluctuating background-a challenge common and sometimes inevitable in liquid-phase TEM videos. We expect our framework to push the potency of liquid-phase TEM to its full quantitative level and to shed insights, in high-throughput and statistically significant fashion, on the nanoscale dynamics of synthetic and biological nanomaterials.

Feedback-controlled active brownian colloids with space-dependent rotational dynamics

The non-thermal nature of self-propelling colloids offers new insights into non-equilibrium physics. The central mathematical model to describe their trajectories is active Brownian motion, where a particle moves with a constant speed, while randomly changing direction due to rotational diffusion. While several feedback strategies exist to achieve position-dependent velocity, the possibility of spatial and temporal control over rotational diffusion, which is inherently dictated by thermal fluctuations, remains untapped. Here, we decouple rotational diffusion from thermal fluctuations. Using external magnetic fields and discrete-time feedback loops, we tune the rotational diffusivity of active colloids above and below its thermal value at will and explore a rich range of phenomena including anomalous diffusion, directed transport, and localization. These findings add a new dimension to the control of active matter, with implications for a broad range of disciplines, from optimal transport to smart materials.

Active colloidal systems can serve as an enabling platform to study complex out-of-equilibrium physical phenomena. Using a magnetic control with a feedback loop, here the authors program the dynamics of active Brownian particles by updating their rotational diffusion coefficient depending on their locations.

Microscopic theory of onset of decaging and bond-breaking activated dynamics in ultradense fluids with strong short-range attractions

We theoretically study thermally activated "in cage" elementary dynamical processes that precede full structural relaxation in ultradense particle liquids interacting via strong short-range attractive forces. The analysis is based on a microscopic theory formulated at the particle trajectory level built on the dynamic free energy concept and an explicit treatment of how attractive forces control the formation and lifetime of physical bonds. Mean time scales for bond breaking, the early stage of cage escape, and non-Fickian displacement by a fixed amount are analyzed in the repulsive glass, bonded repulsive (attractive) glass, fluid, and dense gel regimes. The theory predicts a strong length-scale-dependent growth of these time scales with attractive force strength at fixed packing fraction, a much weaker slowing down with density at fixed attraction strength, and a strong decoupling of the shorter bond-breaking time with the other two time scales that are controlled mainly by perturbed steric caging. All results are in good accord with simulations, and additional testable predictions are made. The classic statistical mechanical projection approximation of replacing all bare attractive and repulsive forces with a single effective force determined by pair structure incurs major errors for describing processes associated with thermally activated escape from transiently localized states.

Heterogeneity of dynamics in a modulated colloidal liquid

We study the dynamics of a system of two dimensional colloidal particles subjected to a spatially periodic external potential using Brownian dynamics simulations. We characterize the dynamics in the system by the mean square displacements and the self-van Hove function. The static density plots suggest that system gets into modulated liquid phase in presence of the external potential. We find that diffusion coefficients, obtained from long time mean sqaure displacements, decay exponentially with increasing potential strength. The self-van Hove functions computed from the distribution of particle displacemets in a given time interval show non-gaussian behaviour in directions both parallel and transverse to the external modulation. This suggests heterogeneous dynamics and is supported by particle mobilities and residence times.

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T_g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T_g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface.

A comparative study on the dynamic heterogeneity of supercooled polymers under nanoconfinement

Dynamic heterogeneity (DH) is a universal property of glass transition phenomena. In this work, we perform a comparative analysis of DH for pure polymer and polymer/nanoparticle composite systems in both film and bulk states via molecular dynamics simulations. We find that the dynamic gradient and the faster average dynamics due to the presence of a free surface are two leading factors, resulting from a nanoconfinement effect, which influence different parts of DH in a film system. The dynamic gradient results from differences in dynamics at different distances from the mobile surface, which induces a large deviation from the Gaussian distribution for the displacement distribution in the film. At the same time, the maximum string size which describes the region size for cooperative motion (dynamic correlation) can also be influenced by the dynamic gradient, although this influence is much weaker than that on the displacement distribution. On the other hand, reflecting temporal fluctuations of dynamics or temporal parts of DH, characteristic peak times of the non-Gaussian parameter and string size, and the ratio between persistent times and exchange times which describe the dynamic exchange properties, are mainly influenced by the faster dynamics on average. Our results demonstrate that measuring different properties (dynamic distribution, dynamic correlation or dynamic exchange) place an emphasis on distinct temporal and spatial parts of DH. It is necessary to use combinational measurements of these properties to give a complete picture of DH in nanoconfinement environments.

Continuous-time random-walk approach to supercooled liquids: Self-part of the van Hove function and related quantities

From equilibrium molecular dynamics (MD) simulations of a bead-spring model for short-chain glass-forming polymer melts we calculate several quantities characterizing the single-monomer dynamics near the (extrapolated) critical temperature [Formula: see text] of mode-coupling theory: the mean-square displacement g₀(t), the non-Gaussian parameter [Formula: see text] and the self-part of the van Hove function [Formula: see text] which measures the distribution of monomer displacements r in time t. We also determine these quantities from a continuous-time random walk (CTRW) approach. The CTRW is defined in terms of various probability distributions which we know from previous analysis. Utilizing these distributions the CTRW can be solved numerically and compared to the MD data with no adjustable parameter. The MD results reveal the heterogeneous and non-Gaussian single-particle dynamics of the supercooled melt near [Formula: see text]. In the time window of the early [Formula: see text] relaxation [Formula: see text] is large and [Formula: see text] is broad, reflecting the coexistence of monomer displacements that are much smaller ("slow particles") and much larger ("fast particles") than the average at time t, i.e. than [Formula: see text]. For large r the tail of [Formula: see text] is compatible with an exponential decay, as found for many glassy systems. The CTRW can reproduce the spatiotemporal dependence of [Formula: see text] at a qualitative to semiquantitative level. However, it is not quantitatively accurate in the studied temperature regime, although the agreement with the MD data improves upon cooling. In the early [Formula: see text] regime we also analyze the MD results for [Formula: see text] via the space-time factorization theorem predicted by ideal mode-coupling theory. While we find the factorization to be well satisfied for small r, both above and below [Formula: see text] , deviations occur for larger r comprising the tail of [Formula: see text]. The CTRW analysis suggests that single-particle "hops" are a contributing factor for these deviations.

Colloidal probe dynamics in gelatin solution during the sol-gel transition

The dynamics of the colloidal probes in a gelatin solution during the time-dependent sol-gel transition was investigated by multi-particle tracking. The relationship between the relaxation of the medium at the critical gel point and the mean square displacement of the probes was elucidated. Based on this understanding, the critical gel point of gelatin and the corresponding critical exponent n were unambiguously determined by the loss angle criterion and the time-cure superposition. The shift factors of the latter are further used to estimate the time/length-scale evolution of the gelatin during the sol-gel transition. The growth of the medium length scale crossed with the two measuring length scales successively at the pre-gel regime. Coinciding with the length-scale crossovers, the probability density function (PDF) of the probe displacements displayed two transient peaks of non-Gaussianity. In the post-gel regime, the third peak of Gaussianity suggested inhomogeneity in the gel network. The non-Gaussianity results from the bifurcation of diffusivity. The present work showed that the non-Gaussian dynamics of the probes are not the direct equivalence of that of the medium, but an effect of length-scale coupling.

Diffusing diffusivity: Rotational diffusion in two and three dimensions

We consider the problem of calculating the probability distribution function (pdf) of angular displacement for rotational diffusion in a crowded, rearranging medium. We use the diffusing diffusivity model and following our previous work on translational diffusion [R. Jain and K. L. Sebastian, J. Phys. Chem. B 120, 3988 (2016)], we show that the problem can be reduced to that of calculating the survival probability of a particle undergoing Brownian motion, in the presence of a sink. We use the approach to calculate the pdf for the rotational motion in two and three dimensions. We also propose new dimensionless, time dependent parameters, αrot,2D and αrot,3D, which can be used to analyze the experimental/simulation data to find the extent of deviation from the normal behavior, i.e., constant diffusivity, and obtain explicit analytical expressions for them, within our model.

Morphology dictated heterogeneous dynamics in two-dimensional aggregates

Particulate aggregates occur in a variety of non-equilibrium steady-state morphologies ranging from finite-size compact crystalline structures to non-compact string-like conformations. This diversity is due to the competition between pair-wise short range attraction and long range repulsion between particles. We identify different microscopic mechanisms in action by following the simulated particle trajectories for different morphologies in two dimensions at a fixed density and temperature. In particular, we show that the compact clusters are governed by symmetric caging of particles by their nearest neighbors while sidewise asymmetric binding of particles leads to non-compact aggregates. The measured timescales for these two mechanisms are found to be distinctly different providing phenomenological evidence of a relation between microstructure and dynamics of particulate aggregates. Supporting these findings, the time dependent diffusivity is observed to differ across the morphological hierarchy, while the average long-time dynamics is, in general, sub-diffusive at 'low' temperatures. Finally, one generic relation between diffusivity and structural randomness, applicable to simple equilibrium systems, is validated for complex aggregate forming systems through further analysis of the same system at different temperatures.

Active microrheology in a colloidal glass

We study the dynamics of a probe particle driven by a constant force through a colloidal glass of hard spheres. This nonequilibrium and anisotropic problem is investigated using a new implementation of the mode-coupling approximation with multiple relaxation channels and Langevin dynamics simulations. A force threshold is found, below which the probe remains localized, while above it the probe acquires a finite velocity. We focus on the localized regime, comparing theory and simulations concerning the dynamics in the length scale of the cage and the properties of the transition to the delocalized regime, such as the critical power-law decay of the probe correlation function. Probe van Hove functions predicted by the theory show exponential tails reminiscent of an intermittent dynamics of the probe. This scenario is microscopically supported by simulations.

Diffusion in a Crowded, Rearranging Environment

It has been found in many experiments that the mean square displacement of a Brownian particle x(T) diffusing in a rearranging environment is strictly Fickian, obeying ⟨(x(T))(2)⟩ ∝ T, but the probability distribution function for the displacement is not Gaussian. An explanation of this is that the diffusivity of the particle itself is changing as a function of time. Models for this diffusing diffusivity have been solved analytically in the limit of small time, but simulations were necessary for intermediate and large times. We show that one of the diffusing diffusivity models is equivalent to Brownian motion in the presence of a sink and introduce a class of models for which it is possible to find analytical solutions. Our solution gives ⟨(x(T))(2)⟩ ∝ T for all times and at short times the probability distribution function of the displacement is exponential which crosses over to a Gaussian in the limit of long times and large displacements.

A molecular dynamics simulations study on the relations between dynamical heterogeneity, structural relaxation, and self-diffusion in viscous liquids

We perform molecular dynamics simulations for viscous liquids to study the relations between dynamical heterogeneity, structural (α) relaxation, and self-diffusion. For atomistic models of supercooled water, polymer melts, and an ionic liquid, we characterize the space-time characteristics of dynamical heterogeneity by the degree of deviations from Gaussian displacement statistics (α2), the size of clusters comprising highly mobile particles (S(w)), and the length of strings consisting of cooperatively moving particles (L(w)). Comparison of our findings with previous simulation results for a large variety of viscous liquids, ranging from monoatomic liquids to silica melt, reveals a nearly universal decoupling between the time scales of maximum non-Gaussian parameter (τ(α2)) and the time constant of the α relaxation (τ(α)) upon cooling, explicitly, τ(α2) ∝τ(α)(3/4). Such uniform relation was not observed between the peak times of S(w) or L(w) and τ(α). On the other hand, the temperature-dependent time scale of maximum string length (τ(L)) follows the inverse of the self-diffusion coefficient (D) for various systems at sufficiently low temperatures, i.e., τ(L) ∝ D(-1). These observations are discussed in view of a breakdown of the Stokes-Einstein relation for the studied systems. It is found that the degree of deviation from this relation is correlated with the stretching of the α relaxation.

Potential energy landscape and inherent dynamics of a hard-sphere fluid

Hard-sphere models exhibit many of the same kinds of supercooled-liquid behavior as more realistic models of liquids, but the highly nonanalytic character of their potentials makes it a challenge to think of that behavior in potential energy landscape terms. We show here that it is possible to calculate an important topological property of hard-sphere landscapes, the geodesic pathways through those landscapes, and to do so without artificially coarse-graining or softening the potential. We show, moreover, that the rapid growth of the lengths of those pathways with increasing packing fraction quantitatively predicts the precipitous decline in diffusion constants in a glass-forming hard-sphere mixture model. The geodesic paths themselves can be considered as defining the intrinsic dynamics of hard spheres, so it is also revealing to find that they (and therefore the features of the underlying potential energy landscape) correctly predict the occurrence of dynamic heterogeneity and nonzero values of the non-Gaussian parameter. The success of these landscape predictions for the dynamics of such a singular model emphasizes that there is more to potential energy landscapes than is revealed by looking at the minima and saddle points.

Continuous-time random-walk approach to supercooled liquids. II. Mean-square displacements in polymer melts

The continuous-time random walk (CTRW) describes the single-particle dynamics as a series of jumps separated by random waiting times. This description is applied to analyze trajectories from molecular dynamics (MD) simulations of a supercooled polymer melt. Based on the algorithm presented by Helfferich et al. [Phys. Rev. E 89, 042603 (2014)], we detect jump events of the monomers. As a function of temperature and chain length, we examine key distributions of the CTRW: the jump-length distribution (JLD), the waiting-time distribution (WTD), and the persistence-time distribution (PTD), i.e., the distribution of waiting times for the first jump. For the equilibrium (polymer) liquid under consideration, we verify that the PTD is determined by the WTD. For the mean-square displacement (MSD) of a monomer, the results for the CTRW model are compared with the underlying MD data. The MD data exhibit two regimes of subdiffusive behavior, one for the early α process and another at later times due to chain connectivity. By contrast, the analytical solution of the CTRW yields diffusive behavior for the MSD at all times. Empirically, we can account for the effect of chain connectivity in Monte Carlo simulations of the CTRW. The results of these simulations are then in good agreement with the MD data in the connectivity-dominated regime, but not in the early α regime where they systematically underestimate the MSD from the MD.

Continuous-time random-walk approach to supercooled liquids. I. Different definitions of particle jumps and their consequences

Single-particle trajectories in supercooled liquids display long periods of localization interrupted by "fast moves." This observation suggests a modeling by a continuous-time random walk (CTRW). We perform molecular dynamics simulations of equilibrated short-chain polymer melts near the critical temperature of mode-coupling theory Tc and extract "moves" from the monomer trajectories. We show that not all moves comply with the conditions of a CTRW. Strong forward-backward correlations are found in the supercooled state. A refinement procedure is suggested to exclude these moves from the analysis. We discuss the repercussions of the refinement on the jump-length and waiting-time distributions as well as on characteristic time scales, such as the average waiting time ("exchange time") and the average time for the first move ("persistence time"). The refinement modifies the temperature (T) dependence of these time scales. For instance, the average waiting time changes from an Arrhenius-type to a Vogel-Fulcher-type T dependence. We discuss this observation in the context of the bifurcation of the α process and (Johari) β process found in many glass-forming materials to occur near Tc. Our analysis lays the foundation for a study of the jump-length and waiting-time distributions, their temperature and chain-length dependencies, and the modeling of the monomer dynamics by a CTRW approach in the companion paper [J. Helfferich et al., Phys. Rev. E 89, 042604 (2014)].

Avalanches mediate crystallization in a hard-sphere glass

By molecular-dynamics simulations, we have studied the devitrification (or crystallization) of aged hard-sphere glasses. First, we find that the dynamics of the particles are intermittent: Quiescent periods, when the particles simply "rattle" in their nearest-neighbor cages, are interrupted by abrupt "avalanches," where a subset of particles undergo large rearrangements. Second, we find that crystallization is associated with these avalanches but that the connection is not straightforward. The amount of crystal in the system increases during an avalanche, but most of the particles that become crystalline are different from those involved in the avalanche. Third, the occurrence of the avalanches is a largely stochastic process. Randomizing the velocities of the particles at any time during the simulation leads to a different subsequent series of avalanches. The spatial distribution of avalanching particles appears random, although correlations are found among avalanche initiation events. By contrast, we find that crystallization tends to take place in regions that already show incipient local order.

Coherent neutron scattering and collective dynamics on mesoscale

By combining, and modestly extending, a variety of theoretical concepts for the dynamics of liquids in the supercooled regime, we formulate a simple analytic model for the temperature and wavevector dependent collective density fluctuation relaxation time that is measurable using coherent dynamic neutron scattering. Comparison with experiments on the ionic glass-forming liquid Ca-K-NO3 in the lightly supercooled regime suggests the model captures the key physics in both the local cage and mesoscopic regimes, including the unusual wavevector dependence of the collective structural relaxation time. The model is consistent with the idea that the decoupling between diffusion and viscosity is reflected in a different temperature dependence of the collective relaxation time at intermediate wavevectors and near the main (cage) peak of the static structure factor. More generally, our analysis provides support for the ideas that decoupling information and growing dynamic length scales can be at least qualitatively deduced by analyzing the collective relaxation time as a function of temperature and wavevector, and that there is a strong link between dynamic heterogeneity phenomena at the single and many particle level. Though very simple, the model can be applied to other systems, such as molecular liquids.

Slow relaxations and stringlike jump motions in fragile glass-forming liquids: breakdown of the Stokes-Einstein relation

We perform molecular dynamics simulation on a glass-forming liquid binary mixture with the soft-core potential in three dimensions. We investigate crossover of the configuration changes caused by stringlike jump motions. With lowering the temperature T, the motions of the particles composing strings become larger in sizes and displacements, while those of the particles surrounding strings become smaller. Then the contribution of the latter to time-correlation functions tends to be long-lived as T is lowered. As a result, the relaxation time τ(α) and the viscosity η grow more steeply than the inverse diffusion constant D(-1) at low T, leading to breakdown of the Stokes-Einstein relation. At low T, the diffusion occurs as activation processes and may well be described by short-time analysis of rare jump motions with broken bonds and large displacements. Some characteristic features of the Van Hove self-correlation function arise from escape jumps over high potential barriers. We also visualize the particle motions at string formation taking place in a very short time.

Decoupling of rotational and translational diffusion in supercooled colloidal fluids

We use confocal microscopy to directly observe 3D translational and rotational diffusion of tetrahedral clusters, which serve as tracers in colloidal supercooled fluids. We find that as the colloidal glass transition is approached, translational and rotational diffusion decouple from each other: Rotational diffusion remains inversely proportional to the growing viscosity whereas translational diffusion does not, decreasing by a much lesser extent. We quantify the rotational motion with two distinct methods, finding agreement between these methods, in contrast with recent simulation results. The decoupling coincides with the emergence of non-Gaussian displacement distributions for translation whereas rotational displacement distributions remain Gaussian. Ultimately, our work demonstrates that as the glass transition is approached, the sample can no longer be approximated as a continuum fluid when considering diffusion.

Space-time correlated two-particle hopping in glassy fluids: structural relaxation, irreversibility, decoupling, and facilitation

The microscopic nonlinear Langevin equation (NLE) theory of correlated two-particle dynamics in dense fluids of spherical particles is extended to construct a predictive model of multiple correlated hopping and recaging events of a pair of tagged particles as a function of their initial separation. Modest coarse graining over the liquid structural disorder allows contact to be made with various definitions of irreversible particle motion within the context of a multistate Markov model. The correlated space-time hopping process that underlies structural relaxation can also be analyzed in the context of kinetically constrained models. The dependence of microscopically defined mean persistence and exchange times, their distributions, and relaxation-diffusion decoupling on hard-sphere fluid volume fraction is derived from a model in which irreversible jumps serve as the nucleating persistence event. For a subset of questions, the predictions of the two-particle theory are compared with results from the earlier single-particle NLE approach.

Analysis of a growing dynamic length scale in a glass-forming binary hard-sphere mixture

We examine a length scale that characterizes the spatial extent of heterogeneous dynamics in a glass-forming binary hard-sphere mixture up to the mode-coupling volume fraction ϕ(c). First, we characterize the system's dynamics. Then, we utilize a method [Phys. Rev. Lett. 105, 217801 (2010)] to extract and analyze the ensemble-independent dynamic susceptibility χ(4)(t) and the dynamic correlation length ξ(t) for a range of times between the β and α relaxation times. We find that in this time range the dynamic correlation length follows a volume fraction-independent master curve ξ(t)~ln(t). For longer times, ξ(t) departs from this master curve and remains constant up to the largest time at which we can determine the length accurately. In addition to the previously established correlation τ(α)~exp[ξ(τ(α))] between the α relaxation time, τ(α), and the dynamic correlation length at this time, ξ(τ(α)), we also find a similar correlation for the diffusion coefficient D~exp[ξ(τ(α))(θ)] with θ≈0.6. We discuss the relevance of these findings for different theories of the glass transition.

Activated dynamics in dense fluids of attractive nonspherical particles. II. Elasticity, barriers, relaxation, fragility, and self-diffusion

In paper II of this series we apply the center-of-mass version of Nonlinear Langevin Equation theory to study how short-range attractive interactions influence the elastic shear modulus, transient localization length, activated dynamics, and kinetic arrest of a variety of nonspherical particle dense fluids (and the spherical analog) as a function of volume fraction and attraction strength. The activation barrier (roughly the natural logarithm of the dimensionless relaxation time) is predicted to be a rich function of particle shape, volume fraction, and attraction strength, and the dynamic fragility varies significantly with particle shape. At fixed volume fraction, the barrier grows in a parabolic manner with inverse temperature nondimensionalized by an onset value, analogous to what has been established for thermal glass-forming liquids. Kinetic arrest boundaries lie at significantly higher volume fractions and attraction strengths relative to their dynamic crossover analogs, but their particle shape dependence remains the same. A limited universality of barrier heights is found based on the concept of an effective mean-square confining force. The mean hopping time and self-diffusion constant in the attractive glass region of the nonequilibrium phase diagram is predicted to vary nonmonotonically with attraction strength or inverse temperature, qualitatively consistent with recent computer simulations and colloid experiments.

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.

Theory of aging, rejuvenation, and the nonequilibrium steady state in deformed polymer glasses

The nonlinear Langevin equation theory of segmental relaxation, elasticity, and mechanical response of polymer glasses is extended to describe the coupled effects of physical aging, mechanical rejuvenation, and thermal history. The key structural variable is the amplitude of density fluctuations, and segmental dynamics proceeds via stress-modified activated barrier hopping on a dynamic free-energy profile. Mechanically generated disorder (rejuvenation) is quantified by a dissipative work argument and increases the amplitude of density fluctuations, thereby speeding up relaxation beyond that induced by the landscape tilting mechanism. The theory makes testable predictions for the time evolution and nonequilibrium steady state of the alpha relaxation time, density fluctuation amplitude, elastic modulus, and other properties. Model calculations reveal a rich dependence of these quantities on preaging time, applied stress, and temperature that reflects the highly nonlinear competition between physical aging and mechanical disordering. Thermal history is "erased" in the long-time limit, although the nonequilibrium steady state is not the literal "fully rejuvenated" freshly quenched glass. The present work provides the conceptual foundation for a quantitative treatment of the nonlinear mechanical response of polymer glasses under a variety of deformation protocols.

Segmental dynamics in polymers: from cold melts to ageing and stressed glasses

Recent progress in developing statistical mechanical theories of supercooled polymer melts and glasses is reviewed. The focus is on those approaches that are either explicitly formulated for polymers, or are applications of more generic theories to interpret polymeric phenomena. These include two configurational entropy theories, a percolated free volume distribution model, and the activated barrier hopping nonlinear Langevin theory. Both chemically-specific and universal aspects are discussed. After a brief summary of classic phenomenological approaches, a discussion of the relevant length scales and key experimental phenomena in both the supercooled liquid and glassy solid state is presented including ageing and nonlinear mechanical response. The central concepts that underlie the theories in the molten state are then summarized and key predictions discussed, including the glass transition in oriented polymer liquids and deformed rubber networks. Physical ageing occurs in the nonequilibrium glass, and theories for its consequences on the alpha relaxation are discussed. Very recent progress in developing a segment scale theory for the dramatic effects of external stress on polymer glasses, including acceleration of relaxation, yielding, plastic flow and strain hardening, is summarized. The article concludes with a discussion of outstanding theoretical challenges.

Anomalous yet Brownian

We describe experiments using single-particle tracking in which mean-square displacement is simply proportional to time (Fickian), yet the distribution of displacement probability is not Gaussian as should be expected of a classical random walk but, instead, is decidedly exponential for large displacements, the decay length of the exponential being proportional to the square root of time. The first example is when colloidal beads diffuse along linear phospholipid bilayer tubes whose radius is the same as that of the beads. The second is when beads diffuse through entangled F-actin networks, bead radius being less than one-fifth of the actin network mesh size. We explore the relevance to dynamic heterogeneity in trajectory space, which has been extensively discussed regarding glassy systems. Data for the second system might suggest activated diffusion between pores in the entangled F-actin networks, in the same spirit as activated diffusion and exponential tails observed in glassy systems. But the first system shows exceptionally rapid diffusion, nearly as rapid as for identical colloids in free suspension, yet still displaying an exponential probability distribution as in the second system. Thus, although the exponential tail is reminiscent of glassy systems, in fact, these dynamics are exceptionally rapid. We also compare with particle trajectories that are at first subdiffusive but Fickian at the longest measurement times, finding that displacement probability distributions fall onto the same master curve in both regimes. The need is emphasized for experiments, theory, and computer simulation to allow definitive interpretation of this simple and clean exponential probability distribution.

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.