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

The distinguishable-particle lattice model of glasses in three dimensions

The nature of glassy states in realistic finite dimensions is still under fierce debate. Lattice models can offer valuable insights and facilitate deeper theoretical understanding. Recently, a disordered-interacting lattice model with distinguishable particles in two dimensions (2D) has been shown to produce a wide range of dynamical properties of structural glasses, including the slow and heterogeneous characteristics of the glassy dynamics, various fragility behaviors of glasses, and so on. These findings support the usefulness of this model for modeling structural glasses. An important question is whether such properties still hold in the more realistic three dimensions. In this study, we aim to extend the distinguishable-particle lattice model (DPLM) to three dimensions (3D) and explore the corresponding glassy dynamics. Through extensive kinetic Monte Carlo simulations, we found that the 3D DPLM exhibits many typical glassy behaviors, such as plateaus in the mean square displacement of particles and the self-intermediate scattering function, dynamic heterogeneity, variability of glass fragilities, and so on, validating the effectiveness of the DPLM in a broader realistic setting. The observed glassy behaviors of the 3D DPLM appear similar to those of its 2D counterpart, in accordance with recent findings in molecular models of glasses. We further investigate the role of void-induced motions in dynamical relaxations and discuss their relation to dynamic facilitation. As lattice models tend to keep the minimal but important modeling elements, they are typically much more amenable to analysis. Therefore, we envisage that the DPLM will benefit future theoretical developments, such as the configuration tree theory, towards a more comprehensive understanding of structural glasses.

Unraveling the Glass-like Dynamic Heterogeneity in Ring Polymer Melts: From Semiflexible to Stiff Chain

Ring polymers are an intriguing class of polymers with unique physical properties, and understanding their behavior is important for developing accurate theoretical models. In this study, we investigate the effect of chain stiffness and monomer density on the static and dynamic behaviors of ring polymer melts using molecular dynamics simulations. Our first focus is on the non-Gaussian parameter of center-of-mass displacement as a measure of dynamic heterogeneity, which is commonly observed in glass-forming liquids. We find that the non-Gaussianity in the displacement distribution increases with the monomer density and stiffness of the polymer chains, suggesting that excluded volume interactions between centers of mass have a strong effect on the dynamics of ring polymers. We then analyze the relationship between the radius of gyration and monomer density for semiflexible and stiff ring polymers. Our results indicate that the relationship between the two varies with chain stiffness, which can be attributed to the competition between repulsive forces inside the ring and from adjacent rings. Finally, we study the dynamics of bond-breakage virtually connected between the centers of mass of rings to analyze the exchanges of intermolecular networks of bonds. Our results demonstrate that the dynamic heterogeneity of bond-breakage is coupled with the non-Gaussianity in ring polymer melts, highlighting the importance of the bond-breaking method in determining the intermolecular dynamics of ring polymer melts. Overall, our study sheds light on the factors that govern the dynamic behaviors of ring polymers.

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Violation of Stokes-Einstein and Stokes-Einstein-Debye relations in polymers at the gas-supercooled liquid coexistence

Molecular dynamics simulations are performed on a system of model linear polymers to look at the violations of Stokes-Einstein (SE) and Stokes-Einstein-Debye (SED) relations near the mode coupling theory transition temperatureT_cat three (one higher and two lower) densities. At low temperatures, both lower density systems show stable gas-supercooled-liquid coexistence whereas the higher density system is homogeneous. We show that monomer density relaxation exhibits SE violation for all three densities, whereas molecular density relaxation shows a weak violation of the SE relation nearT_cin both lower density systems. This study identifies disparity in monomer mobility and observation of jumplike motion in the typical monomer trajectories resulting in the SE violations. In addition to the SE violation, a weak SED violation is observed in the gas-supercooled-liquid coexisting domains of the lower densities. Both lower density systems also show a decoupling of translational and rotational dynamics in this polymer system.

Acoustic resonance in periodically sheared glass: damping due to plastic events

Using molecular dynamics simulation, we study acoustic resonance in a low-temperature model glass by applying a small periodic shear at a boundary wall. Shear wave resonance occurs as the frequency ω approaches ωl = πc⊥l/L (l = 1, 2…). Here, c⊥ is the transverse sound speed and L is the cell width. At resonance, large-amplitude sound waves appear after many cycles even if the applied strain γ0 is very small. They then induce plastic events, which are heterogeneous on the mesoscopic scale and intermittent on timescales longer than the oscillation period tp = 2π/ω. We visualize them together with the extended elastic strains around them. These plastic events serve to damp sounds. We obtain the nonlinear damping Q-1 = tan δ due to the plastic events near the first resonance at ω ≅ ω1, which is linear in γ0 and independent of ω. After many resonant cycles, we observe an increase in the shear modulus (measured after switching-off the oscillation). We also observe catastrophic plastic events after a very long time (∼103tp), which induce system-size elastic strains and cause a transition from resonant to off-resonant states. At resonance, stroboscopic diffusion becomes detectable.

Local Density Fluctuation Governs the Divergence of Viscosity Underlying Elastic and Hydrodynamic Anomalies in a 2D Glass-Forming Liquid

If a liquid is cooled rapidly to form a glass, its structural relaxation becomes retarded, producing a drastic increase in viscosity. In two dimensions, strong long-wavelength fluctuations persist, even at low temperature, making it difficult to evaluate the microscopic structural relaxation time. This Letter shows that, in a 2D glass-forming liquid, relative displacement between neighbor particles yields a relaxation time that grows in proportion to the viscosity. In addition to thermal elastic vibrations, hydrodynamic fluctuations are found to affect the long-wavelength dynamics, yielding a logarithmically diverging diffusivity in the long-time limit.

What does the instantaneous normal mode spectrum tell us about dynamical heterogeneity in glass-forming fluids?

We examine the instantaneous normal mode spectrum of model metallic and polymeric glass-forming liquids. We focus on the localized modes in the unstable part of the spectrum [unstable localized (UL) modes] and find that the particles making the dominant contribution to the participation ratio form clusters that grow upon cooling in a fashion similar to the dynamical heterogeneity in glass-forming fluids, i.e., highly mobile (or immobile) particles form clusters that grow upon cooling; however, a comparison of the UL mode clusters to the mobile and immobile particle clusters indicates that they are distinct entities. We also show that the cluster size provides an alternate method to distinguish localized and delocalized modes, offering a significant practical advantage over the finite-size scaling approach. We examine the trajectories of particles contributing most to the UL modes and find that they have a slightly enhanced mobility compared to the average, and we determine a characteristic time quantifying the persistence time of this excess mobility. This time scale is proportional to the structural relaxation time τ_α of the fluid, consistent with a prediction by Zwanzig [Phys. Rev. 156, 190 (1967)] for the lifetime of collective excitations in cooled liquids. Evidently, these collective excitations serve to facilitate relaxation but do not actually participate in the motion associated with barrier crossing events governing activated transport. They also serve as a possible concrete realization of the "facilitation" clusters postulated in previous modeling of glass-forming liquids.

Spurious violation of the Stokes-Einstein-Debye relation in supercooled water

The theories of Brownian motion, the Debye rotational diffusion model, and hydrodynamics together provide us with the Stokes–Einstein–Debye (SED) relation between the rotational relaxation time of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\boldsymbol{\ell }}$$\end{document}ℓ-th degree Legendre polynomials \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ, and viscosity divided by temperature, η/T. Experiments on supercooled liquids are frequently performed to measure the SED relations, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓk_BT/η and D_t\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ, where D_t is the translational diffusion constant. However, the SED relations break down, and its molecular origin remains elusive. Here, we assess the validity of the SED relations in TIP4P/2005 supercooled water using molecular dynamics simulations. Specifically, we demonstrate that the higher-order \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ values exhibit a temperature dependence similar to that of η/T, whereas the lowest-order \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\boldsymbol{\tau }}}_{{\boldsymbol{\ell }}}$$\end{document}τℓ values are decoupled with η/T, but are coupled with the translational diffusion constant D_t. We reveal that the SED relations are so spurious that they significantly depend on the degree of Legendre polynomials.

Theory of applying shear strains from boundary walls: Linear response in glasses

We construct a linear response theory of applying shear deformations from boundary walls in the film geometry in Kubo's theoretical scheme. Our method is applicable to any solids and fluids. For glasses, we assume quasi-equilibrium around a fixed inherent state. Then, we obtain linear-response expressions for any variables including the stress and the particle displacements, even though the glass interior is elastically inhomogeneous. In particular, the shear modulus can be expressed in terms of the correlations between the interior stress and the forces from the walls. It can also be expressed in terms of the inter-particle correlations, as has been shown in the previous literature. Our stress relaxation function includes the effect of the boundary walls and can be used for inhomogeneous flow response. We show the presence of long-ranged, long-lived correlations among the fluctuations of the forces from the walls and the displacements of all the particles in the cell. We confirm these theoretical results numerically in a two-dimensional model glass. As an application, we describe emission and propagation of transverse sounds after boundary wall motions using these time-correlation functions. We also find resonant sound amplification when the frequency of an oscillatory shear approaches that of the first transverse sound mode.

Two-tag correlations and nonequilibrium fluctuation-response relation in ageing single-file diffusion

Spatiotemporally correlated motions of interacting Brownian particles, confined in a narrow channel of infinite length, are studied in terms of statistical quantities involving two particles. A theoretical framework that allows analytical calculation of two-tag correlations is presented on the basis of the Dean-Kawasaki equation describing density fluctuations in colloidal systems. In the equilibrium case, the time-dependent Einstein relation holds between the two-tag displacement correlation and the response function corresponding to it, which is a manifestation of the fluctuation-dissipation theorem for the correlation of density fluctuations. While the standard procedure of closure approximation for nonlinear density fluctuations is known to be obstructed by inconsistency with the fluctuation-dissipation theorem, this difficulty is naturally avoided by switching from the standard Fourier representation of the density field to the label-based Fourier representation of the vacancy field. In the case of ageing dynamics started from equidistant lattice configuration, the time-dependent Einstein relation is violated, as the two-tag correlation depends on the waiting time for equilibration while the response function is not sensitive to it. Within linear approximation, however, there is a simple relation between the density (or vacancy) fluctuations and the corresponding response function, which is valid even if the system is out of equilibrium. This non-equilibrium fluctuation-response relation can be extended to the case of nonlinear fluctuations by means of closure approximation for the vacancy field.

Isolating long-wavelength fluctuation from structural relaxation in two-dimensional glass: cage-relative displacement

It has recently been revealed that long-wavelength fluctuation exists in two-dimensional (2D) glassy systems, having the same origin as that given by the Mermin-Wagner theorem for 2D crystalline solids. In this paper, we discuss how to characterise quantitatively the long-wavelength fluctuation in a molecular dynamics simulation of a lightly supercooled liquid. We employ the cage-relative mean-square displacement (MSD), defined on relative displacement to its cage, to quantitatively separate the long-wavelength fluctuation from the original MSD. For increasing system size the amplitude of acoustic long wavelength fluctuations not only increases but shifts to later times causing a crossover with structural relaxation of caging particles. We further analyse the dynamic correlation length using the cage-relative quantities. It grows as the structural relaxation becomes slower with decreasing temperature, uncovering an overestimation by the four-point correlation function due to the long-wavelength fluctuation. These findings motivate the usage of cage-relative MSD as a starting point for analysis of 2D glassy dynamics.

Brownian versus Newtonian devitrification of hard-sphere glasses

In a recent molecular dynamics simulation work it has been shown that glasses composed of hard spheres crystallize via cooperative, stochastic particle displacements called avalanches [E. Sanz et al., Proc. Natl. Acad. Sci. USA 111, 75 (2014)PNASA60027-842410.1073/pnas.1308338110]. In this Rapid Communication we investigate if such a devitrification mechanism is also present when the dynamics is Brownian rather than Newtonian. The research is motivated in part by the fact that colloidal suspensions, an experimental realization of hard-sphere systems, undergo Brownian motion. We find that Brownian hard-sphere glasses do crystallize via avalanches with very similar characteristics to those found in the Newtonian case. We briefly discuss the implications of these findings for experiments on colloids.

Identifying time scales for violation/preservation of Stokes-Einstein relation in supercooled water

The violation of the Stokes-Einstein (SE) relation D ~ (η/T)^-1 between the shear viscosity η and the translational diffusion constant D at temperature T is of great importance for characterizing anomalous dynamics of supercooled water. Determining which time scales play key roles in the SE violation remains elusive without the measurement of η. We provide comprehensive simulation results of the dynamic properties involving η and D in the TIP4P/2005 supercooled water. This enabled the thorough identification of the appropriate time scales for the SE relation Dη/T. In particular, it is demonstrated that the temperature dependence of various time scales associated with structural relaxation, hydrogen bond breakage, stress relaxation, and dynamic heterogeneities can be definitely classified into only two classes. That is, we propose the generalized SE relations that exhibit "violation" or "preservation." The classification depends on the examined time scales that are coupled or decoupled with the diffusion. On the basis of the classification, we explain the physical origins of the violation in terms of the increase in the plateau modulus and the nonexponentiality of stress relaxation. This implies that the mechanism of SE violation is attributed to the attained solidity upon supercooling, which is in accord with the growth of non-Gaussianity and spatially heterogeneous dynamics.

Repetition and pair-interaction of string-like hopping motions in glassy polymers

The dynamics of many glassy systems are known to exhibit string-like hopping motions each consisting of a line of particles displacing one another. By using the molecular dynamics simulations of glassy polymers, we show that these motions become highly repetitive back-and-forth motions as temperature decreases and do not necessarily contribute to net displacements. Particle hops which constitute string-like motions are reversed with a high probability, reaching 73% and beyond at low temperature. The structural relaxation rate is then dictated not by a simple particle hopping rate but instead by the rate at which particles break away from hopping repetitions. We propose that disruption of string repetitions and hence also structural relaxations are brought about by pair-interactions between strings.

Structural origin of fractional Stokes-Einstein relation in glass-forming liquids

In many glass-forming liquids, fractional Stokes-Einstein relation (SER) is observed above the glass transition temperature. However, the origin of such phenomenon remains elusive. Using molecular dynamics simulations, we investigate the break- down of SER and the onset of fractional SER in a model of metallic glass-forming liquid. We find that SER breaks down when the size of the largest cluster consisting of trapped atoms starts to increase sharply at which the largest cluster spans half of the simulations box along one direction, and the fractional SER starts to follows when the largest cluster percolates the entire system and forms 3-dimentional network structures. Further analysis based on the percolation theory also confirms that percolation occurs at the onset of the fractional SER. Our results directly link the breakdown of the SER with structure inhomogeneity and onset of the fraction SER with percolation of largest clusters, thus provide a possible picture for the break- down of SER and onset of fractional SER in glass-forming liquids, which is is important for the understanding of the dynamic properties in glass-forming liquids.

A structural signature of the breakdown of the Stokes–Einstein relation in metallic liquids

The study provides a possible structural origin for the breakdown of the Stokes–Einstein relation in metallic liquids.

Unveiling Dimensionality Dependence of Glassy Dynamics: 2D Infinite Fluctuation Eclipses Inherent Structural Relaxation

By using large-scale molecular dynamics simulations, the dynamics of two-dimensional (2D) supercooled liquids turns out to be dependent on the system size, while the size dependence is not pronounced in three-dimensional (3D) systems. It is demonstrated that the strong system-size effect in 2D amorphous systems originates from the enhanced fluctuations at long wavelengths which are similar to those of 2D crystal phonons. This observation is further supported by the frequency dependence of the vibrational density of states, consisting of the Debye approximation in the low-wave-number limit. However, the system-size effect in the intermediate scattering function becomes negligible when the length scale is larger than the vibrational amplitude. This suggests that the finite-size effect in a 2D system is transient and also that the structural relaxation itself is not fundamentally different from that in a 3D system. In fact, the dynamic correlation lengths estimated from the bond-breakage function, which do not suffer from those enhanced fluctuations, are not size dependent in either 2D or 3D systems.

Decoupled length scales for diffusivity and viscosity in glass-forming liquids

The growth of the characteristic length scales both for diffusion and viscosity is investigated by molecular dynamics utilizing the finite-size effect in a binary Lennard-Jones mixture. For those quantities relevant to the diffusion process (e.g., the hydrodynamic value and the spatial correlation function), a strong system-size dependence is found. In contrast, it is weak or absent for the shear relaxation process. Correlation lengths are estimated from the decay of the spatial correlation functions. We find the length scale for viscosity decouples from the one of diffusivity, featured by a saturated length even in high supercooling. This temperature-independent behavior of the length scale is reminiscent of the unapparent structure change upon supercooling, implying the manifestation of configuration entropy. Whereas for the diffusion process, it is manifested by relaxation dynamics and dynamic heterogeneity. The Stokes-Einstein relation is found to break down at the temperature where the decoupling of these lengths happens.

Calculation of displacement correlation tensor indicating vortical cooperative motion in two-dimensional colloidal liquids

As an indicator of cooperative motion in a system of Brownian particles that models two-dimensional colloidal liquids, a displacement correlation tensor is calculated analytically and compared with numerical results. The key idea for the analytical calculation is to relate the displacement correlation tensor, which is a kind of four-point space-time correlation, to the Lagrangian two-time correlation of the deformation gradient tensor. Tensorial treatment of the statistical quantities, including the displacement correlation itself, allows capturing the vortical structure of the cooperative motion. The calculated displacement correlation also implies a negative long-time tail in the velocity autocorrelation, which is a manifestation of the cage effect. Both the longitudinal and transverse components of the displacement correlation are found to be expressible in terms of a similarity variable, suggesting that the cages are nested to form a self-similar structure in the space-time.

Sucrose diffusion in aqueous solution

The diffusion of sugar in aqueous solution is important both in nature and in technological applications, yet measurements of diffusion coefficients at low water content are scarce. We report directly measured sucrose diffusion coefficients in aqueous solution. Our technique utilises a Raman isotope tracer method to monitor the diffusion of non-deuterated and deuterated sucrose across a boundary between the two aqueous solutions. At a water activity of 0.4 (equivalent to 90 wt% sucrose) at room temperature, the diffusion coefficient of sucrose was determined to be approximately four orders of magnitude smaller than that of water in the same material. Using literature viscosity data, we show that, although inappropriate for the prediction of water diffusion, the Stokes-Einstein equation works well for predicting sucrose diffusion under the conditions studied. As well as providing information of importance to the fundamental understanding of diffusion in binary solutions, these data have technological, pharmaceutical and medical implications, for example in cryopreservation. Moreover, in the atmosphere, slow organic diffusion may have important implications for aerosol growth, chemistry and evaporation, where processes may be limited by the inability of a molecule to diffuse between the bulk and the surface of a particle.

Insights from Single-File Diffusion into Cooperativity in Higher Dimensions

Diffusion in colloidal suspensions can be very slow due to the cage effect, which confines each particle within a short radius on one hand, and involves large-scale cooperative motions on the other. In search of insight into this cooperativity, here the authors develop a formalism to calculate the displacement correlation in colloidal systems, mainly in the two-dimensional (2D) case. To clarify the idea for it, studies are reviewed on cooperativity among the particles in the one-dimensional (1D) case, i.e. the single-file diffusion (SFD). As an improvement over the celebrated formula by Alexander and Pincus on the mean-square displacement (MSD) in SFD, it is shown that the displacement correlation in SFD can be calculated from Lagrangian correlation of the particle interval in the one-dimensional case, and also that the formula can be extended to higher dimensions. The improved formula becomes exact for large systems. By combining the formula with a nonlinear theory for correlation, a correction to the asymptotic law for the MSD in SFD is obtained. In the 2D case, the linear theory gives description of vortical cooperative motion.

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.

Spatial correlations of elementary relaxation events in glass-forming liquids

The dynamical facilitation scenario, by which localized relaxation events promote nearby relaxation events in an avalanche process, has been suggested as the key mechanism connecting the microscopic and the macroscopic dynamics of structural glasses. Here we investigate the statistical features of this process via numerical simulations of a model structural glass. First we show that the relaxation dynamics of the system occurs through particle jumps that are irreversible, and that cannot be decomposed in smaller irreversible events. Then we show that each jump does actually trigger an avalanche. The characteristics of this avalanche change upon cooling, suggesting that the relaxation dynamics crossovers from a noise dominated regime, where jumps do not trigger other relaxation events, to a regime dominated by the facilitation process, where a jump triggers more relaxation events.

Correlation between crystalline order and vitrification in colloidal monolayers

We investigate experimentally the relationship between local structure and dynamical arrest in a quasi-2d colloidal model system which approximates hard discs. We introduce polydispersity to the system to suppress crystallisation. Upon compression, the increase in structural relaxation time is accompanied by the emergence of local hexagonal symmetry. Examining the dynamical heterogeneity of the system, we identify three types of motion: 'zero-dimensional' corresponding to β-relaxation, 'one-dimensional' or stringlike motion and '2D' motion. The dynamic heterogeneity is correlated with the local order, that is to say locally hexagonal regions are more likely to be dynamically slow. However, we find that lengthscales corresponding to dynamic heterogeneity and local structure do not appear to scale together approaching the glass transition.

Dynamics in a tetrahedral network glassformer: vibrations, network rearrangements, and diffusion

We perform molecular dynamics simulation on a tetrahedral network glassformer using a model for viscous SiO2 by Coslovich and Pastore [J. Phys.: Condens. Matter 21, 285107 (2009)]. In this system, Si and O particles form a random network at low temperature T. We attach an ellipsoid to each particle to represent its time-averaged vibration tensor. We then examine the anisotropic vibrations of Si and O, where the ellipsoid orientations are correlated with the network. The ellipsoids exhibit marked vibrational heterogeneity. The configuration changes occur as breakage and reorganization of the network, where only one or two particles undergo large jumps at each rearrangement leading to diffusion. To the time-correlation functions, however, the particles surrounding these largely displaced ones yield significantly T-dependent contributions, resulting in a weak violation of the Stokes-Einstein relation. This crossover is mild in silica due to the small Si-O bond numbers per particle, while it is strong in fragile glassformers with large coordination numbers. On long timescales, jump events tend to occur in the same regions forming marked dynamic heterogeneity. We also calculate the diffusion constants and the viscosity. The diffusion obeys activation dynamics and may be studied by short-time analysis of irreversible jumps.

Connecting short and long time dynamics in hard-sphere-like colloidal glasses

Glass-forming materials are characterized by an intermittent motion at the microscopic scale. Particles spend most of their time rattling within the cages formed by their neighbors, and seldom jump to a different cage. In molecular glass formers the temperature dependence of the jump features, such as the average caging time and jump length, characterizes the relaxation processes and allows for a short-time prediction of the diffusivity. Here we experimentally investigate the cage-jump motion of a two-dimensional hard-sphere-like colloidal suspension, where the volume fraction is the relevant parameter controlling the slowing down of the dynamics. We characterize the volume fraction dependence of the cage-jump features and show that, as in molecular systems, they allow for a short time prediction of the diffusivity.

From cage-jump motion to macroscopic diffusion in supercooled liquids

The evaluation of the long term stability of a material requires the estimation of its long-time dynamics. For amorphous materials such as structural glasses, it has proven difficult to predict the long-time dynamics starting from static measurements. Here we consider how long one needs to monitor the dynamics of a structural glass to predict its long-time features. We present a detailed characterization of the statistical features of the single-particle intermittent motion, and show that single-particle jumps are the irreversible events leading to the relaxation of the system. This allows us to evaluate the diffusion constant on the time-scale of the jump duration, which is small and temperature independent, i.e. well before the system enters the diffusive regime. The prediction is obtained by analyzing the particle trajectories via a parameter-free algorithm.

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.

Analytical calculation of four-point correlations for a simple model of cages involving numerous particles

Dynamics of a one-dimensional system of Brownian particles with short-range repulsive interaction (diameter σ) is studied with a liquid-theoretical approach. The mean square displacement, the two-particle displacement correlation, and the overlap-density-based generalized susceptibility are calculated analytically by way of the Lagrangian correlation of the interparticulate space, instead of the Eulerian correlation of density that is commonly used in the standard mode-coupling theory. In regard to the mean square displacement, the linear analysis reproduces the established result on the asymptotic subdiffusive behavior of the system. A finite-time correction is given by incorporating the effect of entropic nonlinearity with a Lagrangian version of mode-coupling theory. The notorious difficulty in derivation of the mode-coupling theory concerning violation of the fluctuation-dissipation theorem is found to disappear by virtue of the Lagrangian description. The Lagrangian description also facilitates analytical calculation of four-point correlations in the space-time, such as the two-particle displacement correlation. The two-particle displacement correlation, which is asymptotically self-similar in the space-time, illustrates how the cage effect confines each particle within a short radius on one hand and creates collective motion of numerous particles on the other hand. As the time elapses, the correlation length grows unlimitedly, and the generalized susceptibility based on the overlap density converges to a finite value which is an increasing function of the density. The distribution function behind these dynamical four-point correlations and its extension to three-dimensional cases, respecting the tensorial character of the two-particle displacement correlation, are also discussed.

Spatiotemporal heterogeneity of local free volumes in highly supercooled liquid

We discuss the spatiotemporal behavior of local density and its relation to dynamical heterogeneity in a highly supercooled liquid by using molecular dynamics simulations of a binary mixture with different particle sizes in two dimensions. To trace voids heterogeneously existing with lower local densities, which move along with the structural relaxation, we employ the minimum local density for each particle in a time window whose width is set along with the structural relaxation time. Particles subject to free volumes correspond well to the configuration rearranging region of dynamical heterogeneity. While the correlation length for dynamical heterogeneity grows with temperature decrease, no growth in the correlation length of heterogeneity in the minimum local density distribution takes place. A comparison of these results with those of normal mode analysis reveals that superpositions of lower-frequency soft modes extending over the free volumes exhibit spatial correlation with the broken bonds. This observation suggests a possibility that long-ranged vibration modes facilitate the interactions between fragile regions represented by free volumes, to induce dynamical correlations at a large scale.

Formation of double glass in binary mixtures of anisotropic particles: dynamic heterogeneities in rotations and displacements

We study glass behavior in a mixture of elliptic and circular particles in two dimensions at low temperatures using an orientation-dependent Lennard-Jones potential. The ellipses have a mild aspect ratio (∼1.2) and tend to align at low temperatures, while the circular particles play the role of impurities disturbing the ellipse orientations at a concentration of 20%. These impurities have a size smaller than that of the ellipses and attract them in the homeotropic alignment. As a result, the coordination number around each impurity is mostly 5 or 4 in glassy states. We realize double glass, where both the orientations and the positions are disordered but still hold mesoscopic order. We find a strong heterogeneity in the flip motions of the ellipses, which sensitively depends on the impurity clustering. In our model, a small fraction of the ellipses still undergo flip motions relatively rapidly even at low temperatures. In contrast, the nonflip rotations (with angle changes not close to ±π) are mainly caused by the cooperative configuration changes involving many particles. Then, there arises a long-time heterogeneity in the nonflip rotations closely correlated with the dynamic heterogeneity in displacements.