Relationship between bond-breakage correlations and four-point correlations in heterogeneous glassy dynamics: configuration changes and vibration modes

Emerging Mesoscale Flows and Chaotic Advection in Dense Active Matter

We study two models of overdamped self-propelled disks in two dimensions, with and without aligning interactions. Both models support active mesoscale flows, leading to chaotic advection and transport over large length scales in their homogeneous dense fluid states, away from dynamical arrest. They form streams and vortices reminiscent of multiscale flow patterns in turbulence. We show that the characteristics of these flows do not depend on the specific details of the active fluids, and result from the competition between crowding effects and persistent propulsions. This observation suggests that dense active suspensions of self-propelled particles present a type of "active turbulence" distinct from collective flows reported in other types of active systems.

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.

Johari-Goldstein β relaxation in glassy dynamics originates from two-scale energy landscape

Supercooled liquids undergo complicated structural relaxation processes, which have been a long-standing problem in both experimental and theoretical aspects of condensed matter physics. In particular, past experiments widely observed for many types of molecular liquids that relaxation dynamics separated into two distinct processes at low temperatures. One of the possible interpretations is that this separation originates from the two-scale hierarchical topography of the potential energy landscape; however, it has never been verified. Molecular dynamics simulations are a promising approach to tackle this issue, but we must overcome laborious difficulties. First, we must handle a model of molecular liquids that is computationally demanding compared to simple spherical models, which have been intensively studied but show only a slower process: α relaxation. Second, we must reach a sufficiently low-temperature regime where the two processes become well-separated. Here, we handle an asymmetric dimer system that exhibits a faster process: Johari-Goldstein β relaxation. Then, we employ the parallel tempering method to access the low-temperature regime. These laborious efforts enable us to investigate the potential energy landscape in detail and unveil the first direct evidence of the topographic hierarchy that induces the β relaxation. We also successfully characterize the microscopic motions of particles during each relaxation process. Finally, we study the correlation between low-frequency modes and two relaxation processes. Our results establish a fundamental and comprehensive understanding of experimentally observed relaxation dynamics in supercooled liquids.

BOTAN: BOnd TArgeting Network for prediction of slow glassy dynamics by machine learning relative motion

Recent developments in machine learning have enabled accurate predictions of the dynamics of slow structural relaxation in glass-forming systems. However, existing machine learning models for these tasks are mostly designed such that they learn a single dynamic quantity and relate it to the structural features of glassy liquids. In this study, we propose a graph neural network model, "BOnd TArgeting Network," that learns relative motion between neighboring pairs of particles, in addition to the self-motion of particles. By relating the structural features to these two different dynamical variables, the model autonomously acquires the ability to discern how the self motion of particles undergoing slow relaxation is affected by different dynamical processes, strain fluctuations and particle rearrangements, and thus can predict with high precision how slow structural relaxation develops in space and time.

Disordered Collective Motion in Dense Assemblies of Persistent Particles

We explore the emergence of nonequilibrium collective motion in disordered nonthermal active matter when persistent motion and crowding effects compete, using simulations of a two-dimensional model of size polydisperse self-propelled particles. In stark contrast with monodisperse systems, we find that polydispersity stabilizes a homogeneous active liquid at arbitrary large persistence times, characterized by remarkable velocity correlations and irregular turbulent flows. For all persistence values, the active fluid undergoes a nonequilibrium glass transition at large density. This is accompanied by collective motion, whose nature evolves from near-equilibrium spatially heterogeneous dynamics at small persistence, to a qualitatively different intermittent dynamics when persistence is large. This latter regime involves a complex time evolution of the correlated displacement field.

Comparing Microscopic and Macroscopic Dynamics in a Paradigmatic Model of Glass-Forming Molecular Liquid

Glass transition is a most intriguing and long-standing open issue in the field of molecular liquids. From a macroscopic perspective, glass-forming systems display a dramatic slowing-down of the dynamics, with the inverse diffusion coefficient and the structural relaxation times increasing by orders of magnitude upon even modest supercooling. At the microscopic level, single-molecule motion becomes strongly intermittent, and can be conveniently described in terms of “cage-jump” events. In this work, we investigate a paradigmatic glass-forming liquid, the Kob–Andersen Lennard–Jones model, by means of Molecular Dynamics simulations, and compare the macroscopic and microscopic descriptions of its dynamics on approaching the glass-transition. We find that clear changes in the relations between macroscopic timescales and cage-jump quantities occur at the crossover temperature where Mode Coupling-like description starts failing. In fact, Continuous Time Random Walk and lattice model predictions based on cage-jump statistics are also violated below the crossover temperature, suggesting the onset of a qualitative change in cage-jump motion. Interestingly, we show that a fully microscopic relation linking cage-jump time- and length-scales instead holds throughout the investigated temperature range.

Breakdown of the Stokes-Einstein relation in supercooled liquids: A cage-jump perspective

The breakdown of the Stokes-Einstein relation in supercooled liquids, which is the increase in the ratio τ_ατ_D between the two macroscopic times for structural relaxation and diffusion on decreasing the temperature, is commonly ascribed to dynamic heterogeneities, but a clear-cut microscopic interpretation is still lacking. Here, we tackle this issue exploiting the single-particle cage-jump framework to analyze molecular dynamics simulations of soft disk assemblies and supercooled water. We find that τ_ατ_D∝⟨t_p⟩⟨t_c⟩, where ⟨t_p⟩ and ⟨t_c⟩ are the cage-jump times characterizing slow and fast particles, respectively. We further clarify that this scaling does not arise from a simple term-by-term proportionality; rather, the relations τ_α∝⟨t_p⟩⟨Δr_J ²⟩ and τ_D∝⟨t_c⟩⟨Δr_J ²⟩ effectively connect the macroscopic and microscopic timescales, with the mean square jump length ⟨Δr_J ²⟩ shrinking on cooling. Our work provides a microscopic perspective on the Stokes-Einstein breakdown and generalizes previous results on lattice models to the case of more realistic glass-formers.

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.

Viscoelastic shear stress relaxation in two-dimensional glass-forming liquids

Translational dynamics of 2D glass-forming fluids is strongly influenced by soft, long-wavelength fluctuations first recognized by D. Mermin and H. Wagner. As a result of these fluctuations, characteristic features of glassy dynamics, such as plateaus in the mean-squared displacement and the self-intermediate scattering function, are absent in two dimensions. In contrast, Mermin-Wagner fluctuations do not influence orientational relaxation, and well-developed plateaus are observed in orientational correlation functions. It has been suggested that, by monitoring translational motion of particles relative to that of their neighbors, one can recover characteristic features of glassy dynamics and thus disentangle the Mermin-Wagner fluctuations from the 2D glass transition. Here we use molecular dynamics simulations to study viscoelastic relaxation in two and three dimensions. We find different behavior of the dynamic modulus below the onset of slow dynamics (determined by the orientational or cage-relative correlation functions) in two and three dimensions. The dynamic modulus for 2D supercooled fluids is more stretched than for 3D supercooled fluids and does not exhibit a plateau, which implies the absence of glassy viscoelastic relaxation. At lower temperatures, the 2D dynamic modulus starts exhibiting an intermediate time plateau and decays similarly to the 2D dynamic modulus. The differences in the glassy behavior of 2D and 3D glass-forming fluids parallel differences in the ordering scenarios in two and three dimensions.

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.

Local structure in deeply supercooled liquids exhibits growing lengthscales and dynamical correlations

Glasses are among the most widely used of everyday materials, yet the process by which a liquid's viscosity increases by 14 decades to become a glass remains unclear, as often contradictory theories provide equally good descriptions of the available data. Knowledge of emergent lengthscales and higher-order structure could help resolve this, but this requires time-resolved measurements of dense particle coordinates-previously only obtained over a limited time interval. Here we present an experimental study of a model colloidal system over a dynamic window significantly larger than previous measurements, revealing structural ordering more strongly linked to dynamics than previously found. Furthermore we find that immobile regions and domains of local structure grow concurrently with density, and that these regions have low configurational entropy. We thus show that local structure plays an important role at deep supercooling, consistent with a thermodynamic interpretation of the glass transition rather than a principally dynamic description.

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.

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.

Long-wavelength fluctuations and the glass transition in two dimensions and three dimensions

Phase transitions significantly differ between 2D and 3D systems, but the influence of dimensionality on the glass transition is unresolved. We use microscopy to study colloidal systems as they approach their glass transitions at high concentrations and find differences between two dimensions and three dimensions. We find that, in two dimensions, particles can undergo large displacements without changing their position relative to their neighbors, in contrast with three dimensions. This is related to Mermin-Wagner long-wavelength fluctuations that influence phase transitions in two dimensions. However, when measuring particle motion only relative to their neighbors, two dimensions and three dimensions have similar behavior as the glass transition is approached, showing that the long-wavelength fluctuations do not cause a fundamental distinction between 2D and 3D glass transitions.

Mermin-Wagner fluctuations in 2D amorphous solids

In a recent commentary, J. M. Kosterlitz described how D. Thouless and he got motivated to investigate melting and suprafluidity in two dimensions [Kosterlitz JM (2016) J Phys Condens Matter 28:481001]. It was due to the lack of broken translational symmetry in two dimensions-doubting the existence of 2D crystals-and the first computer simulations foretelling 2D crystals (at least in tiny systems). The lack of broken symmetries proposed by D. Mermin and H. Wagner is caused by long wavelength density fluctuations. Those fluctuations do not only have structural impact, but additionally a dynamical one: They cause the Lindemann criterion to fail in 2D in the sense that the mean squared displacement of atoms is not limited. Comparing experimental data from 3D and 2D amorphous solids with 2D crystals, we disentangle Mermin-Wagner fluctuations from glassy structural relaxations. Furthermore, we demonstrate with computer simulations the logarithmic increase of displacements with system size: Periodicity is not a requirement for Mermin-Wagner fluctuations, which conserve the homogeneity of space on long scales.

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.

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.

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.

Particle jumps in structural glasses

Particles in structural glasses rattle around temporary equilibrium positions, that seldom change through a process which is much faster than the relaxation time, known as particle jump. Since the relaxation of the system is due to the accumulation of many such jumps, it could be possible to connect the single particle short time motion to the macroscopic relaxation by understanding the features of the jump dynamics. Here we review recent results in this research direction, clarifying the features of particle jumps that have been understood and those that are still under investigation, and examining the role of particle jumps in different theories of the glass transition.

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.

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.

Mutual information reveals multiple structural relaxation mechanisms in a model glass former

Among the key challenges to our understanding of solidification in the glass transition is that it is accompanied by little apparent change in structure. Recently, geometric motifs have been identified in glassy liquids, but a causal link between these motifs and solidification remains elusive. One 'smoking gun' for such a link would be identical scaling of structural and dynamic lengthscales on approaching the glass transition, but this is highly controversial. Here we introduce an information theoretic approach to determine correlations in displacement for particle relaxation encoded in the initial configuration of a glass-forming liquid. We uncover two populations of particles, one inclined to relax quickly, the other slowly. Each population is correlated with local density and geometric motifs. Our analysis further reveals a dynamic lengthscale similar to that associated with structural properties, which may resolve the discrepancy between structural and dynamic lengthscales.

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.

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.

Simple picture of supercooled liquid dynamics: dynamic scaling and phenomenology based on clusters

Although it is now well established that in glassy liquids, slow structural relaxation accompanies a correlated structural rearrangement, the role of such a correlation in the transport anomaly, and thus in the slow dynamics, remains unclear. In this paper, we argue from a hydrodynamic viewpoint that a correlated structure (cluster) with a characteristic size ξ sustains the long-lived stress and dynamically couples with the hydrodynamic fluctuations; therefore, the dynamics of this cluster is the origin of the mesoscopic nature of anomalous hydrodynamic transport. Based on this argument, we derive a dynamic scaling law for τ(α) (or η, where η is the macroscopic shear viscosity) as a function of ξ: τ(α)([proportionality]η)[proportionality]ξ(4). We provide a simple explanation for basic features of anomalous transport, such as the breakdown of the Stokes-Einstein relation and the length-scale-dependent decoupling between viscosity and diffusion. The present study further suggests a different physical picture: Through the coarse graining of smaller-scale fluctuations (</~ξ), the supercooled liquid dynamics can be regarded as the dynamics of normal (cluster) liquids composed of units with a typical size of ξ. Although the correlation length of hydrodynamic transport ξ and the dynamic heterogeneity size ξ(DH), which is determined by the usual four-point correlation function, reflect some aspects of the cooperative effects, the correspondence between ξ and ξ(DH) is not one to one. We highlight the possibility that ξ(DH) overestimates the actual collective transport range at a low degree of supercooling.

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.