基于分離粒子隨機擴散理論的色譜模擬

郭 彥 傅應(yīng)強 孫寅璐 陳天南 趙健偉,*

(1南京信息工程大學環(huán)境科學與工程學院,江蘇省大氣環(huán)境監(jiān)測與污染控制高技術(shù)研究重點實驗室,南京210044; 2南京大學化學化工學院,生命分析化學國家重點實驗室,南京210093)

1 引言

分子擴散,尤其是在微納米尺度受限空間內(nèi)的分子擴散是化學學科中的核心問題之一,微化學反應(yīng)器、分子篩催化、晶體生長、分離富集等基本化學過程都與之密切相關(guān).1-6分子擴散行為對于這些化學過程中產(chǎn)物的組成和分布具有決定作用.了解該過程對提升器件性能、優(yōu)化操作以及新型反應(yīng)器件的設(shè)計等都尤為關(guān)鍵.對分子擴散行為的研究,除了利用一些先進的實驗方法外,隨著計算技術(shù)的進步,模擬仿真也成為有力的研究手段.計算機模擬在研究分子擴散過程中有顯著的優(yōu)勢,它可以從微觀模擬每個粒子的擴散行為,動態(tài)追蹤并顯示擴散過程,分析各種環(huán)境因素對擴散行為的影響,了解不同操作條件下的擴散機制,進而解釋實驗中所觀察到的現(xiàn)象,為實驗化學提供微觀理論基礎(chǔ).

分子動力學在描述微觀擴散行為方面具有獨特的優(yōu)勢,7,8但其計算工作量巨大,現(xiàn)有的計算技術(shù)水平還不能為模擬較大的擴散體系提供軟硬件支撐.而隨機行走的優(yōu)勢在于提取了制約擴散行為的關(guān)鍵信息,對該行為進行了合理近似,從而有效簡化了計算工作量,使從微觀角度模擬分離擴散成為可能.隨機行走模型已在一些擴散問題的研究中得到較為廣泛的應(yīng)用,如土壤介質(zhì)、生物組織、藥物輸運、神經(jīng)元傳導(dǎo)和晶體生長等.9-13

氣相色譜以其高靈敏度、高分辨率和高選擇性已成為最為常用和重要的分離分析方法.基于宏觀的擴散理論,人們對于色譜的保留值規(guī)律、峰形規(guī)律、色譜條件的優(yōu)化等已有大量的實驗研究.14-16色譜理論的發(fā)展,亟需從微觀角度對色譜的分離規(guī)律以及條件控制開展模擬研究,Siepmann等17應(yīng)用蒙特卡羅方法模擬了烴類在氣液色譜界面上的吸附和分配情況.O?Connell等18應(yīng)用格子-偶極模型在介觀尺度上模擬了肽鏈在反相色譜上的吸附行為. Thornton等19應(yīng)用MM+力場和AM1半經(jīng)驗方法計算了脂肪族和芳香族溶質(zhì)的范德華半徑、與固定相和流動相的相互作用,及這些作用情況在反相色譜分離中的應(yīng)用.Sun等20應(yīng)用分子動力學研究了配體密度對于蛋白質(zhì)在疏水電荷誘導(dǎo)色譜上的吸附、脫附、構(gòu)象變化的影響.Sarti等21應(yīng)用數(shù)學模型研究了蛋白質(zhì)在膜親和相色譜中的吸附平衡和流出曲線.值得注意的是,這些模擬研究盡管對色譜的吸附熱力學、待分離物質(zhì)的物性以及溶劑化環(huán)境等作了較為深入的研究,但尚未涉及到色譜動力學的內(nèi)容,如追蹤每個分離粒子的運動軌跡,驗證速率理論方程和流出曲線方程,研究粒子在色譜柱內(nèi)分布平衡和遷移過程等.這些動力學方面的研究將有助于進一步了解色譜的分離本質(zhì),提升分離器件性能和優(yōu)化分離條件.

氣相色譜中分子擴散既有隨機運動特征,又受限于微尺度空間的具體環(huán)境.基于這兩個特點,本文采用粒子在排布周期性勢壘的二維空間中隨機行走的方法模擬了色譜中粒子的擴散行為,利用粒子與勢壘的彈性碰撞簡化了色譜中待分離粒子與固定相之間的相互作用.重點考察了固定相的填充率、固定相的形狀和柱長因素對色譜動力學行為的影響.

2 模擬與計算方法

采用自行開發(fā)的隨機行走軟件,22,23對色譜填充柱中粒子擴散過程進行模擬(如圖1所示):粒子在給定的初速度條件下在二維受限空間內(nèi)按照隨機的方式運動,在預(yù)設(shè)時間后粒子運動停止,根據(jù)給定的條件計算粒子的運動軌跡并分析粒子的運動行為.

2.1 分離柱

圖1 算法模型示意圖Fig.1 Schematic diagram of the simulated modelThe tetragonal region,the blank area,the dashed box,the dot,and the curve represent the stationary phase,the cavity,the cell, the particle,and its trajectory,respectively.

分離柱定義為二維平面空間,該空間在x方向上無限延伸,y方向上有上下兩個邊界.空間內(nèi)周期排布著單元格,如圖1b中的虛線框所示.單元格中的多邊形為勢壘,用來模擬色譜中的固定相,如圖1b中的正方形所示,勢壘粒徑為正方形對角線長度.用勢壘面積除以單元格面積可以求出勢壘填充率.二維空間內(nèi)粒子可以在空腔內(nèi)(無勢壘區(qū)域)自由運動但不能進入這些勢壘.基于色譜分離的特點,其熱力學和動力學過程都發(fā)生在曲線流出方向,也就是色譜可以近似看成為一維體系,因而本模擬主要關(guān)注粒子擴散的x方向.如無特別說明,分離柱均是長度為1 m的填充柱.

2.2 分離物

以忽略質(zhì)量和體積的粒子作為分離物,如圖1b中的圓點所示.粒子的位置用其所在單元格當中的相對于單元格中心點的相對位置來描述:

其中,Parrel是粒子的相對位置對應(yīng)的位矢,Parabs是粒子的絕對位置對應(yīng)的位矢,Cellabs是單元格中心點的絕對位置對應(yīng)的位矢.在預(yù)設(shè)時間后,粒子在二維平面空間內(nèi)的運動完全停止,即相對于所在單元格的絕對位置不再發(fā)生變化.

粒子的速度由兩部分組成.一是粒子的擴散速度,該速度為溫度相關(guān)的一個變量.模擬中溫度采用室溫298 K,根據(jù)氣體布朗運動估算,粒子的擴散速度設(shè)定為2×10-3m·s-1.粒子擴散的運動方向隨機生成,在以粒子為中心的0≤α＜2π范圍內(nèi)等概率分布,α為方向角.當粒子達到色譜邊界上時,將按照彈性碰撞模型,即粒子的擴散速度大小不變,方向按照反射定律進行反射后繼續(xù)在色譜柱中運動.

另一個速度為給待分離粒子施加的橫向速度,用來模擬載氣流速,在程序中定義為沿x軸方向的線速度(Vx,單位為mm·s-1).

2.3 擴散過程

粒子和勢壘的相互作用可以用多種模型來描述,本模擬中采用彈性碰撞模型簡化了粒子與固定相的作用.假定粒子在某段時間之內(nèi)速度的方向和大小不發(fā)生變化,即速度和時間步(即粒子的實際行走物理單位與在模擬程序中的時間單位之間的比值,程序中采用0.005 s·step-1)決定了粒子一步隨機行走的步長;在不受到勢壘阻礙的情況下,根據(jù)速度、每步時間和粒子的當前位置即可確定粒子完成一步行走之后所處的位置.若粒子與勢壘發(fā)生碰撞,則可計算出粒子運動到碰撞點所需要的運動時間,將當前時間減去運動時間更新當前時間,將速度按照反射定律進行反射,據(jù)此可以確定碰撞之后的粒子運動軌跡.根據(jù)粒子若干個時間段內(nèi)位置參數(shù)的時間積分,即可以確定整個時間段內(nèi)粒子運動的軌跡,如圖1c所示.

3 結(jié)果與討論

應(yīng)用如上描述的模擬程序,研究了氣相色譜中待分離粒子在固定相中的擴散行為.重點考察了固定相的填充率、固定相的形狀和柱長等結(jié)構(gòu)特性因素對于色譜動力學行為的影響.

3.1 固定相的填充率

固定相填充質(zhì)量的優(yōu)劣直接影響到色譜填充柱的分離效果,在模擬程序中勢壘填充率大小可以反映出色譜中固定相的填充情況.在周期單元一定的條件下,勢壘粒徑設(shè)置小,則填充率小,相當于一般實驗經(jīng)驗中采用固定相填充的較為疏松的色譜柱;勢壘粒徑大,填充率大,相當于固定相填料粒度小,填充的緊實.模擬中勢壘粒徑大小分別設(shè)置為0.80、0.84、0.88、0.92、0.96和0.98 mm,相對應(yīng)的填充率為32%、35%、39%、42%、46%和48%,考察粒子在500 s時間內(nèi)的擴散行為.圖2給出了相同的載氣流速(Vx=2 mm·s-1)條件下,隨著勢壘粒徑的變化,粒子沿x軸向擴散的統(tǒng)計分布圖.每個粒徑條件下的柱狀圖都分別為1000個粒子擴散的統(tǒng)計結(jié)果.圖中曲線為對統(tǒng)計數(shù)據(jù)點的擬合曲線,由圖可見,在所設(shè)定的填充率下,粒子在柱內(nèi)x軸向上的濃度分布符合高斯正態(tài)分布,峰形對稱性好,與實際色譜中粒子的擴散行為非常吻合.當勢壘粒徑從0.80增大到0.98 mm,相當于固定相填料粒度由大變小,圖2顯示峰形逐漸變窄變尖銳,柱效增加,同時峰的中心位置即粒子在柱內(nèi)x軸向上濃度分布的平均距離減小.

圖2 相同載氣流速(Vx=2 mm·s-1)條件下,粒子在預(yù)定時間(500 s)內(nèi)沿x軸向擴散隨勢壘粒徑(0.80、0.84、0.88、0.92、0.96和0.98 mm)變化的統(tǒng)計分布圖Fig.2 Statistical distribution graph of the x axial distances passed by the particles(500 s)with the barrier sizes(0.80, 0.84,0.88,0.92,0.96,and 0.98 mm)under the fixed carrier gas flow rate(Vx=2 mm·s-1)

半峰寬是色譜分離流出曲線的一項重要參數(shù),涉及到色譜動力學和組分分離方程.圖2中統(tǒng)計峰的半峰寬隨著勢壘粒徑的增大,總體變化趨勢是在減小,說明柱效有增大趨勢.我們也系統(tǒng)考察了在其它載氣流速下,相應(yīng)的半峰寬特征,均表現(xiàn)為隨勢壘粒徑的增大而減小,而在相同粒徑條件下,不同載氣流速對半峰寬基本沒有影響.色譜柱參數(shù)中,固定相的填充情況可以用渦流擴散項A來描述:24

上式中,λ為填充不均勻性因子,隨機行走模擬中,勢壘都是以相同的密堆方式堆積,因而勢壘粒徑的改變對λ的影響可以忽略.dp為填料粒徑,當設(shè)置的勢壘粒徑變大時,對應(yīng)于一般實際實驗時使用了小顆粒的填料,填充得緊實密集,使得粒子在擴散中遷移路徑較為接近.由經(jīng)典的van Deemter方程可知dp的減小將導(dǎo)致塔板高度變小,柱效提高,半峰寬減小.在固定相填充率與色譜柱效研究方面,該模擬結(jié)果與一般實驗經(jīng)驗一致.

圖3 在不同的載氣流速條件(Vx=4,2,1 mm·s-1)下,粒子在預(yù)定時間(500 s)內(nèi)x軸向擴散平均距離隨勢壘粒徑(0.80、0.84、0.88、0.92、0.96和0.98 mm)變化的關(guān)系圖Fig.3 Mean x axial distances passed(500 s)versus the barrier particle sizes(0.80,0.84,0.88,0.92,0.96,and 0.98 mm)under fixed carrier gas flow rates(Vx=4,2,1 mm·s-1)

勢壘粒徑對粒子擴散行為的影響還體現(xiàn)在粒子在x軸向上的擴散,也就是粒子在x軸向上的濃度分布,實質(zhì)上描述了粒子的色譜流出曲線.圖3為在不同的壓力條件,即不同載氣流速(Vx分別為4、2和1 mm·s-1)下粒子在預(yù)定時間(500 s)內(nèi)在x軸向擴散距離的平均值與勢壘粒徑大小的關(guān)系圖.由圖可知,在相同壓力條件下,勢壘粒徑越小,粒子沿x軸向擴散距離越大,對應(yīng)于色譜中使用大粒徑的填料,填充率減小,粒子在擴散時與固定相的碰撞作用減少,粒子可以沿相對較直的路徑較快地通過色譜柱.在相同勢壘粒徑條件下,載氣流速越大,粒子的軸向擴散距離越大,相當于在實驗過程中增加了軸向驅(qū)動力,粒子的流動行為更加越明顯.Darcy法則給出對于粒徑較小的填料,可以觀察到流動相線速度與粒徑平方的線性關(guān)系,由于模擬條件設(shè)定粒徑相對較大,相同時間下的軸向距離與勢壘粒徑呈一元回歸響應(yīng).

粒子位移平方的平均值稱為均方位移(MSD),是用來表征粒子運動性質(zhì)的微觀物理量.在模擬程序中定義為

其中N是總步數(shù),dk表示粒子在k步后的位置矢量(d0是初始位置),|dk+t-dk|表示k步到k+t步的位移矢量.

圖4為在相同的載氣流速(Vx=2 mm·s-1)條件下,填充物粒徑大小從0.80增大到0.98 mm時MSD變化圖.由圖可見,當填充物粒徑較大時,MSD隨時間線性變化,表明粒子的運動行為為自由擴散過程,由斜率可以求出擴散系數(shù).如勢壘粒徑為0.96 mm時,擴散系數(shù)為0.049 mm2·s-1,這個值比實際氣體擴散系數(shù)略大.25,26當勢壘粒徑減小,MSD隨時間呈曲線變化,表明粒子除了在固定相中有擴散行為,還有明顯的流動行為,27該流動行為的驅(qū)動力來自于外加載氣壓力和固定相對粒子擴散阻力的綜合作用.因而勢壘粒徑越小,流動行為越明顯,這也與圖3的結(jié)果相一致.

圖4 相同載氣流速(Vx=2 mm·s-1)條件下,不同勢壘粒徑大小(0.80、0.84、0.88、0.92、0.96和0.98 mm)的均方位移(MSD)Fig.4 Mean square displacement(MSD)of several barrier particle sizes(0.80,0.84,0.88,0.92,0.96,and 0.98 mm)with the fixed carrier gas flow rate(Vx=2 mm·s-1)

3.2 固定相的形狀

本文也考察了固定相填料的形狀(勢壘從四邊形到接近圓形的二十邊形)對色譜的影響,模擬結(jié)果如圖5插圖所示.圖中盡管固定相形狀不同,但粒子擴散的峰形和濃度分布的平均距離都非常接近,表明在模擬中固定相形狀對粒子的擴散行為并沒有明顯的影響.在一般色譜實驗經(jīng)驗中,色譜柱采用圓形填料,是由于圓形填料更容易填充緊實,在模擬中,由于固定相都是以預(yù)設(shè)的密集方式排布,因而形狀影響沒有實際中重要.

不同固定相形狀對應(yīng)的半峰寬也都在0.2 m左右的范圍內(nèi)波動,如圖5所示,也說明密集排布的固定相形貌對粒子的擴散行為影響不大,但要注意到,當固定相填料形狀開始接近于圓形如十二、十六和二十邊形時,半峰寬略有變小,相應(yīng)于在固定相密集排布的基礎(chǔ)上,填充的不均勻性因子λ減小也起到了一定作用,使得柱效略有增加.

盡管本文對色譜模擬條件進行了適當?shù)暮喕?這種固定相排布相同但形狀不同對分子擴散微弱的影響,為有些如載體相同固定液不同且分離效果相近的實際實驗,提供了部分理論支持.28,29

3.3 柱 長

圖5 半峰寬與勢壘形狀的關(guān)系圖Fig.5 Relation between FWHM and the barrier shapesThe inset is statistical distribution graph of the x axial distances (500 s)with the barrier shapes(polygon n=4,6,8,12,16,20)under the fixed carrier gas flow rate(Vx=2 mm·s-1)and barrier particle size(0.88 mm).FWHM:full width at half maximum

圖6 給定載氣流速(Vx=2 mm·s-1)和固定相粒徑大小(0.88 mm),固定相形狀(四邊形),粒子沿x軸向擴散隨分離時間(500、1500、3000和5000 s)變化的統(tǒng)計分布圖Fig.6 Statistical distribution graph of the x axial distances with the different of calculation times(500,1500,3000 and 5000 s)under the fixed carrier gas flow rate(Vx=2 mm·s-1), barrier size(0.88 mm)and barrier shape(tetragon) Inset is the FWHM versus square root of column length.

進一步考察了柱長對粒子分離擴散的影響.圖6為在相同的載氣流速下(Vx=2 mm·s-1),相同的固定相填充率及固定相形狀條件下,通過計算步長的增加(1×105、3×105、6×105和10×105步),即在分離時間分別為500、1500、3000和5000 s時考察粒子的擴散行為.在設(shè)定的分離時間條件下,相當于色譜柱長分別為1、3、6和10 m.由圖6可以看出,分離時間設(shè)定越長,即柱長越長,粒子的濃度分布峰越鈍越寬.理想的氣相色譜通常滿足以下條件:色譜柱內(nèi)任何點上兩相的比例恒定;載氣流在色譜柱內(nèi)任何部位均一樣;在兩相中均不發(fā)生分子的縱向擴散作用;組分分子的分配平衡能瞬間完成.30一般的實驗經(jīng)驗中,通常分子的縱向擴散和分配平衡是峰展寬的決定因素.本模擬中,為了能夠追蹤粒子的擴散過程,簡化了粒子與固定相的作用,因而峰展寬的主要因素為粒子的縱向擴散.圖6中柱長越長,在相同起點上的粒子由于縱向擴散作用,使得峰展寬嚴重.圖6插圖為圖6中粒子濃度分布峰的半峰寬隨柱長平方根的變化圖,兩者呈線性變化,符合色譜流出曲線及質(zhì)量平衡規(guī)律.隨著柱長的增加,半峰寬指數(shù)級的增大,相應(yīng)于柱效變差,這也為在達到完全分離前提下盡量用短柱的實驗經(jīng)驗,提供了可靠的模擬基礎(chǔ).

4 結(jié)論

利用隨機行走的算法對氣相色譜中粒子在固定相中的擴散行為進行了模擬,模擬了粒子擴散的運動軌跡和流出曲線,流出曲線均呈高斯分布,與實際色譜粒子的擴散過程一致.文中重點考察了固定相的填充情況、固定相的形狀和柱長等性能指標對粒子擴散行為的影響,模擬結(jié)果表明固定相填充情況和柱長對色譜出峰的影響與實際經(jīng)驗相符,而在相同的密堆排列下,固定相形狀基本對色譜動力學無明顯影響.

(1) Grill,L.;Dyer,M.;Lafferentz,L.;Persson,M.;Peters,M.V.; Hecht,S.Nat.Nanotechnol.2007,2,687.doi:10.1038/ nnano.2007.346

(2)Yariv,E.;Ben-Dov,G.;Dorfman,K.D.Europhys.Lett.2005, 71(6),1008.doi:10.1209/epl/i2005-10171-6

(3) Corma,A.;Díaz-Caba?as,M.J.;Jordá,J.L.;Martínez,C.; Moliner,M.Nature 2006,443,842.doi:10.1038/nature05238

(4) Ridgway,D.;Broderick,G.;Lopez-Campistrous,A.;Ruaini, M.;Winter,P.;Hamilton,M.;Boulanger,P.;Kovalenko,A.; Ellison,M.J.Biophys.J.2008,94(10),3748.doi:10.1529/ biophysj.107.116053

(5) Persson,A.I.;Larsson,M.W.;Stenstr?m,S.;Ohlsson,B.J.; Samuelson,L.;Wallenberg,L.R.Nat.Mater.2004,3,677.doi: 10.1038/nmat1220

(6) Chou,C.F.;Bakajin,O.;Turner,S.W.P.;Duke,T.A.J.;Chan, S.S.;Cox,E.C.;Craighead,H.G.;Austin,R.H.Proc.Natl. Acad.Sci.U.S.A.1999,96(23),13762.

(7)Wang,F.Y.;Liu,Y.H.;Yin,X.;Wang,N.;Wang,D.X.;Gao,Y. J.;Zhao,J.W.J.Appl.Phys.2010,108,074311.doi:10.1063/ 1.3477323

(8)Liu,Y.H.;Zhao,J.W.;Wang,F.Y.Phys.Rev.B 2009,80, 1154117.

(9) Rao,P.S.C.;Jessup,R.E.;Addiscott,T.M.Soil Sci.1982,133, 342.doi:10.1097/00010694-198206000-00002

(10) Safford,R.E.;Bassingthwaighte,E.A.;Bassingthwaighte,J.B. J.Gen.Physiol.1978,72,513.doi:10.1085/jgp.72.4.513

(11) Santamaria,F.;Wils,S.;De Schutter,E.;Augustine,G.J. Neuron 2006,52,635.doi:10.1016/j.neuron.2006.10.025

(12) Saltzman,W.M.Drug Delivery:Engineering Principles for Drug Therapy;Oxford University Press:New York,2001;pp 1-30.

(13) Zubarev,E.R.;Pralle,M.U.;Li,L.;Stupp,S.I.Science 1999, 283,523.doi:10.1126/science.283.5401.523

(14)McNair,H.M.;Miller,J.M.Basic Gas Chromatography;John Wiley&Sons:Hoboken,New Jersey,2009;pp 1-267.

(15) Poole,C.F.The Essence of Chromatography;Elsevier Science B.V.:Amsterdam,the Netherlands,2003;pp 80-86.

(16) Heftmann,E.Fundamentals and Applications of Chromatography and Related Differential Migration Methods, in Chromatography;Elsevier Science B.V.:Amsterdam,the Netherlands,2004;pp 319-364.

(17)Wick,C.D.;Siepmann,J.I.;Schure,M.R.Anal.Chem.2002, 74,3518.doi:10.1021/ac0200116

(18) Makrodimitris,K.;Fernandez,E.J.;Woolf,T.B.;O?Connell,J. P.Anal.Chem.2005,77,1243.doi:10.1021/ac048812r

(19)Turowski,M.;Yamakawa,N.;Meller,J.;Kimata,K.;Ikegami, T.;Hosoya,K.;Tanaka,N.;Thornton,E.R.J.Am.Chem.Soc 2003,125,13836.doi:10.1021/ja036006g

(20) Zhang,L.;Zhao,G.;Sun,Y.J.Phys.Chem.B 2010,114,2203. doi:10.1021/jp903852c

(21) Dimartino,S.;Boi,C.;Sarti,G.C.J.Chromatogr.A 2011,1218, 1677.doi:10.1016/j.chroma.2010.11.056

(22) Fu,Y.;Chen,L.;Ke,J.;Gao,Y.;Zhang,S.;Li,S.;Chen,T.; Zhao,J.Mol.Simul.2012,38(6),491.doi:10.1080/ 08927022.2011.649427

(23) Ke,J.Y.;Fu,Y.Q.;Chen,L.L.;Chen,T.N.;Zhang,S.J.; Wang,H.B.;Zhao,J.W.J.Fudan Univ.2012,51(2),251. [柯佳穎,傅應(yīng)強,陳莉莉,陳天南,張世界,王洪波,趙健偉.復(fù)旦大學學報,2012,51(2),251.]

(24) van Deemter,J.J.;Zuiderweg,F.J.;Klinkenberg,A.Chem. Eng.Sci.1956,5,271.doi:10.1016/0009-2509(56)80003-1

(25) Schwander,J.;Stauffer,B.;Sigg,A.Ann.Glac.1988,10,141.

(26) McAfee,M.;Lindstr?m,J.;Johansson,W.J.Soil.Sci.1989,40 (4),707.doi:10.1111/ejs.1989.40.issue-4

(27) Saxton,M.J.;Jacobson,K.Annu.Rev.Biophys.Biomol.Struct. 1997,26,373.doi:10.1146/annurev.biophys.26.1.373

(28) Sander,L.C.;Wise,S.A.Anal.Chem.1995,67(18),3284.doi: 10.1021/ac00114a027

(29) Welsch,T.;Michalke,D.J.Chromatogr.A 2003,1000,935.doi: 10.1016/S0021-9673(03)00503-X

(30)Keulemans,A.L.M.;Kwantes,A.Vapour Phase Chromatography;Academic Press:New York,1957;p 15.