Universal collapse of the viscosity of supercooled fluids

Universal collapse of the viscosity of supercooled fluids

N. B. Weingartner, C. Pueblo, F. S. Nogueira,, K. F. Kelton, and Z. Nussinov

Department of Physics and Institute of Materials Science and Engineering
Washington University
St. Louis, Missouri 63130, USA
Institute for Theoretical Solid State Physics,
IFW Dresden, PF 270116, 01171 Dresden, Germany
Institut fur Theoretische Physik III
Ruhr-Universitat Bochum, Universitatsstrasse 150
DE-44801 Bochum, Germany

Corresponding author: zohar@wuphys.wustl.edu

All liquids in nature can be supercooled to form a glass. Surprisingly, although this phenomenon has been employed for millennia [1, 2, 3, 4, 5, 6], it still remains ill-understood. Perhaps the most puzzling feature of supercooled liquids is the dramatic increase in their viscosity as the temperature () is lowered. This precipitous rise has long posed a fundamental theoretical challenge. Numerous approaches currently attempt to explain this phenomenon. When present, data collapse points to an underlying simplicity in various branches of science. In this Letter, we report on a 16 decade data collapse of the viscosity of 45 different liquids of all known types. Specifically, the viscosity of supercooled liquids scaled by their value at their respective equilibrium melting temperature ( is, for all temperatures , a universal function of where is a constant that does not change significantly from one liquid to another. This exceptionally plain behavior hints at a link between glassy dynamics and the conventional equilibrium melting transition in all known supercooled fluids.

A long standing endeavor of science is to understand one of the most prevalent and enigmatic states of matter- the “glass” [1, 7, 8, 9]. Any system may be made glassy by supercooling (although the ease in which this can be achieved varies greatly). “Supercooling” refers to a rapid cooling of the liquid below its melting temperature so that crystallization cannot occur. As the temperature of supercooled liquids is depressed, they progressively become increasingly sluggish. At sufficiently low temperatures, when the viscosity () exceeds a threshold value of Pascal second, the resulting system is termed a “glass”. At temperatures below the “glass transition temperature” at which the above threshold is reached, the relaxation times are too long to be readily measured. Given their high viscosity, glasses exhibit a solid like rigidity. However, by comparison to ordered solids, glasses are highly amorphous. The “transition” of liquids into the glassy state is very odd. Most (non-glassy) systems exhibit marked changes in all measurable quantities as they undergo transitions from one phase to another. By stark contrast, while the dynamics of supercooled liquids may change dramatically (the viscosity and relaxation times of supercooled liquids may increase by many orders of magnitude en route to forming glass) [8, 9, 10], their static observables (such as structure and thermodynamic observables) display far milder changes.

Figure 1: The scaling relation of Eq. (1). The“General Curve” is given by Eq. (2).

The most common empirical (Vogel-Fulcher-Tamann-Hesse (VFTH)) [11] fit for the viscosity of supercooled liquids is, by now, nearly a century old. According to this fit, the viscosity at a temperature is where and are liquid specific constants. This functional form asserts the existence of a singular temperature () at which the viscosity diverges; has often been termed an “ideal glass transition temperature”. Since viscosity data at temperatures below have been scarce, it has been hard to examine this putative divergence. To date, the VFTH and all other fits and suggested theories argue for special temperature dependences that are not associated with the simple equilibrium melting transition temperature, e.g., [8, 9, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26]. In essence, the existence of the glassy state has always been ascribed to particular processes that are notably different from those present in the standard equilibrium melting (or freezing) transition. A recent theoretical approach to the glass transition [27] prompted us to critically revisit this deeply rooted dogma and to explore a simple possibility- that of a universal collapse of the viscosity data in which the sole important temperature is that of the standard equilibrium melting transition (). Here, is the temperature at which slow cooling leads to the standard liquid to solid transition (or, more precisely, the “liquidus” temperature above which the system is a fluid). That is, we ask whether the glass transition might be simply related to the standard equilibrium melting transition without the need to invoke, from the outset, assumptions about conjectured novel temperatures (such as ), appealing to constrained dynamics/activation processes in a complex energy landscape, or invoke other numerous considerations. Towards this end, we analyzed 45 different supercooled fluids of all known types (including organic, silicate, metallic, and chalcogenide glassformers) to see whether the temperature dependence of the viscosity is dominated by the melting temperature. Specifically, we tested if a collapse of the form


occurs. Here, is a dimensionless material dependent constant and is a universal function. As seen in Fig. 1 , this relation is indeed satisfied. The values of do not vary widely across different glassformers (see Tables 1 and 2 in the Supporting Information). Much further detail and analysis (at temperatures both below and above the melting transition temperature) appear in [28]. The theory of [27] predicts Eq. (1) with a particularly succinct functional form,


This prediction, highlighted in the collapse of Fig. 1, conforms well to the data. An assumption was made in deriving Eq. (2) that we will briefly comment on below and is elaborated on in great detail in [27]. Regardless of theoretical bias, the raw experimental data illustrate that Eq. (1) holds very well (with a function that is equal to or is very close to Eq. (2) in the broad range of values).

In conclusion, we demonstrated that the conventional equilibrium transition temperature governs the dynamics of a supercooled fluid at temperatures . This finding lies contrary to approaches in which numerous other special temperatures have been suggested to dominate the phenomenology of the glass transition. The work of [27] first proposed Eq. (1) by relying on the observation that the same many body Hamiltonian governs both the equilibrium liquid as well as the supercooled liquid. The core observation of [27] was that the spectral problem posed by this single Hamiltonian has singularities only at the energy densities associated with the equilibrium melting transition (this is the only transition that occurs as the system’s energy density (or, equivalently, the temperature) is lowered). This, in turn, implied that the only singularity governing the behavior of supercooled liquids should also be associated with these energies or, equivalently, with the equilibrium melting temperature. A (scale free) Gaussian distribution in the eigenstate decomposition (in which the temperature was the only dimensional parameter governing the width of this distribution) led to Eq. (2). Albeit suggestive, our results should not be taken as a proof of this theory. Rather, we wish to underscore that the universal collapse of Fig. (1) is, on its own, a feature that underlies all glassformers. Our finding constitutes a challenge to the common lore that the standard equilibrium melting temperature does not play a central role in determining the dynamics of supercooled liquids.

Acknowledgements: NW and ZN were supported by the NSF DMR-1411229. ZN thanks the Feinberg foundation visiting faculty program at Weizmann Institute. CP and KFK were supported by the NSF DMR 15-06553 This work was completed while ZN was at the Aspen Center for Physics, which is supported by the National Science Foundation grant PHY-1066293.

References and Notes

  • [1] L. Berthier and M. D. Ediger, Physics Today 69(1), 40 (2016).
  • [2] R. Zallen, “The Physics of Amorphous Solids”, John Wiley & Sons, Inc., pages 23-32 (1983).
  • [3] A. L. Greer and E. Ma, MRS Bulletin 32, 611 (2007).
  • [4] B. C. Hancock and M. Parks, Pharmaceutical Research 17, 397 (2000).
  • [5] M. Telford, Materials Today 7, 36 (2004).
  • [6] M. Wuttig, and N. Yamada, Nature Materials 6, 824 (2007).
  • [7] P. W. Anderson, Science 267, 1615 (1995).
  • [8] L. Berthier and G. Biroli, Rev. Mod. Phys. 83, 587 (2011).
  • [9] F. H. Stillinger and P. G. Debenedetti, Annu. Rev. Condens. Matter Phys. 4, 263 (2013).
  • [10] C. A. Angell, Science 267, 1924 (1995).
  • [11] H. Vogel, Z. Phys. 22, 645 (1921); G. S. Fulcher, J. Am. Ceram. Soc. 8, 339 (1925); G. Tammann and W. Z. Hesse, Anorg. Allgem. Chem. 156, 245 (1926).
  • [12] G. Adam and J. H. Gibbs, J. Chem. Phys. 43, 139 (1965).
  • [13] T. R. Kirkpatrick, D. Thirumalai, and P. G. Wolynes, Phys. Rev. A 40, 1045 (1989).
  • [14] G. Parisi and M. Mezard, J. Chem. Phys. 111 (1999).
  • [15] V. Lubchenko and P. G. Wolynes, Ann. Rev. of Phys. Chemistry 58, 235 (2007).
  • [16] M. H. Cohen and G. S. Grest, Phys. Rev. B 20, 1077-1098 (1979).
  • [17] E. Leutheusser, Phys. Rev. A 29, 2765 (1984).
  • [18] U. Bengtzelius, W. Gotze, and A. Sjoilander, J. Phys. C 17, 5915 (1984).
  • [19] W. Gotze, Complex dynamics of glass-forming liquids: A mode-coupling theory (Oxford University Press, Oxford, 2008).
  • [20] D. Kivelson, S. A. Kivelson, X. Zhao, Z. Nussinov, and G. Tarjus, Physica A, 219, 27-38 (1995).
  • [21] G. Tarjus, S. A. Kivelson, Z. Nussinov, and P. Viot, J. Phys: Condens Matter, 17, R1143-R1182 (2005).
  • [22] Z. Nussinov, Phys. Rev. B, 69, 014208 (2004).
  • [23] Y. S. Elmatad, D. Chandler, and J. P. Garrahan, J. Phys. Chem. B 113, 5563-5567 (2009).
  • [24] Y. S. Elmatad, R. L. Jack, D. Chandler, Jamd . P. Garrahan, PNAS 107, 12793-12798 (2010).
  • [25] M. Blodgett, T. Egami, Z. Nussinov, and K. F. Kelton, Scientific Reports, Volume 5, id. 13837 (2015).
  • [26] J. C. Mauro, Y. Yue, A. J. Ellison, P. K. Gupta, and D. C. Allan, PNAS, 106, 19780-19784 (2009).
  • [27] Z. Nussinov, https://arxiv.org/pdf/1510.03875.pdf (2015).
  • [28] Nicholas B. Weingartner, Chris Pueblo, F. S. Nogueira, K. F. Kelton, and Zohar Nussinov, https://arxiv.org/pdf/1512.04565.pdf (2015).

Supporting Information

In the tables that follow, the parameter in the collapse of Eq. 1 and Fig. 1 is provided. Further detail concerning the analysis leading to these values and many further aspects (both empirical and theoretical) appear in [28].

          Composition [K] [Pa*s]
          BS2 0.157129 1699 5.570596
          Diopside 0.134328 1664 1.5068
          LS2 0.170384 1307 22.198
          OTP 0.069685 329.35 0.02954
          Salol 0.087192 315 0.008884
          Anorthite 0.131345 1823 39.81072
           0.234171 1450 0.01564
           0.154701 1030 0.030197
           0.187851 1350 0.03643
           0.124879 1058 0.0446
          Albite 0.103344 1393 24154952.8
           0.142960 1230 0.021
           0.209359 1283 0.0269
           0.167270 1220 0.0233
           0.109488 1483 0.042
          Vit106a 0.133724 1125 0.131
           0.144521 1193 0.0266
           0.133069 273.15 0.001794
          Glucose 0.079455 419 0.53
          Glycerol 0.108834 290.9 1.9953
           0.185389 1279.226 0.01652
           0.21073 1350.789 0.02288
           0.169362 1363.789 0.04805
          NS2 0.134626 1147 992.274716
           0.103380 1125.409 0.04516
           0.157480 1313 0.01155
           0.156955 1198 0.02044535
Table 1: Values of Relevant Parameters for all liquids studied
          Composition [K] [Pa*s]
           0.244979 1233 0.02625234
           0.148249 1152 0.0268
          D Fructose 0.050124 418 7.31553376
          TNB1 0.07567 472 0.03999447
          Selenium 0.130819 494 2.9512
          CN60.40 0.149085 1170 186.2087
          CN60.20 0.161171 1450 12.5887052
           0.137623 1071 0.03615283
           0.118631 1173 0.03797
           0.137753 1185 0.0256
           0.166699 1226 0.02162
          Isopropylbenzene 0.073845 177 0.086
          ButylBenzene 0.085066 185 0.0992
           0.131969 1199 0.02526
          Vit 1 0.111185 937 36.59823
          Trehalose 0.071056 473 2.71828
          Sec-Butylbenzene 0.080088 190.3 0.071
           0.090948 1873 1.196x
Table 2: Values of Relevant Parameters for all liquids studied (continued)
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description