Stochastic model for the glass transition of simple classical liquids

Microscopic dynamics of enhanced glass-forming ability with minor oxygen addition in bulk metallic glasses

Minor oxygen addition has been proposed as a promising strategy to enhance the performance of metallic glasses, particularly their glass-forming ability. In this work, we investigate the microscopic dynamics of a CuZr glass former with oxygen content up to 2 at. % using molecular dynamics simulations based on specially developed neural network interatomic potentials. Our findings indicate a gradual increase in the glass transition temperature with oxygen addition, with an anomalous peak at 0.4 at. % O. We reveal an anti-correlation of kinetic fragility and dynamic heterogeneity behind this unusual rise, where the system exhibits reduced kinetic fragility alongside more significant dynamic heterogeneity. Using the continuous time random walk method, we show that at 0.4 at. % O, a highly mobile Cu atomic layer forms around O-Zr clusters, resulting in notable dynamic heterogeneity. This dynamic behavior is closely linked to the bonding pattern within the O-Zr network, particularly favoring the configuration with edge and surface sharing. In addition, such structures contribute to a more compact O-Zr network, leading to lower kinetic fragility. These findings provide detailed insights into the microscopic dynamics behind the effects of minor oxygen additions.

Waiting time dependence of aging

Aging phenomena have been observed in many non-equilibrium systems such as polymers and glasses, where physical properties depend on the waiting time between the starting time of observation and the time when the temperature is changed. The aging is classified into two types on the basis of the waiting time dependence of an instantaneous relaxation time: When the relaxation time is always an increasing function of the waiting time, the aging is called Type I and when it depends on the protocol of the temperature change, the aging is called Type II. Aging of a random walk in three dimensions is investigated when the free energy landscape controlling the jump rate responds to temperature change with a delay. It is shown that the intermediate scattering function of the random walk model exhibits Type II aging. It is also shown that the relaxation time of the free energy landscape can be deduced from the waiting time dependence of the instantaneous relaxation time.

Fickian Non-Gaussian Diffusion in Glass-Forming Liquids

Fickian yet non-Gaussian diffusion (FnGD), a most intriguing open issue in soft matter, is generically associated with some dynamical and/or structural heterogeneity of the environment. Here we investigate the features of FnGD in glass-forming liquids, the epitome of dynamical heterogeneity, drawing on experiments on hard-sphere colloidal suspensions and simulations of a simple model of molecular liquid. We demonstrate that FnGD strengthens on approaching the glass transition, by identifying distinct timescales for Fickianity, τ_{F}, and for restoring of Gaussianity, τ_{G}>τ_{F}, as well as their associated length scales, ξ_{F} and ξ_{G}. We find τ_{G}∝τ_{F}^{γ} with γ≃1.8 for both systems. In the deep FnGD regime, the displacement distributions display exponential tails. We show that, in simulations, the time-dependent decay lengths l(t) at different temperatures all collapse onto a power-law master curve [l(t)/(ξ_{G})]∝(t/τ_{G})^{α}, with α=0.33. A similar collapse, if less sharp, is also found in experiments, seemingly with the same exponent α. We further discuss the connections of the timescales and length scales characterizing FnGD with structural relaxation and dynamic heterogeneity.

Brownian motion with alternately fluctuating diffusivity: Stretched-exponential and power-law relaxation

We investigate Brownian motion with diffusivity alternately fluctuating between fast and slow states. We assume that sojourn-time distributions of these two states are given by exponential or power-law distributions. We develop a theory of alternating renewal processes to study a relaxation function which is expressed with an integral of the diffusivity over time. This relaxation function can be related to a position correlation function if the particle is in a harmonic potential and to the self-intermediate scattering function if the potential force is absent. It is theoretically shown that, at short times, the exponential relaxation or the stretched-exponential relaxation are observed depending on the power-law index of the sojourn-time distributions. In contrast, at long times, a power-law decay with an exponential cutoff is observed. The dependencies on the initial ensembles (i.e., equilibrium or nonequilibrium initial ensembles) are also elucidated. These theoretical results are consistent with numerical simulations.

Transport dynamics of complex fluids

Thermal motion in complex fluids is a complicated stochastic process but ubiquitously exhibits initial ballistic, intermediate subdiffusive, and long-time diffusive motion, unless interrupted. Despite its relevance to numerous dynamical processes of interest in modern science, a unified, quantitative understanding of thermal motion in complex fluids remains a challenging problem. Here, we present a transport equation and its solutions, which yield a unified quantitative explanation of the mean-square displacement (MSD), the non-Gaussian parameter (NGP), and the displacement distribution of complex fluids. In our approach, the environment-coupled diffusion kernel and its time correlation function (TCF) are the essential quantities that determine transport dynamics and characterize mobility fluctuation of complex fluids; their time profiles are directly extractable from a model-free analysis of the MSD and NGP or, with greater computational expense, from the two-point and four-point velocity autocorrelation functions. We construct a general, explicit model of the diffusion kernel, comprising one unbound-mode and multiple bound-mode components, which provides an excellent approximate description of transport dynamics of various complex fluidic systems such as supercooled water, colloidal beads diffusing on lipid tubes, and dense hard disk fluid. We also introduce the concepts of intrinsic disorder and extrinsic disorder that have distinct effects on transport dynamics and different dependencies on temperature and density. This work presents an unexplored direction for quantitative understanding of transport and transport-coupled processes in complex disordered media.

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.

Renewal events in glass-forming liquids

On cooling toward the glass transition temperature, glass-forming liquids display long periods of localized motion interrupted by fast "jumps" in the single-particle trajectories. Several theoretical models based on these single-particle jumps have been proposed, most prominently the continuous-time random walk (CTRW). The central assumption of the CTRW is that jumps are renewal events, i.e. that the internal clock of a particle can be reset upon a jump. In this paper, I present an easy-to-implement method to test whether jumps detected in a supercooled liquid or glass are renewal events or not. The test was applied to molecular dynamics simulations of a short-chain polymer melt, demonstrating that the jumps can in fact be treated as renewal events. The test further revealed that additional relaxation processes are present which are not accounted for in the CTRW picture, highlighting the limitations of this approach. The notion of renewal events in glass-forming systems could be a very important building block for the interpretation of aging and the glass transition. Furthermore, it could have practical implications for the study of non-equilibrium dynamics in glasses as well as mechanical rejuvenation.

Continuous-time random-walk approach to supercooled liquids. II. Mean-square displacements in polymer melts

The continuous-time random walk (CTRW) describes the single-particle dynamics as a series of jumps separated by random waiting times. This description is applied to analyze trajectories from molecular dynamics (MD) simulations of a supercooled polymer melt. Based on the algorithm presented by Helfferich et al. [Phys. Rev. E 89, 042603 (2014)], we detect jump events of the monomers. As a function of temperature and chain length, we examine key distributions of the CTRW: the jump-length distribution (JLD), the waiting-time distribution (WTD), and the persistence-time distribution (PTD), i.e., the distribution of waiting times for the first jump. For the equilibrium (polymer) liquid under consideration, we verify that the PTD is determined by the WTD. For the mean-square displacement (MSD) of a monomer, the results for the CTRW model are compared with the underlying MD data. The MD data exhibit two regimes of subdiffusive behavior, one for the early α process and another at later times due to chain connectivity. By contrast, the analytical solution of the CTRW yields diffusive behavior for the MSD at all times. Empirically, we can account for the effect of chain connectivity in Monte Carlo simulations of the CTRW. The results of these simulations are then in good agreement with the MD data in the connectivity-dominated regime, but not in the early α regime where they systematically underestimate the MSD from the MD.

Continuous-time random-walk approach to supercooled liquids. I. Different definitions of particle jumps and their consequences

Single-particle trajectories in supercooled liquids display long periods of localization interrupted by "fast moves." This observation suggests a modeling by a continuous-time random walk (CTRW). We perform molecular dynamics simulations of equilibrated short-chain polymer melts near the critical temperature of mode-coupling theory Tc and extract "moves" from the monomer trajectories. We show that not all moves comply with the conditions of a CTRW. Strong forward-backward correlations are found in the supercooled state. A refinement procedure is suggested to exclude these moves from the analysis. We discuss the repercussions of the refinement on the jump-length and waiting-time distributions as well as on characteristic time scales, such as the average waiting time ("exchange time") and the average time for the first move ("persistence time"). The refinement modifies the temperature (T) dependence of these time scales. For instance, the average waiting time changes from an Arrhenius-type to a Vogel-Fulcher-type T dependence. We discuss this observation in the context of the bifurcation of the α process and (Johari) β process found in many glass-forming materials to occur near Tc. Our analysis lays the foundation for a study of the jump-length and waiting-time distributions, their temperature and chain-length dependencies, and the modeling of the monomer dynamics by a CTRW approach in the companion paper [J. Helfferich et al., Phys. Rev. E 89, 042604 (2014)].

Minimal model for short-time diffusion in periodic potentials

We investigate the dynamics of a single, overdamped colloidal particle, which is driven by a constant force through a one-dimensional periodic potential. We focus on systems with large barrier heights where the lowest-order cumulants of the density field, that is, average position and the mean-squared displacement, show nontrivial (nondiffusive) short-time behavior characterized by the appearance of plateaus. We demonstrate that this "cage-like" dynamics can be well described by a discretized master equation model involving two states (related to two positions) within each potential valley. Nontrivial predictions of our approach include analytic expressions for the plateau heights and an estimate of the "de-caging time" obtained from the study of deviations from Gaussian behavior. The simplicity of our approach means that it offers a minimal model to describe the short-time behavior of systems with hindered dynamics.

Metastable-state dynamics of a liquid: a free-energy landscape study

Using the time dependence of density fluctuations in a supercooled liquid obtained from the solutions of the equations of nonlinear fluctuating hydrodynamics (NFH), the evolution of the system in the free energy landscape is studied. A crossover from a continuous fluid type dynamics to that of hopping between different free energy minima is observed as the liquid is increasingly supercooled. We demonstrate that our results are also in agreement with equilibrium density functional analysis of the same system. The density field obtained in the numerical solution of the NFH equations are further analyzed to introduce complimentary density of voids in the supercooled liquid state and its static and dynamic correlations are computed. The nature of the relaxation of vacancy correlations are observed to be similar to that of the density fluctuations.

Aging and nonergodicity beyond the Khinchin theorem

The Khinchin theorem provides the condition that a stationary process is ergodic, in terms of the behavior of the corresponding correlation function. Many physical systems are governed by nonstationary processes in which correlation functions exhibit aging. We classify the ergodic behavior of such systems and suggest a possible generalization of Khinchin's theorem. Our work also quantifies deviations from ergodicity in terms of aging correlation functions. Using the framework of the fractional Fokker-Planck equation, we obtain a simple analytical expression for the two-time correlation function of the particle displacement in a general binding potential, revealing universality in the sense that the binding potential only enters into the prefactor through the first two moments of the corresponding Boltzmann distribution. We discuss applications to experimental data from systems exhibiting anomalous dynamics.

Subdiffusive motion in kinetically constrained models

We discuss a kinetically constrained model in which real-valued local densities fluctuate in time, as introduced recently by Bertin, Bouchaud, and Lequeux. We show how the phenomenology of this model can be reproduced by an effective theory of mobility excitations propagating in a disordered environment. Both excitations and probe particles have subdiffusive motion, characterized by different exponents and operating on different time scales. We derive these exponents, showing that they depend continuously on one of the parameters of the model.

Exploring the potential energy landscape of glass-forming systems: from inherent structures via metabasins to macroscopic transport

In this review a systematic analysis of the potential energy landscape (PEL) of glass-forming systems is presented. Starting from the thermodynamics, the route towards the dynamics is elucidated. A key step in this endeavor is the concept of metabasins. The relevant energy scales of the PEL can be characterized. Based on the simulation results for some glass-forming systems one can formulate a relevant model system (ideal Gaussian glass-former) which can be treated analytically. The macroscopic transport can be related to the microscopic hopping processes, using either the strong relation between energy (thermodynamics) and waiting times (dynamics) or, alternatively, the concepts of the continuous-time random walk. The relation to the geometric properties of the PEL is stressed. The emergence of length scales within the PEL approach as well as the nature of finite-size effects is discussed. Furthermore, the PEL view is compared to other approaches describing the glass transition.

Negative differential mobility of weakly driven particles in models of glass formers

We study the response of probe particles to weak constant driving in kinetically constrained models of glassy systems, and show that the probe's response can be nonmonotonic and give rise to negative differential mobility: increasing the applied force can reduce the probe's drift velocity in the force direction. Other significant nonlinear effects are also demonstrated, such as the enhancement with increasing force of the probe's fluctuations away from the average path, a phenomenon known in other contexts as giant diffusivity. We show that these results can be explained analytically by a continuous-time random walk approximation where there is decoupling between persistence and exchange times for local displacements of the probe. This decoupling is due to dynamic heterogeneity in the glassy system, which also leads to bimodal distributions of probe particle displacements. We discuss the relevance of our results to experiments.

Free-energy landscape for a tagged particle in a dense hard-sphere fluid

Exploiting the thermodynamic potential functional provided by density functional theory, we determine analytically the free-energy landscape (FEL) in a hard-sphere fluid. The FEL is represented in the three-dimensional coordinate space of the tagged particle. We also analyze the distribution of the free-energy barrier between adjacent basins and show that the most provable value and the average of the free-energy barrier are increasing functions of the density. Since the size of the cooperatively rearranging region (CRR) is also increased as the density is raised [Yoshidome, Phys. Rev. E 76, 021506 (2007)], the present result is consistent with the Adam-Gibbs theory in which the increase of the activation energy is due to the increase of the size of the CRR.

Large-amplitude jumps and non-Gaussian dynamics in highly concentrated hard sphere fluids

Our microscopic stochastic nonlinear Langevin equation theory of activated dynamics has been employed to study the real-space van Hove function of dense hard sphere fluids and suspensions. At very short times, the van Hove function is a narrow Gaussian. At sufficiently high volume fractions, such that the entropic barrier to relaxation is greater than the thermal energy, its functional form evolves with time to include a rapidly decaying component at small displacements and a long-range exponential tail. The "jump" or decay length scale associated with the tail increases with time (or particle root-mean-square displacement) at fixed volume fraction, and with volume fraction at the mean alpha relaxation time. The jump length at the alpha relaxation time is predicted to be proportional to a measure of the decoupling of self-diffusion and structural relaxation. At long times corresponding to mean displacements of order a particle diameter, the volume fraction dependence of the decay length disappears. A good superposition of the exponential tail feature based on the jump length as a scaling variable is predicted at high volume fractions. Overall, the theoretical results are in good accord with recent simulations and experiments. The basic aspects of the theory are also compared with a classic jump model and a dynamically facilitated continuous time random-walk model. Decoupling of the time scales of different parts of the relaxation process predicted by the theory is qualitatively similar to facilitated dynamics models based on the concept of persistence and exchange times if the elementary event is assumed to be associated with transport on a length scale significantly smaller than the particle size.

Universal nature of particle displacements close to glass and jamming transitions

We examine the structure of the distribution of single particle displacements (van Hove function) in a broad class of materials close to glass and jamming transitions. In a wide time window comprising structural relaxation, van Hove functions reflect the coexistence of slow and fast particles (dynamic heterogeneity). The tails of the distributions exhibit exponential, rather than Gaussian, decay. We argue that this behavior is universal in glassy materials and should be considered the analog, in space, of the stretched exponential decay of time correlation functions. We introduce a dynamical model that describes quantitatively numerical and experimental data in supercooled liquids, colloidal hard spheres, and granular materials. The tails of the distributions directly explain the decoupling between translational diffusion and structural relaxation observed in glassy materials.

Free energy landscape and cooperatively rearranging region in a hard sphere glass

Exploiting the density functional theory, we calculate the free energy landscape (FEL) of the hard sphere glass in three dimensions. From the FEL, we estimate the number of the particles in the cooperatively rearranging region (CRR). We find that the density dependence of the number of the particles in the CRR is expressed as a power law function of the density. Analyzing the relaxation process in the CRR, we also find that the string motion is the elementary process for the structural relaxation, which leads to the natural definition of the simultaneously rearranging region as the particles displaced in the string motion.

On the role of frustration on the glass transition and polyamorphism of mesoscopically heterophase liquids

The model of heterophase fluctuations is developed accounting frustration of the mesoscopic solidlike fluctuons. Within the framework of this model, the glass transition and polyamorphous transformations are considered. It is shown that the frustration increases the temperature range in which the heterophase liquid state exists. the upper and lower boundaries of this temperature range are determined. These boundaries separate different phase states-amorphous solid, heterophase liquid, and fluid phases. Polyamorphous liquid-liquid transitions in the liquid are investigated. Frustration can call forth continuous fluid-solid phase transformation avoiding the first- or second-order phase transition. Conditions under which the first-order phase transition fraction takes place are formulated. Two scenarios of the first-order liquid-liquid polyamorphous transformation are described. As an example the glacial phase formation and the first-order liquid-liquid phase transition in triphenyl phosphate are considered and discussed. Impact of frustration on the liquid crystallization and crystallinity of the glassy state is studied.

Direct observation of spatiotemporal dependence of anomalous diffusion in inhomogeneous fluid by sampling-volume-controlled fluorescence correlation spectroscopy

The direct observation of a spatiotemporal behavior of anomalous diffusion in aqueous polymer [hyaluronan (HA)] solution was achieved by fluorescence correlation spectroscopy (FCS) using a modified instrument, enabling continuous change of the confocal volume of a microscope, namely, sampling-volume-controlled (SVC) FCS (SVC-FCS). Since HA chains form a mesh structure with a pore size of about 10-40 nm, the observed diffusion coefficient (Dobs) is markedly dependent on the diffusion distance (L). By SVC-FCS, the curve of the distance dependence of diffusion coefficient was directly obtained as a continuous profile in L = 245-600 nm showing evidence of anomalous diffusion. On plotting Dobs against either of the sampling time (tauobs) or the diffusion distance (L), Dobs turnover was observed near the anomalous diffusion area. The appearance of this turnover is attributed to the nonuniform mesh structure that can be observed only by a fast observation and that should be dynamically averaged by polymer motions with large tauobs. This behavior is similar to that revealed in glass, colloidal systems, and gel solutions using dynamic light scattering, neutron scattering, and other techniques.

Cooling rate dependence of specific heat in systems out of equilibrium

The anomaly of specific heat in systems out of equilibrium, especially the measurement procedure dependence of specific heat, is investigated by means of free energy landscape. Introducing measurement procedure which is based on experimental method, we propose a calculation method of specific heat in systems out of equilibrium and find an abrupt change in specific heat between annealed and quenched states. For longer observation time the change in specific heat occurs at lower temperature and becomes sharper. For slower cooling of a system the transition temperature becomes lower. This cooling rate dependence of the transition temperature is consistent with experiments and thus the abrupt change in specific heat can be regarded as the glass transition which is thermally identified.

Dynamical fluctuations in ion conducting glasses: slow and fast components in lithium metasilicate

Molecular dynamics simulations of lithium metasilicate (Li2SiO3) glass have been performed. Dynamic heterogeneity of lithium ions has been examined in detail over 4 ns at 700 K. Type A particles show slow dynamics in accordance with a long tail of waiting time distribution of jump motion and localized jumps within neighboring sites (fracton), while type B particles show fast dynamics related to cooperative jumps and a strong forward correlated motion (Lévy flight). Mutual changes of two kinds of dynamics with the relatively long time scale have been observed. The changes cause an extremely large fluctuation of the mean squared displacements as well as the squared displacement of each particle, which depends on the time window of observation. Localized jump motion (fracton) cannot contribute to the long-time-translational diffusion but it can contribute to the rotational diffusion. On the other hand, forward correlated jump motion mainly contributes to long-time-translational diffusion and not to the short-time-rotational diffusion, although this can be a slower part of the rotational diffusion. These results have been compared with those of simple glass-forming liquids exhibiting the dynamic heterogeneity near T(g). The translation-rotation paradox can be explained by the characteristics of slow and fast dynamics.

Dynamic entropy as a measure of caging and persistent particle motion in supercooled liquids

The length-scale dependence of the dynamic entropy is studied in a molecular dynamics simulation of a binary Lennard-Jones liquid above the mode-coupling critical temperature T(c). A number of methods exist for estimating the entropy of dynamical systems, and we utilize an approximation based on calculating the mean first-passage time (MFPT) for particle displacement because of its tractability and its accessibility in real and simulation measurements. The MFPT dynamic entropy S(epsilon) is defined as equal to the inverse of the average first-passage time for a particle to exit a sphere of radius epsilon. This measure of the degree of chaotic motion allows us to identify characteristic time and space scales and to quantify the increasingly correlated particle motion and intermittency occurring in supercooled liquids. In particular, we identify a "cage" size defining the scale at which the particles are transiently localized, and we observe persistent particle motion at intermediate length scales beyond the scale where caging occurs. Furthermore, we find that the dynamic entropy at the scale of one interparticle spacing extrapolates to zero as the mode-coupling temperature T(c) is approached.

Characteristics of slow and fast ion dynamics in a lithium metasilicate glass

Molecular dynamics simulations of lithium metasilicate (Li(2)SiO(3)) glass have been performed. The motion of lithium ions is divided into slow (A) and fast (B) categories in the glassy state. The waiting time distribution of the jump motion of each component shows power law behavior with different exponents. Slow dynamics are caused by localized jump motions and the long waiting time. On the other hand, the fast dynamics of the lithium ions in Li(2)SiO(3) are characterized as Lévy flight caused by cooperative jumps. Short intervals of jump events also occur in the fast dynamics in the short time region. Both the temporal and spatial terms contribute to the dynamics acceleration and the heterogeneity caused by these two kinds of dynamics is illustrated. The slow dynamics characteristics of the "glass transition" and in the "mixed alkali effect" are discussed.