Multiple length and time scales of dynamic heterogeneities in model glass-forming liquids: A systematic analysis of multi-point and multi-time correlations

Selecting relevant structural features for glassy dynamics by information imbalance

We numerically investigate the identification of relevant structural features that contribute to the dynamical heterogeneity in a model glass-forming liquid. By employing the recently proposed information imbalance technique, we select these features from a range of physically motivated descriptors. This selection process is performed in a supervised manner (using both dynamical and structural data) and an unsupervised manner (using only structural data). We then apply the selected features to predict future dynamics using a machine learning technique. One of the advantages of the information imbalance technique is that it does not assume any model a priori, i.e., it is a non-parametric method. Finally, we discuss the potential applications of this approach in identifying the dominant mechanisms governing the glassy slow dynamics.

Glassy dynamics in a liquid of anisotropic molecules: Bifurcation of relaxation spectrum

In experimental and theoretical studies of glass transition phenomena, one often finds a sharp crossover in dynamical properties at a temperature Tcr. A bifurcation of a relaxation spectrum is also observed at a temperature TB≈Tcr; both lie significantly above the glass transition temperature. In order to better understand these phenomena, we introduce a new model of glass-forming liquids, a binary mixture of prolate and oblate ellipsoids. This model system exhibits sharp thermodynamic and dynamic anomalies, such as the specific heat jump during heating and a sharp variation in the thermal expansion coefficient around a temperature identified as the glass transition temperature, Tg. The same temperature is obtained from the fit of the calculated relaxation times to the Vogel-Fulcher-Tammann (VFT) form. As the temperature is lowered, the calculated single peak rotational relaxation spectrum splits into two peaks at TB above the estimated Tg. Similar bifurcation is also observed in the distribution of short-to-intermediate time translational diffusion. Interrogation of the two peaks reveals a lower extent of dynamic heterogeneity in the population of the faster mode. We observe an unexpected appearance of a sharp peak in the product of rotational relaxation time τ2 and diffusion constant D at a temperature Tcr, close to TB, but above the glass transition temperature. Additionally, we coarse-grain the system into cubic boxes, each containing, on average, ∼62 particles, to study the average dynamical properties. Clear evidence of large-scale sudden changes in the diffusion coefficient and rotational correlation time signals first-order transitions between low and high-mobility domains.

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

Spectroscopic and Chemical Properties of Ionic Liquids: Computational Study

A brief account is given of highlights of our computational efforts - often in collaboration with experimental groups - to understand spectroscopic and chemical properties of ionic liquids (ILs). Molecular dynamics, including their inhomogeneous character, responsible for key spectral features observed in dielectric absorption, infra-red (IR) and fluorescence correlation spectroscopy (FCS) measurements are elucidated. Mechanisms of chemical processes involving imidazolium-based ILs are illustrated for CO2 capture and related reactions, transesterification of cellulose, and Au nanocluster-catalyzed Suzuki cross-coupling reaction with attention paid to differing roles of IL ions. A comparison with experiments is also made.

Dynamical susceptibilities near ideal glass transitions

Building on the recently derived inhomogeneous mode-coupling theory, we extend the generalized mode-coupling theory of supercooled liquids to inhomogeneous environments. This provides a first-principles-based, systematic, and rigorous way of deriving high-point dynamical susceptibilities from variations of the many-body dynamic structure factors with respect to their conjugate field. This framework allows for a fully microscopic possibility to probe for collective relaxation mechanisms in supercooled liquids near the mode-coupling glass transition. The behavior of these dynamical susceptibilities is then studied in the context of simplified self-consistent relaxation models.

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.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Observation of Collective Molecular Dynamics in a Chalcogenide Glass: Results from X-ray Photon Correlation Spectroscopy

The structural relaxation processes in a Ge₃As₅₂S₄₅ molecular chalcogenide glass sample were directly studied by X-ray photon correlation spectroscopy (XPCS). XPCS was conducted at the first sharp diffraction peak at q = 1.16 Å^-1 at temperatures ranging from 123 K to above the glass transition at 328 K, and the results showed two different dynamical regimes. At a low temperature, the observed glass dynamics are slow and dominated by X-ray-photon-induced effects, which are temperature independent. At a higher temperature, we observed a dramatic decrease in the fluctuation timescales, indicating that the dynamics were mainly due to the intermolecular correlation of the As₄S₃ molecule in the glass. The timescales in this high-temperature range agree well with those determined from measurements of the Newtonian viscosity. Our XPCS studies suggest an extended length scale of the relaxation process in glassy Ge₃As₅₂S₄₅ from the single molecule to the intermolecular range across the glass transition, providing a unique direct probe of the dynamics beyond the length scales of the individual molecule.

Thermodynamics and its correlation with dynamics in a mean-field model and pinned systems: A comparative study using two different methods of entropy calculation

A recent study introduced a novel mean-field model system where each particle over and above the interaction with its regular neighbors interacts with k extra pseudo-neighbors. Here, we present an extensive study of thermodynamics and its correlation with the dynamics of this system. We surprisingly find that the well-known thermodynamic integration (TI) method of calculating the entropy provides unphysical results. It predicts vanishing of the configurational entropy at temperatures close to the onset temperature of the system and negative values of the configurational entropy at lower temperatures. Interestingly, well below the temperature at which the configurational entropy vanishes, both the collective and the single-particle dynamics of the system show complete relaxation. Negative values of the configurational entropy are unphysical, and complete relaxation when the configurational entropy is zero violates the prediction of the random first-order transition theory (RFOT). However, the entropy calculated using the two-phase thermodynamics (2PT) method remains positive at all temperatures for which we can equilibrate the system, and its values are consistent with RFOT predictions. We find that with an increase in k, the difference in the entropy computed using the two methods increases. A similar effect is also observed for a system where a randomly selected fraction of the particles are pinned in their positions in the equilibrated liquid. We show that the difference in entropy calculated via the 2PT and TI methods increases with pinning density.

Nonlinear measurements of kinetics and generalized dynamical modes. I. Extracting the one-dimensional Green's function from a time series

Often, a single correlation function is used to measure the kinetics of a complex system. In contrast, a large set of k-vector modes and their correlation functions are commonly defined for motion in free space. This set can be transformed to the van Hove correlation function, which is the Green's function for molecular diffusion. Here, these ideas are generalized to other observables. A set of correlation functions of nonlinear functions of an observable is used to extract the corresponding Green's function. Although this paper focuses on nonlinear correlation functions of an equilibrium time series, the results are directly connected to other types of nonlinear kinetics, including perturbation-response experiments with strong fields. Generalized modes are defined as the orthogonal polynomials associated with the equilibrium distribution. A matrix of mode-correlation functions can be transformed to the complete, single-time-interval (1D) Green's function. Diagonalizing this matrix finds the eigendecays. To understand the advantages and limitation of this approach, Green's functions are calculated for a number of models of complex dynamics within a Gaussian probability distribution. Examples of non-diffusive motion, rate heterogeneity, and range heterogeneity are examined. General arguments are made that a full set of nonlinear 1D measurements is necessary to extract all the information available in a time series. However, when a process is neither dynamically Gaussian nor Markovian, they are not sufficient. In those cases, additional multidimensional measurements are needed.

Connecting real glasses to mean-field models

We propose a novel model for a glass-forming liquid, which allows us to switch in a continuous manner from a standard three-dimensional liquid to a fully connected mean-field model. This is achieved by introducing k additional particle-particle interactions, which thus augments the effective number of neighbors of each particle. Our computer simulations of this system show that the structure of the liquid does not change with the introduction of these pseudo-neighbors and by means of analytical calculations, and we determine the structural properties related to these additional neighbors. We show that the relaxation dynamics of the system slows down very quickly with the increase in k and that the onset and the mode-coupling temperatures increase. The systems with high values of k follow the mode-coupling theory power law behavior for a larger temperature range compared to the ones with lower values of k. The dynamic susceptibility indicates that the dynamic heterogeneity decreases with the increase in k, whereas the non-Gaussian parameter is independent of it. Thus, we conclude that with the increase in the number of pseudo-neighbors, the system becomes more mean-field-like. By comparing our results with previous studies on mean-field-like systems, we come to the conclusion that the details of how the mean-field limit is approached are important since they can lead to different dynamical behavior in this limit.

Translational and rotational critical-like behaviors in the glass transition of colloidal ellipsoid monolayers

Critical-like behaviors have been found in translational degrees of freedom near the glass transition of spherical particle systems mainly with local polycrystalline structures, but it is not clear if criticality exists in more general glassy systems composed of nonspherical particles without crystalline structures. Here, through experiments and simulations, we show critical-like behaviors in both translational and rotational degrees of freedom in monolayers of monodisperse colloidal ellipsoids in the absence of crystalline orders. We find rich features of the Ising-like criticality in structure and slow dynamics at the ideal glass transition point ϕ₀, showing the thermodynamic nature of glass transition at ϕ₀ A dynamic criticality is found at the mode-coupling critical point ϕ_c for the fast-moving clusters whose critical exponents increase linearly with fragility, reflecting a dynamic glass transition. These results cast light on the glass transition and explain the mystery that the dynamic correlation lengths diverge at two different temperatures.

Coarse-graining strategy for modeling effective, highly diffusive fluids with reduced polydispersity: A dynamical study

We present a coarse-graining strategy for reducing the number of particle species in mixtures to achieve a simpler system with higher diffusion while preserving the total particle number and characteristic dynamic features. As a system of application, we chose the bidisperse Lennard-Jones-like mixture, discovered by Kob and Andersen [Phys. Rev. Lett. 73, 1376 (1994)], possessing a slow dynamics due to the fluid's multi-component character with its apparently unconventional choice for the pair potential of the type-A-type-B arrangement. We further established in a so-formed coarse-grained and temperature-independent monodisperse system an equilibrium structure with a radial distribution function resembling its mixture counterpart. This one-component system further possesses similar dynamic features such as glass transition temperature and critical exponents while subjected to Newtonian mechanics. This strategy may finally lead to the manufacturing of new nanoparticle/colloidal fluids by experimentally modeling only the outcoming effective pair potential(s) and no other macroscopic quantity.

Micelle Heterogeneity from the 2D Kinetics of Solute Rotation

The chemical and physical properties of microstructured materials vary with position. The photophysics of solute molecules can measure these local properties, but they often show multiple rates (rate dispersion), which complicates the interpretation. In the case of micelles, rate dispersion in a solute's anisotropy decay has been assigned to either local anisotropy or heterogeneity in the local viscosity. To resolve this conflict, the rotation of PM597 molecules in SDS micelles has been measured by polarized MUPPETS (multiple population-period transient spectroscopy). This 2D technique shows that heterogeneity is strong and that local anisotropy is minimal. The results suggest that on a subnanosecond time scale, the solute sees only one strong fluctuation of the micelle structure. The anisotropic, average structure emerges on longer time scales.

Surface Fluctuations Dominate the Slow Glassy Dynamics of Polymer-Grafted Colloid Assemblies

Colloids grafted with a corona layer of polymers show glassy behavior that covers a wide range of fragilities, with this behavior being tunable through variations in grafting density and grafting chain length. We find that the corona roughness, which is maximized for long chain lengths and sparse grafting, is directly correlated to the concentration-dependence of the system relaxation time (fragility). Relatively rougher colloids result in stronger liquids because their rotational motions become orientationally correlated across the whole system even at low particle loadings leading to an essentially Arrhenius-like concentration-dependence of the relaxation times near the glass transition. The smoother colloids do not show as much orientational correlation except at higher densities leading to fragile behavior. We therefore propose that these materials are an ideal model to study the physical properties of the glass transition.

Single molecule studies reveal temperature independence of lifetime of dynamic heterogeneity in polystyrene

Polymeric systems close to their glass transition temperature are known to exhibit heterogeneous dynamics that evolve both over time and space, comparable to the dynamics of small molecule glass formers. It remains unclear how temperature influences the degree of heterogeneous dynamics in such systems. In the following report, a fluorescent perylene dicarboximide probe molecule that reflects the full breadth of heterogeneity of the host was used to examine the temperature dependence of the dynamic heterogeneity lifetime in polystyrene at several temperatures ranging from the glass transition to 10 K above this temperature via single molecule microscopy. Contrary to prior reports, no apparent temperature dependence of time scales associated with dynamic heterogeneity was detected; indeed, the probe molecules report characteristic dynamic heterogeneity lifetimes 100-300 times the average alpha-relaxation time (τ_α) of the polystyrene host at all temperatures studied.

Structural-dynamical transition in the Wahnström mixture

In trajectory space, dynamical heterogeneities in glass-forming liquids correspond to the emergence of a dynamical phase transition between an active phase poor in local structure and an inactive phase which is rich in local structure. We support this scenario with the study of a model additive mixture of Lennard-Jones particles, quantifying how the choice of the relevant structural and dynamical observable affects the transition in trajectory space. We find that the low mobility, structure-rich phase is dominated by icosahedral order. Applying a non-equilibrium rheological protocol, we connect local order to the emergence of mechanical rigidity.

Revealing the Link between Structural Relaxation and Dynamic Heterogeneity in Glass-Forming Liquids

Despite the use of glasses for thousands of years, the nature of the glass transition is still mysterious. On approaching the glass transition, the growth of dynamic heterogeneity has long been thought to play a key role in explaining the abrupt slowdown of structural relaxation. However, it still remains elusive whether there is an underlying link between structural relaxation and dynamic heterogeneity. Here, we unravel the link by introducing a characteristic time scale hiding behind an identical dynamic heterogeneity for various model glass-forming liquids. We find that the time scale corresponds to the kinetic fragility of liquids. Moreover, it leads to scaling collapse of both the structural relaxation time and dynamic heterogeneity for all liquids studied, together with a characteristic temperature associated with the same dynamic heterogeneity. Our findings imply that studying the glass transition from the viewpoint of dynamic heterogeneity is more informative than expected.

Growing length scale accompanying vitrification: A perspective based on nonsingular density fluctuations

A model for describing growing length scale accompanying the vitrification is introduced, in which we assume that in a subsystem whose density is above a certain threshold value, ρ_{c}, due to steric constraints, particle rearrangements are highly suppressed for a sufficiently long time period (∼structural relaxation time). We regard such a subsystem as a glassy cluster. With this assumption and without introducing any complicated thermodynamic arguments, we predict that with compression (increasing average density ρ) at a fixed temperature T in supercooled states, the characteristic length of the clusters, ξ, diverges as ξ∼(ρ_{c}-ρ)^{-2/d}, where d is the spatial dimensionality. This ξ measures the average persistence length of the steric constraints in blocking the rearrangement motions and is determined by the subsystem density. Additionally, with decreasing T at a fixed ρ, the length scale diverges in the same manner as ξ∼(T-T_{c})^{-2/d}, for which ρ is identical to ρ_{c} at T=T_{c}. The exponent describing the diverging length scale is the same as the one predicted by some theoretical models and indeed has been observed in some simulations and experiments. However, the basic mechanism for this divergence is different; that is, we do not invoke thermodynamic anomalies associated with the thermodynamic phase transition as the origin of the growing length scale. We further present arguements for the cooperative properties of the structural relaxation based on the clusters.

Role of the Pair Correlation Function in the Dynamical Transition Predicted by Mode Coupling Theory

In a recent study, we have found that for a large number of systems the configurational entropy at the pair level S_{c2}, which is primarily determined by the pair correlation function, vanishes at the dynamical transition temperature T_{c}. Thus, it appears that the information of the transition temperature is embedded in the structure of the liquid. In order to investigate this, we describe the dynamics of the system at the mean field level and, using the concepts of the dynamical density functional theory, show that the dynamical transition temperature depends only on the pair correlation function. Thus, this theory is similar in spirit to the microscopic mode coupling theory (MCT). However, unlike microscopic MCT, which predicts a very high transition temperature, the present theory predicts a transition temperature that is similar to T_{c}. This implies that the information of the dynamical transition temperature is embedded in the pair correlation function.

Fractional Viscosity Dependence of Reaction Kinetics in Glass-Forming Liquids

The diffusion of molecules in complex systems such as glasses and cell cytoplasm is slow, heterogeneous, and sometimes nonergodic. The effects of such intriguing diffusion on the kinetics of chemical and biological reactions remain elusive. In this Letter, we report that the kinetics of the polymer loop formation reaction in a Kob-Andersen (KA) glass forming liquid is influenced significantly by the dynamic heterogeneity. The diffusion coefficient D of a KA liquid deviates from the Stokes-Einstein relation at low temperatures and D shows a fractional dependence on the solvent viscosity η_{s}, i.e., D∼η_{s}^{-ξ_{D}} with ξ_{D}=0.85. The dynamic heterogeneity of a KA liquid affects the rate constant k_{rxn} of the loop formation and leads to the identical fractional dependence of k_{rxn} on η_{s} with k_{rxn}∼η_{s}^{-ξ} and ξ=ξ_{D}, contrary to reactions in dynamically homogeneous solutions where k_{rxn}∼η_{s}^{-1}.

Probing heterogeneous dynamics from spatial density correlation in glass-forming liquids

We numerically investigate the connection between spatial density correlation and dynamical heterogeneity in glass-forming liquids. We demonstrate that the cluster size defined by the spatial aggregation of densely packed particles (DPPs) can better capture the difference between the dynamics of the Lennard-Jones glass model and the Weeks-Chandler-Andersen truncation model than the commonly used pair correlation functions. More interestingly, we compare the mobility of DPPs and loosely packed particles, and we find that high local density correlates well with slow dynamics in systems with relatively hard repulsive interactions but links to mobile ones in the system with soft repulsive interactions at one relaxation time scale. Our results show clear evidence that the above model dependence behavior stems from the hopping motion of DPPs at the end of the caging stage due to the compressive nature of soft repulsive spheres, which activates the dynamics of DPPs in the α relaxation stage.

Significant difference in the dynamics between strong and fragile glass formers

Glass-forming liquids are often classified into strong glass formers with nearly Arrhenius behavior and fragile ones with super-Arrhenius behavior. We reveal a significant difference in the dynamics between these two types of glass formers through molecular dynamics simulations: In strong glass formers, the relaxation dynamics of density fluctuations is nondiffusive, whereas in fragile glass formers it exhibits diffusive behavior. We demonstrate that this distinction is a direct consequence of the fundamental difference in the underlying elementary relaxation process between these two dynamical classes of glass formers. For fragile glass formers, a density-exchange process proceeds the density relaxation, which takes place locally at the particle level in normal states but is increasingly cooperative and nonlocal as the temperature is lowered in supercooled states. On the other hand, in strong glass formers, such an exchange process is not necessary for density relaxation due to the presence of other local relaxation channels. Our finding provides a novel insight into Angell's classification scheme from a hydrodynamic perspective.

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.

How Heterogeneous Are Trehalose/Glycerol Cryoprotectant Mixtures? A Combined Time-Resolved Fluorescence and Computer Simulation Investigation

Heterogeneity and molecular motions in representative cryoprotectant mixtures made of trehalose and glycerol are investigated in the temperature range 298 ≤ T (K) ≤ 353, via time-resolved fluorescence Stokes shift and anisotropy measurements, and molecular dynamics simulations of four-point density-time correlations and H-bond relaxations. Mixtures containing 5 and 20 wt % of trehalose along with neat glycerol are studied. Viscosity coefficients for these systems lie in the range 0.30 < η (P) < 23. Measured solute (Coumarin 153) rotation and solvation times reveal a substantial departure from the hydrodynamic viscosity dependence, suggesting the strong microheterogeneous nature of these systems. Fluorescence anisotropy decays are highly nonexponential, reflecting a non-Markovian character of the medium friction. A complete missing of the Stokes shift dynamics in these systems at 298 K but partial detection of it at other higher temperatures (shift magnitude being ∼400-600 cm^-1) indicates rigid solute environments. An amorphous solid-like feature emerges in the simulated radial distribution functions at these temperatures. Analyses of mean squared displacements reveal rattling-in-a-cage motion, non-Gaussian displacement distributions, and strong dynamic heterogeneity features. Simulated dynamic structure factors and four-point correlations hint, respectively, at very long α-relaxation and correlated time scales at 298 K. This explains the long solute rotation times (∼80-200 ns) measured at 298 K. Stretched exponential decay of the simulated H-bond relaxations with long time scales further highlights the strong temporal heterogeneity and slow dynamics inherent to these systems. In summary, this work provides the first insight into the molecular motions and interspecies interaction in a representative cryoprotectant mixture, and stimulates further study to investigate the interconnection between cryoprotection and dynamic heterogeneity.

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.

Heterogeneous dynamics and its length scale in simple ionic liquid models: a computational study

We numerically investigate the dynamic heterogeneity and its length scale found in coarse-grained ionic liquid model systems.

Microscopic origin of temporal heterogeneities in translational dynamics of liquid water

Liquid water is known to reorient via a combination of large angular jumps (due to exchange of hydrogen bonding (H-bond) partners) and diffusive orientations. Translation of the molecule undergoing the orientational jump and its initial and final H-bond acceptors plays a key role in the microscopic reorientation process. Here, we partition the translational dynamics into those occurring during intervals when rotating water molecules (and their initial and final H-bonding partners) undergo orientational jump and those arising when molecules wait between consecutive orientational jumps. These intervals are chosen in such a way that none of the four possible H-bonds involving the chosen water molecule undergo an exchange process within its duration. Translational dynamics is analysed in terms of the distribution of particle displacements, van Hove functions, and its moments. We observe that the translational dynamics, calculated from molecular dynamics simulations of liquid water, is fastest during the orientational jumps and slowest during periods of waiting. The translational dynamics during all temporal intervals shows an intermediate behaviour. This is the microscopic origin of temporal dynamic heterogeneity in liquid water, which is mild at 300 K and systematically increases with supercooling. Study of such partitioned dynamics in supercooled water shows increased disparity in dynamics occurring in the two different types of intervals. Nature of the distribution of particle displacements in supercooled water is investigated and it reveals signatures non-Gaussian behaviour.

Universal rescaling of flow curves for yield-stress fluids close to jamming

The experimental flow curves of four different yield-stress fluids with different interparticle interactions are studied near the jamming concentration. By appropriate scaling with the distance to jamming all rheology data can be collapsed onto master curves below and above jamming that meet in the shear-thinning regime and satisfy the Herschel-Bulkley and Cross equations, respectively. In spite of differing interactions in the different systems, master curves characterized by universal scaling exponents are found for the four systems. A two-state microscopic theory of heterogeneous dynamics is presented to rationalize the observed transition from Herschel-Bulkley to Cross behavior and to connect the rheological exponents to microscopic exponents for the divergence of the length and time scales of the heterogeneous dynamics. The experimental data and the microscopic theory are compared with much of the available literature data for yield-stress systems.

Dynamic heterogeneity in crossover spin facilitated model of supercooled liquid and fractional Stokes-Einstein relation

Kinetically constrained models have gained much interest as models that assign the origins of interesting dynamic properties of supercooled liquids to dynamical facilitation mechanisms that have been revealed in many experiments and numerical simulations. In this work, we investigate the dynamic heterogeneity in the fragile-to-strong liquid via Monte Carlo method using the model that linearly interpolates between the strong liquid-like behavior and the fragile liquid-like behavior by an asymmetry parameter b. When the asymmetry parameter is sufficiently small, smooth fragile-to-strong transition is observed both in the relaxation time and the diffusion constant. Using these physical quantities, we investigate fractional Stokes-Einstein relations observed in this model. When b is fixed, the system shows constant power law exponent under the temperature change, and the exponent has the value between that of the Frederickson-Andersen model and the East model. Furthermore, we investigate the dynamic length scale of our systems and also find the crossover relation between the relaxation time. We ascribe the competition between energetically favored symmetric relaxation mechanism and entropically favored asymmetric relaxation mechanism to the fragile-to-strong crossover behavior.

Structural heterogeneities at the origin of acoustic and transport anomalies in glycerol glass-former

By means of large scale molecular dynamics simulations, we explore mesoscopic properties of prototypical glycerol glass-former above and below the glass transition. The model used, in excellent agreement with various experimental techniques, permits to carefully study the structure and the vibrational dynamics. We find that a medium range order is present in glycerol glass-former and arises from hydrogen bond network extension. The characteristic size of the structural heterogeneities is related to the anomalous properties of acoustic vibrations (Rayleigh scattering, "mode softening," and Boson Peak) in the glassy state. Finally the characteristic size of these heterogeneities, nearly constant in temperature, is also connected to the cross-over between structural relaxation and diffusion in liquid glycerol.

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.

Time scale of dynamic heterogeneity in model ionic liquids and its relation to static length scale and charge distribution

We find a general power-law behavior: , where ζ_dh ≈ 1.2 for all the ionic liquid models, regardless of charges and the length scale of structural relaxation.

Correlation of local order with particle mobility in supercooled liquids is highly system dependent

We investigate the connection between local structure and dynamical heterogeneity in supercooled liquids. Through the study of four different models, we show that the correlation between a particle's mobility and the degree of local order in nearby regions is highly system dependent. Our results suggest that the correlation between local structure and dynamics is weak or absent in systems that conform well to the mean-field picture of glassy dynamics and strong in those that deviate from this paradigm. Finally, we investigate the role of order-agnostic point-to-set correlations and reveal that they provide similar information content to local structure measures, at least in the system where local order is most pronounced.