孫艷麗 王華光 張澤新

(蘇州大學(xué),軟凝聚態(tài)物理及交叉研究中心,蘇州 215006)

1 引 言

玻璃是自然界和生產(chǎn)生活中一種常見(jiàn)的物質(zhì)形態(tài)[1],在食品[2]、醫(yī)學(xué)[3]和材料等[4]領(lǐng)域應(yīng)用廣泛.雖然玻璃材料已經(jīng)使用了數(shù)千年,但是人們對(duì)于液態(tài)如何轉(zhuǎn)化成玻璃態(tài)的物理機(jī)理尚未很好地理解.半個(gè)世紀(jì)以來(lái),各國(guó)科學(xué)家提出了很多關(guān)于玻璃化轉(zhuǎn)變的理論模型[5],其中比較重要的有自由體積模型[6]、Adam-Gibbs熱力學(xué)理論[7]以及模耦合理論(model coupling theory,MCT)[8?10],但是,這些理論也只能解釋轉(zhuǎn)變過(guò)程中的部分實(shí)驗(yàn)現(xiàn)象,玻璃化轉(zhuǎn)變的物理機(jī)理一直是世界性的科學(xué)難題之一.

相較于大量的理論研究,玻璃化轉(zhuǎn)變的實(shí)驗(yàn)研究發(fā)展相對(duì)緩慢,而且傳統(tǒng)的實(shí)驗(yàn)研究都是采用各向同性的球形粒子作為模型來(lái)研究玻璃[8,11?14].眾所周知,真實(shí)的原子和分子體系通常為各向異性的非球形粒子,因此,設(shè)計(jì)各向異性粒子作為模型體系是近年來(lái)玻璃化轉(zhuǎn)變的研究熱點(diǎn)[10,15?21],其中最簡(jiǎn)單的各向異性粒子模型就是橢球膠體體系.例如,Han等[22]通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn),在二維(2D)的橢球膠體懸浮液中,可以觀測(cè)到平動(dòng)玻璃和轉(zhuǎn)動(dòng)玻璃兩步的玻璃化轉(zhuǎn)變.Latz等[23]通過(guò)計(jì)算機(jī)模擬研究發(fā)現(xiàn),在橢球組成的三維(3D)體系中,當(dāng)長(zhǎng)徑比大于2.5時(shí),可以形成一種取向的玻璃.此外,最近的研究表明,各向異性的橢球粒子的確可以引發(fā)新的動(dòng)力學(xué)現(xiàn)象,這為玻璃化轉(zhuǎn)變過(guò)程中研究平動(dòng)和轉(zhuǎn)動(dòng)的相互關(guān)系提供了一個(gè)新的思路[24,25].

研究玻璃化轉(zhuǎn)變的傳統(tǒng)體系都是由各向同性粒子或者各向異性粒子組成的單一體系.最近,各向同性和各向異性粒子混合體系下的玻璃態(tài)行為受到了計(jì)算機(jī)模擬研究者的關(guān)注[26].計(jì)算機(jī)模擬的研究結(jié)果表明,隨著面積分?jǐn)?shù)?的增加,這種各向同性和各向異性粒子混合體系也會(huì)形成玻璃,體系的動(dòng)力學(xué)和結(jié)構(gòu)行為與兩種類型粒子的組分以及體系的?密切相關(guān)[27].然而,在實(shí)驗(yàn)上研究這種混合體系的玻璃化轉(zhuǎn)變的工作還鮮有報(bào)道.本文通過(guò)設(shè)計(jì)物理實(shí)驗(yàn),考察了膠體橢球與圓球混合體系隨著?增加時(shí)發(fā)生玻璃化轉(zhuǎn)變過(guò)程中的結(jié)構(gòu)和動(dòng)力學(xué)的演變規(guī)律,發(fā)現(xiàn)隨著?的增加,體系結(jié)構(gòu)上仍然保持著無(wú)序的狀態(tài),橢球可以有效地抑制圓球結(jié)晶;同時(shí),體系的動(dòng)力學(xué)明顯減緩,弛豫時(shí)間在接近MCT預(yù)測(cè)的玻璃化轉(zhuǎn)變點(diǎn)迅速增大并發(fā)散.通過(guò)對(duì)體系內(nèi)的快速粒子參與的協(xié)同重排區(qū)域(cooperative rearrangements regions,CRRs)進(jìn)行分析,發(fā)現(xiàn)CRRs形狀、大小和位置與橢球密切相關(guān).

2 實(shí)驗(yàn)過(guò)程與方法

研究玻璃化轉(zhuǎn)變的經(jīng)典體系是Kob-Anderson模型,該模型由數(shù)量比例為4:1的大小圓球組成,這樣的體系可以有效抑制結(jié)晶而形成玻璃[28].因此,本文的實(shí)驗(yàn)中,也采用圓球和橢球比例為4:1的混合體系來(lái)研究玻璃化轉(zhuǎn)變,其中圓球是尺寸為1.6μm的聚苯乙烯微球,購(gòu)買自Thermo Scientific公司;橢球是2.5μm的聚苯乙烯微球(購(gòu)買自Thermo Scientific公司),通過(guò)物理拉伸方法制備得到[29],制備所得的橢球長(zhǎng)軸約為2.0μm,短軸約為1.2μm.實(shí)驗(yàn)時(shí),首先使用移液槍吸取適量混合液置于兩蓋玻片之間使其形成單層樣品,利用紫外光膠對(duì)其密封.密封后的樣品經(jīng)過(guò)10 h靜置,然后在顯微鏡載物臺(tái)上靜置平衡1 h后進(jìn)行觀察和數(shù)據(jù)采集(圖1).為了減小可能的樣品槽的邊界效應(yīng),我們選擇了遠(yuǎn)離邊界的區(qū)域進(jìn)行數(shù)據(jù)采集.視野中體系面積分?jǐn)?shù)?的計(jì)算公式為

其中nS(nE)是圓球(橢球)的數(shù)目,σS是圓球粒子的直徑,a是橢球的長(zhǎng)軸,b是橢球的短軸,AT表示觀察的整個(gè)視野的面積.本文采用10組不同?的體系來(lái)進(jìn)行實(shí)驗(yàn)研究,分別為?=0.38,0.50,0.53,0.61,0.63,0.65,0.67,0.72,0.75,0.77.數(shù)據(jù)采集時(shí)先進(jìn)行快拍,其速度為100幀/s,連續(xù)拍攝10 s,然后慢拍,速度為3幀/s,連續(xù)拍攝1000 s.最后以通用的粒子追蹤程序來(lái)確定每個(gè)粒子在不同時(shí)刻的位置信息[30],從而對(duì)體系的結(jié)構(gòu)和動(dòng)力學(xué)行為進(jìn)行定量的分析.

圖1 (a)樣品在物鏡上的示意圖;(b)?=0.38的樣品顯微鏡照片F(xiàn)ig.1.(a)Schematic of experimental setup:a monolayer of colloidal suspension(spheres and ellipsoids mixture)confined between two coverslips,making a quasi-2D sample;(b)bright-field microscope image of the colloidal suspension for ?=0.38.

3 結(jié)果與討論

3.1 體系的結(jié)構(gòu)特征

玻璃化轉(zhuǎn)變過(guò)程中的一個(gè)顯著特征是靜態(tài)結(jié)構(gòu)上保持無(wú)序的狀態(tài).為了研究橢球與圓球混合體系隨?增加過(guò)程中的結(jié)構(gòu)性質(zhì),我們首先分析了體系的徑向分布函數(shù):

其中n表示粒子數(shù)密度,表示區(qū)域內(nèi)的粒子分布,??表示對(duì)時(shí)間和空間上的平均.徑向分布函數(shù)表征了體系內(nèi)距離單個(gè)粒子的一個(gè)給定距離r處找到其他粒子的概率[13,31].圖2給出了體系在不同?下徑向分布函數(shù)的變化規(guī)律.可以看出,隨著探測(cè)距離r的增大,g(r)的峰很快消失,說(shuō)明體系處于一個(gè)非長(zhǎng)程有序的狀態(tài).此外,隨著體系濃度的增加,g(r)在短程距離關(guān)聯(lián)性會(huì)增強(qiáng),峰值變高,插圖所示的體系在不同?下g(r)第一峰峰高g1很好地說(shuō)明了這一點(diǎn).

圖2 徑向分布函數(shù)g(r),從下到上體系面積分?jǐn)?shù)?依次為0.38,0.50,0.53,0.61,0.63,0.65,0.67,0.72,0.75,0.77;插圖為各面積分?jǐn)?shù)?下對(duì)應(yīng)的徑向分布函數(shù)的第一峰峰高g1Fig.2.Pair-correlation function g(r)for spheres at different area fraction ?(from bottom to top:0.38,0.50,0.53,0.61,0.63,0.65,0.67,0.72,0.75,0.77).The inset shows g1,the height of the first peak of g(r),as a function of area fraction ?.

為了直觀地表示體系的有序/無(wú)序的結(jié)構(gòu),計(jì)算了體系的泰森多邊形分布圖[32].由圖3可以直觀地識(shí)別出每個(gè)粒子周圍的近鄰粒子數(shù),不同顏色的多邊形代表粒子周圍的近鄰粒子數(shù)不同.對(duì)于完美的2D晶體體系,每個(gè)粒子周圍有6個(gè)近鄰粒子,對(duì)應(yīng)的是六邊形,表明體系結(jié)構(gòu)有序;而如果粒子周圍不是6個(gè)近鄰粒子,對(duì)應(yīng)的是非六邊形的多邊形,表明體系結(jié)構(gòu)無(wú)序.因此,泰森多邊形可以方便地分析體系內(nèi)缺陷產(chǎn)生、聚集和變化等情況.對(duì)于單分散的圓球體系,隨著?增加,體系中近鄰粒子數(shù)為6的粒子會(huì)不斷增加,最后形成2 D晶體.對(duì)于圓球和橢球混合體系,發(fā)現(xiàn)當(dāng)?較低時(shí),近鄰粒子數(shù)隨機(jī)分布,呈現(xiàn)無(wú)序狀態(tài)(圖3(a));當(dāng)?較大時(shí),粒子的近鄰粒子數(shù)并沒(méi)有發(fā)生本質(zhì)的變化(圖3(b)),這說(shuō)明體系依然保持無(wú)序的結(jié)構(gòu).這主要是由于體系中橢球粒子的存在導(dǎo)致了缺陷的產(chǎn)生,從而有效地抑制了圓球粒子的結(jié)晶,使得體系在高濃度時(shí)仍然保持著無(wú)序的結(jié)構(gòu),有利于形成玻璃態(tài).

為了研究橢球?qū)w系結(jié)構(gòu)的影響,我們計(jì)算了粒子的局部取向序參量:

其中ni是粒子i的最近鄰粒子數(shù),θij是粒子i與其最相連的第j個(gè)粒子的連線與X軸的夾角.當(dāng)ψ6i=1時(shí),表示粒子呈現(xiàn)完美的六角形排列;當(dāng)ψ6i=0時(shí),表示粒子呈現(xiàn)隨機(jī)排列的無(wú)序結(jié)構(gòu)[33].圖4給出了?=0.38(圖4(a))和?=0.63(圖4(b))時(shí),粒子的取向序參量分布圖.由圖4(a)可以看出,當(dāng)體系?比較低時(shí),粒子的取向序參量值總體比較低,說(shuō)明體系為無(wú)序狀態(tài).當(dāng)?增大,發(fā)現(xiàn)橢球周圍的圓球粒子的取向序參量值比其他區(qū)域低,這表明橢球粒子的存在抑制了周圍粒子的結(jié)晶,引起體系的缺陷,使得體系保持無(wú)序結(jié)構(gòu)(圖4(b)).

圖3 泰森多邊形分布圖 (a)?=0.38;(b)?=0.63;圖中的色標(biāo)表示粒子的最近鄰粒子的個(gè)數(shù)(參見(jiàn)圖標(biāo))Fig.3.Snapshots of voronoi diagram corresponding to the systems of ? =0.38(a)and ? =0.63(b)respectively.The color code indicates the number of nearest neighbors for particles(see the legend).

圖4 體系取向序參量分布圖 (a)?=0.38;(b)?=0.63;顏色代表粒子局部取向序參量ψ6i值(參見(jiàn)色帶圖)Fig.4.Snapshot of local order parameter, ψ6i,for ? =0.38(a),and ? =0.63(b).The color of particles represents the value of ψ6i(see the color bar).

3.2 體系的動(dòng)力學(xué)特征

在膠體體系中,體系的動(dòng)力學(xué)強(qiáng)烈地依賴于粒子的濃度.為了量化橢球與圓球混合體系在不同?的動(dòng)力學(xué)特征,我們計(jì)算了體系的均方位移(mean square displacement,MSD):

其中ri(t)是粒子i在時(shí)刻t的位置,ri(0)是該粒子最初的位置[11].圖5給出了隨著?的增加,體系MSD的變化規(guī)律.發(fā)現(xiàn)當(dāng)?較低時(shí),體系中粒子的動(dòng)力學(xué)較快,隨時(shí)間增加,MSD呈現(xiàn)類似自由擴(kuò)散的線性增長(zhǎng);當(dāng)?增加時(shí),MSD增長(zhǎng)趨勢(shì)逐漸變緩.當(dāng)?達(dá)到0.75或者更高時(shí),發(fā)現(xiàn)MSD出現(xiàn)了平臺(tái)期,并且?越高,平臺(tái)期對(duì)應(yīng)的時(shí)間越長(zhǎng).這種行為是因?yàn)樵诟邼舛润w系下,粒子的運(yùn)動(dòng)受到周圍近鄰粒子的限制.普遍認(rèn)為的物理圖景是:粒子被束縛在近鄰粒子所構(gòu)成的“籠子”(cage)里,導(dǎo)致MSD出現(xiàn)平臺(tái)期;隨著時(shí)間的增加,有些粒子會(huì)逃出“籠子”,從而在長(zhǎng)時(shí)間尺度上MSD曲線又會(huì)上升;而當(dāng)?非常高時(shí),粒子很難從“籠子”中逃離,MSD的平臺(tái)期延長(zhǎng),這種MSD隨?的變化行為是玻璃化轉(zhuǎn)變中最顯著的動(dòng)力學(xué)特征[11].

圖5 不同面積分?jǐn)?shù)下粒子均方位移的變化Fig.5.Mean squared displacement of particles for different area fractions.

在玻璃動(dòng)力學(xué)研究中,另一個(gè)經(jīng)典的物理學(xué)參量就是自散射函數(shù)(self-intermediate scattering function),

其中N為粒子總數(shù),q為自散射矢量,它的大小由體系靜態(tài)結(jié)構(gòu)因子決定[22].根據(jù)體系的MSD平臺(tái)期的動(dòng)力學(xué)行為,我們選定的自散射矢量q=4.0.圖6給出了不同?下自散射函數(shù)的變化趨勢(shì),可以看出,當(dāng)體系濃度較低時(shí),該函數(shù)在所研究的時(shí)間尺度上衰減得十分迅速,表明較快的動(dòng)力學(xué)弛豫.隨著體系?的增加,衰減過(guò)程逐漸減緩,這說(shuō)明粒子運(yùn)動(dòng)弛豫趨緩.當(dāng)體系濃度達(dá)到0.75及以上時(shí),自散射函數(shù)出現(xiàn)了類似MSD的平臺(tái)期.這同樣表明,體系濃度增加后,出現(xiàn)的“籠子”效應(yīng)抑制了體系的運(yùn)動(dòng)弛豫,這是玻璃化轉(zhuǎn)變中最典型的動(dòng)力學(xué)行為.

圖6 不同面積分?jǐn)?shù)下自散射函數(shù)FS(q,t)的衰減變化Fig.6.Self-intermediate scattering function,FS(q,t)for different area fractions.

圖7 不同面積分?jǐn)?shù)對(duì)應(yīng)的弛豫時(shí)間,實(shí)線表示MCT進(jìn)行的參數(shù)擬合,虛線表示擬合得到的玻璃化轉(zhuǎn)變點(diǎn)Fig.7.The relaxation time corresponding to different area fractions.The solid line indicates the fitting curve according to the MCT,and the position of the dotted line represents the fitted glass transition point.

根據(jù)MCT,自散射函數(shù)FS(q,t)長(zhǎng)時(shí)間的弛豫過(guò)程滿足e?(t/τ)β,體系的弛豫時(shí)間可以通過(guò)擬合FS(q,t)的長(zhǎng)時(shí)間弛豫過(guò)程得到[10].為了研究和預(yù)測(cè)體系的玻璃化轉(zhuǎn)變點(diǎn),我們對(duì)體系的弛豫時(shí)間和?的關(guān)系進(jìn)行擬合[8]:τ～(?c??)?γ,其中?c是MCT預(yù)測(cè)的玻璃化轉(zhuǎn)變點(diǎn);γ為體系擬合得到的參數(shù),γ=1.擬合結(jié)果如圖7所示,實(shí)線表示擬合曲線,虛線是擬合得到的玻璃化轉(zhuǎn)變點(diǎn),?c≈0.76.我們發(fā)現(xiàn),隨著?的增加,尤其在接近玻璃化轉(zhuǎn)變點(diǎn)時(shí),體系的弛豫時(shí)間迅速增大并發(fā)散,符合MCT.

3.3 協(xié)同重排

協(xié)同重排是玻璃體系的特征物理現(xiàn)象.為了研究橢球與圓球混合體系中玻璃化轉(zhuǎn)變過(guò)程中的協(xié)同重排行為,我們分析了快速粒子參與協(xié)同重排的空間分布[34].根據(jù)前人研究的方法,將一個(gè)弛豫時(shí)間內(nèi)運(yùn)動(dòng)距離為前10%的圓球粒子定義為快速粒子.圖8(a)和圖8(b)展示了玻璃體系?=0.77中CRRs的典型圖片.其中快速粒子用紅色圓點(diǎn)表示,其余圓球粒子顯示為黑色圓點(diǎn),橢球粒子顯示為藍(lán)色橢圓,箭頭指示粒子運(yùn)動(dòng)方向.可以看出,快速粒子呈現(xiàn)出團(tuán)簇狀和鏈狀結(jié)構(gòu),并且集中出現(xiàn)在橢球粒子的周圍.這表明體系的CRRs與橢球的位置有很大的關(guān)聯(lián)性,橢球的存在改變了圓球粒子局部的動(dòng)力學(xué)行為.

圖8 快速粒子參與協(xié)同重排的分布圖(?=0.77)Fig.8.Snapshots of the particles that take part in cooperatively rearrangement for ?=0.77.

4 總 結(jié)

研究了橢球與圓球混合體系隨著?增加時(shí)發(fā)生玻璃化轉(zhuǎn)變的行為.通過(guò)結(jié)構(gòu)參量徑向分布函數(shù)、取向序參量以及泰森多邊形的分析,發(fā)現(xiàn)隨著?的增加,體系結(jié)構(gòu)上仍然保持為無(wú)序的狀態(tài),表明體系進(jìn)入玻璃態(tài).在動(dòng)力學(xué)方面,隨著體系濃度的增加,均方位移和自散射函數(shù)均出現(xiàn)了平臺(tái)期,體系的動(dòng)力學(xué)明顯減緩,粒子的運(yùn)動(dòng)受到了周圍粒子的“籠子”效應(yīng)限制,呈現(xiàn)典型的玻璃化轉(zhuǎn)變的動(dòng)力學(xué)演化行為.體系弛豫時(shí)間隨著?的變化符合模耦合理論,隨?的增加而增長(zhǎng),并在玻璃化轉(zhuǎn)變點(diǎn)發(fā)散.體系中粒子的協(xié)同重排區(qū)域出現(xiàn)了團(tuán)簇和和鏈狀結(jié)構(gòu),這些重排區(qū)域與橢球的存在密切相關(guān).

此外,本文通過(guò)改變橢球與圓球混合體系中面積分?jǐn)?shù)?,研究了體系玻璃化轉(zhuǎn)變過(guò)程的結(jié)構(gòu)和動(dòng)力學(xué).對(duì)于兩種類型粒子的混合比例,按照經(jīng)典的圓球混合體系Kob-Anderson模型來(lái)確定.未來(lái)的研究工作可以系統(tǒng)地改變混合比例,考察混合比對(duì)體系的結(jié)構(gòu)和動(dòng)力學(xué)產(chǎn)生的影響.

[1]Angell C A 1995Science267 1924

[2]Slade L,Levine H,Ievolella J,Wang M 1993J.Sci.Food Agric.63 133

[3]Zahn K,Lenke R,Maret G 1999Phys.Rev.Lett.82 2721

[4]Debenedetti P G,Truskett T M,Lewis C P,Stillinger F H 2001Adv.Chem.Eng.28 21

[5]Wen P 2017Acta Phys.Sin.66 176407(in Chinese)[聞平2017物理學(xué)報(bào)66 176407]

[6]Fox T G,Flory P J 1950J.Appl.Phys.21 581

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

[8]van Megen W,Underwood S M 1993Phys.Rev.Lett.70 2766

[9]Debenedetti P G,Stillinger F H 2001Nature410 259

[10]Gotze W,Sjogren L 1992Rep.Prog.Phys.55 241

[11]Weeks E R,Crocker J C,Levitt A C,Schofield A,Weitz D A 2000Science287 627

[12]Kegel W K,Van B A 2000Science287 290

[13]Zhang Z,Xu N,Chen D T N,Yunker P,Alsayed A M,Aptowicz K B,Habdas P,Liu A J,Nagel S R,Yodh A G 2009Nature459 230

[14]Yunker P,Zhang Z,Yodh A G 2010Phys.Rev.Lett.104 015701

[15]Chong S H,Moreno A J,Sciortino F,Kob W 2005Phys.Rev.Lett.94 215701

[16]Yatsenko G,Schweizer K S 2007J.Chem.Phys.126 014505

[17]Tripathy M,Schweizer K S 2009J.Chem.Phys.130 244906

[18]Jadrich R,Schweizer K S 2012Phys.Rev.E86 061503

[19]Kramb R C,Zhang R,Schweizer K S,Zukoski C F 2010Phys.Rev.Lett.105 055702

[20]Kramb R C,Zhang R,Schweizer K S,Zukoski C F 2011J.Chem.Phys.134 014503

[21]Kang K,Dhont J K G 2013Phys.Rev.Lett.110 015901

[22]Zheng Z,Wang F,Han Y 2011Phys.Rev.Lett.107 065702

[23]Letz M,Schilling R,Latz A 2000Phys.Rev.E62 5173

[24]Jadrich R,Schweizer K S 2012Phys.Rev.E86 061503

[25]Kramb R C,Zhang R,Schweizer K S,Zukoski C F 2011J.Chem.Phys.134 014503

[26]Xu W S,Duan X,Sun Z Y,An L J 2015J.Chem.Phys.142 224506

[27]Takae K,Onuki A 2013Phys.Rev.E88 042317

[28]Toxvaerd S,Schr?der T B,Dyre J C 2009J.Chem.Phys.130 224501

[29]Champion J A,Katare Y K,Mitragotri S 2007PNAS104 11901

[30]Liu H X,Chen K,Hou M Y 2015Acta Phys.Sin.64 116302(in Chinese)[劉海霞,陳科,厚美瑛2015物理學(xué)報(bào)64 116302]

[31]Chen K 2017Acta Phys.Sin.66 178201(in Chinese)[陳科2017物理學(xué)報(bào) 66 178201]

[32]Gasser U 2009J.Phys.:Condens.Matter21 203101

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

[34]Zhang Z,Yunker P J,Habdas P,Yodh A G 2011Phys.Rev.Lett.107 208303