Particle jumps in structural glasses

Unraveling the dynamic slowdown in supercooled water: The role of dynamic disorder in jump motions

When a liquid is rapidly cooled below its melting point without inducing crystallization, its dynamics slow down significantly without noticeable structural changes. Elucidating the origin of this slowdown has been a long-standing challenge. Here, we report a theoretical investigation into the mechanism of the dynamic slowdown in supercooled water, a ubiquitous yet extraordinary substance characterized by various anomalous properties arising from local density fluctuations. Using molecular dynamics simulations, we found that the jump dynamics, which are elementary structural change processes, deviate from Poisson statistics with decreasing temperature. This deviation is attributed to slow variables competing with the jump motions, i.e., dynamic disorder. The present analysis of the dynamic disorder showed that the primary slow variable is the displacement of the fourth nearest oxygen atom of a jumping molecule, which occurs in an environment created by the fluctuations of molecules outside the first hydration shell. As the temperature decreases, the jump dynamics become slow and intermittent. These intermittent dynamics are attributed to the prolonged trapping of jumping molecules within extended and stable low-density domains. As the temperature continues to decrease, the number of slow variables increases due to the increased cooperative motions. Consequently, the jump dynamics proceed in a higher-dimensional space consisting of multiple slow variables, becoming slower and more intermittent. It is then conceivable that with further decreasing temperature, the slowing and intermittency of the jump dynamics intensify, eventually culminating in a glass transition.

The influence of molecular shape on glass-forming behavior in a minimalist trimer model

In this study, we employed molecular dynamics simulations to probe the influence of molecular morphological changes on the dynamic behavior of a model consisting of trimer molecules. This model, comprising a chain of three particles, facilitates the exploration of variations in the internal angle between these particles. Our findings highlight the significant impact of molecular conformation: systems with more linear conformations, characterized by larger internal angles, exhibit relaxation times several orders of magnitude greater than their counterparts with smaller internal angles. Furthermore, we delve into the role of angular interaction rigidity, uncovering a pronounced deceleration in dynamics and an increase in dynamic heterogeneity as rigidity escalates. This model not only provides insights into azobenzene-type systems but also sets the stage for subsequent research into the microscopic nuances of related systems, with potential extensions to composite systems.

Glassy phases of the Gaussian core model

We present results from molecular dynamics simulations exploring the supercooled dynamics of the Gaussian Core Model in the low- and intermediate-density regimes. In particular, we analyse the transition from the low-density hard-sphere-like glassy dynamics to the high-density one. The dynamics at low densities is well described by the caging mechanism, giving rise to intermittent dynamics. At high densities, the particles undergo a more continuous motion in which the concept of cage loses its meaning. We elaborate on the idea that these different supercooled dynamics are in fact the precursors of two different glass states.

Noncanonical Relationship between Heterogeneity and the Stokes-Einstein Breakdown in Deep Eutectic Solvents

The relationship between Stokes-Einstein breakdown (SEB) and dynamical heterogeneity (DH) is of paramount importance in the physical chemistry of complex fluids. In this work, we employ neutron scattering to probe the DH and SEB in a series of deep eutectic solvents (DESs) composed of acetamide and lithium salts. Quasielastic neutron scattering experiments reveal SEB in the jump diffusion of acetamide, represented by a fractional Stokes-Einstein relationship. Among these DESs, lithium perchlorate exhibits the most pronounced SEB while lithium bromide displays the weakest. Concurrently, elastic incoherent neutron scans identify that bromide DES is the most heterogeneous and perchlorate is the least. For the first time, our study unveils a counterintuitive incommensurate relationship between DH and SEB. Further, it reveals the intricate contrasting nature of the SEB-DH relationship when investigated in proximity to the glass-transition temperature and further away from it.

Local Plastic Response and Slow Heterogeneous Dynamics of Supercooled Liquids

We demonstrate, via numerical simulations, that the relaxation dynamics of supercooled liquids correlates well with a plastic length scale measuring a particle's response to impulsive localized perturbations and weakly to measures of local elasticity. We find that the particle averaged plastic length scale vanishes linearly in temperature and controls the super-Arrhenius temperature dependence of the relaxation time. Furthermore, we show that the plastic length scale of individual particles correlates with their typical displacement at the relaxation time. In contrast, the local elastic response only correlates with the dynamics on the vibrational timescale.

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.

Finite-size scaling and thermodynamics of model supercooled liquids: long-range concentration fluctuations and the role of attractive interactions

We compute partial structure factors, Kirkwood-Buff integrals (KBIs) and chemical potentials of model supercooled liquids with and without attractive interactions. We aim at investigating whether relatively small differences in the tail of the radial distribution functions result in contrasting thermodynamic properties. Our results suggest that the attractive potential favours the nucleation of long-range structures. Indeed, upon decreasing temperature, Bathia-Thornton structure factors display anomalous behaviour in the k→0 limit. KBIs extrapolated to the thermodynamic limit confirm this picture, and excess coordination numbers identify the anomaly with long-range concentration fluctuations. By contrast, the purely repulsive system remains perfectly miscible for the same temperature interval and only reveals qualitatively similar concentration fluctuations in the crystalline state. Furthermore, differences in both isothermal compressibilities and chemical potentials show that thermodynamics is not entirely governed by the short-range repulsive part of the interaction potential, emphasising the nonperturbative role of attractive interactions. Finally, at higher density, where both systems display nearly identical dynamical properties and repulsive interactions become dominant, the anomaly disappears, and both systems also exhibit similar thermodynamic properties.

Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective

Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.

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.

Nonlinear measurements of kinetics and generalized dynamical modes. II. Application to a simulation of solvation dynamics in an ionic liquid

Solvation dynamics in ionic liquids show features that are often associated with supercooled liquids, including "stretched" nonexponential relaxation. To better understand the mechanism behind the stretching, the nonlinear mode-correlation methods proposed in Paper I [S. R. Hodge and M. A. Berg, J. Chem. Phys. 155, 024122 (2021)] are applied to a simulation of a prototypical ionic liquid. A full Green's function is recovered. In addition, specific tests for non-Gaussian dynamics are made. No deviations from Gaussian dynamics are found. This finding is incompatible with rate heterogeneity as a cause of the nonexponential relaxation and appears to be in conflict with an earlier multidimensional analysis of the same data. Although this conflict is not resolved here, this work does demonstrate the practicality of mode-correlation analysis in the face of finite datasets and calculations.

Mean-Field Caging in a Random Lorentz Gas

The random Lorentz gas (RLG) is a minimal model of both percolation and glassiness, which leads to a paradox in the infinite-dimensional, d → ∞ limit: the localization transition is then expected to be continuous for the former and discontinuous for the latter. As a putative resolution, we have recently suggested that, as d increases, the behavior of the RLG converges to the glassy description and that percolation physics is recovered thanks to finite-d perturbative and nonperturbative (instantonic) corrections [Biroli et al. Phys. Rev. E 2021, 103, L030104]. Here, we expand on the d → ∞ physics by considering a simpler static solution as well as the dynamical solution of the RLG. Comparing the 1/d correction of this solution with numerical results reveals that even perturbative corrections fall out of reach of existing theoretical descriptions. Comparing the dynamical solution with the mode-coupling theory (MCT) results further reveals that, although key quantitative features of MCT are far off the mark, it does properly capture the discontinuous nature of the d → ∞ RLG. These insights help chart a path toward a complete description of finite-dimensional glasses.

Jump-precursor state emerges below the crossover temperature in supercooled o-terphenyl

In a supercooled liquid, the crossover temperature T_{c} separates a high-temperature region of diffusive dynamics from a low-temperature region of activated dynamics. A molecular-dynamics simulation of all-atom, flexible o-terphenyl [Eastwood et al., J. Phys. Chem. B 117, 12898 (2013)10.1021/jp402102w] is analyzed with advanced statistical methods to reveal the molecular features associated with this crossover. The simulations extend to an α-relaxation time of 14 μs (272.5 K), two orders of magnitude slower than at T_{c} (290 K). At T_{c} and below, a distinct state emerges that immediately precedes an orientational jump. Compared to the initial, tightly caged state, this jump-precursor state has a looser cage, with solid-angular excursions of 0.054-0.0125 × 4π sr. At T_{c} (290 K), rate heterogeneity is already the dominant cause of stretched relaxation. Exchange within the distribution of rates is faster than α relaxation at T_{c}, but becomes equal to it at the lowest temperature simulated (272.5 K). The results trend toward a recent experimental observation near the glass transition (243 K) [Kaur et al., Phys. Rev. E 98, 040603(R) (2018)10.1103/PhysRevE.98.040603], which saw exchange substantially slower than α relaxation. Overall, the dynamic crossover comprises multiple phenomena: the development of heterogeneity, an increasing jump size, an emerging jump-precursor state, and a lengthening exchange time. The crossover is neither sharp, nor a simple superposition of the high- and low-temperature regimes; it is a broad region that contains unique and complex phenomena.

Understanding charged vesicle suspensions as Wigner glasses: dynamical aspects

Suspensions of charged vesicles in water with added salt are widespread in nature and industrial production. Here we investigate, via Brownian dynamics simulations, a model that grasps the key features of these systems, with bidisperse colloidal beads interacting via a hard-core and an electrostatic double layer potential. Our goal is to focus on a set of interaction parameters that is not generic but measured in recent experiments, and relevant for a class of consumer products, such as liquid fabric softeners. On increasing the volume fraction in a range relevant to real formulation, we show that the dynamics become progressively slower and heterogeneous, displaying the typical signatures of an approaching glass transition. On lowering the salt concentration, which corresponds to increasing the strength and range of the electrostatic repulsion, the emergence of glassy dynamics becomes significantly steeper, and, remarkably, occurs at volume fractions well below the hard-sphere glass transition. The volume fraction dependence of the structural relaxation time at different salt concentration is well described through a functional law inspired by a recently proposed model (Krausser et al 2015 Proc. Natl Acad. Sci. USA 112 13762). According to our results, the investigated system may be thought of as a Wigner glass, i.e. a low-density glassy state stabilized by long-range repulsive interactions. Overall, our study suggests that glassy dynamics plays an important role in controlling the stability of these suspensions.

Connection of large amplitude angular jump motions with temporal heterogeneity in aqueous solutions

It has now been established that large angular jumps do take place when a rotating water molecule exchanges its hydrogen bond (H-bond) identity. This motion differs from the small angular diffusional steps occurring within short time intervals which define the 'Debye diffusion model' of water dynamics. We intend to investigate whether these two processes do eventually complement each other. In this present investigation the orientational dynamics of water in its mixture with a small hydrophobic molecule 1,2-dimethoxy ethane (DME) is studied microscopically using the all-atom classical molecular dynamics (MD) simulation technique. We found that the reorientational motions of water molecules are governed by continuous making and breaking of intermolecular H-bonds with their partners. We characterise these H-bond reorientation motions with the description of the "large amplitude angular jump model" and explore the coupling between the rotational and translational motions. By following the trajectories of each molecule in the solutions we describe the orientational dynamics of liquid water with a 'continuous time random walk' (CTRW) approach. Finally, we explore the diffusivity distribution through the jump properties of the water molecules, which successfully leads to the inherent transient heterogeneity of the solutions. We observe that the heterogeneity increases with increasing DME content in the mixtures. Our study correlates the coupling between rotational and translational motions of water molecules in the mixtures.

Packets of Diffusing Particles Exhibit Universal Exponential Tails

Brownian motion is a Gaussian process described by the central limit theorem. However, exponential decays of the positional probability density function P(X,t) of packets of spreading random walkers, were observed in numerous situations that include glasses, live cells, and bacteria suspensions. We show that such exponential behavior is generally valid in a large class of problems of transport in random media. By extending the large deviations approach for a continuous time random walk, we uncover a general universal behavior for the decay of the density. It is found that fluctuations in the number of steps of the random walker, performed at finite time, lead to exponential decay (with logarithmic corrections) of P(X,t). This universal behavior also holds for short times, a fact that makes experimental observations readily achievable.

Understanding the Origin of the Breakdown of the Stokes-Einstein Relation in Supercooled Water at Different Temperature-Pressure Conditions

A recent experiment has measured the viscosity of water down to approximately 244 K and up to 300 MPa. The correct viscosity and translational diffusivity data at various temperature-pressure (T-P) state points allowed for checking the validity of the Stokes-Einstein (SE) relation, which accounts for the coupling between translational self-diffusion and medium viscosity. The diffusion-viscosity decoupling increases with decreasing temperature, but the increasing pressure reduces the extent of the decoupling. Earlier simulation studies explained the breakdown of the SE relation in terms of the location of the Widom line, emanating from the liquid-liquid critical point (LLCP). Although these studies made a significant contribution to the current understanding of the above phenomena, a detailed molecular picture is still lacking. Recently, our group has explained the diffusion-viscosity decoupling from a jump-diffusion perspective. The jump-diffusion coefficient, emanating from the jump translation of water molecules, is calculated using a quantitative approach for different temperatures at ambient pressure. It has been observed that jump-diffusion is the key factor for diffusion-viscosity decoupling in supercooled water. The same method is adopted in the present work to estimate the jump-diffusion coefficient for different T-P state points and, thereby, explains the role of jump-diffusion for the different extents of the SE relation breakdown at different pressures. The residual diffusion coefficient, the other component of the total diffusion that originates from small step displacement and that is calculated by subtracting the jump-diffusion coefficient from the total diffusion, is seen to be fairly coupled to the viscosity at the entire range of temperature and pressure. Furthermore, we have calculated the average number of H-bonds per water molecule and the tetrahedral order for different T-P state points and investigated an approximate correlation between the average local structure and the contribution of the jump-diffusion to the total diffusion of water. This study, therefore, puts forward a new perspective for explaining the SE relation breakdown in supercooled water under different pressure conditions.

Glassy dynamics predicted by mutual role of free and activation volumes

Broadband Dielectric Spectroscopy (BDS) at elevated pressures and Positron Annihilation Lifetime Spectroscopy (PALS) are employed to elucidate the importance of the ratio of activation and free volumes during vitrification. We show that this ratio has a linear correlation with the structural relaxation of glass forming liquids in a wide temperature range hence engendering it as a vital input in the description of the dynamic glass transition.

Dynamics of dense hard sphere colloidal systems: A numerical analysis

The applicability to dense hard sphere colloidal suspensions of a general coarse-graining approach called Record Dynamics (RD) is tested by extensive molecular dynamics simulations. We reproduce known results as logarithmic diffusion and the logarithmic decay of the average potential energy per particle. We provide quantitative measures for the cage size and identify the displacements of single particles corresponding to intermittent cage breakings. We then partition the system into spatial domains and show that, within each domain, a subset of such intermittent events called quakes constitutes a log-Poisson process, as predicted by RD. Specifically, quakes are shown to be statistically independent and Poisson distributed with an average depending on the logarithm of time. Finally, we discuss the nature of the dynamical barriers surmounted by quakes and link RD to the phenomenology of aging hard sphere colloids.

Automated detection of many-particle solvation states for accurate characterizations of diffusion kinetics

Discrete-space kinetic models, i.e., Markov state models, have emerged as powerful tools for reducing the complexity of trajectories generated from molecular dynamics simulations. These models require configuration-space representations that accurately characterize the relevant dynamics. Well-established, low-dimensional order parameters for constructing this representation have led to widespread application of Markov state models to study conformational dynamics in biomolecular systems. On the contrary, applications to characterize single-molecule diffusion processes have been scarce and typically employ system-specific, higher-dimensional order parameters to characterize the local solvation state of the molecule. In this work, we propose an automated method for generating a coarse configuration-space representation, using generic features of the solvation structure-the coordination numbers about each particle. To overcome the inherent noisy behavior of these low-dimensional observables, we treat the features as indicators of an underlying, latent Markov process. The resulting hidden Markov models filter the trajectories of each feature into the most likely latent solvation state at each time step. The filtered trajectories are then used to construct a configuration-space discretization, which accurately describes the diffusion kinetics. The method is validated on a standard model for glassy liquids, where particle jumps between local cages determine the diffusion properties of the system. Not only do the resulting models provide quantitatively accurate characterizations of the diffusion constant, but they also reveal a mechanistic description of diffusive jumps, quantifying the heterogeneity of local diffusion.

Comparison of single particle dynamics at the center and on the surface of equilibrium glassy films

Glasses prepared by vapor depositing molecules onto a properly prepared substrate can have enhanced kinetic stability when compared with glasses prepared by cooling from the liquid state. The enhanced stability is due to the high mobility of particles at the surface, which allows them to find lower energy configurations than for liquid cooled glasses. Here we use molecular dynamics simulations to examine the temperature dependence of the single particle dynamics in the bulk of the film and at the surface of the film. First, we examine the temperature dependence of the self-intermediate scattering functions for particles in the bulk and at the surface. We then examine the temperature dependence of the probability of the logarithm of single particle displacements for bulk and surface particles. Both bulk and surface particle displacements indicate populations of slow and fast particles, i.e., heterogeneous dynamics. We find that the temperature dependence of the surface dynamics mirrors the bulk despite being several orders of magnitude faster.

Single molecule translocation in smectics illustrates the challenge for time-mapping in simulations on multiple scales

Understanding the connections between the characteristic dynamical time scales associated with a coarse-grained (CG) and a detailed representation is central to the applicability of the coarse-graining methods to understand molecular processes. The process of coarse graining leads to an accelerated dynamics, owing to the smoothening of the underlying free-energy landscapes. Often a single time-mapping factor is used to relate the time scales associated with the two representations. We critically examine this idea using a model system ideally suited for this purpose. Single molecular transport properties are studied via molecular dynamics simulations of the CG and atomistic representations of a liquid crystalline, azobenzene containing mesogen, simulated in the smectic and the isotropic phases. The out-of-plane dynamics in the smectic phase occurs via molecular hops from one smectic layer to the next. Hopping can occur via two mechanisms, with and without significant reorientation. The out-of-plane transport can be understood as a superposition of two (one associated with each mode of transport) independent continuous time random walks for which a single time-mapping factor would be rather inadequate. A comparison of the free-energy surfaces, relevant to the out-of-plane transport, qualitatively supports the above observations. Thus, this work underlines the need for building CG models that exhibit both structural and dynamical consistency to the underlying atomistic model.

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

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

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.

Elastically driven intermittent microscopic dynamics in soft solids

Soft solids with tunable mechanical response are at the core of new material technologies, but a crucial limit for applications is their progressive aging over time, which dramatically affects their functionalities. The generally accepted paradigm is that such aging is gradual and its origin is in slower than exponential microscopic dynamics, akin to the ones in supercooled liquids or glasses. Nevertheless, time- and space-resolved measurements have provided contrasting evidence: dynamics faster than exponential, intermittency and abrupt structural changes. Here we use 3D computer simulations of a microscopic model to reveal that the timescales governing stress relaxation, respectively, through thermal fluctuations and elastic recovery are key for the aging dynamics. When thermal fluctuations are too weak, stress heterogeneities frozen-in upon solidification can still partially relax through elastically driven fluctuations. Such fluctuations are intermittent, because of strong correlations that persist over the timescale of experiments or simulations, leading to faster than exponential dynamics.

Cage Size and Jump Precursors in Glass-Forming Liquids: Experiment and Simulations

Glassy dynamics is intermittent, as particles suddenly jump out of the cage formed by their neighbors, and heterogeneous, as these jumps are not uniformly distributed across the system. Relating these features of the dynamics to the diverse local environments explored by the particles is essential to rationalize the relaxation process. Here we investigate this issue characterizing the local environment of a particle with the amplitude of its short time vibrational motion, as determined by segmenting in cages and jumps the particle trajectories. Both simulations of supercooled liquids and experiments on colloidal suspensions show that particles in large cages are likely to jump after a small time-lag, and that, on average, the cage enlarges shortly before the particle jumps. At large time-lags, the cage has essentially a constant size, which is smaller for longer-lasting cages. Finally, we clarify how this coupling between cage size and duration controls the average behavior and opens the way to a better understanding of the relaxation process in glass-forming liquids.

Differential Variance Analysis: a direct method to quantify and visualize dynamic heterogeneities

Many amorphous materials show spatially heterogenous dynamics, as different regions of the same system relax at different rates. Such a signature, known as Dynamic Heterogeneity, has been crucial to understand the nature of the jamming transition in simple model systems and is currently considered very promising to characterize more complex fluids of industrial and biological relevance. Unfortunately, measurements of dynamic heterogeneities typically require sophisticated experimental set-ups and are performed by few specialized groups. It is now possible to quantitatively characterize the relaxation process and the emergence of dynamic heterogeneities using a straightforward method, here validated on video microscopy data of hard-sphere colloidal glasses. We call this method Differential Variance Analysis (DVA), since it focuses on the variance of the differential frames, obtained subtracting images at different time-lags. Moreover, direct visualization of dynamic heterogeneities naturally appears in the differential frames, when the time-lag is set to the one corresponding to the maximum dynamic susceptibility. This approach opens the way to effectively characterize and tailor a wide variety of soft materials, from complex formulated products to biological tissues.

V-T theory for the self-intermediate scattering function in a monatomic liquid

In V-T theory the atomic motion is harmonic vibrations in a liquid-specific potential energy valley, plus transits, which move the system rapidly among the multitude of such valleys. In its first application to the self intermediate scattering function (SISF), V-T theory produced an accurate account of molecular dynamics (MD) data at all wave numbers q and time t. Recently, analysis of the mean square displacement (MSD) resolved a crossover behavior that was not observed in the SISF study. Our purpose here is to apply the more accurate MSD calibration to the SISF, and assess the results. We derive and discuss the theoretical equations for vibrational and transit contributions to the SISF. The time evolution is divided into three successive intervals: the vibrational interval when the vibrational contribution alone accurately accounts for the MD data; the crossover when the vibrational contribution saturates and the transit contribution becomes resolved; and the diffusive interval when the transit contribution alone accurately accounts for the MD data. The resulting theoretical error is extremely small at all q and t. V-T theory is compared to mode-coupling theories for the MSD and SISF, and to recent developments in Brownian motion experiments and theory.