玻璃與玻璃化轉(zhuǎn)變

李艷偉　孫昭艷　安立佳

(中國科學(xué)院長春應(yīng)用化學(xué)研究所，高分子物理與化學(xué)國家重點實驗室，長春130022)

?

玻璃與玻璃化轉(zhuǎn)變

李艷偉孫昭艷*安立佳

(中國科學(xué)院長春應(yīng)用化學(xué)研究所，高分子物理與化學(xué)國家重點實驗室，長春130022)

摘要：玻璃與玻璃態(tài)的應(yīng)用極其廣泛。玻璃化轉(zhuǎn)變是一種典型的非晶液-固轉(zhuǎn)變，當(dāng)轉(zhuǎn)變發(fā)生時體系的結(jié)構(gòu)并沒有明顯變化，因而我們無法將其歸類于已有的任何相變類型。作為凝聚態(tài)物理和軟物質(zhì)領(lǐng)域的核心問題，玻璃化轉(zhuǎn)變的研究已有近70年的歷史。然而，時至今日，人們還是無法回答玻璃態(tài)的本質(zhì)是什么這一基本問題。本文簡述了玻璃態(tài)的性質(zhì)以及伴隨玻璃化轉(zhuǎn)變發(fā)生的一些基本物理現(xiàn)象，并總結(jié)了半個世紀(jì)以來一些與玻璃化轉(zhuǎn)變相關(guān)的理論，以期加深讀者對玻璃及玻璃化轉(zhuǎn)變的認(rèn)識。

關(guān)鍵詞：玻璃；玻璃化轉(zhuǎn)變；過冷液體

www.dxhx.pku.edu.cn

玻璃是人類最早使用的材料之一[1]。大約公元前3500年，古巴比倫人已經(jīng)可以在簡單的作坊里制作玻璃。戰(zhàn)國時期，羅馬人制造的玻璃經(jīng)由印度傳入中國。最初人們誤認(rèn)為玻璃是極其珍貴的天然寶藏，后由東晉著名煉丹家葛洪指出玻璃為人造品，他在其《抱樸子內(nèi)篇——論仙》中說：“外國作水精碗，實是合五種灰以作之。”水精指的就是我們現(xiàn)在說的玻璃。根據(jù)葛洪的記載，當(dāng)時在今兩廣一帶已能制作玻璃。中國古代對玻璃的稱呼很多，有記載的就有十多種，如琉璃、瑯玕、水玉、琉琳、繆琳、陸離、罐玉、玻黎、藥玉、水精等。當(dāng)代考古學(xué)中又把不透明玻璃稱之為“料器”，把半透明的稱之為“琉璃”，把透明度高的稱之為“玻璃”。我國著名玻璃研究專家中國科學(xué)院院士干福熹教授建議“統(tǒng)一稱為玻璃”[2]。

幾乎所有的物質(zhì)在一定條件下都可以形成玻璃[3－6]。液體在快速加壓或快速降溫時會避開結(jié)晶并進(jìn)入過冷態(tài)，當(dāng)過冷度繼續(xù)增加時，體系的黏度不斷升高，直至體系變得過于黏稠而失去流動性，此時材料的性質(zhì)已與固體無異。這種處于非晶態(tài)的固體被稱為玻璃。玻璃以其優(yōu)越的力學(xué)及光學(xué)性能，在材料領(lǐng)域一直占據(jù)著重要的地位。另外，處于玻璃態(tài)的生命系統(tǒng)，其細(xì)胞和組織的結(jié)構(gòu)并未遭到破壞，因此，在未來醫(yī)學(xué)(如復(fù)活冰凍人技術(shù))及微生物學(xué)(如研究古老的南極孢子)中都將繼續(xù)發(fā)揮重要作用。

玻璃態(tài)作為對立于晶態(tài)的另一種固體形式，既具有類似于晶體的力學(xué)性質(zhì)又具有類似于液體的無序結(jié)構(gòu)[7－9]。這種奇特的混合性質(zhì)使得科學(xué)家們著迷于對其物理本質(zhì)的研究。關(guān)于玻璃態(tài)本質(zhì)的研究一直是凝聚態(tài)物理及軟物質(zhì)領(lǐng)域的重要內(nèi)容，也是至今懸而未決的難題。早在2005年，國際著名雜志Science在其創(chuàng)刊125周年時，將“玻璃態(tài)的本質(zhì)是什么”這一問題列為21世紀(jì)亟待解決的重大科學(xué)問題之一[10]。然而，研究發(fā)現(xiàn)，我們既不能用已有的液體理論來認(rèn)識玻璃，也不能用已有的晶體理論來理解玻璃，因此，發(fā)展完善的理論解釋玻璃的物理本質(zhì)一直是科學(xué)家們的不懈追求[5]。自上個世紀(jì)40年代以來，關(guān)于玻璃化的理論可謂層出不窮，美國科學(xué)院院士哈佛大學(xué)Weitz教授曾開玩笑說，玻璃化轉(zhuǎn)變的理論比提出它們的科學(xué)家還多。遺憾的是，迄今為止沒有一個理論能解釋玻璃化轉(zhuǎn)變過程中的所有現(xiàn)象，已有的理論也只是在某些特定的過冷區(qū)間和特定的體系中才與實驗或模擬結(jié)果吻合。

由于玻璃化轉(zhuǎn)變本身的復(fù)雜性加上原子、分子體系尺度太小，運(yùn)動過快，從而難以研究其微觀機(jī)理。膠體是一個很好的替代原子和小分子體系的模型體系，它是由尺寸介于幾納米到幾微米的顆粒物分散在一定溶劑中形成的混合物。膠體的尺度適中，既不會因尺度過大而沉淀，也不會因尺度太小而難于觀察。膠體作為原子、分子的模型體系已經(jīng)被廣泛認(rèn)可。從物質(zhì)的分類上來說，膠體屬于軟物質(zhì)，它具有軟物質(zhì)的典型特點，即體系是由熵效應(yīng)和熱運(yùn)動驅(qū)動的體系，量子效應(yīng)對體系影響不大。實驗上常用的膠體顆粒通常是由高分子制備的，這些膠體顆粒的性質(zhì)可以通過一些手段進(jìn)行有效的調(diào)控，如改變粒子直徑、粒子間相互作用等。除了膠體體系還有很多體系可以用來研究玻璃化轉(zhuǎn)變行為，如直徑在1微米以上的顆粒物質(zhì)體系[11]、高分子體系[12,13]、非晶合金體系[9]等。

研究玻璃化轉(zhuǎn)變的課題組通常被劃分在軟凝聚態(tài)物理或者高分子物理方向。如哈佛大學(xué)Weitz教授領(lǐng)導(dǎo)的軟凝聚態(tài)物質(zhì)研究組報道了許多重要的關(guān)于玻璃體系的實驗研究成果，杜克大學(xué)的物理系及芝加哥大學(xué)的James Franck研究所長期以來從事相關(guān)研究工作。另外，由于de Gennes的影響，法國在研究玻璃化轉(zhuǎn)變方面有很強(qiáng)的研究力量，如蒙彼利埃(Université de Montpellier)大學(xué)、巴黎高師(Ecole Normale Superieure)等研究機(jī)構(gòu)在玻璃化轉(zhuǎn)變理論和模擬方面的研究成果非常引人注目。

關(guān)于玻璃化轉(zhuǎn)變的研究通常是從熱力學(xué)和動力學(xué)兩個角度開展的。人們期望借此明晰玻璃化轉(zhuǎn)變本質(zhì)上是熱力學(xué)控制的還是動力學(xué)控制的。目前，人們已經(jīng)有了大量的關(guān)于玻璃體系結(jié)構(gòu)與動力學(xué)的知識積累，并且已然確定了二者之間存在一定程度的關(guān)聯(lián)[14,15]，但是，動力學(xué)上的急劇變慢與結(jié)構(gòu)上哪些性質(zhì)的變化相對應(yīng)還未確定。雖然玻璃化轉(zhuǎn)變問題的復(fù)雜性超乎想象，但是，我們依然堅信，目前遇到的瓶頸只是黎明前的黑暗，玻璃化轉(zhuǎn)變問題的解決必將推動相關(guān)學(xué)科的發(fā)展。

以下我們簡要總結(jié)了玻璃體系的一些基本物理現(xiàn)象和玻璃化轉(zhuǎn)變的一些基本理論，重點討論玻璃體系結(jié)構(gòu)與動力學(xué)的關(guān)系。

1　玻璃體系的基本物理現(xiàn)象

1.1強(qiáng)性和脆性

不同材料黏度的對數(shù)對玻璃化轉(zhuǎn)變溫度與溫度的比值(Tg/T)作圖，被稱為Angell圖[16]。如圖1所示，不同的體系在Angell圖中表現(xiàn)出不同的特征。有些體系的黏度符合Arrhenius關(guān)系，即：

圖1　黏度η與玻璃化轉(zhuǎn)變溫度Tg與溫度T的比值關(guān)系示意圖，也稱為Angell圖[16]

這些體系在Angell圖中表現(xiàn)為直線的特征，這一類玻璃體系被稱為強(qiáng)性(strong)玻璃，我們平時見到的窗戶玻璃屬于這一類；另外一些體系較大地偏離了Arrhenius關(guān)系，在Angell圖中呈現(xiàn)出較大的曲率，這一類玻璃體系被稱為脆性(fragile)玻璃。顯然，我們可以用Angell圖中各曲線在Tg(玻璃化轉(zhuǎn)變點)處的斜率來表示體系脆性的強(qiáng)弱，此斜率值越大，體系越脆。

需要注意的是，以上提及的Tg是人為定義的(η= 1013時的溫度)，并非是某個熱力學(xué)上的轉(zhuǎn)變點，因此，強(qiáng)性和脆性的區(qū)分也帶有一定的人為性。我們可以想象，如果人類生命的時間尺度是1030而不是102的話，我們做一個更長時間的實驗，給Tg一個更深過冷度的定義的話，那么，目前對體系強(qiáng)性和脆性的區(qū)分可能就需要有較大的修正。總之，強(qiáng)性和脆性的區(qū)分僅僅是為了方便研究做的定量上的區(qū)別[2]。

從物理的角度講，如果認(rèn)為強(qiáng)性玻璃是一類黏性非常大的液體的話也未嘗不可，然而，這種說法對脆性玻璃來說并不合適，因為脆性玻璃體系黏度的劇變僅發(fā)生在接近Tg的某個溫度區(qū)間附近，這個劇變背后隱含哪些物理機(jī)制顯然需要更深入的研究。

1.2動力學(xué)變慢及兩步松弛

圖2　三維Kob-Andersen模型體系的自中間散射函數(shù)Fs(qp,t)

體系從初始的平衡狀態(tài)過渡到另一個平衡狀態(tài)的過程稱為松弛過程，恢復(fù)到新平衡態(tài)所需要的時間稱為松弛時間τ。在光散射和中子散射實驗中，通常用中間散射函數(shù)(intermediate scattering function)來表征體系的松弛行為。在計算機(jī)模擬中，自關(guān)聯(lián)部分的中間散射函數(shù)(圖2)也經(jīng)常用來表征玻璃體系的松弛行為。其定義如下：

在高溫下，體系動力學(xué)的衰減呈指數(shù)形式；低溫下，出現(xiàn)了衰減的平臺區(qū)，這時的松弛行為表現(xiàn)為兩步松弛模式。隨著溫度的降低，體系需要更長的時間才能使中間散射函數(shù)松弛到0，即體系的動力學(xué)變得越來越慢。通常定義體系的松弛時間τ為中間散射函數(shù)衰減到1/e的時間。

過冷液體的松弛時間與溫度之間滿足Vogel-Fulcher-Tamman(VFT)方程[17－19]：

由此方式獲得的T0被稱為VFT玻璃化轉(zhuǎn)變溫度，其數(shù)值上約等于理想玻璃化轉(zhuǎn)變溫度或Kauzmann溫度TK。

1.3動力學(xué)異質(zhì)性

動力學(xué)異質(zhì)性是過去十幾年發(fā)現(xiàn)的玻璃體系動力學(xué)上的另一個重要特征。圖3給出了動力學(xué)異質(zhì)性的典型實例，粒子在空間中的運(yùn)動表現(xiàn)為有些區(qū)域運(yùn)動得慢，而有些區(qū)域運(yùn)動得快。動力學(xué)異質(zhì)性被認(rèn)為是造成過冷液體不同于簡單流體的諸多行為(如斯托克斯-愛因斯坦(Stokes-Einstein)關(guān)系)的原因[20－25]。另外，探索動力學(xué)異質(zhì)性的結(jié)構(gòu)起源也是近年來研究的焦點問題之一[26－30]。

圖3　二元Lennard-Jones(LJ)過冷流體中單粒子位移的空間軌跡圖[31]

2　玻璃化轉(zhuǎn)變的基本理論

自20世紀(jì)40年代以來，玻璃化轉(zhuǎn)變的理論不斷涌現(xiàn)，如：自由體積理論(free volume theory)[32,33]、Adam-Gibbs理論[34]、Goldstein勢能圖景理論(Goldstein′s energy landscape theory)[35]、模耦合理論(mode coupling theory，MCT)[36]、動力學(xué)約束模型(kinetically constrained model，KCM)[37]、自旋玻璃理論(spin glass theory)[38]、無規(guī)一階轉(zhuǎn)變理論(random first order transition theory，RFOT)[39,40]、受挫有限疇理論(frustration-limited domain theory)[41]、介觀平均場理論(mesoscopic mean-field theory)[42]、副本計算方法(replica calculation)[43]等。下面介紹其中幾種常見的理論。

2.1自由體積理論

自由體積理論是早期針對玻璃化轉(zhuǎn)變現(xiàn)象而提出來的諸多維象理論之一[44]。它的理論假設(shè)為：在過冷液體中存在一定量的局域自由體積，這些自由體積為鄰近粒子的運(yùn)動提供了空間。該理論最終得出了體系黏度的VFT變化規(guī)律。后來Cohen等[32,33]擴(kuò)展了自由體積理論，他們認(rèn)為液體體系包含類固區(qū)和類液區(qū)，類液區(qū)內(nèi)體積可以自由交換，并以此保證液體的流動。低溫下，類固區(qū)的增多而相應(yīng)類液區(qū)的減少使得體系內(nèi)進(jìn)行體積交換變得困難，這導(dǎo)致了體系動力學(xué)上的變慢。在玻璃化轉(zhuǎn)變點時，體系的類固區(qū)跨越了整個體系，即體系發(fā)生了逾滲。這時類液區(qū)的通道被堵塞，體系不能發(fā)生流動。

自由體積理論為我們描述了一幅清晰而具體的物理圖像，并解釋了一部分的實驗現(xiàn)象。然而，其不足之處在于，自由體積作為理論的唯一參量，在實驗中幾乎無法測量。另外，自由體積并非是一個固定值，它與冷卻速率有關(guān)，即便在Tg以下，體系并沒有完全被凍結(jié)，體系的自由體積仍然會改變。

2.2模耦合理論

在玻璃化轉(zhuǎn)變的諸多理論中，模耦合理論(MCT)是少數(shù)可以定量計算的理論之一[36]。計算中只需知道體系的密度漲落(如靜態(tài)結(jié)構(gòu)因子)，通過Zwanzig-Mori投影算符技術(shù)，最終得出體系的動態(tài)關(guān)聯(lián)函數(shù)。我們知道，過冷態(tài)與其相應(yīng)非過冷態(tài)的靜態(tài)結(jié)構(gòu)因子差別并不顯著，然而，為什么能通過沒有顯著差別的量推導(dǎo)出具有顯著差別的中間散射函數(shù)呢？這是因為MCT方程組中引入了非線性反饋機(jī)制(為人熟知的蝴蝶效應(yīng)便是此機(jī)制)，靜態(tài)結(jié)構(gòu)因子上小的區(qū)別便會在中間散射函數(shù)上有大的反饋。MCT理論成功預(yù)測了動力學(xué)的兩步松弛模式以及中間散射函數(shù)上出現(xiàn)平臺區(qū)的現(xiàn)象。另外，MCT預(yù)測體系松弛時間隨溫度的變化規(guī)律滿足MCT公式：

其中τ0，γ為常數(shù)，TC為模耦合玻璃化轉(zhuǎn)變溫度。顯然，體系的松弛時間在T = TC處發(fā)散。然而，以往報道的實驗和模擬結(jié)果卻并未觀察到τ在TC處發(fā)散，并且TC總是明顯大于Tg。另外，近期的研究結(jié)果表明：MCT理論在更高的維度下并不成立[45,46]，也無法解釋吸引相互作用對體系的非微擾作用[47,48]。

2.3熱力學(xué)理論

液體的熱容(heat capacity)比晶體大，因此，降溫過程中，液體的熵下降幅度比晶體快。如圖4(a)所示，當(dāng)液體在低于熔點并避開結(jié)晶時便進(jìn)入過冷液體區(qū)，若在過冷液體區(qū)保持其原有的熵下降速率的話，在某個溫度TK(又被稱為理想玻璃化轉(zhuǎn)變溫度)，過冷液體的熵將等于晶體的熵。如果定義液體的熵與相同溫度下晶體的熵的差值為過剩熵(excess entropy)，那么此時的過剩熵為0。當(dāng)繼續(xù)降溫時，過冷液體的熵繼續(xù)降低，并在非零溫度時就達(dá)到負(fù)值(由于完美晶體在T = 0時的熵為0)，顯然這違反了熱力學(xué)第三定律。這就是所謂的“Kauzmann佯謬”(Kauzmann paradox)[49]。事實上，在高于TK的某個溫度點Tg時，體系就已發(fā)生了玻璃化轉(zhuǎn)變。

玻璃化轉(zhuǎn)變中的許多熱力學(xué)理論都是基于Kauzmann的思想。Adam-Gibbs理論就是其中一個[34]，該理論認(rèn)為玻璃體系的松弛行為源于體系中的一系列子分區(qū)松弛到各自的能量極小點。這些相互獨立的子分區(qū)被稱為協(xié)同重整區(qū)(cooperative rearranging regions，簡稱CRR)。由于溫度的降低導(dǎo)致了CRR可以歷經(jīng)的能量極小點的減小，這意味著體系的構(gòu)型熵(configurational entropy)的減小。Adam-Gibbs理論運(yùn)用統(tǒng)計力學(xué)將各個CRR作為相互獨立的多體體系處理，最終導(dǎo)出了聯(lián)系體系松弛動力學(xué)和熱力學(xué)的公式：

其中τ0和C為常數(shù)，SC為體系的構(gòu)型熵。模擬中通過直接計算體系的構(gòu)型熵(體系的總熵減去體系的振動熵)與松弛時間證明Adam-Gibbs理論公式在高于二維的體系中符合良好，如圖4(b)所示。

以上Adam-Gibbs理論并沒有給出關(guān)于CRR尺寸大小及CRR中可以遍歷的構(gòu)型數(shù)的預(yù)測。上世紀(jì)80年代末提出的RFOT理論或者被稱為“馬賽克”理論(mosaic theory)[39]給出了一定的解釋。RFOT理論認(rèn)為玻璃體系存在很多的亞穩(wěn)態(tài)。不同區(qū)域的亞穩(wěn)態(tài)(被稱為mosaic態(tài))之間必然有一個界面，此不相容界面的界面能使得體系之間不同亞穩(wěn)態(tài)的轉(zhuǎn)換變得困難。另外，體系的熵效應(yīng)(或者熱漲落)又促使體系能夠發(fā)生這樣的轉(zhuǎn)變。當(dāng)體系熵效應(yīng)與界面效應(yīng)達(dá)到平衡時，可定義出體系亞穩(wěn)態(tài)區(qū)域的最大長度尺度ξ(也可理解為Adam-Gibbs理論中CRR區(qū)最大值)。

圖4　(a)液體、過冷液體、玻璃和晶體的熵與溫度之間的關(guān)系[54]；(b)三維KA模型的Adam-Gibbs關(guān)系圖[55]

其中d為體系的維度，θ為與維度有關(guān)的量，在三維體系建議值為1.5(初始計算值為2.0)。Y和SC分別為體系的界面能和構(gòu)型熵，它們是與溫度有關(guān)的量。RFOT理論認(rèn)為，玻璃體系包含許多尺寸為R (R <ξ)的亞穩(wěn)區(qū)域，這些區(qū)域就像給體系打了一個個的馬賽克一樣，每個區(qū)域可以遍歷的狀態(tài)數(shù)在exp(RdSC)的量級。顯然，當(dāng)R≥ξ時，亞穩(wěn)態(tài)的概念將失去意義，體系進(jìn)入熱力學(xué)上的理想玻璃態(tài)，因此，RFOT理論預(yù)測了玻璃體系低于TC的行為，溫度范圍在TK< T < TC的區(qū)也被稱為RFOT主導(dǎo)區(qū)。近年來發(fā)展的點對集關(guān)聯(lián)長度即是RFOT理論的一個重要應(yīng)用[28－30,50－53]。

3　玻璃體系的結(jié)構(gòu)與動力學(xué)

玻璃化轉(zhuǎn)變伴隨著體系動力學(xué)上的顯著變化，但是對關(guān)聯(lián)函數(shù)卻不發(fā)生明顯變化，因而使得很多人開始質(zhì)疑玻璃化轉(zhuǎn)變過程中結(jié)構(gòu)與動力學(xué)是否一一對應(yīng)。目前，人們還不知道玻璃化轉(zhuǎn)變過程中到底是哪種結(jié)構(gòu)參量在起作用。通過同構(gòu)型系綜平均(isoconfigurational ensemble)可以確定的是，玻璃體系結(jié)構(gòu)與動力學(xué)必然存在某種聯(lián)系[14,15,56]。近年來，人們發(fā)展了許多方法來描述玻璃體系的結(jié)構(gòu)特征，如鍵取向有序[57－60]、點對集關(guān)聯(lián)長度[28－30,50－53]、局部有序結(jié)構(gòu)(如二十面體團(tuán)簇結(jié)構(gòu)[61,62]、面心立方團(tuán)簇結(jié)構(gòu)[63]和正四面體團(tuán)簇結(jié)構(gòu)[64,65])等。目前，這些方法在解釋玻璃體系結(jié)構(gòu)與動力學(xué)這一基礎(chǔ)問題上至少有以下兩方面的問題需要解決：一方面需要解決該結(jié)構(gòu)特征是否對眾多玻璃體系都普適的問題。基于有序結(jié)構(gòu)或者有序度出發(fā)來描述玻璃體系的結(jié)構(gòu)，往往對弱受挫的多分散硬球體系(這些體系中，粒子的相對位置漲落并不是很大，往往能形成較好的局部有序結(jié)構(gòu))能給出很好的結(jié)果[66,67]，但是在二分散的體系中很難給出合理的結(jié)果。另一方面是這種結(jié)構(gòu)特征是否與動力學(xué)有直接的對應(yīng)關(guān)系，也就是說，體系的靜態(tài)關(guān)聯(lián)長度與動態(tài)關(guān)聯(lián)長度和松弛時間尺度是否有直接的對應(yīng)關(guān)系。但是已有的研究工作表明，隨著溫度的降低，靜態(tài)關(guān)聯(lián)長度比動態(tài)關(guān)聯(lián)長度增加的幅度小得多。接下來我們再介紹幾種描述玻璃體系結(jié)構(gòu)與動力學(xué)關(guān)系的方法[24,28,68－71]。

3.1鍵取向有序

Tanaka等[57－60]認(rèn)為鍵取向有序可以作為玻璃體系動力學(xué)異質(zhì)性的結(jié)構(gòu)起源，他們把體系中高鍵取向有序度的區(qū)域稱為中度結(jié)晶有序(medium-range crystalline order，簡稱MRCO)區(qū)。處于MRCO區(qū)的粒子運(yùn)動較慢，并主導(dǎo)了體系慢松弛過程。通過對比多分散性為11.98%(此多分散度既可以避免結(jié)晶又允許局部有序結(jié)構(gòu)在低溫下能較穩(wěn)定的存在)的Lennard-Jones(LJ)體系和Weeks-Chandler-Andersen (WCA)體系，可以看出LJ和WCA體系的動力學(xué)并不能與MRCO定義的靜態(tài)關(guān)聯(lián)長度完美對應(yīng)[70]。如圖5(a)所示，LJ和WCA體系在相同溫度下的動態(tài)關(guān)聯(lián)長度ξ4存在明顯差別(尤其是低溫下)，這與兩個體系在相同溫度下的動力學(xué)存在明顯差異的結(jié)果相符。但是，兩個體系的靜態(tài)關(guān)聯(lián)長度ξ6差別不大(圖5(b))。圖5(c)和圖5(d)分別給出了兩個體系的ξ4和ξ6隨松弛時間τα的變化規(guī)律。顯然，LJ和WCA體系ξ4隨τα的變化有相同的標(biāo)度關(guān)系，而兩個體系ξ6隨τα的變化有明顯不同的標(biāo)度關(guān)系。這說明由MRCO定義的靜態(tài)關(guān)聯(lián)長度并不能體現(xiàn)LJ和WCA體系動力學(xué)的差別。需要注意的這里所描述的體系是二維下多分散的體系，受挫度較小，這說明即便在弱受挫的玻璃體系里，以有序度表征的靜態(tài)關(guān)聯(lián)長度ξ6也不能與體系的動力學(xué)性質(zhì)相對應(yīng)。而在三維的強(qiáng)受挫玻璃體系里，ξ6更難給出讓人信服的結(jié)果。

圖5　動態(tài)關(guān)聯(lián)長度ξ4(a)與通過有序度定義的靜態(tài)關(guān)聯(lián)長度ξ6(b)隨溫度的變化以及動態(tài)關(guān)聯(lián)長度ξ4(c)與靜態(tài)關(guān)聯(lián)長度ξ6(d)隨松弛時間的變化[70]

近期報道的其他結(jié)果也表明，通過局部有序來理解玻璃體系的結(jié)構(gòu)與動力學(xué)有很強(qiáng)的體系依賴性，在粒子位置漲落不是太大的弱受挫玻璃體系里有可能給出較好的結(jié)果，而在強(qiáng)受挫的玻璃體系里幾乎沒有什么對應(yīng)性[67]。

3.2點對集關(guān)聯(lián)

基于RFOT理論[39,40]，Biroli等[30]提出了通過點對集重疊參量獲得體系靜態(tài)長度的方法。后來，Kob等[50]發(fā)展了此方法，實現(xiàn)了一次模擬即可獲得點對集關(guān)聯(lián)長度的方法。具體做法為：首先任意選取一個平衡構(gòu)型，然后固定某個區(qū)域(可以是不同的幾何形狀)之外的粒子(圖6)，未被固定的粒子在這種無定形邊界(amorphous boundary)條件下繼續(xù)演化。這種受限形式也被稱為固定場(pinning field)。倘若點對集關(guān)聯(lián)能夠有效表征玻璃體系的結(jié)構(gòu)性質(zhì)，那么，即使采取不同幾何形狀的固定場，以此獲得的點對集靜態(tài)關(guān)聯(lián)長度與動力學(xué)性質(zhì)之間的關(guān)系應(yīng)不隨幾何形狀的改變而改變[28]。

圖7給出了由不同幾何形狀的固定場獲得的約化靜態(tài)關(guān)聯(lián)長度與松弛時間之間的關(guān)系。結(jié)果暗示這種線性關(guān)系只存在于轉(zhuǎn)變溫度T0附近，而在深過冷區(qū)，靜態(tài)關(guān)聯(lián)長度與松弛時間的關(guān)系顯然與幾何形狀密切相關(guān)。這可能預(yù)示著在深過冷區(qū)，玻璃體系存在多種松弛模式。同時，這一結(jié)果也從另一方面說明用點對集靜態(tài)關(guān)聯(lián)長度來描述玻璃體系的動力學(xué)行為，尤其是在深過冷區(qū)仍然不能令人滿意。

圖6　不同固定場幾何形狀的示意圖[28]

圖7　不同幾何固定場條件下松弛時間與約化點對集關(guān)聯(lián)長度的關(guān)系[28]

3.3局域密堆積粒子的空間關(guān)聯(lián)

前述幾種描述玻璃化轉(zhuǎn)變過程中的結(jié)構(gòu)變化均未能給出與動力學(xué)呈現(xiàn)一一對應(yīng)的理想結(jié)果。事實上，當(dāng)我們從直覺上思考玻璃化轉(zhuǎn)變時，似乎玻璃化轉(zhuǎn)變過程中的一個重要信息就是液體變得更加稠密。而密度(或體積的倒數(shù))作為表征粒子堆砌緊密程度的重要表征量，也許就是解開結(jié)構(gòu)與動力學(xué)對應(yīng)之謎的金鑰匙。在數(shù)學(xué)上，有一種可以準(zhǔn)確計算粒子局域體積的方法，即Voronoi方法。在一個過冷流體中，將所有流體粒子的質(zhì)心點標(biāo)出，并作所有粒子與其近鄰粒子質(zhì)心的垂直平分線(二維體系)或垂直平分面(三維體系)，這些垂直平分線或垂直平分面所圍成的多邊形或多面體即稱為Voronoi多邊形或Voronoi多面體[72]。將三維體系進(jìn)行Voronoi切分后，每個粒子都對應(yīng)一個Voronoi多面體，這些多面體體積(用ν表示)的大小反映了粒子所處的局域環(huán)境的疏密情況，越小，粒子的局域密度(用ρl表示，ρl= 1/ν)越大。然而，目前尚未發(fā)現(xiàn)通過局域密堆積粒子的空間關(guān)聯(lián)找到與動力學(xué)之間存在一一對應(yīng)關(guān)系的直接證據(jù)。但仍有許多科學(xué)家堅信，局域密度的空間關(guān)聯(lián)是最能夠代表玻璃結(jié)構(gòu)與液體結(jié)構(gòu)本質(zhì)區(qū)別的重要物理量。因此，發(fā)展合適的結(jié)構(gòu)表征手段并尋找結(jié)構(gòu)與動力學(xué)之間的關(guān)系仍然屬于未解的科學(xué)難題。

4　結(jié)束語

玻璃，從極度落后的遠(yuǎn)古時代到科技發(fā)達(dá)的今天，已有超過四千年的歷史。但是人們對玻璃及玻璃化轉(zhuǎn)變的物理本質(zhì)這一基本問題仍未得到確切的答案。但這并不意味著人們對玻璃仍舊一無所知，科學(xué)家們已經(jīng)能從不同角度闡述玻璃或者玻璃化轉(zhuǎn)變過程中的一些現(xiàn)象。正如偉大的科學(xué)家愛因斯坦說過的那樣，科學(xué)絕不是也永遠(yuǎn)不會是一本寫完了的書。每一項重大成就都會帶來新的問題；任何一個發(fā)展隨著時間的推移都會出現(xiàn)新的嚴(yán)重的困難。正是隨著人們對玻璃化轉(zhuǎn)變認(rèn)識的不斷深入，才會發(fā)掘出更多的問題。每年美國物理學(xué)會的三月份年會都會有一個與玻璃及玻璃化轉(zhuǎn)變相關(guān)的分會場，專門討論研究人員在玻璃領(lǐng)域取得的最新成果。在物理學(xué)領(lǐng)域，許多人堅信如果能找到物理學(xué)的精髓，那么只需一頁紙就夠了。但是，自然界本身就是紛繁復(fù)雜的，也許本身就沒有一個統(tǒng)一的規(guī)律可循。即使對于玻璃這個在龐大自然界中顯得如此渺小的問題，其實也有其多面性和復(fù)雜性，因而，對于玻璃及玻璃化轉(zhuǎn)變的認(rèn)識更應(yīng)該是由許多知識一點一滴積累起來的。契科夫說過，想喝水時，仿佛能喝下整個海洋似的，這是信仰；等到真的喝起來，一共也只能喝兩杯罷了，這就是科學(xué)。但是，把每個人喝下的水都匯集起來，那里就會是一片汪洋大海。相信經(jīng)過科學(xué)家們持續(xù)不斷的努力，玻璃及玻璃化轉(zhuǎn)變的物理本質(zhì)之謎最終一定會揭開！

參考文獻(xiàn)

[1]李青會,干福熹,顧冬紅.自然科學(xué)史研究, 2007, 26 (2), 234.

[2]干福熹.自然雜志, 2006, 28 (4), 187.

[3] Ediger, M. D.; Harrowell, P. J. Chem. Phys. 2012, 137, 080901.

[4] Berthier, L.; Biroli, G. Rev. Mod. Phys. 2011, 83, 587.

[5] Lubchenko, V.; Wolynes, P. G.Annu. Rev. Phys. Chem. 2007, 58, 235.

[6] Cavagna,A. Physics Reports 2009, 476, 51.

[7] Binder, K.; Kob, W. Glassy Materials and Disordered Solids; World Scientific: Singapore, 2005.

[8] Berthier, L.; Biroli, G.; Bouchaud, J.-P.; Cipelletti, L.; van Saarloos, W. Dynamical Heterogeneities in Glasses, Colloids, and Granular Media; Oxford University Press: Oxford, 2011.

[9]汪衛(wèi)華.物理學(xué)進(jìn)展, 2013, 33 (5), 177.

[10] Couzin, J. Science 2005, 309, 83.

[11] Jaeger, H. M.; Nagel, S. R.; Behringer, R. P. Rev. Mod. Phys. 1996, 68, 1259.

[12] Ngai, K. L.; Roland, C. M. Macromolecules 1993, 26, 6824.

[13] Siretanu, I.; Saadaoui, H.; Chapel, J.-P.; Drummond, C. Macromolecules 2015, 48, 2787.

[14] Widmer-Cooper,A.; Harrowell, P. Phys. Rev. Lett. 2006, 96, 185701.

[15] Widmer-Cooper,A.; Perry, H.; Harrowell, P.; Reichman, D. R. Nat. Phys. 2008, 4, 711.

[16] Angell, C.A. Science 1995, 267, 1924.

[17] Vogel, H. Phys. Z. 1921, 22, 645.

[18] Fulcher, G. S. J. Am. Ceram. Soc. 1925, 8, 339.

[19] Tammann, G.; Hesse, W. Z.Anorg. Allg. Chem. 1926, 156, 245.

[20] Sengupta, S.; Karmakar, S.; Dasgupta, C.; Sastry, S. J. Chem. Phys. 2013, 138, 12A548.

[21] Griffin, P. J.; Sangoro, J. R.; Wang, Y.; Holt,A. P.; Novikov, V. N.; Sokolov,A. P.; Wojnarowska, Z.; Paluch, M.; Kremer, F. Soft Matter 2013, 9, 10373.

[22] Xu, L.; Mallamace, F.; Yan, Z.; Starr, F. W.; Buldyrev, S. V.; Eugene Stanley, H. Nat. Phys. 2009, 5, 565.

[23] Becker, S. R.; Poole, P. H.; Starr, F. W. Phys. Rev. Lett. 2006, 97, 055901.

[24] Li, Y. W.; Zhu, Y. L.; Sun, Z. Y. J. Chem. Phys. 2015, 142, 124507.

[25] Charbonneau, P.; Jin, Y.; Parisi, G.; Zamponi, F. Proc. Natl. Acad. Sci. U. S. A. 2014, 111, 15025.

[26] Smarajit, K.; Chandan, D.; Srikanth, S. Rep. Prog. Phys. 2016, 79, 016601.

[27] Russo, J.; Tanaka, H. Proc. Natl. Acad. Sci. U. S. A. 2015, 112, 6920.

[28] Li, Y. W.; Xu, W. S.; Sun, Z. Y. J. Chem. Phys. 2014, 140, 124502.

[29] Hocky, G. M.; Markland, T. E.; Reichman, D. R. Phys. Rev. Lett. 2012, 108, 225506.

[30] Biroli, G.; Bouchaud, J. P.; Cavagna,A.; Grigera, T. S.; Verrocchio, P. Nat. Phys. 2008, 4, 771.

[31] Berthier, L. Physics 2011, 4, 42.

[32] Cohen, M. H.; Turnbull, D. J. Chem. Phys. 1959, 31, 1164.

[33] Turnbull, D.; Cohen, M. H. J. Chem. Phys. 1961, 34, 120.

[34] Adam, G.; Gibbs, J. H. J. Chem. Phys. 1965, 43, 139.

[35] Goldstein, M. J. Chem. Phys. 1969, 51, 3728.

[36] G?tze, W. Complex Dynamics of Glass-Forming Liquids: A Mode-Coupling Theory; Oxford University Press: Oxford, 2008.

[37] Fredrickson, G. H.;Andersen, H. C. Phys. Rev. Lett. 1984, 53, 1244.

[38] Mézard, M.; Parisi, G.; Virasoro, M. Spin Glass Theory and Beyond; World Scientific: Singapore, 1988.

[39] Kirkpatrick, T. R.; Thirumalai, D.; Wolynes, P. G. Phys. Rev. A 1989, 40, 1045.

[40] Bouchaud, J. P.; Biroli, G. J. Chem. Phys. 2004, 121, 7347.

[41] Kivelson,A.; Kivelson, S.A.; Zhao, X.; Nussinov, Z.; Tarjus, G. Physica A 1995, 218, 27.

[42] Chamberlin, R. V. Phys. Rev. Lett. 1999, 82, 2520.

[43] Parisi, G.; Zamponi, F. J. Chem. Phys. 2005, 123, 144501.

[44] Eyring, H.; Wynne-Jones, W. F. K. J. Chem. Phys. 1936, 4, 293.

[45] Ikeda,A.; Miyazaki, K. Phys. Rev. Lett. 2010, 104, 25570.

[46] Charbonneau, P.; Ikeda,A.; Parisi, G.; Zamponi, F. Proc. Natl. Acad. Sci. U. S. A. 2012, 109, 13939.

[47] Berthier, L.; Tarjus, G. Phys. Rev. Lett. 2009, 103, 170601.

[48] Berthier, L.; Tarjus, G. J. Chem. Phys. 2011, 134, 214503.

[49] Kauzmann, W. Chem. Rev. 1948, 43, 219.

[50] Kob, W.; Roldán-Vargas, S.; Berthier, L. Nat. Phys. 2012, 8, 164.

[51] Yaida, S.; Berthier, L.; Charbonneau, P.; Tarjus, G. arxiv 2015, 1511.03573.

[52] Berthier, L.; Kob, W. Phys. Rev. E 2012, 85, 011102.

[53] Charbonneau, P.; Tarjus, G. Phys. Rev. E 2013, 87, 042305.

[54] Chang, S. S.; Bestul,A. B. J. Chem. Phys. 1972, 56, 503.

[55] Sengupta, S.; Karmakar, S.; Dasgupta, C.; Sastry, S. Phys. Rev. Lett. 2012, 109, 095705.

[56] Widmer-Cooper,A.; Harrowell, P.; Fynewever, H. Phys. Rev. Lett. 2004, 93, 135701.

[57] Shintani, H.; Tanaka, H. Nat. Phys. 2006, 2, 200.

[58] Kawasaki, T.;Araki, T.; Tanaka, H. Phys. Rev. Lett. 2007, 99, 215701.

[59] Tanaka, H.; Kawasaki, T.; Shintani, H.; Watanabe, K. Nat. mater. 2010, 9, 324.

[60] Watanabe, K.; Kawasaki, T.; Tanaka, H. Nat. Mater. 2011, 10, 512.

[61] Tarjus, G.; Kivelson, S.A.; Nussinov, Z.; Viot, P. J. Phys.-Condens. Mat. 2005, 17, R1143.

[62] Spaepen, F. Nature 2000, 408, 781.

[63] Leocmach, M.; Tanaka, H. Nat. Commun. 2012, 3, 974.

[64] Anikeenko,A. V.; Medvedev, N. N. Phys. Rev. Lett. 2007, 98, 235504.

[65] Xia, C.; Li, J.; Cao, Y.; Kou, B.; Xiao, X.; Fezzaa, K.; Xiao, T.; Wang, Y. Nat. Commun. 2015, 6, 8409.

[66] Bruning, R.; St-Onge, D.A.; Patterson, S.; Kob, W. J. Phys.-Condens Mat. 2009, 21, 035117.

[67] Hocky, G. M.; Coslovich, D.; Ikeda,A.; Reichman, D. R. Phys. Rev. Lett. 2014, 113, 157801.

[68] Xu, W. S.; Sun, Z. Y.;An, L. J. J. Phys.-Condens Mat. 2012, 24, 325101.

[69] Xu, W. S.; Sun, Z. Y.;An, L. J. J. Chem. Phys. 2012, 137, 104509.

[70] Xu, W. S.; Sun, Z. Y.;An, L. J. Phys. Rev. E 2012, 86, 041506.

[71] Li, Y. W.; Zhu, Y. L.; Sun, Z. Y. arxiv 2015, 1509.07046.

[72] Finney, J. Proc. R. Soc. London, Ser. A 1970, 319, 479.

?教學(xué)研究與改革?

Glass and Glass Transition

LI Yan-WeiSUN Zhao-Yan*AN Li-Jia
(State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun 130022, P. R. China)

Abstract:Glass and glassy state cover a wide range of applications. Different from liquid-crystal transition, glass transition is non-crystalline liquid-solid transition, which is not accompanied with obvious changes in structure. Therefore, we cannot classify it into any known type of phase transition. The research on the glass transition has become a key topic in condensed matter physics for several decades. However, the nature of glass state is still far away from being understood until now. Here we provide some basic physical phenomena of glass-forming systems during the glass transition, review some theoretical perspectives of the glass transition in the past few decades, and introduce some commonly used methods to investigate the structure and dynamics of glass formers, with an attempt to deepen the understanding of disordered solids and glass transition of readers.

Key Words:Glass; Glass transition; Supercooled liquids

中圖分類號：O6；G64

doi:10.3866/PKU.DXHX20160301

*通訊作者，Email: zysun@ciac.ac.cn

基金資助：科技部973項目(2012CB821500)