Dynamics Of Supercooled Liquids Research Articles

Emergent facilitation and glassy dynamics in supercooled liquids

In supercooled liquids, dynamical facilitation refers to a phenomenon where microscopic motion begets further motion nearby, resulting in spatially heterogeneous dynamics. This is central to the glassy relaxation dynamics of such liquids, which show super-Arrhenius growth of relaxation timescales with decreasing temperature. Despite the importance of dynamical facilitation, there is no theoretical understanding of how facilitation emerges and impacts relaxation dynamics. Here, we present a theory that explains the microscopic origins of dynamical facilitation. We show that dynamics proceeds by localized bond-exchange events, also known as excitations, resulting in the accumulation of elastic stresses with which new excitations can interact. At low temperatures, these elastic interactions dominate and facilitate the creation of new excitations near prior excitations. Using the theory of linear elasticity and Markov processes, we simulate a model, which reproduces multiple aspects of glassy dynamics observed in experiments and molecular simulations, including the stretched exponential decay of relaxation functions, the super-Arrhenius behavior of relaxation timescales as well as their two-dimensional finite-size effects. The model also predicts the subdiffusive behavior of the mean squared displacement (MSD) on short, intermediate timescales. Furthermore, we derive the phonon contributions to diffusion and relaxation, which when combined with the excitation contributions produce the two-step relaxation processes, and the ballistic-subdiffusive-diffusive crossover MSD behaviors commonly found in supercooled liquids.

Local vs. cooperative: Unraveling glass transition mechanisms with SEER

Which phenomenon slows down the dynamics in supercooled liquids and turns them into glasses is a long-standing question of condensed matter. Most popular theories posit that as the temperature decreases, many events must occur in a coordinated fashion on a growing length scale for relaxation to occur. Instead, other approaches consider that local barriers associated with the elementary rearrangement of a few particles or "excitations" govern the dynamics. To resolve this conundrum, our central result is to introduce an algorithm, Systematic Excitation ExtRaction, which can systematically extract hundreds of excitations and their energy from any given configuration. We also provide a measurement of the activation energy, characterizing the liquid dynamics, based on fast quenching and reheating. We use these two methods in a popular liquid model of polydisperse particles. Such polydisperse models are known to capture the hallmarks of the glass transition and can be equilibrated efficiently up to millisecond time scales. The analysis reveals that cooperative effects do not control the fragility of such liquids: the change of energy of local barriers determines the change of activation energy. More generally, these methods can now be used to measure the degree of cooperativity of any liquid model.

Atomic-level mechanisms of short-circuit diffusion in materials

Abstract This paper reviews the recent progress in understanding the atomic mechanisms of short-circuit diffusion along materials interfaces, such as grain and interphase boundaries, as well as lattice and interfacial dislocations/disconnections. Recent atomistic computer simulations have shown that short-circuit diffusion is dominated by collective atomic rearrangements in the form of strings and rings of mobile atoms. The process is dynamically heterogeneous in space and time and has many features in common with atomic dynamics in supercooled glass-forming liquids. We discuss examples of grain boundary, interphase boundary, and dislocation diffusion in metals and alloys, including the solute effect on the diffusion rates and mechanisms. Interphase boundaries are exemplified by Al–Si interfaces with diverse orientation relationships and atomic structures. The hierarchy of short-circuit diffusion paths in materials is reviewed by comparing the rates of grain boundary, interphase boundary, and dislocation diffusion. Future directions in the field of short-circuit diffusion in defect core regions are discussed.

Observation of an isothermal glass transition in metallic glasses.

Glass transition, commonly manifested upon cooling a liquid, is continuous and cooling rate dependent. For decades, the thermodynamic basis in liquid-glass transition has been at the center of debate. Here, long-time isothermal annealing was conducted via molecular dynamics simulations for metallic glasses to explore the connection of physical aging in supercooled liquid and glassy states. An anomalous two-step aging is observed in various metallic glasses, exhibiting features of supercooled liquid dynamics in the first step and glassy dynamics in the second step, respectively. Furthermore, the transition potential energy is independent of initial states, proving that it is intrinsic for a metallic glass at a given temperature. We propose that the observed dynamic transition from supercooled liquid dynamics to glassy dynamics could be glass transition manifested isothermally. On this basis, glass transition is no longer cooling rate dependent, but is shown as a clear phase boundary in the temperature-energy phase diagram. Hence, a modified out-of-equilibrium phase diagram is proposed, providing new insights into the nature of glass transition.

How Probes Affect the Local Dynamics in Supercooled Liquid Revealed by Single Molecule Fluorescence Microscopy

AbstractDue to its high sensitivity and spatial resolution, single molecule fluorescence microscopy is widely used to investigate the dynamic heterogeneity of supercooled liquids around the glass transition temperature, powerful to access microscopic information rather than the merely ensemble averaging. Plenty of works are done, provided a deeper insight into the dynamic heterogeneity of supercooled liquids. However, controversies exist when different kinds of probes are adopted for investigation, which report the local dynamics with discrepancies as a result of their diverse photophysical properties. In this review, a guidance is provided for the optical setup in wide‐field configuration and sample preparation, together with the recent progress on the photophysical properties of probes influencing the reported local dynamics, such as the probe photostability, the intermolecular interaction between probe and host matrix, covalent bonding of probe to host molecule as well as the probe size.

Predicting self-diffusion coefficients in semi-crystalline and amorphous solid dispersions using free volume theory

The self-diffusion coefficient of active ingredients (AI) in polymeric solid dispersions is one of the essential parameters for the rational formulation design in life sciences. Measuring this parameter for products in their application temperature range can, however, be difficult to realise and time-consuming (due to the slow kinetics of diffusion). The aim of this study is to present a simple and time-saving platform for predicting the AI self-diffusivity in amorphous and semi-crystalline polymers on the basis of a modified version of Vrentas’ and Duda’s free volume theory (FVT) [A. Mansuri, M. Völkel, T. Feuerbach, J. Winck, A.W.P. Vermeer, W. Hoheisel, M. Thommes, Modified free volume theory for self-diffusion of small molecules in amorphous polymers, Macromolecules. (2023)]. The predictive model discussed in this work requires pure-component properties as its input and covers the approximate temperature range of T < 1.2 Tg, the whole compositional range of the binary mixtures (as long as a molecular mixture is present), and the whole crystallinity range of the polymer. In this context, the self-diffusion coefficients of the AIs imidacloprid, indomethacin, and deltamethrin were predicted in polyvinylpyrrolidone, polyvinylpyrrolidone/vinyl acetate, polystyrene, polyethylene, and polypropylene. The results highlight the profound importance of the kinetic fragility of the solid dispersion on the molecular migration; a property which in some cases might entail higher self-diffusion coefficients despite an increase in the molecular weight of the polymer. We interpret this observation within the context of the theory of heterogeneous dynamics in glass-formers [M.D. Ediger, Spatially heterogeneous dynamics in supercooled liquids, Annu. Rev. Phys. Chem. 51 (2000) 99–128] by attributing it to the stronger presence of “fluid-like” mobile regions in fragile polymers offering facilitated routes for the AI diffusion within the dispersion. The modified FVT further allows for identifying the influence of some structural and thermophysical material properties on the translational mobility of AIs in binary dispersions with polymers. In addition, estimates of self-diffusivity in semi-crystalline polymers are provided by further accounting for the tortuosity of the diffusion paths and the chain immobilisation at the interface of the amorphous and crystalline phases.

From Single-Particle to Collective Dynamics in Supercooled Liquids.

It has been recognized recently that the considerable difference between photon correlation (PCS) and dielectric (BDS) susceptibility spectra arises from their respective association with single-particle and collective dynamics. This work presents a model that captures the narrower width and shifted peak position of collective dynamics (BDS), given the single-particle susceptibility derived from PCS studies. Only one adjustable parameter is required to connect the spectra of collective and single-particle dynamics. This constant accounts for cross-correlations between molecular angular velocities and the ratio of the first- and second-rank single-particle relaxation times. The model is tested for three supercooled liquids, glycerol, propylene glycol, and tributyl phosphate, and is shown to provide a good account of the difference between BDS and PCS spectra. Because PCS spectra appear to be rather universal across a range of supercooled liquids, this model provides a first step toward rationalizing the more material-specific dielectric loss profiles.

Dynamics of supercooled liquids from static averaged quantities using machine learning

We introduce a machine-learning approach to predict the complex non-Markovian dynamics of supercooled liquids from static averaged quantities. Compared to techniques based on particle propensity, our method is built upon a theoretical framework that uses as input and output system-averaged quantities, thus being easier to apply in an experimental context where particle resolved information is not available. In this work, we train a deep neural network to predict the self intermediate scattering function of binary mixtures using their static structure factor as input. While its performance is excellent for the temperature range of the training data, the model also retains some transferability in making decent predictions at temperatures lower than the ones it was trained for, or when we use it for similar systems. We also develop an evolutionary strategy that is able to construct a realistic memory function underlying the observed non-Markovian dynamics. This method lets us conclude that the memory function of supercooled liquids can be effectively parameterized as the sum of two stretched exponentials, which physically corresponds to two dominant relaxation modes.

Inherent-state melting and the onset of glassy dynamics in two-dimensional supercooled liquids

Below the onset temperature To, the equilibrium relaxation time of most glass-forming liquids exhibits glassy dynamics characterized by a super-Arrhenius temperature dependence. In this supercooled regime, the relaxation dynamics also proceeds through localized elastic excitations corresponding to hopping events between inherent states, i.e., potential-energy-minimizing configurations of the liquid. Despite its importance in distinguishing the supercooled regime from the high-temperature regime, the microscopic origin of To is not yet known. Here, we construct a theory for the onset temperature in two dimensions and find that an inherent-state melting transition, described by the binding-unbinding transition of dipolar elastic excitations, delineates the supercooled regime from the high-temperature regime. The corresponding melting transition temperature is in good agreement with the onset temperature found in various two-dimensional (2D) atomistic models of glass formers and an experimental binary colloidal system confined to a water-air interface. Additionally, we find the predictions for the renormalized elastic moduli to agree with the experimentally observed values for the latter 2D colloidal system. We further discuss the predictions of our theory on the displacement and density correlations at supercooled conditions, which are consistent with observations of the Mermin-Wagner fluctuations in experiments and molecular simulations.

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.

Long-Wavelength Fluctuations in Quasi-2D Supercooled Liquids.

We use molecular simulation to characterize the dynamics of supercooled liquids confined in quasi-2D slit geometries. Similar to bulk supercooled liquids, the confined systems exhibit subdiffusive dynamics on intermediate time scales arising from particle localization inside their neighbor cages, followed by an eventual crossover to diffusive behavior as cage rearrangement occurs. The quasi-2D confined liquids also exhibit signatures of long-wavelength fluctuations (LWFs) in the lateral directions parallel to the confining walls, reminiscent of the collective displacements observed in 2D but not 3D systems. The magnitude of the LWFs increases with the lateral dimensions of systems with the same particle volume fraction and confinement length scale, consistent with the logarithmic scaling predicted for 2D Mermin-Wagner fluctuations. The amplitude of the fluctuations is a nonmonotonic function of the confinement length scale because of a competition between caging and strengthening LWFs upon approaching the 2D limit. Our findings suggest that LWFs may play an important role in understanding the behavior of confined supercooled liquids due to their prevalence over a surprisingly broad range of particle densities and confinement length scales.

Understanding dynamics in coarse-grained models. II. Coarse-grained diffusion modeled using hard sphere theory.

The first paper of this series [J. Chem. Phys. 158, 034103 (2023)] demonstrated that excess entropy scaling holds for both fine-grained and corresponding coarse-grained (CG) systems. Despite its universality, a more exact determination of the scaling relationship was not possible due to the semi-empirical nature. In this second paper, an analytical excess entropy scaling relation is derived for bottom-up CG systems. At the single-site CG resolution, effective hard sphere systems are constructed that yield near-identical dynamical properties as the target CG systems by taking advantage of how hard sphere dynamics and excess entropy can be analytically expressed in terms of the liquid packing fraction. Inspired by classical equilibrium perturbation theories and recent advances in constructing hard sphere models for predicting activated dynamics of supercooled liquids, we propose a new approach for understanding the diffusion of molecular liquids in the normal regime using hard sphere reference fluids. The proposed "fluctuation matching" is designed to have the same amplitude of long wavelength density fluctuations (dimensionless compressibility) as the CG system. Utilizing the Enskog theory to derive an expression for hard sphere diffusion coefficients, a bridge between the CG dynamics and excess entropy is then established. The CG diffusion coefficient can be roughly estimated using various equations of the state, and an accurate prediction of accelerated CG dynamics at different temperatures is also possible in advance of running any CG simulation. By introducing another layer of coarsening, these findings provide a more rigorous method to assess excess entropy scaling and understand the accelerated CG dynamics of molecular fluids.

Cooperative activated hopping dynamics in binary glass-forming liquids: effects of the size ratio, composition, and interparticle interactions.

Slow dynamics in supercooled and glassy liquids is an important research topic in soft matter physics. Compared to the traditionally focused one-component systems, glassy dynamics in mixture systems adds in a rich set of new complexities, which are fundamentally interesting and also relevant for many technological applications. In this paper, we apply the recently developed self-consistent cooperative hopping theory (SCCHT) to systematically investigate the effects of the size ratio, composition and interparticle interactions on the cooperative activated hopping dynamics of matrix (in larger size) and penetrant (in smaller size) particles in varied binary sphere mixture model systems, with a specific focus on ultrahigh mixture packing fractions that mimic the deeply supercooled glass transition conditions for molecular/polymeric mixture materials. Analysis shows that in these high activation barrier cases, the long-range elastic distortion associated with a matrix particle hopping over its cage confinement always generates an elastic barrier of a nonnegligible magnitude, although the ratio between the elastic barrier and local barrier contribution is sensitively dependent on all three mixture-specific system factors considered in this work. SCCHT predicts two general scenarios of penetrant-matrix cooperative activated hopping dynamics: matrix/penetrant co-hopping (regime 1) or the penetrant mean barrier hopping time shorter than that of the matrix (regime 2). Increasing the penetrant-to-matrix size ratio or the penetrant-matrix cross-attraction strength is found to universally enlarge the composition window of regime 1. Diverse dynamical properties characterising different aspects of the cooperative activated hopping process, including the penetrant and matrix transient localization lengths, penetrant and matrix hopping jump distances, different types of local and elastic activated barriers, and matrix long-time diffusivity, relaxation time and dynamic fragility are quantitatively studied against a wide range of variations over the three system factors. Of particular interest is the universal "anti-plasticization" phenomenon achievable for sufficiently strong cross-attractive interactions. The prospects this work opens for the exploration of a wide variety of polymer-based mixture materials are briefly discussed at the end.

Quantum uncertainty effects in the dynamics of supercooled liquids: A molecular dynamics study.

Dynamics of density fluctuations in quantum supercooled liquids is analyzed using molecular dynamics simulations. In contrast to the classical case, the uncertainty in the particle position (delocalization of quantum particle in space) leads to significant differences in the dynamics of quantum liquids, both in the short- and long-time limits. The effect of uncertainty is found to be significant for length scales smaller than the uncertainty itself, and diminishes as the length scale grows. The dynamic heterogeneity of the system at short times is enhanced due to uncertainty. In the intermediate (β-relaxation) time regime, the heterogeneity tends to get suppressed due to quantum uncertainty. The probability distribution of particle displacements shows highly nonclassical behavior with double-peak structure at short timescales.

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 machine learning techniques for predicting glassy dynamics.

In the quest to understand how structure and dynamics are connected in glasses, a number of machine learning based methods have been developed that predict dynamics in supercooled liquids. These methods include both increasingly complex machine learning techniques and increasingly sophisticated descriptors used to describe the environment around particles. In many cases, both the chosen machine learning technique and choice of structural descriptors are varied simultaneously, making it hard to quantitatively compare the performance of different machine learning approaches. Here, we use three different machine learning algorithms-linear regression, neural networks, and graph neural networks-to predict the dynamic propensity of a glassy binary hard-sphere mixture using as structural input a recursive set of order parameters recently introduced by Boattini et al. [Phys. Rev. Lett. 127, 088007 (2021)]. As we show, when these advanced descriptors are used, all three methods predict the dynamics with nearly equal accuracy. However, the linear regression is orders of magnitude faster to train, making it by far the method of choice.

Tagged particle dynamics in supercooled quantum liquids

We analyze the dynamics of quantum supercooled liquids in terms of tagged particle dynamics. Unlike the classical case, uncertainty in the position of a particle in quantum liquid leads to qualitative changes. We demonstrate these effects in the dynamics of the first two moments of displacements, namely, the mean-squared displacement, 〈Δr^{2}(t)〉, and 〈Δr^{4}(t)〉. Results are presented for a hard-sphere liquid using mode-coupling theory formulation and simulation on a binary Lennard-Jones liquid. As the quantumness (controlled by the de Broglie thermal wavelength) is increased, a nonzero value of the moments at zero time leads to significant deviations from the classical behavior in the initial dynamics. Initial displacement shows ballistic behavior 〈Δr^{2}(t)〉∼t^{2}, but, as a result of large uncertainty in the position, the dynamical effects become weaker with increasing quantumness over this timescale.

Microscopic origin of excess wings in relaxation spectra of supercooled liquids

Glass formation is encountered in diverse materials. Experiments have revealed that dynamic relaxation spectra of supercooled liquids generically become asymmetric near the glass transition temperature, $T_g$, where an extended power law emerges at high frequencies. The microscopic origin of this "wing" remains unknown, and was so far inaccessible to simulations. Here, we develop a novel computational approach and study the equilibrium dynamics of model supercooled liquids near $T_g$. We demonstrate the emergence of a power law wing in numerical spectra, which originates from relaxation at rare, localised regions over broadly-distributed timescales. We rationalise the asymmetric shape of relaxation spectra by constructing an empirical model associating heterogeneous activated dynamics with dynamic facilitation, which are the two minimal physical ingredients revealed by our simulations. Our work offers a glimpse of the molecular motion responsible for glass formation at relevant experimental conditions.

How the cation size impacts on the relaxational and diffusional dynamics of supercooled butylammonium-based ionic liquids: DPEBA–TFSI versus BTMA–TFSI

Abstract Li-bis(trifluoromethylsulfonyl)imide based ionic liquids with either butyl-trimethylammonium or N,N-dimethyl-N-(2-(propionyloxy)-ethyl)butan-1-ammonium as the anion were studied using proton and fluorine relaxometry as well as using field-gradient diffusometry to gain separate access to cation and anion dynamics in these compounds. The transport parameters obtained for these ionic liquids are compared with the estimates based on the conductivity data from literature and from the present work. The impact of cation size on correlation effects, the latter parameterized in terms of various Haven ratios, is mapped out.

Revisiting impulsive stimulated thermal scattering in supercooled liquids: Relaxation of specific heat and thermal expansion.

Impulsive stimulated thermal scattering (ISTS) allows one to access the structural relaxation dynamics in supercooled molecular liquids on a time scale ranging from nanoseconds to milliseconds. Till now, a heuristic semi-empirical model has been commonly adopted to account for the ISTS signals. This model implicitly assumes that the relaxation of specific heat, C, and thermal expansion coefficient, γ, occur on the same time scale and accounts for them via a single stretched exponential. This work proposes two models that assume disentangled relaxations, respectively, based on the Debye and Havriliak-Negami assumptions for the relaxation spectrum and explicitly accounting for the relaxation of C and γ separately in the ISTS response. A theoretical analysis was conducted to test and compare the disentangled relaxation models against the stretched exponential. The former models were applied to rationalize the experimental ISTS signals acquired on supercooled glycerol. This allows us to simultaneously retrieve the frequency-dependent specific heat and thermal expansion up to the sub-100MHz frequency range and further to compare the fragility and time scale probed by thermal, mechanical, and dielectric susceptibilities.