Dynamical heterogeneity in a highly supercooled liquid: consistent calculations of correlation length, intensity, and lifetime

Structural Origin of Dynamic Heterogeneity in Supercooled Liquids

As a liquid is supercooled toward the glass transition point, its dynamics slow significantly, provided that crystallization is avoided. With increased supercooling, the particle dynamics become more spatially heterogeneous, a phenomenon known as dynamic heterogeneity. Since its discovery, this characteristic of metastable supercooled liquids has garnered considerable attention in glass science. However, the precise physical origins of dynamic heterogeneity remain elusive and widely debated. In this perspective, we examine the relationship between dynamic heterogeneity and structural order, based on numerical simulations of fragile liquids with isotropic potentials and strong liquids with directional interactions. We demonstrate that angular ordering, arising from many-body steric interactions, plays a crucial role in the slow dynamics and dynamic cooperativity of fragile liquids. Additionally, we explore how the growth of static order correlates with slower dynamics. In fragile liquids exhibiting super-Arrhenius behavior, the spatial extent of regions with high angular order grows upon cooling, and the sequential propagation of particle rearrangements within these ordered regions increases the activation energy for particle motion. In contrast, strong liquids with spatially constrained local ordering display a distinct "two-state" dynamic characteristic, marked by a transition between two Arrhenius-type behaviors. We argue that dynamic heterogeneity, irrespective of a liquid's fragility, arises from underlying structural order, with its spatial extent determined by static ordering. This perspective aims to deepen our understanding of the interplay between structural and dynamic properties in metastable supercooled liquids.

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.

Single-particle dynamics near the glass transition of a metallic glass

The single-particle dynamics of the glass-forming Cu_{50}Zr_{50} alloy, from the supercooled liquid well above the glass-transition temperature, T_{g} to the glassy state, is studied by using the molecular dynamics simulations. When the liquid is cooled below 1.2T_{g}, the dynamics heterogeneity characterized by the cage-jump motion becomes increasingly pronounced. The analyses based on the continuous time random walk method indicate that the liquid falls out of equilibrium in the present simulation time scale when it is cooled into the regime below 1.02T_{g}. However, we find that the jump length and the jump rate do not display the non-equilibrium behaviors even in the glassy state below T_{g}, which allows us to study the intrinsic dynamic characteristics through T_{g}. The mean waiting time between two successive jumps has a rapid growth following the Vogel-Fulcher-Tammann law as the non-equilibrium regime is approached, in analogy with the temperature behaviors of transport properties for fragile supercooled liquids. In contrast, the jump rate maintains the Arrhenius decay and the jump length has even a weaker temperature dependence when the liquid is cooled into glassy state. We find that a pronounced enhancement of the spatial correlation of jumps occurs accompanied by the glass transition: the string-like cooperative jumps dominate the fast motion instead of the uncorrelated and individual jumps. Our work offers an insight into the equilibrium effect of the single-particle dynamics in glass transition.

Analysis of structural correlations in a model binary 3D liquid through the eigenvalues and eigenvectors of the atomic stress tensors

It is possible to associate with every atom or molecule in a liquid its own atomic stress tensor. These atomic stress tensors can be used to describe liquids' structures and to investigate the connection between structural and dynamic properties. In particular, atomic stresses allow to address atomic scale correlations relevant to the Green-Kubo expression for viscosity. Previously correlations between the atomic stresses of different atoms were studied using the Cartesian representation of the stress tensors or the representation based on spherical harmonics. In this paper we address structural correlations in a 3D model binary liquid using the eigenvalues and eigenvectors of the atomic stress tensors. This approach allows to interpret correlations relevant to the Green-Kubo expression for viscosity in a simple geometric way. On decrease of temperature the changes in the relevant stress correlation function between different atoms are significantly more pronounced than the changes in the pair density function. We demonstrate that this behaviour originates from the orientational correlations between the eigenvectors of the atomic stress tensors. We also found correlations between the eigenvalues of the same atomic stress tensor. For the studied system, with purely repulsive interactions between the particles, the eigenvalues of every atomic stress tensor are positive and they can be ordered: λ1 ≥ λ2 ≥ λ3 ≥ 0. We found that, for the particles of a given type, the probability distributions of the ratios (λ2/λ1) and (λ3/λ2) are essentially identical to each other in the liquids state. We also found that λ2 tends to be equal to the geometric average of λ1 and λ3. In our view, correlations between the eigenvalues may represent "the Poisson ratio effect" at the atomic scale.

Fluctuation analysis of time-averaged mean-square displacement for the Langevin equation with time-dependent and fluctuating diffusivity

The mean-square displacement (MSD) is widely utilized to study the dynamical properties of stochastic processes. The time-averaged MSD (TAMSD) provides some information on the dynamics which cannot be extracted from the ensemble-averaged MSD. In particular, the relative standard deviation (RSD) of the TAMSD can be utilized to study the long-time relaxation behavior. In this work, we consider a class of Langevin equations which are multiplicatively coupled to time-dependent and fluctuating diffusivities. Various interesting dynamics models such as entangled polymers and supercooled liquids can be interpreted as the Langevin equations with time-dependent and fluctuating diffusivities. We derive a general formula for the RSD of the TAMSD for the Langevin equation with the time-dependent and fluctuating diffusivity. We show that the RSD can be expressed in terms of the correlation function of the diffusivity. The RSD exhibits the crossover at the long time region. The crossover time is related to a weighted average relaxation time for the diffusivity. Thus the crossover time gives some information on the relaxation time of fluctuating diffusivity which cannot be extracted from the ensemble-averaged MSD. We discuss the universality and possible applications of the formula via some simple examples.

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.

Dynamic length scales in glass-forming liquids: an inhomogeneous molecular dynamics simulation approach

In this work, we numerically investigate a new method for the characterization of growing length scales associated with spatially heterogeneous dynamics of glass-forming liquids. This approach, motivated by the formulation of the inhomogeneous mode-coupling theory (IMCT) [Biroli, G.; et al. Phys. Rev. Lett. 2006 97, 195701], utilizes inhomogeneous molecular dynamics simulations in which the system is perturbed by a spatially modulated external potential. We show that the response of the two-point correlation function to the external field allows one to probe dynamic correlations. We examine the critical properties shown by this function, in particular, the associated dynamic correlation length, that is found to be comparable to the one extracted from standardly employed four-point correlation functions. Our numerical results are in qualitative agreement with IMCT predictions but suggest that one has to take into account fluctuations not included in this mean-field approach to reach quantitative agreement. Advantages of our approach over the more conventional one based on four-point correlation functions are discussed.

Computer simulation of diffusion in silica liquid under temperature and pressure

We have studied the diffusion mechanism in silica liquid following a new approach where the diffusion rate is estimated via the rate of SiO(x) → SiO(x±1) and the mean square displacement of Si particles per SiO(x) → SiO(x±1). Molecular dynamics simulation has been conducted for a model consisting of 1998 particles over a wide range of temperatures (3000-4500 K) and pressure (from 0 to 25.75 GPa). Our results show that the rate of SiO(x) → SiO(x±1) increases either with increasing the temperature or pressure. Further, we find that SiO(x) → SiO(x±1) is heterogeneously distributed through the network structure of the liquid. In particular, it is concentrated on a small section of Si particles in a low-temperature regime and at ambient pressure. The spatial localisation of SiO(x) → SiO(x±1) originates from the fact that the stable unit in low- and high-pressure regime is SiO4 and SiO6, respectively. The major change in the diffusion mechanism under pressure or temperature concerns the change in the distribution of SiO(x) → SiO(x±1) through the network structure. It is finally shown that the spatial localisation of SiO(x) → SiO(x±1) is responsible for the dynamics heterogeneity and the diffusion anomaly for silica liquid. This finding supports the concept that as the temperature approaches the glass transition point, SiO(x) → SiO(x±1) spatially localises such that the diffusivity drops and the dynamics are anomalously slow.

Correlation effect for dynamics in silica liquid

We study numerically the diffusion mechanism in silica liquid via molecular dynamics simulation. For this purpose we examine the evolution of structural units SiO(x) (x=4-6) for different times and at temperatures from 3000 to 4500 K. Simulation shows that the diffusivity of the silicon particle is performed through the transition SiO(x)→SiO(x±1), i.e., the bond-breaking and bond-reformation events. As a SiO(x)→SiO(x±1) transition occurs, one oxygen particle moves out of or into the coordination shell, leading to a collective movement of Si particles. Other types of transitions, for instance, SiO(4)→SiO(6) or SiO(6)→SiO(4), are negligible. We establish an expression for the diffusion coefficient that shows that the diffusivity is not proportional to the rate of SiO(x)→SiO(x±1) because it is strongly localized in the network structure. A high degree of localization of SiO(x)→SiO(x±1) leads to a heterogeneous dynamics. We find that the dynamics slowdown is determined by two terms: The first one concerns the change in the statistic property related to the fraction of non-four-coordinated units (SiO(3), SiO(5), SiO(6), and SiO(7)) and the second term concerns the correlation effect. Furthermore, we show that the correlation coefficient depends on both the fraction of the back-forth SiO(x)→SiO(x±1) transition and the degree of localization of SiO(x)→SiO(x±1). Our finding qualitatively supports the ideal that anomalously slow dynamics near the glass-transition point is caused by a strong localization of SiO(x)→SiO(x±1).

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

We investigate the dynamic heterogeneities of glassy particle systems in the theoretical schemes of bond breakage and four-point correlation functions. In the bond-breakage scheme, we introduce the structure factor S(b)(q,t) and the susceptibility χ(b)(t) to detect the spatial correlations of configuration changes. Here χ(b)(t) attains a maximum at t=t(b)(max) as a function of time t, where the fraction of the particles with broken bonds φ(b)(t) is about 1/2. In the four-point scheme, treating the structure factor S(4)(q,t) and the susceptibility χ(4)(t), we detect superpositions of the heterogeneity of bond breakage and that of thermal low-frequency vibration modes. While the former grows slowly, the latter emerges quickly to exhibit complex space-time behavior. In two dimensions, the vibration modes extending over the system yield significant contributions to the four-point correlations, which depend on the system size logarithmically. A maximum of χ(4)(t) is attained at t=t(4)(max), where these two contributions become of the same order. As a result, t(4)(max) is considerably shorter than t(b)(max).