Abstract
The dynamics of glass is of importance in materials science but its nature has not yet been fully understood. Here we report that a verification of the temperature dependencies of the primary relaxation time or viscosity in the ultraslowing/ultraviscous domain of glass-forming systems can be carried out via the analysis of the inverse of the Dyre–Olsen temperature index. The subsequent analysis of experimental data indicates the possibility of the self-consistent description of glass-forming low-molecular-weight liquids, polymers, liquid crystals, orientationally disordered crystals and Ising spin-glass-like systems, as well as the prevalence of equations associated with the ‘finite temperature divergence’. All these lead to a new formula for the configurational entropy in glass-forming systems. Furthermore, a link to the dominated local symmetry for a given glass former is identified here. Results obtained show a new relationship between the glass transition and critical phenomena.
Similar content being viewed by others
Introduction
The glass transition constitutes one of the outstanding problems of condensed matter physics1, with fundamental importance for many areas ranging from geophysics to modern food science and material engineering2,3,4,5. Over the last decades, an important progress in this field has been achieved, primarily because of the progress in the state-of-the-art experimental, theoretical and numerical insights. Hence, a conceptual breakthrough seems to emerge6,7,8,9,10,11,12,13,14,15,16,17,18.
Recently, however, serious doubts regarding the most fundamental experimental artifact, namely the parameterization of viscosity η(T) or primary relaxation time τ(T) on approaching the glass temperature (Tg), appeared15,19,20,21. The importance of this issue is related to the fact that it is often associated with a hypothetical link between the glass transition and critical phenomena physics, that is, the thrilling possibility of a new great unification within condensed matter physics14,18. Moreover, relations portraying the evolution of η(T) or τ(T) are often considered as a checkpoint for theoretical models3,13. Decades of studies have led to the prevailing conviction that the ultimate parameterization is possible via the Vogel–Fulcher–Tammann (VFT) equation3,13,22,23,24:
where DT denotes the fragility strength coefficient and T0 is the VFT singular temperature ( and ).
However, by following the compelling analysis of 43 sets of τ(T) data for low-molecular-weight ultraviscous liquids, Hecksher et al.15 recently concluded: ‘There is no compelling evidence for the VFT prediction that the relaxation time diverges at a finite temperature […] with a dynamic divergence of the VFT form lack a direct experimental basis…. we suggest that in the search for the correct theory for ultraviscous liquid dynamics, theories not predicting a dynamic divergence of the VFT form should be focused on’. McKenna25 extended this outlook to the case of polymer glass formers. These findings were strengthened by the derivation by Mauro et al.19 of the equation introduced empirically by Waterton26:
The subsequent analysis of experimental data confirmed the advantage of the Waterton-Mauro (WM, called also MYEGA) equation without a ‘finite temperature divergence’ over the VFT one11,19,20,21,27. Thus, the question regarding τ(T) or η(T) portrayal seemed to be clarified. However, the problem becomes more complex if glass-forming liquid crystals (LCs) and orientationally disordered crystals (ODICs; plastic crystals) are also considered. For these materials, the prevalence of the ‘critical-like’ equation with the clear ‘finite temperature divergence’ was clearly shown28,29:
where and : most often Tg−T0≈30–40 K and Tg−TC≈10–15 K.
Surprisingly, the latter description seems to be applicable also for selected low-molecular-weight liquids and polymers28,29. Worth mentioning is that such dynamics is also commonly accepted for spin-glass-like (SGL) systems5,28. Remarkable is furthermore the approximate coincidence of power exponents in the mentioned systems, namely φ≈9 for LC28,29, φ=9–15 in ODIC28,29 and φ=9–12 in SGL3,13,28.
In this respect, it is important to point out the link between anomalous dynamics in supercooled glass-forming systems and the cooperatively rearranged regions, introduced by the Adam and Gibbs (AG) model3,13,30, leading to the following output relation30:
where Δμ defines that the free-energy barrier between cooperatively rearranged region and SC is the configurational entropy related to the difference between the entropy of the metastable disordered supercooled system and the stable crystal state.
Angell et al.31,32,33 noted that the AG equation can be converted into the VFT one by assuming:
for T0=TK, where the latter denotes the Kauzmann temperature34.
Although equation (5) became popular in the semi-phenomenological analysis of supercooled glass formers3,13, some puzzling issues remain. Mauro11,19 has shown that the ‘divergence-less’ equation (2) can also be derived from the AG dependence by omitting equation (5) for the configurational entropy. Moreover, the analysis of experimental data indicated that the coincidence between T0 and TK may be limited35. Eckmann and Procaccia36 showed that for a soft potential based model, the configurational entropy can stay positive for , concluding that the Kauzmann temperature may not exist.
Particularly noteworthy is that the parameterization of τ(T) or η(T) via ‘finite temperature divergence’ relations was used as an important argument to support the hypothesis, which correlates the glass transition with a hidden phase transition below Tg. For instance, Colby16 suggested that equation (3), with exponent φ=9, may be ‘universal’ for glass-forming low–molecular-weight liquids and polymers, although for some of them, a multiplicative factor exp(E/kT) should be included. The ‘critical-like’ divergence was associated with a hypothetical phase transition at TC<Tg, linked to the correlation length ξ(T)=ξ0(T−TC)−v with (refs 16, 17). However, a subsequent thorough analysis of τ(T) experimental data did not confirm the suggested superior-fitting features28,29,37. Notwithstanding, there are several low-molecular-weight and polymeric glass formers for which the behaviour suggested by Colby holds surprisingly well28,29.
It is worth mentioning that in 2010, Tanaka et al.18 carried out an in-depth simulation analysis of six glass-forming systems and stated: ‘Our results suggest a far more direct link than thought before between glass transition and critical phenomena. Indeed, the glass transition may be a new type of critical phenomenon where a structural order parameter is directly linked to slowness’. Tanaka et al.18 linked the ‘pretransitional’ behaviour of the relaxation time to the Ising-like critical power-law divergence of the correlation length towards the ideal glass transition with the exponent , where d defines the spatial dimensionality. This led to the VFT-type equation . Tanaka et al.18, further suggested that glass transition may be encountered in the ‘Ising universality class for systems with a non-conserved order parameter’.
At this point, some basic concepts of critical phenomena physics are worth recalling38,39. The pretransitional behaviour for different physical magnitudes (X(T)) is described by the relation X(T)|T−TC|λ, where the values of critical exponents (λ) are universal and depend solely on the dimensionality of the order parameter, the space dimensionality and the range of interactions. The symmetry is playing an essential role, whose changes are quantified by the order parameter. Microscopically, different systems can be ascribed in universality classes. It is also worth to bear in mind that the physics of critical phenomena is considered one of the greatest ‘unification’ successes in condensed matter physics of the 20th century2,38,39. To the best of the authors’ knowledge, no clear manifestation of the above fundamental features for the glass-transition phenomenon has been reported so far. This paper aims to fill this gap and resolve the Gordian knot of the τ(T) or η(T) evolution on approaching the glass transition.
Results
Fragility and the Dyre–Olsen temperature index
The Super-Arrhenius (SA) evolution of relaxation time or viscosity in the ultraslowing/ultraviscous domain can be generally expressed as , where stands for the normalized apparent activation energy (AAE) and R denotes the gas constant3,13. The departure from the simple Arrhenius behaviour with Ea=const is characterized by the fragility coefficient (m), introduced by Angell et al.40,41 as the metric at the normalized plots log10τ or log10η versus Tg/T:
Although the fragility coefficient m provides a metric for glass-forming liquids within a wide diversity of materials, it is not the only parameter of this type3,13. A decade ago, Dyre–Olsen (DO)42 noted the formal importance of the Grüneisen parameter (γ)43:
It indicates that if the restoring force acting on an atom is non-linear in displacement, the individual phonon frequency ωi changes with the volume of crystal V, and parameter γ reaches a transition from zero (linear displacement) to γ>0 (ref. 43). In the SA domain, the AAE increases with a decreasing temperature, allowing a dramatic increase from constant values to , hence giving rise to a non-linear behaviour in the log10τ versus 1/T representation. This formal similarity allowed to propose the metric of the SA behaviour by the ‘Grüneisen-style’ activation energy temperature index42:
The DO temperature index can serve as the alternative metric of fragility, which can be also presented as , where and denote the apparent activation enthalpy44.
Assuming the hypothetical ‘universal’ value for the prefactor τ0=10−14 s (refs 3, 13, 15), one obtains: , for τ(Tg)=102 s and τ0=10−14 s. Subsequently, taking into account the recent suggestion that 16<m<175 (ref. 45), one can estimate: . It is noteworthy that in 1997, Simon et al.46 indicated the upper fragility limit m=214, yielding IDO(Tg)=12.375. Very recently, McKenna47 indicated that if T0 is an ideal glass transition, the value of m evaluated at T=T0=Tg would tend towards infinity for all materials. The same can be expected for IDO(Tg).
Heksher et al.15 showed that the analysis of IDO(T) may represent a decisive tool for testing the validity of the given equation for fitting τ(T) experimental data in low-molecular-weight glass-forming ultraviscous liquids. To quantify how exactly the VFT and the popular Bässler–Avramov48,49 (BA) equations describe the data, they compared experimental and ‘theoretical’ values of IDO(T), namely15:
Derivations of IDO(T) for equations (2) and (3) are presented in the Methods section. In the ultraslowing domain, near Tg, they are as follows:
A glimpse on equations (9) and (10) shows that reciprocals of DO indexes yield a linear dependence , namely:
The discussed evolutions of DO indexes are shown graphically in Fig. 1. Two characteristic features are worth mentioning, namely: and . It is worth stressing that while the physics associated with the fragility index m has been very broadly discussed3,13,15,20,31,40,41,42,43,44,45,50,51, the meaning and properties of the DO index is poorly explored so far15,42.
For the estimation of IDO(T), the evolution of Ea(T) has to be determined. This, however, constitutes a notable issue. Hecksher et al.15 calculated Ea(T) from the SA relation, assuming τ0=10−14 s as an ‘universal averaged’ value of the prefactor. However, in the systems analysed in this paper τ0=10−11 s÷10−16 s (refs 28, 29, 44). In the Methods section, we present a novel way of determining Ea(T), avoiding the biasing influence of poorly estimated values of the prefactor τ0 for the given system.
Experimental data analysis
This paper focuses on the ultraslow/ultraviscous domain as in Hecksher et al.15, for or , which usually occurs for (refs 50, 51, 52, 53). By heating above TB, the so-called dynamic cross-over occurs. The discussion on recent advances regarding this phenomenon can be found in Martinez-Garcia et al.52 The analysis of experimental data presented below does not only cover the most often discussed supercooled low-molecular-weight liquids (L) and multimeric and polymeric (P) glass formers as in Hecksher et al.15, but also LCs, plastic crystals (ODICs) and SGL systems. The analysis is based on 16 selected sets of τ(T) experimental data representing these types of glass-forming systems. However, to clarify the discussion, the key results are presented for a smaller number of characteristic systems. They are L including diethyl phthalate (DEP), glycerol, sorbitol, ethanol28,54, P including EPON 828 and polystyrene (MW=700) (refs 28, 55), ODIC (plastic crystal) including C7-OH, C8-OH (ref. 53), CNc6 (ref. 57), CNadm58,59 and NPaNPG0.3 (ref. 56), and SGL including ferrofluid liquid, with 5% solution of single-domain particles of the amorphous alloy Fe1−xCx (x=0.2–0.3) (ref. 57). Details regarding all tested systems are collected in the Supplementary Table S1. The presented analysis is focused on the primary relaxation time τ(T), because the modern broad-band dielectric spectroscopy can yield sets of high precision and consistent experimental data3,13. For the presentation of experimental data, we apply the linearized derivative-based analysis introduced in Drozd-Rzoska and Rzoska44. The VFT equation allows the following transformation of τ(T) data:
For the critical-like equation (3):
where stands for the normalized apparent enthalpy.
The linear behaviour of equations (12) and (13) show the valid temperature domain of the VFT and ‘critical-like’ equations. The subsequent linear regression fit can yield optimal values of relevant parameters44. Results of such analysis for selected systems are shown in Fig. 2a–e. It is visible that a clear ‘critical-like’ parameterization is preferred for all LC rod-like compounds, as well as for sorbitol, ethanol and polystyrene (Fig. 2b). This behaviour can be linked to the uniaxial structure of molecules, which may lead to the local orientational symmetry (OS) arrangements28. The ‘critical-like’ behaviour also dominates the dynamics of plastic crystals (ODICs)28,29, where a positional symmetry (PS) is a fundamental feature (Fig. 2a). Regarding EPON 828, glycerol and DEP, the critical-like parameterization definitively fails (Fig. 2c). Because these compounds cannot be linked neither to OS nor to PS cases, we assigned them as no-symmetry (NS) compounds.
For ‘NS’ glass formers, the VFT equation fits the experimental data fairly well as shown in Fig. 2e. It is noteworthy that for ethanol and sorbitol, both equations (1) and (3) seem to offer the acceptable parameterization (Fig. 2b).
For a further insight into the optimal τ(T) parameterization, an analysis exploring behaviour has been carried out (Fig. 3 and equation (11)). It is visible that for none of the tested compounds, the evolution expected for BA49 () or WM19 () relations take place. Solely, patterns expected for the VFT or the ‘critical-like’ description occurs. Experimental results show that such patterns are associated with three ‘terminal, convergence’ cases: (i) , for ODICs and SGL (PS case); (ii) n≈1.53, for all LCs and also low-molecular-weight sorbitol, propanol and polystyrene (OS case); (iii) n≈1, for DEP, EPON 828 and glycerol (NS case). Worth to mention is that the points (i) and (ii) are related to approximately the same value of the ‘critical’ exponent φ=8.5~12 in equation (3), being remarkable that the plot allows to separate OS and PS cases.
A summary of the analysis covering 16 glass-forming systems is presented in Fig. 4 and in the Supplementary Table S1. Apart from the already-mentioned manifestation of three ‘characteristic’ values of the coefficient n, which can be related to the local symmetry, two more issues are notable. First, the clear access of the SGL system to the PS group, together with ODICs. Second, the ‘intermediate’ cases between three characteristic values of the coefficient n. For ethanol, the coefficient is n≈1.2, whereas for propanol, being a ‘longer’ and thus more uniaxial compound and from the same homologous series, the coefficient is n≈1.5, as for LCs and polystyrene. One of the ‘intermediate’ cases is cyclooctanol (plastic crystal), for which n≈0.8. In this case, the discrepancy may be linked to some distortion from the optimal positional ordering.
Consequences for the configurational entropy
Assuming the validity of equation (4), recalling equation (18) from Martinez-Garcia et al.52 and the definition of the DO index by equation (8), the following relation can be formulated:
By applying experimental findings from the previous section, the above equation can be written as:
The latter directly yields the following solution:
where the ‘singular’ temperature can be defined as and the power exponent as .
Equation (16) represents a new generalized dependence for the evolution of the configurational entropy of glass-forming liquids, allowing to construct the normalized/scaling plot shown in Fig. 5.
The scaling Fig. 5 covers the variety of ultraslowing glass-forming systems ranging from low-molecular liquids and polymers to plastic crystals, LCs and SGLs. Clearly visible are three characteristic (universal?) values of the power exponent: n≈1.53, n=1 and n≈0.17, which can be linked to different dominating molecular symmetries, as shown in the right part of Fig. 5. Worth to mention is also that equation (16) has been derived by the solution of the first-order differential equation (15), without underlying assumptions that were used in the derivation of equation (5). The latter is recovered for n=1. It is also noteworthy that equation (16) indicates a possible link between the ‘critical’ temperature (TC) in equation (3) and the VFT singular temperature (T0) in equation (1)).
Discussion
In the year 2008, Hecksher et al.15 boosted the discussion regarding one of the most fundamental experimental artifacts within glass-transition physics, namely the validity of the VFT equation, the existence of the ‘final temperature divergence’ <Tg, and the ultimate way of τ(T) or η(T) parameterization. This discussion and subsequent papers28,29,52 have indicated that parameterizations not associated with a ‘finite temperature divergence’ may be crucial. All these findings indirectly questioned the connection between the glass transition and a hidden phase transition below Tg. However, parallel studies exploring the linearized derivative and distortion sensitive analysis led to a clear proof that the ‘critical-like’ description with a finite temperature divergence may be optimal for liquid crystalline glass formers, plastic crystals (not tested in Erwin and Colby17) and even for selected low-molecular liquids and polymers27,28.
This paper shows that the ‘finite temperature divergence’ in ultraslowing glass-forming systems on approaching Tg may be feasible not only for low-molecular-weight liquids and polymers but also for LCs, plastic crystals and even SGL systems. Surprisingly, a consistency seems to be possible, despite the fact that a single equation for portraying the dynamics in the ultraslowing domain for an arbitrary glass former may not exist. It is worth stressing that the optimal parameterization via the VFT equation should be considered solely for selected compounds, and other systems such as parameterization is inherently non-optimal. This fact ought to be taken into account when validating different equations for portraying the τ(T) or η(T) experimental behaviour.
We would like to stress the results of the analysis focused on BA48,49 and WM (MYEGA) equations11,19,26 without the finite temperature divergence.
Worth to mention is also a group of further equations of such type2,3,13,15, which should be tested in the future.
In conclusion, the obtained key forms for τ(T) or η(T) parameterization can be linked to predefined values of the parameter/exponent , as well as to the dominating local symmetry. They are clearly associated with a finite temperature divergence below the glass temperature. All these results indicate qualitatively new links between glass transition and critical phenomena.
Methods
Analysis of experimental data
The analysis of ‘dynamic’ experimental data employing higher-order derivatives of τ(T) or η(T) experimental data is most often avoided due to scatter associated with the inherent experimental error. In this paper, this basic issue has been overcome due to the novel implementation of Savitzky–Golay (SG) filtering procedure introduced in Martinez-Garcia et al.52
Determination of the AAE from τ(T) experimental data
For estimating the DO temperature index (IDO(T)), the experimental evolution of the AAE is needed (equation (8)). So far, it was calculated from the ‘general’ apparent SA equation via , assuming for τ0=10−14 s, which is most often considered as the ‘typical’ or ‘universal’ value, at least for ultraviscous liquids. Notwithstanding, there is clear evidence that the value of τ0 ranges from 10−11–10−16 s (ref. 49). The erroneous value τ0 for a given experimental system can notably influence estimations of Ea(T). In this paper, experimental values of Ea(T) were calculated in a novel way, avoiding τ0-related problems. The AAE was calculated via . Ultimate estimation was supported by the SG filtering procedure. Subsequently, the following ‘general’ first-order differential equation was solved:
The ‘noise-related’ scatter of the calculated Ea(T) values was minimized by the SG filtering procedure for each tested glass-forming system. In a next step, the DO was calculated via the equation (8).
Derivations of the DO temperature indexes for selected equations describing τ(T) or η(T) evolutions
The form of IDO(T) for VFT equation (1) and BA relation can be found in Erwin and Colby17. The DO temperature index for the WM (MYEGA) equation can be obtained by linking equation (2), the SA equation and IDO(T) definition (equation (8)):
Forms of IDO(T) for the critical-like case in the ultraslowing/ultraviscous dynamic domain can be obtained by linking equation (3) and the general SA dependence :
This work focuses on the ultraviscous/ultraslowing dynamic domain close to Tg, that is, for T→TC, hence:
The latter equation can be re-written as:
The function can be expanded into the Taylor series because of its inherent features () if |δ(T)−1|≤1. This condition is valid for ultraslowing ODICs, LCs and SGL systems tested in this paper. Taking this into account, one obtains:
Then, recalling IDO(T) definition (equation (8)):
The evolutions of and are presented graphically in Fig. 6.
The analytic form of the index is given in the bracket in equation (24) and shown graphically in Figure 6, together with the approximation . This figure is supplemented by the experimental DO index for C7-OH (ODIC; plastic crystal) where the critical-like parameterization of τ(T) is optimal (see Fig. 2). The validity of for portraying experimental data in the ultraslowing domain (below TB) is clearly visible.
Additional information
How to cite this article: Martinez-Garcia, J. C. et al. A universal description of ultraslow glass dynamics. Nat. Commun. 4:1823 doi: 10.1038/ncomms2797 (2013).
References
Kennedy, D. 125th Anniversary Issue: 125 outstanding problems in all of science: what is the nature of the glassy state. Science 309, 83 (2005).
Jones, R. A. L. Introduction to Soft Matter Oxford University Press: NY, (2005).
Donth, E. The Glass Transition. Relaxation Dynamics in Liquids and Disordered Materials. Vol. 48 (, Springer Series in Material Sci. II: Springer, Berlin, (1998).
Poirier, J.-P. Introduction to the Physics of the Earth’s Interior Cambridge University: Cambridge, (2000).
Mezzenga, R., Schurtenberger, P., Burbidge, A. & Michel, M. Understanding foods as soft materials. Nat. Mater. 4, 729–740 (2005).
Andersen, H. C. Molecular dynamics studies of heterogeneous dynamics and dynamic crossover in supercooled atomic liquids. Proc. Natl Acad. Sci. USA 102, 6686–6691 (2005).
Lubchenko, V. & Wolyness, P. G. Theory of Structural Glasses and Supercooled Liquids. Ann. Rev. Phys. Chem. 58, 235–266 (2007).
Kivelson, S. A. & Tarjus, G. In search of a theory of supercooled liquids. Nat. Mater. 7, 831–833 (2008).
Roland, C. M. Characteristic relaxation times and their invariance to thermodynamic conditions. Soft Matter 4, 2316–2322 (2008).
Hedges, L. O., Jack, R. L., Garrahan, J. & Chandler, D. A. Dynamic order–disorder in atomistic models of structural glass formers. Science 6, 1309–1313 (2009).
Mauro, J. C. Through a glass, darkly: dispelling three common misconceptions in glass science. Int. J. Appl. Glass Sci. 2, 245–261 (2011).
Angell, A. & Klein, I. S. Glass physics: prigogine and defay say relax. Nat. Phys. 7, 750–751 (2011).
Ngai, K. L. Relaxation and Diffusion in Complex Systems Springer: Berlin, (2011).
Tanaka, H. Bond orientational order in liquids: Towards a unified description of water-like anomalies, liquid-liquid transition, glass transition, and crystallization. Europ. Phys. J. E 35, colloquium 1–84 (2012).
Hecksher, T., Nielsen, A. I., Olsen, N. B. & Dyre, J. C. Little evidence for dynamic divergences in ultraviscous molecular liquids. Nat. Phys. 4, 737–741 (2008).
Colby, R. H. Dynamic scaling approach to glass formation. Phys. Rev. E 61, 1783–1792 (2000).
Erwin, B. M. & Colby, R. H. Temperature dependence of relaxation times and the length scale of cooperative motion for glass-forming liquids. J. Non-Cryst. Solids 307, 225–231 (2002).
Tanaka, H., Kawasaki, T., Shintani, H. & Watanabe, K. Critical-like behaviour of glass-forming liquids. Nat. Mater. 9, 324–331 (2010).
Mauro, J. C., Yue, Y., Ellison, A. J., Gupta, P. K. & Allan, D. C. Viscosity of glass-forming liquids. Proc. Natl Acad. Sci. USA 106, 19780–19784 (2009).
Lunkenheimer, P., Kastner, S., Köhler, M. & Loidl, A. Temperature development of glassy α-relaxation dynamics determined by broadband dielectric spectroscopy. Phys. Rev. E 81, 051504 (2010).
Martinez-Garcia, J. C., Tamarit, J. Ll. & Rzoska, S. J. Enthalpy space analysis of the evolution of the primary relaxation time in ultraslowing systems. J. Chem. Phys. 134, 024512–024517 (2011).
Vogel, H. Temperaturabhängigkeitsgesetz der Viskosität von Flüssigkeiten. Phys. Zeit 22, 645–646 (1921).
Fulcher, G. S. Analysis of recent measurements of the viscosity of glasses. J. Am. Ceram. Soc. 8, 339–355 (1925).
Tammann, G. Glasses as supercooled liquids. J. Soc. Glass Technol. 9, 166–185 (1925).
McKenna, G. B. Diverging views on glass transition. Nat. Phys. 4, 673–674 (2008).
Waterton, S. C. The viscosity-temperature relationship and some inferences on the nature of molten and of plastic glass. J. Soc. Glass Technol. 16, 244–247 (1932).
Rouxel, T. Thermodynamics of viscous flow and elasticity of glass forming liquids in the glass transition range. J. Chem. Phys. 135, 184501–184516 (2011).
Drozd-Rzoska, A., Rzoska, S. J. & Paluch, M. Universal critical-like scaling of dynamic properties in symmetry-selected glass formers. J. Chem. Phys. 129, 184509–184515 (2008).
Drozd-Rzoska, A., Rzoska, S. J., Pawlus, S., Martinez-Garcia, J. C. & Tamarit, J. -L. Evidence for critical-like behavior in ultraslowing glass-forming systems. Phys. Rev. E 82, 031501 (2010).
Adam, G. & Gibbs, J. H. On the temperature dependence of cooperative relaxation properties in glass-forming liquids. J. Chem. Phys. 43, 139–146 (1965).
Angell, C. A. Entropy and fragility in supercooling liquids. J. Res. Nat. Inst. Stand. Technol. 102, 171–185 (1997).
Richert, R. & Angell, C. A. Dynamics of glass-forming liquids. V. On the link between molecular dynamics and configurational entropy. J. Chem. Phys. 108, 9016–9027 (1998).
Martinez, L.-M. & Angell, C. A. A thermodynamic connection to the fragility of glass-forming liquids. Nature 410, 663–667 (2001).
Kauzmann, W. The nature of the glassy state and the behavior of liquids at low temperatures. Chem. Rev. 43, 219–256 (1948).
Tanaka, H. Relation between thermodynamics and kinetics of glass-forming liquids. Phys. Rev. Lett. 90, 05570 (2003).
Eckmann, J. P. & Procaccia, I. Ergodicity and slowing down in glass-forming systems with soft potentials: no finite-temperature singularities. Phys. Rev. E 78, 011503 (2008).
Richert, R. Scaling vs. Vogel-Fulcher-type structural relaxation in deeply supercooled materials. Physica A 287, 26–36 (2000).
Stanley, H. E. Introduction to Phase Transitions and Critical Phenomena Oxford Univ. Press.: NY, (1987).
Anisimov, M. A. Critical Phenomena in Liquids and in Liquid Crystals Gordon and Breach: Reading, (1992).
Angell, C. A. in Relaxations in Complex Systems (ed. Ngai K. L. 3NRL: Washington, DC, (1985).
Böhmer, R., Ngai, K. L., Angell, C. A. & Plazek, J. D. Nonexponential relaxations in strong and fragile glass formers. J. Chem. Phys. 99, 4201–4209 (1993).
Dyre, J. C. & Olsen, N. B. Landscape equivalent of the shoving model. Phys. Rev. E 69, 042501 (2004).
Grüneisen, E. Theorie des festen Zustandes einatomiger elemente. Ann. Phys. 12, 257–306 (1912).
Drozd-Rzoska, A. & Rzoska, S. J. On the derivative-based analysis for temperature and pressure evolution of dielectric relaxation times in vitrifying liquids. Phys. Rev. E 73, 041502 (2006).
Wang, L.-M. & Mauro, J. C. An upper limit to kinetic fragility in glass-forming liquids. J. Chem. Phys. 134, 044522 (2011).
Simon, S. L., Plazek, D. J., Sobieski, J. W. & McGregor, E. T. Physical aging of a polyetherimide: Volume recovery and its comparison to creep and enthalpy measurements. J. Polym. Sci. Part B. Polym. Phys. 35, 929–936 (1997).
McKenna, G. B. Deformation and flow of matter: Interrogating the physics of materials using rheological methods. J. Rheol. 56, 113–158 (2012).
Bässler, H. Viscous flow in supercooled liquids analysed in terms of transport theory for random media with energetic disorder. Phys. Rev. Lett. 58, 767–770 (1987).
Avramov, I. Viscosity in disordered media. J. Non-Cryst. Solids 351, 3163–3173 (2005).
Stickel, F., Fischer, E. W. & Richert, R. Dynamics of glass‐forming liquids. I. Temperature‐derivative analysis of dielectric relaxation data. J. Chem. Phys. 102, 6251 (1995).
Richert, R. & Angell, C. A. Dynamics of glass-forming liquids. V. On the link between molecular dynamics and configurational entropy. J. Chem. Phys. 108, 9016–9026 (1998).
Martinez-Garcia, J. C., Martinez-Garcia, J., Rzoska, S. J. & Hulliger, J. The new insight into dynamic crossover in glass forming liquids from the apparent enthalpy analysis. J. Chem. Phys. 137, 064501 (2012).
Martinez-Garcia, J. C. et al. Disentangling the secondary relaxations in the orientationally disordered mixed crystals: cycloheptanol+cyclooctanol two-component system. J. Phys. Chem. B 114, 6099–6106 (2010).
Brand, R., Lunkenheimer, P., Schneider, U. & Loidl, A. Excess wing in the dielectric loss of glass-forming ethanol: a relaxation process. Phys. Rev. B 62, 8878–8883 (2000).
Schwartz, A., Colmenero, J. & Alegria, A. Dielectric study of the segmental relaxation of low and high molecular weight polystyrenes under hydrostatic pressure. J. Non-Cryst. Solids 353, 4298–4302 (2007).
Drozd-Rzoska, A., Rzoska, S. J., Pawlus, S. & Tamarit, J.Ll. Dynamics crossover and dynamic scaling description in vitrification of orientationally disordered crystal. Phys. Rev. B 73, 224205 (2006).
Jönsson, P., García-Palacios, J. L., Hansen, M. F. & Norblad, P. Relaxation in interacting nanoparticle systems. J. Mol. Liq. 114, 131–140 (2004).
Martinez-Garcia, J. C. et al. Universal critical-like scaling of dynamics plastic crystals. J. Non-Cryst. Solids 357, 329–333 (2011).
Martinez-Garcia, J. C. et al. α-relaxation dynamics of orientanionally disordered mixed crystals composed of Cl-adamantane and CN-adamantane. J. Chem. Phys. 132, 164516 (2010).
Acknowledgements
S.J.R. and A.D.-R. were supported by the Grant no. N N202 231737 from National Centre for Science (NCN, Poland). We would like to thank Iris Martinez-Heer for her valuable contribution to improve the final version of the article.
Author information
Authors and Affiliations
Contributions
J.C.M.-G. and S.J.R. conceived the method and created the final shape of the paper. A.D.R. and S.J.R. carried out experiments associated with some sets of experimental data. A.D.R. yielded the linearized derivative-based analysis. J.M.-G. was responsible for the Savitzky–Golay methodology.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Information
Supplementary Table S1 and Supplementary References (PDF 80 kb)
Rights and permissions
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License. To view a copy of this license, visit http://creativecommons.org/licenses/by-nc-sa/3.0/
About this article
Cite this article
Martinez-Garcia, J., Rzoska, S., Drozd-Rzoska, A. et al. A universal description of ultraslow glass dynamics. Nat Commun 4, 1823 (2013). https://doi.org/10.1038/ncomms2797
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms2797
This article is cited by
-
New paradigm for configurational entropy in glass-forming systems
Scientific Reports (2022)
-
Activation volume in superpressed glass-formers
Scientific Reports (2019)
-
The generalized Vogel-Fulcher-Tamman equation for describing the dynamics of relaxor ferroelectrics
Scientific Reports (2019)
-
Fragility and basic process energies in vitrifying systems
Scientific Reports (2015)
-
Divergent dynamics and the Kauzmann temperature in glass forming systems
Scientific Reports (2014)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.