橢球與圓球混合膠體體系的玻璃化轉(zhuǎn)變?

孫艷麗 王華光 張澤新

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

1 引 言

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

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

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

2 實驗過程與方法

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

其中nS(nE)是圓球(橢球)的數(shù)目,σS是圓球粒子的直徑,a是橢球的長軸,b是橢球的短軸,AT表示觀察的整個視野的面積.本文采用10組不同?的體系來進(jìn)行實驗研究,分別為?=0.38,0.50,0.53,0.61,0.63,0.65,0.67,0.72,0.75,0.77.數(shù)據(jù)采集時先進(jìn)行快拍,其速度為100幀/s,連續(xù)拍攝10 s,然后慢拍,速度為3幀/s,連續(xù)拍攝1000 s.最后以通用的粒子追蹤程序來確定每個粒子在不同時刻的位置信息[30],從而對體系的結(jié)構(gòu)和動力學(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)變過程中的一個顯著特征是靜態(tài)結(jié)構(gòu)上保持無序的狀態(tài).為了研究橢球與圓球混合體系隨?增加過程中的結(jié)構(gòu)性質(zhì),我們首先分析了體系的徑向分布函數(shù):

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

圖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ù)?下對應(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 ?.

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

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

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

圖3 泰森多邊形分布圖 (a)?=0.38;(b)?=0.63;圖中的色標(biāo)表示粒子的最近鄰粒子的個數(shù)(參見圖標(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值(參見色帶圖)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 體系的動力學(xué)特征

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

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

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

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

其中N為粒子總數(shù),q為自散射矢量,它的大小由體系靜態(tài)結(jié)構(gòu)因子決定[22].根據(jù)體系的MSD平臺期的動力學(xué)行為,我們選定的自散射矢量q=4.0.圖6給出了不同?下自散射函數(shù)的變化趨勢,可以看出,當(dāng)體系濃度較低時,該函數(shù)在所研究的時間尺度上衰減得十分迅速,表明較快的動力學(xué)弛豫.隨著體系?的增加,衰減過程逐漸減緩,這說明粒子運動弛豫趨緩.當(dāng)體系濃度達(dá)到0.75及以上時,自散射函數(shù)出現(xiàn)了類似MSD的平臺期.這同樣表明,體系濃度增加后,出現(xiàn)的“籠子”效應(yīng)抑制了體系的運動弛豫,這是玻璃化轉(zhuǎn)變中最典型的動力學(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ù)對應(yīng)的弛豫時間,實線表示MCT進(jìn)行的參數(shù)擬合,虛線表示擬合得到的玻璃化轉(zhuǎ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)長時間的弛豫過程滿足e?(t/τ)β,體系的弛豫時間可以通過擬合FS(q,t)的長時間弛豫過程得到[10].為了研究和預(yù)測體系的玻璃化轉(zhuǎn)變點,我們對體系的弛豫時間和?的關(guān)系進(jìn)行擬合[8]:τ～(?c??)?γ,其中?c是MCT預(yù)測的玻璃化轉(zhuǎn)變點;γ為體系擬合得到的參數(shù),γ=1.擬合結(jié)果如圖7所示,實線表示擬合曲線,虛線是擬合得到的玻璃化轉(zhuǎn)變點,?c≈0.76.我們發(fā)現(xiàn),隨著?的增加,尤其在接近玻璃化轉(zhuǎn)變點時,體系的弛豫時間迅速增大并發(fā)散,符合MCT.

3.3 協(xié)同重排

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

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

4 總 結(jié)

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

此外,本文通過改變橢球與圓球混合體系中面積分?jǐn)?shù)?,研究了體系玻璃化轉(zhuǎn)變過程的結(jié)構(gòu)和動力學(xué).對于兩種類型粒子的混合比例,按照經(jīng)典的圓球混合體系Kob-Anderson模型來確定.未來的研究工作可以系統(tǒng)地改變混合比例,考察混合比對體系的結(jié)構(gòu)和動力學(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é)報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é)報64 116302]

[31]Chen K 2017Acta Phys.Sin.66 178201(in Chinese)[陳科2017物理學(xué)報 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