Docsity
Docsity

Prepare for your exams
Prepare for your exams

Study with the several resources on Docsity


Earn points to download
Earn points to download

Earn points by helping other students or get them with a premium plan


Guidelines and tips
Guidelines and tips

Pressure's Impact on Ice's Melting Temp: O:H-O Bond Relaxation & Phase Transitions, Study notes of Dynamics

The relationship between pressure and the melting temperature of ice, focusing on the O:H-O bond relaxation and its impact on quasisolid phase boundaries. the compressive effects on the O:H nonbond and H-O bond, the role of quantum calculations and phonon spectrometrics, and various phenomena such as the Mpemba effect and Hofmeister effect. It also mentions studies on water ice and its phase diagram.

What you will learn

  • How do compressive effects impact the O:H nonbond and H-O bond?
  • How do quantum calculations and phonon spectrometrics contribute to understanding water ice?
  • How does pressure affect the melting temperature of ice?
  • What are some phenomena related to water ice and its phase transitions?
  • What is the role of O:H-O bond relaxation in quasisolid phase boundaries?

Typology: Study notes

2021/2022

Uploaded on 09/12/2022

seshadrinathan_hin
seshadrinathan_hin 🇺🇸

4.5

(15)

231 documents

1 / 14

Toggle sidebar

This page cannot be seen from the preview

Don't miss anything!

bg1
Why does pressure melt ice?
Compression shortens the O:H nonbond and lengthens the H-O bond
simultaneously via O-O Coulomb repulsion. The H-O elongation and
its energy loss lower the melting point.
Hydrogen-Bond Relaxation Dynamics: Resolving Mysteries of Water Ice. Coord. Chem. Rev. 2015, 285, 109-65.
Ice Regelation: Hydrogen-Bond Extraordinary Recoverability and Water Quasisolid-Phase-Boundary Dispersivity.
Scientific Reports 2015, 5, 13655.
Hidden Force Opposing Ice Compression. Chem Sci 2012, 3, 1455-60.
Chang Q Sun, NOVITAS, School of Electrical and Electronic Engineering, Nanyang Technological
University, Singapore 639798
Introduction
Regelation, i.e., ice melts under compression and freezes again when the pressure is relieved, remains
puzzling since its discovery in 1850’s by Faraday, James Thomson and his brother William Thomson
(Later Lord Kevin) in 1850’s 1,2. Here we show that hydrogen bond (O:H-O) cooperativity and its
extraordinary recoverability resolve this anomaly. The H-O bond and the O:H nonbond possesses each a
specific heat
x(T/Dx) whose Debye temperature Dx is proportional to its characteristic phonon
frequency x according to Einstein’s relationship. A superposition of the
x(T/Dx) curves for the H-O
bond (x = H, H ~3200 cm-1) and the O:H nonbond (x = L, L ~ 200 cm-1
, DL = 198 K) yields two
intersecting temperatures that define the liquid/quasi-solid/solid phase boundaries. Compression shortens
the O:H nonbond and stiffens its phonon but does the opposite to the H-O bond through O-O Coulomb
repulsion, which closes up the intersection temperatures and hence depress the melting temperature of ice.
Reproduction of the Tm(P) profile clarifies that the H-O bond energy EH determines the Tm with derivative
of EH = 3.97 eV for bulk water and ice. Oxygen atom always finds bonding partners to retain its sp3-
orbital hybridization once the O:H breaks, which ensures O:H-O bond recoverability to its original state
once the pressure is removed.
Faraday1 noted that :
two pieces of thawing ice, if put together, adhere and become one; at a place where liquefaction was
proceeding, congelation suddenly occurs. The effect will take place in air, in water, or in vacuo. It will
pf3
pf4
pf5
pf8
pf9
pfa
pfd
pfe

Partial preview of the text

Download Pressure's Impact on Ice's Melting Temp: O:H-O Bond Relaxation & Phase Transitions and more Study notes Dynamics in PDF only on Docsity!

Why does pressure melt ice?

Compression shortens the O:H nonbond and lengthens the H-O bond

simultaneously via O-O Coulomb repulsion. The H-O elongation and

its energy loss lower the melting point.

Hydrogen-Bond Relaxation Dynamics: Resolving Mysteries of Water Ice. Coord. Chem. Rev. 2015, 285 , 109-65.Ice Regelation: Hydrogen-Bond Extraordinary Recoverability and Water Quasisolid-Phase-Boundary Dispersivity. Scientific Reports 2015, 5 , 13655.Hidden Force Opposing Ice Compression. Chem Sci 2012, 3 , 1455-60.

Chang Q Sun, NOVITAS, School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798

Introduction Regelation, i.e., ice melts under compression and freezes again when the pressure is relieved, remains puzzling since its discovery in 1850’s by Faraday, James Thomson and his brother William Thomson (Later Lord Kevin) in 1850’s 1,2^. Here we show that hydrogen bond (O:H-O) cooperativity and its extraordinary recoverability resolve this anomaly. The H-O bond and the O:H nonbond possesses each a specific heat x (T/Dx ) whose Debye temperature Dx is proportional to its characteristic phonon frequency x according to Einstein’s relationship. A superposition of the x (T/Dx ) curves for the H-O bond (x = H, H ~3200 cm-1^ ) and the O:H nonbond (x = L, L ~ 200 cm-1, DL = 198 K) yields two intersecting temperatures that define the liquid/quasi-solid/solid phase boundaries. Compression shortens the O:H nonbond and stiffens its phonon but does the opposite to the H-O bond through O-O Coulomb repulsion, which closes up the intersection temperatures and hence depress the melting temperature of ice. Reproduction of the T (^) m(P) profile clarifies that the H-O bond energy EH determines the T (^) m with derivative of E (^) H = 3.97 eV for bulk water and ice. Oxygen atom always finds bonding partners to retain its sp 3 - orbital hybridization once the O:H breaks, which ensures O:H-O bond recoverability to its original state once the pressure is removed.

Faraday 1 noted that :

two pieces of thawing ice, if put together, adhere and become one; at a place where liquefaction was proceeding, congelation suddenly occurs. The effect will take place in air, in water, or in vacuo. It will

occur at every point where the two pieces of ice touch; but not with ice below the freezing-point, i.e., with dry ice, or ice so cold as to be everywhere in the solid state ’.

Faraday suggested that there may be a thin liquid-like layer of nascent ice on the surface, ready to be converted to solid on contact with another layer. James Thomson^2 explained this observations in terms of pressure melting based on equilibrium thermodynamics available in his day, and it was his brother, William, who verified the result experimentally^3. This led to a dispute with Faraday, who observed that blocks of ice would stick together by freezing under mild pressure (which one observes with ice cubes in a basket in modern refrigerators). There is a body of modern literature suggesting that Faraday’s surmise of an anomalous ice layer may be correct.

The Regelation can easily be demonstrated by looping a wire around a block of ice with a heavy weight attached to it. This loaded wire melts the local ice gradually until the wire passing through the entire block. The wire's track will refill as soon as it passes, so the ice block will remain solid even after wire passes completely through. Another example is that a glacier can exert a sufficient amount of pressure on its lower surface to lower the melting point of its ice, allowing liquid water flow from the base of a glacier to lower elevations when the temperature of the air is above the freezing point of water. The regelation is exceedingly interesting, because of its relation to glacial action under nature circumstances^4 , in its bearing upon molecular action 5 , and self-repairing of damaged living cells.

It is usual in ‘normal’ materials that compression raises the critical temperature ( T C) at all phase transitions6-8^ ; however, according to the phase diagram of water and ice, the freezing temperature of liquid water is lowered to -22°C by applying 210 MPa pressure; stretching ice (i.e. tensile, or negative, pressure) has the opposite effect - ice melts at +6.5°C when subjected to -95 MPa pressure^9. Conversely, the T (^) C for ice drops from 280 to 150 K at the transition from ordered ice-VIII to proton-disordered ice-VII phase when P is increased from 1 to 50 GPa10-12^. A molecular-dynamics (MD) study of a nanowire cutting through ice suggests that the transition mode and the cutting rate depend on the wetting properties of the wire - hydrophobic and thicker wires cut ice faster^13.

However, a consistent understanding with numerical reproduction of regelation has yet been achieved despite intensive investigations. It might be true that regelation can occur for substances with the property of expanding upon freezing, but mechanisms for both freezing expansion and regelation remain unclear^14. These issues are beyond the scope of classical thermodynamics in terms of equation of states, which inspires alternative ways of thinking and approaching to unlocking these puzzles.

interaction (EL ~ 0.1 eV, left-handed side), the H-O exchange interaction (EH ~ 4.0 eV, right-handed side), and the Coulomb repulsion (C-repulsion) between electron pairs (paring green dots) on oxygen ions. A combination of these interactions with external stimulus dislocates O atoms in the same direction by different amounts. The relaxation proceeds along the potential paths with respect to the H atom (in grey) coordination origin under compression (linked blue spheres) or tension (linked red spheres further moves left). Springs are analogous the respective interactions. The dH0 and d (^) L0 in (b) are the respective segmental length references at 4 C.

Generally, external stimuli, such as stressing and heating modulate the length d ( T , P ) and energy E ( T , P ) of the representative bond along a path denoted f ( T , P ) 6. For instance, compression stores energy into a substance by shortening and stiffening all bonds with possible plastic deformation while tension does the opposite, as illustrated in Figure 1a, and formulated as follows 15 :

  (^)        

0 0 0 0

, 1 1

, 1

T P b (^) T P T V b T^ V b

d P T d t dt p dp t dt p v dv E P T E (^) E

  

 (^)      (^)    (^)     (^)  (^)     (^)     

   

where T 0 and P 0 are the ambient referential conditions. The (t) is the thermal expansion coefficient.

    v  v p  is the compressibility ( p < 0, compressive stress) or extensibility ( p > 0 tensile stress). The

v is the volume of a bond (cross sectional area times length). The ( t ) is the specific heat of the representative bond in Debye approximation. The integration of the ( t ) from 0 K to the melting point ( T m) approximates the bond energy by omitting experimental conditions as the ( t ) for constant volume deviates only 3% from that of constant pressure 15.

O:H-O bond asymmetric and short-range potentials An extended tetrahedron containing two water molecules and four identical O:H-O bonds has unified the length scale and mass density of molecular packed tetrahedrally in water ice on statistic average 22. This extension has also turned out the O:H-O bond with asymmetric, short-range O:H, H-O and O---O interactions, see Figure 1b 23. The O:H-O bond is segmented into a shorter H-O polar-covalent bond with a stronger exchange interaction uH(r) and a longer O:H nonbond with a weaker nonbond interaction u (^) L (r). The two segments are coupled by Couomb repulsion between electron pairs on adjacent oxygen atoms

u (^) C(r) 18,24^. All interactions are limited to the specific segment without any decay acrossing the respective region. The O:H-O bond links the O---O in both the solid and liquid H 2 O phase, regardless of phase structures or topologic configurations 22.

The O:H-O bond performs as an asymmetrical oscillator pair. Under the O---O Coulomb coupling, external excitation such as cooling 14 , compressing 18 , salting 25 , and clustering 24 always relaxes the O:H and H-O in the same direction but by different amounts. Because of the strength disparity between the two segments, compression shortens and stiffens the O:H nonbond (left hand side of the O:H-O bond) and simultaneously lengthens and softens the H-O bond (right hand side). The COlomb repulsion makes the O:H-O bond recover completely its initial states once the compression is relieved. Conversely, once the O:H nonbond breaks, oxygen atom finds immediately bonding partner to retain its sp^3 -orbital hybridization that occurs at temperatures at 5 K 26 and above even in gaseous phase 27.

With the aid of quantum calculations, Lagrangian oscillating mechanics and Fourier fluid thermo dynamics, and phonon spectrometrics, we have been able to consistently and quantitatively resolve quantitatively a few issues such as: 1) Mpemba effect – hot water freezes quicker than its cold 16 , 2) supersolid skins for the slipperiness of ice and the hydrophobic and tough skin of water liquid 28 , 3) ice expansion and mass density oscillation over full temperatures range 14 , 4) anomalies of water molecules with fewer than four nearest neighbors in clusters and droplets^24 , 5) Hofmeister effect – NaCl mediation of O-O repulsion 25 , 6) density-geometry-dimension correlation of molecules packed in water and ice 22 , 7) low compressibility and proton centralization of ice,^18 and, 7) mapping the local potential paths for the O:H-O bond relaxing with stimulus 23 , etc. Progress made insofar has formed the subject of a recent treatise 17.

Results and discussion

O:H-O bond extraordinary recoverability

Figure 2a shows that a molecular dynamics (MD) decomposition of the measured V-P profile of Ice-VIII at 80 K 29 truns out that the d x asymmetric relaxation proceeds until proton symmetrization occuring at 0.22 nm and 60 GPa. The subscript x = H and L reresnet for the H-O and the O:H,respectively. The d L shortens monotonically by 4.3% from 0.1768 to 0.1692 nm and the d H lengthens by 2.8% from 0.0975 to 0.1003 nm when the pressure is increased from 0 to 20 GPa 18. The d L equals the d H at 0.11 nm and 60 GPa, towards proton centrolization in the O:H-O bond 30-32. Figure 2b shows the x cooperative shift of

eV at 60 GPa; the E H decreases monotonically from 3.97 eV to 1.16 eV at 60 GPa. Different from situation of ‘normal’ substance, compression lowers the total energy of the O:H-O bond rather than raise it. The O:H-O bond will fully recover its initial states once the compression is relieved without any plastic deformation.

Table 1. Pressure-dependence of the O:H-O segmental cohesive energy E (^) x and the net gain at each quasi- equilibrium state under compression. Unlike ‘normal’ substance that gains energy with possible plastic deformation under compression, O:H-O bond always losses energy and tends to recover from its higher energy state to lower initial state without any plastic deformation.

P (GPa) E L (eV)^ E H (eV)^ E (^) H+L (P)-E (^) H+L (0) 0 0.046 3.97 (^0) 5 0.098 3.64 (^) -0. 10 0.141 3.39 (^) -0. 15 0.173 3.19 -0. 20 0.190 3.04 (^) -0. 30 0.247 2.63 (^) -1. 40 0.250 2.13 (^) -1. 50 0.216 1.65 (^) -2. 60 0.160 1.16 (^) -2.

As expected, compression shortens the d L , increases the L and E L of the O:H nonbond; the H-O bond responses oppositely to compression, resulting in d H elongation, H and bond energy E H reduction, which can be formulated in the reduced forms as follows (Ex valids at P < 30 GPa; d x / d x0 = 1 + x1 P + x2 P^2 for instance):

2 2 H^2 L^4

H L

H L

d d P P E E

 

   

  ^   

  ^     

  ^ ^ 

E (^) H dictating the T m

The following proves that E H dictates the T (^) m for melting, T (^) mE H. Acceding to eq (2), The T m changes in the following relationship but x = L or H is yet to be known 6 ,

0 0 0 0 0 0

d ( ) 1 ( ) 1 1 1 ( )

V^ P H C x V x^ H P C x x H

pdv s p d dp T P E P (^) dp T P^ ^ ^ E ^ ^ E ^ ^ E

 .

(4)

Eq (3) defines the slope of d (^) x (P):

0 1 2

d (^) [ 2 ] d

d (^) x dx (^) x xP p ^ ^ ^ . (5) Generally, pressure raises the T (^) m but ice responses to pressure in the opposite – T (^) m drops when the pressure is increased. Reproduction of the measured P -dependent T (^) m for melting (Figure 3a) 34 requires that the integral in eq (4) must be positive. Only the d (^) H in Eq. (3) meets this criterion (  (^) x 1 > 0 and  x2 > 0).

Therefore, the H-O bond E H dominates the T m.

Furthermore, matching the T (^) m ( P ) profile using Eq. (5) yields an EH value of 3.97 eV at 0.1 MPa(1 atm pressure) by taking the H atomic diameter of 0.106 nm as the diameter of the H-O bond 35. This E H value agrees with the energy of 4.66 eV for dissociating the H-O bond of water molecules deposited on a TiO 2 substrate with less than a monolayer coverage, and 5.10 eV for dissociating water monomers in the gaseous phase 36. Molecular undercoordination shortens the H-O bond and raises its cohesive energy from the bulk value of 3.97 to 4.66 and to 5.10 eV when the O:H-O bond is subject to molecular undercoordination^25.

However, in dealing with the representative for all bonds of the entire specimen, it is necessary to consider the specific heat per bond that is obtained by dividing the bulk specific heat by the total number of bonds 6. For a specimen of other usual materials, one bond represents all on average; therefore the thermal response is the same for all the bonds, without any discrimination among all bonds in cooling contraction and thermal expansion 38.

For water ice, however, the representative O:H-O bond is composed of two segments with strong disparity in the specific heat of the Debye approximation,  x (T,  Dx ) 14. These two segments response to a thermal excitation differently. Two parameters characterize the specific heat curves each. One is the Debye temperature Dx and the other is the thermal integral of the  x (T,  Dx ) from 0 K to the T mx. The  Dx determines the rate at which the specific-heat curve reaches its saturation. The  x (T,  Dx ) curve of a segment with a relatively lower Dx value reaches saturation more rapidly than the other segment, since the  Dx , which is lower than T (^) mx , is proportional to the characteristic vibration frequency  x of the respective segment, kDx = ћ x , according to Einstein’s relation, 39 where k and ћ are constants.

Conversely, the integral of  x (T,  Dx ) from 0 K to the T mx determines the cohesive energy per bond E (^) x^6. The T mx is the temperature at which the vibration amplitude of an atom or a molecule expands abruptly to more than 3% of its diameter irrespective of the environment or the size of a molecular cluster 39,40^. Thus we have:

0 0 L

/ /E 4.0/0.1 ~ 40

mH mL

DL DH DH L H T T H dt^ L dt^ EH

   

^ ^ ^ ^ ^ 

^ ^ ^ ^ 

 ^ ^ ^ ^ 

Analysis of the temperature-dependence of water surface tension^37 yielded  DL = 198 K < 273 K (T (^) m) and E (^) L = 0.095 eV compared with E (^) H = 3.97 eV for bulk water ice 25. Hence,  DH (^)  16   DL  3200 K. The

O:H specific heat (^)  L ends at 273 K and the H-O specific heat (^)  (^) H ends at T ≥ 3200 K ( T mH). The area

covered by the  H curve is 40 times greater that covered by the  L curve.

The superposition of these two  x (T,  Dx ) curves implies that the heat capacity of water ice differs from that of other, ‘normal’, materials. Such a  x (T,  Dx ) disparity yields temperature regions with different  L /  H ratios over the full temperature range; see Figure 3b. These regions correspond to phases of liquid and solid (  L /  H < 1), and quasisolid (  L /  H > 1). The intersecting temperatures (  L /  H = 1) correspond to

extreme densities at boundaries of the quasisolid phase (viscose and jelly like). The high-temperature boundary corresponds to the maximal density at 4 C and the lower to the crystallization of bulk water.

Numerical and experimental observations 14,17,22^ confirmed that cooling shortens the O:H nonbond in the liquid phase at temperature above 4 C and in the solid phase below 258 K for bulk at different rates because  L /  H < 1 in both regime. However, Cooling shortens the H-O bond in the quasisolid phase (277- 258 K). The other counterpart in the O:H-O bond responses to cooling in the opposite direction. This observation clarifies that the segment with lower x value follows the general rule of thermal expansion and drives the thermal relaxation of the O:H-O bond, which evidences the essentiality of considering the disparity of the specific heat of water ice 14.

One can imagine what will happen to the crossing temperatures if one depresses the DH( H) and E H,

and meanwhile, elevates the DL ( L ) and E L by compression or the inverse. Compression (P > 0) raises the DL and E L by stiffening L , and meanwhile, lowers the DH and E H by stiffening L ; however, tension (P < 0) does the opposite. Figure 3b illustrates how the positive P squeezes the quasisolid phase boundaries. The E H determines approximately the T m through dispersing the upper phase boundary. The Dx ( x ) always relax simultaneously in opposite direction under a given stimulus, which will disperse the quasisolid phase boundaries resulting in the observed ‘superheating/supercooling’, as one often refers. In fact, external stimulus can raise/depress the melting/freezing point by phonon relaxation, which is different from the effect of superheating/supercooling 41.

Once the O:H bond breaks, oxygen atoms will find new partners to retain the sp 3 -orbital hybridization, which is the same to diamond oxidation and metal corrosion – oxygen atoms penetrate into the bulk when corrosion occurs 15,27^. Therefore, O:H-O bond has the strong recoverability for O:H-O bond relaxation and dissociation without any plastic deformation.

Conclusion

Numerical reproduction of the pressure effect on T m clarifies that O:H-O bond relaxation in length, energy, and phonon frequency disperses the quasisolid phase boundaries defined by the supposition of the x (T) curves. Compression stiffens the O:H nonbond and softens the H-O bond, which closes up the separation between the crossing points and depresses the melting temperature of ice. Negative pressure does the opposite to raise the T (^) m. Numerical duplication of the T (^) m(P) gives rise to the H-O bond cohesive energy

16 Zhang, X. et al. Hydrogen-bond memory and water-skin supersolidity resolving the Mpemba paradox. Phys. Chem. Chem. Phys. 16 , 22995-23002 (2014). 17 Huang, Y. et al. Hydrogen-bond relaxation dynamics: resolving mysteries of water ice. Coord. Chem. Rev. 285 , 109-165 (2015). 18 Sun, C. Q., Zhang, X. & Zheng, W. T. Hidden force opposing ice compression. Chem Sci 3 , 1455-1460 (2012). 19 Zhang, X. et al. Water nanodroplet thermodynamics: quasi-solid phase-boundary dispersivity. J. Phys. Chem. B 119 , 5265-5269 (2015). 20 Desiraju, G. R. A Bond by Any Other Name. Ang Chem Int Ed 50 , 52-59 (2011). 21 Falvello, L. R. The Hydrogen Bond, Front and Center. Ang Chem Int Ed 49 , 10045-10047 (2010). 22 Huang, Y. et al. Size, separation, structure order, and mass density of molecules packing in water and ice. Scientific Reports 3 , 3005 (2013). 23 Huang, Y. et al. Hydrogen-bond asymmetric local potentials in compressed ice. J. Phys. Chem. B 117 , 13639-13645 (2013). 24 Sun, C. Q. et al. Density, elasticity, and stability anomalies of water molecules with fewer than four neighbors. J. Phys. Chem. Lett. 4 , 2565-2570 (2013). 25 Zhang, X. et al. Mediating relaxation and polarization of hydrogen-bonds in water by NaCl salting and heating. Phys. Chem. Chem. Phys. 16 , 24666-24671 (2014). 26 Guo, J. et al. support: Real-space imaging of interfacial water with submolecular resolution-Supp. Nat. Mater. 13 , 184-189 (2014). 27 Sun, C. Q. Oxidation electronics: bond-band-barrier correlation and its applications. Prog. Mater Sci. 48 , 521-685 (2003). 28 Zhang, X. et al. A common supersolid skin covering both water and ice. Phys. Chem. Chem. Phys. 16 , 22987-22994 (2014). 29 Yoshimura, Y. et al. High-pressure x-ray diffraction and Raman spectroscopy of ice VIII. J Chem Phys 124 , 024502 (2006). 30 Kang, D. D. et al. Quantum similation of thermally driven phase transition and O k-edge absorption of high-pressure ice. Scientifc Reports 3 , 3272 (2013). 31 Benoit, M., Marx, D. & Parrinello, M. Tunnelling and zero-point motion in high-pressure ice. Nature 392 , 258-261 (1998). 32 Goncharov, A. F., Struzhkin, V. V., Mao, H.-k. & Hemley, R. J. Raman Spectroscopy of Dense H2O and the Transition to Symmetric Hydrogen Bonds. Phys. Rev. Lett. 83 , 1998-2001 (1999). 33 Yoshimura, Y. et al. Convergent Raman Features in High Density Amorphous Ice, Ice VII, and Ice VIII under Pressure. J. Phys. Chem. B 115 , 3756-3760 (2011).

34 Malenkov, G. Liquid water and ices: understanding the structure and physical properties. J. Phys.-Condes. Matter 21 , 283101 (2009). 35 Sun, C. Q. et al. Dimension, strength, and chemical and thermal stability of a single C-C bond in carbon nanotubes. J. Phys. Chem. B 107 , 7544-7546 (2003). 36 Harich, S. A. et al. Photodissociation of H2O at 121.6 nm: A state-to-state dynamical picture. J Chem Phys 113 , 10073-10090 (2000). 37 Zhao, M. et al. Atomistic origin, temperature dependence, and responsibilities of surface energetics: An extended broken-bond rule. Phys. Rev. B 75 , 085427 (2007). 38 Gu, M. X., Zhou, Y. C. & Sun, C. Q. Local bond average for the thermally induced lattice expansion. J. Phys. Chem. B 112 , 7992-7995 (2008). 39 Omar, M. A. Elementary Solid State Physics: Principles and Applications. (Addison-Wesley, 1993). 40 Lindemann, F. A. The calculation of molecular natural frequencies. Phys. Z. 11 , 609-612 (1910). 41 Debenedetti, P. G. & Stanley, H. E. Supercooled and glassy water. Phys. Today 56 , 40-46 (2003).