From coupled elementary units to the complexity of the glass transition

Excess wings and asymmetric relaxation spectra in a facilitated trap model

In a recent computer study, we have shown that the combination of spatially heterogeneous dynamics and kinetic facilitation provides a microscopic explanation for the emergence of excess wings in deeply supercooled liquids. Motivated by these findings, we construct a minimal empirical model to describe this physics and introduce dynamic facilitation in the trap model, which was initially developed to capture the thermally activated dynamics of glassy systems. We fully characterize the relaxation dynamics of this facilitated trap model varying the functional form of energy distributions and the strength of dynamic facilitation, combining numerical results and analytic arguments. Dynamic facilitation generically accelerates the relaxation of the deepest traps, thus making relaxation spectra strongly asymmetric, with an apparent "excess" signal at high frequencies. For well-chosen values of the parameters, the obtained spectra mimic experimental results for organic liquids displaying an excess wing. Overall, our results identify the minimal physical ingredients needed to describe excess processes in the relaxation spectra of supercooled liquids.

Structural relaxation and crystallization in supercooled water from 170 to 260 K

The origin of water's anomalous properties has been debated for decades. Resolution of the problem is hindered by a lack of experimental data in a crucial region of temperatures, T, and pressures where supercooled water rapidly crystallizes-a region often referred to as "no man's land." A recently developed technique where water is heated and cooled at rates greater than 10⁹ K/s now enables experiments in this region. Here, it is used to investigate the structural relaxation and crystallization of deeply supercooled water for 170 K < T < 260 K. Water's relaxation toward a new equilibrium structure depends on its initial structure with hyperquenched glassy water (HQW) typically relaxing more quickly than low-density amorphous solid water (LDA). For HQW and T > 230 K, simple exponential relaxation kinetics is observed. For HQW at lower temperatures, increasingly nonexponential relaxation is observed, which is consistent with the dynamics expected on a rough potential energy landscape. For LDA, approximately exponential relaxation is observed for T > 230 K and T < 200 K, with nonexponential relaxation only at intermediate temperatures. At all temperatures, water's structure can be reproduced by a linear combination of two, local structural motifs, and we show that a simple model accounts for the complex kinetics within this context. The relaxation time, τ _rel , is always shorter than the crystallization time, τ _xtal For HQW, the ratio, τ _xtal /τ _rel , goes through a minimum at ∼198 K where the ratio is about 60.

Understanding the nonlinear dynamics of driven particles in supercooled liquids in terms of an effective temperature

In active microrheology, the mechanical properties of a material are tested by adding probe particles which are pulled by an external force. In case of supercooled liquids, strong forcing leads to a thinning of the host material which becomes more pronounced as the system approaches the glass transition. In this work, we provide a quantitative theoretical description of this thinning behavior based on the properties of the Potential Energy Landscape (PEL) of a model glass-former. A key role plays the trap-like nature of the PEL. We find that the mechanical properties in the strongly driven system behave the same as in a quiescent system at an enhanced temperature, giving rise to a well-characterized effective temperature. Furthermore, this effective temperature turns out to be independent of the chosen observable and individually shows up in the thermodynamic and dynamic properties of the system. Based on this underlying theoretical understanding, we can estimate its dependence on temperature and force by the PEL-properties of the quiescent system. We furthermore critically discuss the relevance of effective temperatures obtained by scaling relations for the description of out-of-equilibrium situations.

Finite-size scaling for the glass transition: the role of a static length scale

Over the past decade, computer simulations have had an increasing role in shedding light on difficult statistical physical phenomena, and in particular on the ubiquitous problem of the glass transition. Here in a wide variety of materials the viscosity of a supercooled liquid increases by many orders of magnitude upon decreasing the temperature over a modest range. A natural concern in these computer simulations is the very small size of the simulated systems compared to experimental ones, raising the issue of how to assess the thermodynamic limit. Here we turn this limitation to our advantage by performing finite size scaling on the system size dependence of the relaxation time for supercooled liquids to emphasize the importance of a growing static length scale in the theory of glass transition. We demonstrate that the static length scale that was discovered by us in Physica A 391, 1001 (2012) fits the bill extremely well, allowing us to provide a finite-size scaling theory for the α-relaxation time of the glass transition, including predictions for the thermodynamic limit based on simulations in small systems.

How coupled elementary units determine the dynamics of macroscopic glass-forming systems

We investigate the dynamics of a binary mixture Lennard-Jones system of different system sizes with respect to the importance of the properties of the underlying potential energy landscape (PEL). We show that the dynamics of small systems can be very well described within the continuous time random walk formalism, which is determined solely by PEL parameters. Finite size analysis shows that the diffusivity of large and small systems are very similar. This suggests that the PEL parameters of the small system also determine the local dynamics in large systems. The structural relaxation time, however, displays significant finite size effects. Furthermore, using a nonequilibrium configuration of a large system, we find that causal connections exist between nearby regions of the system. These findings can be described by the coupled landscape model for which a macroscopic system is described by a superposition of elementary systems, each described by its PEL. A minimum coupling is introduced which accounts for the finite size behavior. The coupling strength, as the single adjustable parameter, becomes smaller closer to the glass transition.

Manifestations of dynamical facilitation in glassy materials

By characterizing the dynamics of idealized lattice models with a tunable kinetic constraint, we explore the different ways in which dynamical facilitation manifests itself within the local dynamics of glassy materials. Dynamical facilitation is characterized both by a mobility transfer function, the propensity for highly mobile regions to arise near regions that were previously mobile, and by a facilitation volume, the effect of an initial dynamical event on subsequent dynamics within a region surrounding it. Sustained bursts of dynamical activity-avalanches-are shown to occur in kinetically constrained models, but, contrary to recent claims, we find that the decreasing spatiotemporal extent of avalanches with increased supercooling previously observed in granular experiments does not imply diminishing facilitation. Viewed within the context of existing simulation and experimental evidence, our findings show that dynamical facilitation plays a significant role in the dynamics of systems investigated over the range of state points accessible to molecular simulations and granular experiments.

Spatial correlation of the dynamic propensity of a glass-forming liquid

We present computer simulation results on the dynamic propensity (as defined by Widmer-Cooper et al 2004 Phys. Rev. Lett. 93 135701) in a Kob-Andersen binary Lennard-Jones liquid system consisting of 8788 particles. We compute the spatial correlation function for the dynamic propensity as a function of both the reduced temperature T, and the time scale on which the particle displacements are measured. For T ≤ 0.6, we find that non-zero correlations occur at the largest length scale accessible in our system. We also show that a cluster-size analysis of particles with extremal values of the dynamic propensity, as well as 3D visualizations, reveal spatially correlated regions that approach the size of our system as T decreases, consistently with the behavior of the spatial correlation function. Next, we define and examine the 'coordination propensity', the isoconfigurational average of the coordination number of the minority B particles around the majority A particles. We show that a significant correlation exists between the spatial fluctuations of the dynamic and coordination propensities. In addition, we find non-zero correlations of the coordination propensity occurring at the largest length scale accessible in our system for all T in the range 0.466 < T < 1.0. We discuss the implications of these results for understanding the length scales of dynamical heterogeneity in glass-forming liquids.

Viscosity, shear waves, and atomic-level stress-stress correlations

The Green-Kubo equation relates the macroscopic stress-stress correlation function to a liquid's viscosity. The concept of the atomic-level stresses allows the macroscopic stress-stress correlation function in the equation to be expressed in terms of the space-time correlations among the atomic-level stresses. Molecular dynamics studies show surprisingly long spatial extension of stress-stress correlations and also longitudinal and transverse waves propagating in liquids over ranges which could exceed the system size. The results reveal that the range of propagation of shear waves corresponds to the range of distances relevant for viscosity. Thus our results show that viscosity is a fundamentally nonlocal quantity. We also show that the periodic boundary conditions play a nontrivial role in molecular dynamics simulations, effectively masking the long-range nature of viscosity.