Skip to main content

ORIGINAL RESEARCH article

Front. Chem., 29 September 2020
Sec. Chemical Physics and Physical Chemistry
This article is part of the Research Topic Integrating Timescales from Molecules Up View all 9 articles

Why Is the Range of Timescale So Wide in Glass-Forming Liquid?

  • 1Department of Materials Science and Engineering, Shull-Wollan Center – Joint Institute for Neutron Sciences, University of Tennessee, Knoxville, Knoxville, TN, United States
  • 2Department of Physics and Astronomy, University of Tennessee, Knoxville, Knoxville, TN, United States
  • 3Materials Sciences and Technology Division, Oak Ridge National Laboratory, Oak Ridge, TN, United States

The viscosity and the relaxation time of a glass-forming liquid vary over 15 orders of magnitude before the liquid freezes into a glass. The rate of the change with temperature is characterized by liquid fragility. The mechanism of such a spectacular behavior and the origin of fragility have long been discussed, but it remains unresolved because of the difficulty of carrying out experiments and constructing theories that bridge over a wide timescale from atomic (ps) to bulk (minutes). Through the x-ray diffraction measurement and molecular dynamics simulation for metallic liquids we suggest that large changes in viscosity can be caused by relatively small changes in the structural coherence which characterizes the medium-range order. Here the structural coherence does not imply that of atomic-scale structure, but it relates to the coarse-grained density fluctuations represented by the peaks in the pair-distribution function (PDF) beyond the nearest neighbors. The coherence length is related to fragility and increases with decreasing temperature, and it diverges only at a negative temperature. This analysis is compared with several current theories which predict a phase transition near the glass transition temperature.

Introduction

The viscosity of many liquids, such as water, is of the order of 10−2 poise (= 10−3 Pa.s). Its timescale, defined by the Maxwell relaxation time, τM = η/G, where η is viscosity and G is the high-frequency shear modulus, is of the order of pico-second (ps). Upon cooling liquid viscosity rises rather quickly, if crystallization can be avoided for instance by fast cooling. At low enough temperatures τM becomes so long that the supercooled liquid behaves like a solid. This kinetically frozen liquid is a glass. The transition to the glassy state is defined by the value of η reaching 1013 poise (= 1012 Pa.s), when τM becomes of the order of 103 s. Thus, the timescale of liquid dynamics changes by as much as 15 orders of magnitude over a moderate temperature range. Such a rapid change has direct implications on glass-forming ability and other properties of glass-forming liquids, as well as on applications. The origin of this large change has long been debated without wide agreement (Debenedetti and Stillinger, 2001; March and Tosi, 2002; Dyre, 2006; Lubchenko and Wolynes, 2007; Götze, 2009; Donth, 2010; Berthier and Biroli, 2011; Edigar and Harrowell, 2012; Parisi et al., 2020), and remains one of the glass mysteries.

Our recent research results suggest that the medium-range order (MRO) in liquid plays a crucial role in dynamics of metallic and other liquids (Ryu et al., 2019, 2020; Egami, 2020; Ryu and Egami, 2020). In this article we discuss these results and their wider implications in relation to other theories of liquids. In particular, we point out that our results do not suggest the divergence of viscosity just below the glass transition temperature, Tg, as many other theories do, and provide a resolution to the Kauzmann paradox concerning entropy extrapolating to negative values at low temperatures (Kauzmann, 1948). Our results also challenge the idea that defect-like objects control atomic transport and deformation in liquid and glass and raise questions on some prevailing theories.

Viscosity, Fragility and MRO

The temperature dependence of viscosity can be expressed in terms of the temperature dependent activation energy, Ea(T), as

η(T)=ηexp(Ea(T)kBT).    (1)

Above the crossover temperature, TA, viscosity shows the Arrhenius behavior with a constant value of Ea, and below TA it becomes strongly super-Arrhenius (Angell, 1995; Kivelson et al., 1995), resulting in rapid increase in viscosity culminating to the glass transition. It has been shown by simulations (Iwashita et al., 2013) and by experiments (Iwashita et al., 2017; Shinohara et al., 2019; Ashcraft et al., 2020) that above TA viscosity is determined by a bond cutting dynamics, and τM = τLC, where τLC is the timescale for an atom to lose just one neighbor. Below TA, however, the τMLC ratio increases rapidly with decreasing temperature, as liquid dynamics becomes more cooperative (Bellissard and Egami, 2018). This increase in cooperativity is the cause of the rapid increase in viscosity with decreasing temperature and eventual glass transition. The rate of increase in viscosity just above Tg is characterized by fragility,

m=dlogη(T)d(Tg/T)|Tg.    (2)

A liquid with a large value of m is called fragile, whereas the one with a smaller value of m is called strong (Angell, 1995). The origin of the fragility is still in dispute (Angell, 1995; Novikov and Sokolov, 2004).

The structure of liquid and glass is usually described by the atomic pair-distribution function (PDF), g(r), which describes the distribution of distances between atoms by

g(r)=14πr2Nρ0i,jδ(r-|ri-rj|),    (3)

where rwi is the position of the i-th atom, i = 1, …., N, δ(r) is the δ-function, ρ0 is the atomic number density, and <….> denotes thermal average. It is related to the structure function,

S(Q)=14πQ2Ni,jexp(iQ·[ri-rj]),    (4)

which can be determined by x-ray or neutron diffraction, through the Fourier-transformation,

g(r)=1+12π2ρ0r0[S(Q)-1]sin(Qr)QdQ.    (5)

According to Ornstein and Zernike (1914) the medium-range PDF beyond the first peak decays with r as

G(r)=4πrρ0[g(r)-1]=G0(r)exp(-r/ξs),    (6)

where G0(r) is the G(r) of the ideal glass, and ξs is the structural coherence length which characterizes the MRO. The ideal glass state defined by G0(r) has long-range density correlation without periodicity in the structure (Ryu et al., 2019). Because the medium-range PDF mostly accounts for the first peak of S(Q) (Cargill, 1975; Ryu et al., 2020), the height of the first peak, S(Q1)−1, where Q1 is the position of the first peak, is proportional to ξs (Ryu et al., 2019).

In Ryu et al. (2019) G(r), thus ξs, was measured for Pd42.5Ni7.5Cu30P20 liquid by high-energy x-ray diffraction using electrostatic levitation, from 420 to 1,100 K through the glass transition (573 K). Just above Tg Ea(T) was found to be directly related to ξs by,

Ea(T)=E0(ξs(T)a)3,    (7)

where a is the average neighbor distance (Ryu et al., 2019) and E0 is a scaling parameter. Because

nc(T)=ρ0(ξs(T))3    (8)

is the number of atoms in the coherence volume, (ξs)3, it is indicative of the degree of atomic cooperativity of local dynamics in liquid. In other words, Ea(T) is proportional to the number of atoms involved in the activation process for viscous flow;

Ea(T)=nc(T)EB,      EB=E0ρ0a3.    (9)

EB represents the bond energy per atom, which is of the order of a fraction of eV and is significantly larger than kBTg, whereas nc is relatively small even at Tg, typically below ten. The ratio of EB/kBTg being larger than unity allows small changes in cooperativity nc(T) resulting in large changes in Ea, and the rapid increase in viscosity below TA. Moreover, for various liquids examined by experiments as well as by simulations it was found that nc at Tg is directly linked to fragility by,

nc(Tg)=mm0,    (10)

where m0 = 8.7 overall, 10.7 for metallic liquids, 7.4 for organic liquids, and 7.3 for network liquids (Ryu and Egami, 2020). Thus, fragility is related to the cooperativity of liquid dynamics and also to the “ideality” of the liquid structure. The liquid ideality is defined by the shape of the first peak of S(Q) being close to Lorentzian, as implied by the Ornstein-Zernike form, Equation (6), and by long ξs (Ryu et al., 2020). The Equation (9) appears similar to that of the classical Adam-Gibbs theory (Adam and Gibbs, 1965) in which the critical size of the cooperatively rearranging region, z*, determines viscosity. However, the coherence volume defined here refers to the correlation in bulk liquid in equilibrium, whereas the cooperatively rearranging region is a transient defective object. This point will be discussed below.

Temperature Dependence of MRO and Viscosity

In Ryu et al. (2019) we studied the temperature dependence of the structure for various liquids by experiments and simulations, and showed that the height of the first peak of the structure function, S(Q), and the coherence length follow the Curie-Weiss law,

ξs(T)=CT-TIG,    (11)

where TIG is the ideal glass temperature which is negative. The origin of this behavior was briefly discussed in terms of the atomic-level stresses Egami (2011) in Ryu et al. (2019), and will be explained elsewhere (Egami and Ryu, 2020). Then at Tg,

dd(Tg/T)(ξs(T)ξs(Tg))|Tg=11TIG/Tg=mc,    (12)

We found that,

mc2=mm1.    (13)

with m1 = 613, as shown in Figure 1 for various metallic alloy liquids. Therefore,

TIGTg=1-(m1m)1/2.    (14)

Because

ρ0=fpVa=6fpπa3,    (15)

where Va is the atomic volume and fp is the atomic packing fraction,

aξs(T)=(m)16(m1)1/2(π6m0fp)1/3(TTg-TIGTg),    (16)
CaTg=(m1)1/2(π6m0fp)1/3m1/6.    (17)
FIGURE 1
www.frontiersin.org

Figure 1. The plot of mc2 vs. m for various metallic liquids. The dotted line is for m1 = 613 in Equation (13).

Because the value of fp is similar for all metallic glasses (~0.7), the plots of as(T) against T/Tg should be similar, except for vertical shifts, as shown in Figure 2. With vertical shifts they collapse to a near universal curve up, except for weak dependence above Tg on m (Figure 3). The value of C calculated by Equation (17), Ccalc, is compared to the value of C obtained by fit with Equation (11), Cfit, for various liquids in Figure 4, showing good agreement. This near universality must be the reason for the success of the Kivelson scaling (Kivelson et al., 1995).

FIGURE 2
www.frontiersin.org

Figure 2. The plots of as(T) against T/Tg for various metallic liquids.

FIGURE 3
www.frontiersin.org

Figure 3. The plots of as(T) against T/Tg for various metallic liquids with vertical shifts to form a near universal curve.

FIGURE 4
www.frontiersin.org

Figure 4. The value of C given by Equation (17), Ccalc, plotted against the value of C fit by Equation (11), Cfit, for various liquids. The dotted line is a linear fit, which is virtually identical to the line for Ccalc = Cfit.

From Equations (1, 7, 17) we have

η(T)=ηexp(E0kBT(T1T-TIG)3),    (18)

where

T1=(m1)1/2(π6m0fp)1/3m1/6Tg.    (19)

Thus, the viscosity just above Tg can be described in terms of E0 and TIG. At temperatures above TA the value of Ea becomes constant (= E) even though ξs keeps decreasing. The crossover is a purely dynamic phenomenon (Iwashita et al., 2013), and the MRO is irrelevant to dynamics above TA. With a reasonable crossover, for instance,

Ea(T)=E(ξs(T)b)d(T),    (20)

where d(T) = 3 for T < Tg, d(T)=3(TgT-TgTA)/(1-TgTA) for Tg < T < TA, and d = 0 for T > TA, a realistic temperature dependence of viscosity can be reproduced as shown in Figure 5. Here we compare the experimentally determined viscosity data of PdCuNiP liquid (Kato et al., 2006; Mohr et al., 2019) with those calculated with the Equation (20). The pentagonal symbol denotes the viscosity calculated with the ξs(T) determined from the PDF measured by x-ray diffraction (Ryu et al., 2019), whereas the dashed line is calculated using the Curie-Weiss law, Equation (11). We assumed TA/Tg = 2.0 (Blodgett et al., 2015), used the high-temperature data (Mohr et al., 2019) to determine the values of E (= 0.77 eV) and η (= 1.82 × 10−5 Pa.s), and the low-temperature data (Kato et al., 2006) to determine the value of b (= 3.72 Å).

FIGURE 5
www.frontiersin.org

Figure 5. The temperature dependence of viscosity of PdNiCuP liquid: Experimentally determined viscosity data of PdNiCuP liquid (Kato et al., 2006; Mohr et al., 2019) compared to those calculated. The pentagonal symbol denotes the viscosity calculated with the ξs(T) determined from the PDF measured by x-ray diffraction (Ryu et al., 2019) using Equation (20), whereas the dashed line was calculated using the Curie-Weiss law, Equation (11).

Comparison With Other Theories and Models

Absence of Divergence

The divergence of viscosity was first predicted by the Vogel-Fulcher-Tamman (VFT) model (Vogel, 1921; Fulcher, 1925; Tammann and Hesse, 1926),

η(T)=η0exp(BT-T0).    (21)

Models based upon structural coherence, such as the icosahedral correlation models (Steinhardt et al., 1981; Tomida and Egami, 1995; Tanaka et al., 2010), predict the divergence of structural coherence, thus the divergence of viscosity, below Tg in the vicinity of the Kauzmann temperature, TK (Kauzmann, 1948). The mode-coupling theory (Götze, 2009) predicts the divergence at a temperature, Tc, which is even higher than Tg, and only defect hopping provides mobility below Tc (Biroli et al., 2006). For a long time, it has been difficult to measure the viscosity of simple liquids above Tg because of crystallization. In the absence of serious questioning many of the current theories still assume the divergence of viscosity in the vicinity of TK.

However, a more recent measurement of viscosity using liquid levitator (Blodgett et al., 2015) suggests that the VFT model actually shows poor fit to the data. It is likely that the prediction of the viscosity divergence is based upon poor extrapolation of viscosity to infinity. In fact, many other models do not predict divergence at T > 0 (Cohen and Grest, 1979; Nussinov, 2004; Demetriou et al., 2006; Elmatad et al., 2009; Mauro et al., 2009). According to Equation (7) viscosity diverges when the coherence length ξs(T) diverges. For metallic glasses the value of ξs(Tg)/a ranges from 1 to 2.7, with the average around 1.8. Therefore, the structure is quite far from ideal even at Tg. The temperature at which the ideal state is achieved in extrapolation, TIG, is negative. Thus, viscosity never diverges at T > 0 and entropy does not become negative, resolving the Kauzmann paradox (Kauzmann, 1948).

Nature of Structural Order

Many theories attribute the origin of increased viscosity to development of some structural order which is frustrated and cannot achieve long-range ordering. The most prominent example of such order is the icosahedral order (Sadoc, 1981; Steinhardt et al., 1981). The idea is that because the icosahedral order is incompatible with periodicity it never grows into long-range order (Nelson, 1983; Sethna, 1983). However, such structural orders depend on chemical composition and local chemistry (Gaskell, 1979). Also, this is just a sophisticated version of the nano-crystalline theory which Frank (1952) tried to disprove by suggesting the possible presence of local icosahedral configuration. Note that liquid is stabilized by configurational entropy: The development of local order of a particular atomic configuration will reduce the entropy and destabilize liquid.

On the other hand, we postulate the ideal liquid/glass state by extrapolating the coherence length ξs to infinity (Ryu et al., 2019). This state has very diverse local structures, with widely varying local configurations. For instance, the population of the icosahedral local structure is merely 0.7%. The order parameter, ξs, does not describe the structural order, but the MRO of local density fluctuation. Higher-order peaks of the PDF at large distances include many interatomic distances within the peak. The width of the high-order PDF peaks is about 0.1 nm, and this defines the spatial resolution of the structure in the ideal state. Thus, the MRO describes coarse-grained density fluctuations, and not the atomic-level structural correlations, because the spatial resolution needs to be better at least by an order of magnitude to specify the atomic structure. In our view the local icosahedral ordering which occurs in single element liquid is not indicative of glass formation, but it is likely be that of nano-scale crystallization or quasicrystal formation. A single-element metallic liquid is a very poor glass-former and easily crystallizes. For a single-component liquid the second peak of S(Q), which is more sensitive to crystallinity, diverges at a positive temperature below Tg (Ryu et al., 2019). This suggests that the divergence of the local order just below Tg implies nano-scale crystallization.

Idea of Defects

In crystalline solids atomic transport occurs only through the motion of lattice defects, such as vacancies and interstitial defects. Because the structure of liquid and glass is strongly disordered and appears to be full of defect-like structures, it was only natural to assume that more defective parts of the structure allow easier atomic transport. This led to many ideas of defects in liquid and glass, including free-volume (Cohen and Turnbull, 1959), cooperatively rearranging region (Adam and Gibbs, 1965), shear-transformation-zone (Argon, 1979), and others including ours (Egami et al., 1980). However, the results above suggest that the bulk properties, the MRO, control atomic transport, not those of defects. Would the concept of defect be still relevant in elucidating the atomic transport? Our answer is that the concept of defect defined by the specific static structure is not applicable to liquid and glass. We have to consider the “structure” as a dynamic entity.

It has been recognized for a long time that the definition of defect in amorphous system is arbitrary, in the absence of the reference structure. Various attempts have been made to define the defects, by studying the nature of the static structure, including the approaches using machine-learning (Cubuk et al., 2015; Bapst et al., 2020). However, it became apparent recently that what matters is the dynamics, not the static structure before deformation. In crystalline solids the defect retains its structural identity even after motion, because of the translational symmetry of the host lattice. In other words, defects are topologically protected by the lattice. In liquid and glass, however, the topology of atomic connectivity is open, and defects are not topologically protected. The atomic configuration before the motion of a defect is very different from that after the motion. In the picture of the potential energy landscape (PEL), the system moves from one valley to the other through a saddle-point. It was fond that, at the saddle-point, the potential energy of the system is high enough for the system to melt locally for a very short time (~1 ps) (Ding et al., 2020). Consequently, the system loses the memory of prior thermal history (Fan et al., 2017). The saddle-point is known to be a generator of chaos (Mason and Piiroinen, 2012; Párraga et al., 2018). The simple, usually hand-written, schematic picture of the PEL gives an impression that the pathway from one valley to the next is pre-determined. However, in reality, the kinetic momenta of atoms, which vary rapidly in time, give rise to large uncertainty in the directions toward which the system evolves. The major virtue of the PEL concept is that by removing the kinetic energy the underlying PEL is clearly exposed. However, to describe the dynamics of the system we need to add back the kinetic energy which introduces uncertainty, particularly at the saddle-point.

The local melting at the saddle-point decouples the pathway from a valley of the PEL up to the saddle-point and the pathway down to another valley. Therefore, what determines the nature of the saddle-point, thus the system dynamics, is not the initial state in the prior valley but the nature of the molten state which reflects the bulk property, such as the MRO. The propensity to start the activation process depends on the energy of the initial state, which can be described in terms of the fictive, or effective, temperature in the glassy state (Langer, 2004; Fan et al., 2017). However, once the process of activation over the PEL saddle point starts it does not matter where it started initially. The dynamics of the system at the saddle point is totally controlled by the bulk properties, the coherence volume, to be specific.

Mode-Coupling Theory

The mode-coupling theory (MCT) is one of the most widely used theories of liquid dynamics. It describes the dynamics in terms of continuum hydrodynamic variables, such as density and current auto-correlation functions. It is based on the Boltzmann-type equation of motion initially developed for colloids. In the equation of motion, the dynamics at time t is coupled to the dynamics at a prior time t' through the memory function which represents the frictional force. The input to the theory is the snapshot structure function, S(Q), particularly its first peak. Because the height of the first peak of S(Q) is proportional to ξs (Ryu et al., 2019), the MCT focuses on the MRO, similarly to our approach. The dynamic correlations are determined by the equation of motion, and the feedback through the memory function determines dynamics, leading to the glass transition.

In colloids, particles are in touch with solvent, which is in the hydrodynamic steady state, so that the use of the frictional term in the Boltzmann equation is justified. However, in atomic liquids atoms interact each other directly via the potential force, so the application of the MCT becomes more contorted. Viscosity is given in terms of the stress autocorrelation via the Green-Kubo equation. In the MCT this retention of stress correlation is expressed as the memory function which gives rise to the frictional force. Therefore, the feedback from the memory function can produce a runaway leading to the divergence of correlation time.

In our approach the dynamics is governed by discrete local atomic activation processes. The probability of activation is controlled by the activation energy which is directly related to the structural coherence length ξs. The ξs is an equilibrium property, which depends only on temperature and the elastic constants through the atomic-level stresses, without the feedback loop through the memory function. In colloids local dynamics is closely coupled to local density fluctuations, because density plays the role of temperature in the hard-sphere system. Hard jamming at the critical density leads to divergence of viscosity. In atomic liquids, however, hard jamming never occurs, because atoms are compressible and thermal activation is always possible. Even though the MCT explains the glass transition of colloidal systems, its applicability to atomic liquids should be examined more carefully.

Infinite Dimension Models

The spin-glass theories of Edwards and Anderson (1975) and Sherrington and Kirkpatrick (1975) used the replica method (Aharony, 1975; Emery, 1975) and established the presence of the spin-glass ground state, at least in the infinite dimensions. In many spin-glasses spins interact through the long-range RKKY interaction. The large number of interacting neighboring spins justifies the use of the mean-field approximation. The replica method was applied later to the glass problem (Mézard and Parisi, 2000).

In spin-glasses randomness is quenched, because the spin Hamiltonian does not change with temperature. In contrast in real liquids and glasses the Hamiltonian varies with time and temperature. The number of atoms involved in action, nc, is small. Therefore, a similar mean-field approximation is more difficult to justify, and atomic discreteness becomes central to the dynamics. For instance, at TA, ξs(Tg)/a ≈ 1, so in Equation (7) E0 represents the bond energy and nc, ~ 2. The dynamics is totally local, and the action of cutting a bond determines viscosity and diffusivity. Even at Tg, nc ranges from 2 to 12, whereas the ideal state, where ξs(Tg)/a → ∞, is achieved only at a negative temperature. Thus, the liquid above Tg is very far from the ideal state. The infinite dimension theories may be justified in the ideal state, but they may not be appropriate for the real glass and liquid which are far removed from the ideal state. The glass theories based on exact solutions in infinite dimensions (Parisi et al., 2020) are beautiful, but the success of its application to real liquids and glasses needs to be proven.

Conclusions

The study of the structural medium-range-order (MRO) in metallic liquids, represented by the coherence length, ξs, through diffraction experiment and simulation shows that the MRO is intimately related to local dynamics and viscosity. Namely the activation energy of viscosity is directly related to the number of atoms involved in local atomic rearrangement for structural excitation, nc(T), which is proportional to (ξs)3. The magnitude of nc(T) is relatively small, 2 ~ 12 even at Tg, so that the discrete nature of the atomic structure, represented by the topology of atomic connectivity network, is crucial. Conversely, it means that the unit energy for activation per atomic bond, EB, is relatively large. Therefore, a small increase in nc(T) would result in large increase in the activation energy and viscosity. In our view this must be the reason why the timescale of liquid dynamics changes so rapidly over a moderate temperature range. At the same time, the system is quite far from the point of viscosity divergence which occurs when nc(T) diverges to infinity. Actually, we predict nc(T) to diverge at a negative temperature, by extrapolation with the Curie-Weiss law. Therefore, the Kauzmann catastrophe never occurs at T > 0. Even though the results presented here focus on metallic liquids, the same approach was successful in elucidating the ideality and fragility of network and some organic liquids (Ryu and Egami, 2020; Ryu et al., 2020), suggesting that this approach may be applicable beyond metallic liquids. This view is at odds with some of the theories and ideas. This conflict will be resolved by further theoretical and experimental advances in the future.

Data Availability Statement

All datasets generated for this study are included in the article.

Author Contributions

This work was conceived and authored by TE with assistance by CR. The data in this work were generated and analyzed by CR. Both authors contributed to the article and approved the submitted version.

Funding

This work was supported by the U.S. Department of Energy, Office of Science, Basic Energy Sciences, Materials Sciences and Engineering Division.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Adam, G., and Gibbs, J. H. (1965). On the temperature dependence of cooperative relaxation properties in glass forming liquids. J. Chem. Phys. 43, 139–146. doi: 10.1063/1.1696442

CrossRef Full Text | Google Scholar

Aharony, A. (1975). Tetracritical point in mixed magnetic crystals. Phys. Rev. Lett. 34, 590–593. doi: 10.1103/PhysRevLett.34.590

CrossRef Full Text | Google Scholar

Angell, C. A. (1995). Formation of glasses from liquids and biopolymers. Science 267, 1924–1935. doi: 10.1126/science.267.5206.1924

PubMed Abstract | CrossRef Full Text | Google Scholar

Argon, A. S. (1979). Plastic deformation in metallic glasses. Acta Metall. 27, 47–58. doi: 10.1016/0001-6160(79)90055-5

CrossRef Full Text | Google Scholar

Ashcraft, R., Wang, Z., Abernathy, D. L., Egami, T., and Kelton, K. F. (2020). Experimental determination of the temperature-dependent Van Hove function in a Zr80Pt20 liquid-a structural origin for the dynamics at high temperature. J. Chem. Phys. 152:074506. doi: 10.1063/1.5144256

PubMed Abstract | CrossRef Full Text | Google Scholar

Bapst, V., Keck, T., Grabska-Barwińska, A., Donner, C., Cubuk, E. D., Schoenholz, S. S., et al. (2020). Unveiling the predictive power of static structure in glassy systems. Nat. Phys. 16, 448–454. doi: 10.1038/s41567-020-0842-8

CrossRef Full Text | Google Scholar

Bellissard, J., and Egami, T. (2018). Simple theory of viscosity in liquids. Phys. Rev. E 98:063005. doi: 10.1103/PhysRevE.98.063005

CrossRef Full Text | Google Scholar

Berthier, L., and Biroli, G. (2011). Theoretical perspective on the glass transition and amorphous materials. Rev. Mod. Phys. 83:587. doi: 10.1103/RevModPhys.83.587

CrossRef Full Text | Google Scholar

Biroli, G., Bouchaud, J.-P., Miyazaki, K., and Reichman, D. (2006). Inhomogeneous mode-coupling theory and growing dynamic length in supercooled liquids. Phys. Rev. Lett. 97:195701. doi: 10.1103/PhysRevLett.97.195701

PubMed Abstract | CrossRef Full Text | Google Scholar

Blodgett, M. E., Egami, T., Nussinov, Z., and Kelton, K. F. (2015). Proposal for universality in the viscosity of metallic liquids. Sci. Rep. 5:13837. doi: 10.1038/srep13837

PubMed Abstract | CrossRef Full Text | Google Scholar

Cargill, G. S. III. (1975). Structure of metallic alloy glasses. Solid State Phys. 30, 227–320. doi: 10.1016/S0081-1947(08)60337-9

CrossRef Full Text | Google Scholar

Cohen, M. H., and Grest, G. (1979). Liquid-glass transition, a free-volume approach. Phys. Rev. B 20, 1077–1098. doi: 10.1103/PhysRevB.20.1077

CrossRef Full Text | Google Scholar

Cohen, M. H., and Turnbull, D. (1959). Molecular transport in liquids and glasses. J. Chem. Phys. 31, 1164–1169.

Cubuk, E. D., Schoenholz, S. S., Rieser, J. M., Malone, B. D., Rottler, J., Durian, D. J., et al. (2015). Identifying structural flow defects in disordered solids using machine-learning methods. Phys. Rev. Lett. 114:108001. doi: 10.1103/PhysRevLett.114.108001

PubMed Abstract | CrossRef Full Text | Google Scholar

Debenedetti, P. G., and Stillinger, F. H. (2001). Supercooled liquids and the glass transition. Nature 410, 259–267. doi: 10.1038/35065704

PubMed Abstract | CrossRef Full Text | Google Scholar

Demetriou, M. D., Harmon, J. S., Tao, M., Duan, G., Samwer, K., and Johnson, W. L. (2006). Cooperative shear model for the rheology of glass-forming metallic liquids. Phys. Rev. Lett. 97:065502. doi: 10.1103/PhysRevLett.97.065502

PubMed Abstract | CrossRef Full Text | Google Scholar

Ding, J., Li, L., Wang, N., Tian, L., Asta, M., Ritchie, R. O., et al. (2020). Nature of the saddle states of structural excitations in Cu64Zr36 metallic glasses.

Donth, E. (2010). The Glass Transition: Relaxation Dynamics of Liquids and Disordered Materials. Berlin: Springer.

Dyre, J. C. (2006). The glass transition and elastic models of glass-forming liquids. Rev. Mod. Phys. 78, 953–972. doi: 10.1103/RevModPhys.78.953

CrossRef Full Text | Google Scholar

Edigar, M. D., and Harrowell, P. (2012). Perspective: supercooled liquids and glasses. J. Chem. Phys. 137:080901. doi: 10.1063/1.4747326

PubMed Abstract | CrossRef Full Text | Google Scholar

Edwards, S. F., and Anderson, P. W. (1975). Theory of spin glasses. J. Phys. F Met. Phys. 5, 965–974. doi: 10.1088/0305-4608/5/5/017

CrossRef Full Text | Google Scholar

Egami, T. (2011). Atomic level stresses. Prog. Mater. Sci. 56, 637–653. doi: 10.1016/j.pmatsci.2011.01.004

CrossRef Full Text | Google Scholar

Egami, T. (2020). Local density correlations in liquids. Front. Phys. 8:50. doi: 10.3389/fphy.2020.00050

CrossRef Full Text | Google Scholar

Egami, T., Maeda, K., and Vitek, V. (1980). Structural defects in amorphous solids: a computer simulation study. Phil. Mag. A 41, 883–901. doi: 10.1080/01418618008243894

CrossRef Full Text | Google Scholar

Egami, T., and Ryu, C. W. (2020). Origin of the medium-range order in liquid.

Elmatad, Y. S., Chandler, D., and Garrahan, J. P. (2009). Corresponding states of structural glass formers. J. Phys. Chem. B 113, 5563–5567. doi: 10.1021/jp810362g

PubMed Abstract | CrossRef Full Text | Google Scholar

Emery, V. J. (1975). Critical properties of many-component systems. Phys. Rev. B 11, 239–247. doi: 10.1103/PhysRevB.11.239

CrossRef Full Text | Google Scholar

Fan, Y., Iwashita, T., and Egami, T. (2017). Energy landscape-driven nonequilibrium evolution of inherent structure in disordered material. Nature Commun. 8:15417.

PubMed Abstract | Google Scholar

Frank, F. C. (1952). Supercooling of liquids. Proc. Roy. Soc. Lond. A 215, 43–46. doi: 10.1098/rspa.1952.0194

CrossRef Full Text | Google Scholar

Fulcher, G. S. (1925). Analysis of recent measurements of the viscosity of glasses. J. Am. Ceram. Soc. 8, 339–355. doi: 10.1111/j.1151-2916.1925.tb16731.x

CrossRef Full Text | Google Scholar

Gaskell, P. H. (1979). A new structural model for amorphous transition metal silicides, borides, phosphides, and carbides. J. Non-Cryst. Solids 32, 207–224. doi: 10.1016/0022-3093(79)90073-5

CrossRef Full Text | Google Scholar

Götze, W. (2009). Complex Dynamics of Glass-Forming Liquids. Oxford: Oxford University Press. doi: 10.1093/acprof:oso/9780199235346.001.0001

CrossRef Full Text | Google Scholar

Iwashita, T., Nicholson, D. M., and Egami, T. (2013). Elementary excitations and crossover phenomenon in liquids. Phys. Rev. Lett. 110:205504. doi: 10.1103/PhysRevLett.110.205504

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwashita, T., Wu, B., Chen, W.-R., Tsutsui, S., Baron, A. Q. R., and Egami, T. (2017). Seeing real-space dynamics of liquid water through inelastic x-ray scattering. Sci. Adv. 3:1603079. doi: 10.1126/sciadv.1603079

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, H., Wada, T., Hasegawa, M., Saida, J., Inoue, A., and Chen, H. S. (2006). Fragility and thermal stability of Pt- and Pd- based bulk glass forming liquids and their correlation with deformability. Scripta Mat. 54, 2023–2027. doi: 10.1016/j.scriptamat.2006.03.025

CrossRef Full Text | Google Scholar

Kauzmann, W. (1948). The nature of the glassy state and the behavior of liquids at low temperatures. Chem. Rev. 43, 219–256. doi: 10.1021/cr60135a002

CrossRef Full Text | Google Scholar

Kivelson, D., Kivelson, S. A., Zhao, X., Nussinov, Z., and Tarjus, G. (1995). A thermodynamic theory of supercooled liquids. Physica A 219, 27–38. doi: 10.1016/0378-4371(95)00140-3

CrossRef Full Text | Google Scholar

Langer, J. S. (2004). Dynamics of shear-transformation zones in amorphous plasticity: formulation in terms of an effective disorder temperature. Phys. Rev. E 70:041502.

PubMed Abstract | Google Scholar

Lubchenko, V., and Wolynes, P. G. (2007). Theory of structural glasses and supercooled liquids. Ann. Rev. Phys. Chem. 58, 235–266. doi: 10.1146/annurev.physchem.58.032806.104653

PubMed Abstract | CrossRef Full Text | Google Scholar

March, N. H., and Tosi, M. P. (2002). Introduction to Liquid State Physics. Singapore: World Scientific. doi: 10.1142/4717

CrossRef Full Text | Google Scholar

Mason, J. F., and Piiroinen, P. T. (2012). Saddle-point solutions and grazing bifurcations in an impacting system. Chaos 22:013106. doi: 10.1063/1.3673786

PubMed Abstract | CrossRef Full Text | Google Scholar

Mauro, J. C., Yue, Y., Ellison, A. J., Gupta, P. K., and Allan, D. C. (2009). Viscosity of glass-forming liquids. Proc. Natl. Acad. Sci. U.S.A. 106, 19780–19784. doi: 10.1073/pnas.0911705106

PubMed Abstract | CrossRef Full Text | Google Scholar

Mézard, M., and Parisi, G. (2000). Statistical physics of structural glass. J. Phys. Cond. Mat. 12, 6655–6673. doi: 10.1088/0953-8984/12/29/336

CrossRef Full Text | Google Scholar

Mohr, M., Wunderlich, R. K., Zweiacker, K., Prades-Rödel, S., Sauget, R., Blatter, A., et al. (2019). Surface tension and viscosity of liquid Pd43Cu27Ni10P20 measured in a levitation device under microgravity. npj Microgravity 5:4. doi: 10.1038/s41526-019-0065-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelson, D. R. (1983). Order, frustration, and defects in liquids and glasses. Phys. Rev. B 28, 5515–5535.

Novikov, V. N., and Sokolov, A. P. (2004). Poisson's ratio and the fragility of glass-forming liquids. Nature 431, 961–963. doi: 10.1038/nature02947

PubMed Abstract | CrossRef Full Text | Google Scholar

Nussinov, Z. (2004). Avoided phase transitions and glassy dynamics in geometrically frustrated systems and non-Abelian theories. Phys Rev B 69:014208. doi: 10.1103/PhysRevB.69.014208

CrossRef Full Text | Google Scholar

Ornstein, L. S., and Zernike, F. (1914). Accidental deviations of density and opalescence at the critical point of a single substance. R. Netherlands Acad. Arts Sci. 17, 793–806.

Parisi, G., Urbani, P., and Zamponi, F. (2020). Theory of Simple Glasses. Cambridge: Cambridge University Press. doi: 10.1017/9781108120494

CrossRef Full Text

Párraga, H., Arranz, F. J., Benito, R. M., and Borondo, F. (2018). Above saddle-point regions of order in a sea of chaos in the vibrational dynamics of KCN. J. Phys. Chem. A 122, 3433–3441. doi: 10.1021/acs.jpca.8b00113

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, C. W., Dmowski, W., and Egami, T. (2020). Ideality of liquid structure: a case study for metallic alloy liquids. Phys. Rev. E 101:030601. doi: 10.1103/PhysRevE.101.030601

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, C. W., Dmowski, W., Kelton, K. F., Lee, G. W., Park, E. S., Morris, J. R., et al. (2019). Curie-weiss behavior of liquid structure and ideal glass state. Sci. Rep. 9:18579. doi: 10.1038/s41598-019-54758-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Ryu, C. W., and Egami, T. (2020). Origin of Liquid Fragility.

Sadoc, J. F. (1981). Use of regular polytopes for the mathematical description of the order in amorphous structures. J. Non-Cryst. Solids 44, 1–16.

Sethna, J. P. (1983). Frustration and curvature: glasses and the cholesteric blue phase. Phys. Rev. Lett. 51:2198.

PubMed Abstract | Google Scholar

Sherrington, D., and Kirkpatrick, S. (1975). Solvable model of a spin glass. Phys. Rev. Lett. 35, 1792–1796. doi: 10.1103/PhysRevLett.35.1792

CrossRef Full Text | Google Scholar

Shinohara, Y., Dmowski, W., Iwashita, T., Wu, B., Ishikawa, D., Baron, A. Q. R., et al. (2019). Erratum: Viscosity and real space molecular motion of water: observation with inelastic X-ray scattering [Phys. Rev. E 98, 022604 (2018)]. Phys. Rev. E 100:039904. doi: 10.1103/PhysRevE.98.022604

PubMed Abstract | CrossRef Full Text | Google Scholar

Steinhardt, P. J., Nelson, D. R., and Ronchetti, M. (1981). Icosahedral bond order in supercooled liquids. Phys. Rev. Lett. 47, 1297–1300. doi: 10.1103/PhysRevLett.47.1297

CrossRef Full Text | Google Scholar

Tammann, G., and Hesse, G. (1926). The dependency of viscosity on temperature in hypothermic liquids. Z. Anorg. Allg. Chem. 156, 245–257. doi: 10.1002/zaac.19261560121

CrossRef Full Text | Google Scholar

Tanaka, H., Kawasaki, T., Shintani, H., and Watanabe, K. (2010). Critical-like behavior of glass-forming liquids. Nat Mater. 9, 324–331. doi: 10.1038/nmat2634

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomida, T., and Egami, T. (1995). Molecular-dynamics study of orientational order in liquids and glasses and its relation to the glass transition. Phys. Rev. B 52, 3290–3308. doi: 10.1103/PhysRevB.52.3290

PubMed Abstract | CrossRef Full Text | Google Scholar

Vogel, H. (1921). The temperature dependence law of the viscosity of fluids. Z Phys. 22, 645–646.

Keywords: liquid, liquid dynamics, relaxation time, medium-range correlation, fragility

Citation: Egami T and Ryu CW (2020) Why Is the Range of Timescale So Wide in Glass-Forming Liquid? Front. Chem. 8:579169. doi: 10.3389/fchem.2020.579169

Received: 01 July 2020; Accepted: 27 August 2020;
Published: 29 September 2020.

Edited by:

Rene A. Nome, State University of Campinas, Brazil

Reviewed by:

Masaru Aniya, Kumamoto University, Japan
Carmelo Corsaro, University of Messina, Italy

Copyright © 2020 Egami and Ryu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Takeshi Egami, ZWdhbWkmI3gwMDA0MDt1dGsuZWR1

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.