郭 彥， 傅應(yīng)強， 趙健偉

(1.南京信息工程大學(xué)環(huán)境科學(xué)與工程學(xué)院大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心，江蘇南京 210044；2.南京大學(xué)化學(xué)化工學(xué)院，生命分析化學(xué)國家重點實驗室，江蘇南京 210093)

氣相色譜在揮發(fā)性化學(xué)物質(zhì)的分離、鑒定和定量分析上占有重要的地位，它在環(huán)境、工業(yè)、生物、法醫(yī)等方面也有著突出的應(yīng)用[1 - 3]。氣相填充柱色譜是最早發(fā)展起來的色譜技術(shù)，相關(guān)色譜理論，如經(jīng)典的速率理論、板高方程等也是最早被提出并完善的，在此基礎(chǔ)之上色譜理論被進一步拓展到了開管柱以及液相色譜中?；诤暧^的擴散理論，人們對氣相色譜的保留值、半峰寬，色譜操作條件的優(yōu)化已有大量的實驗研究。在色譜理論的推進和深化過程中，從微觀層面上對它的分離及條件控制的模擬也顯得至關(guān)重要。譬如，Dondi等應(yīng)用蒙特卡羅方法模擬了非線性色譜并優(yōu)化了操作條件，并對Langmuir吸附模型相關(guān)的色譜熱力學(xué)進行了研究[4]。Siepmann等應(yīng)用蒙特卡羅方法研究了氣液色譜中烷烴、醇、芳烴等物質(zhì)的界面吸附、分配、保留機理等[5 - 7]。Tanaka等應(yīng)用實驗結(jié)合半經(jīng)驗量化計算的方法研究了氫/氘同位素在疏水作用上的影響[8]。Sun等運用分子動力學(xué)模擬了配體密度對蛋白質(zhì)在疏水電荷誘導(dǎo)色譜上吸附、脫附和構(gòu)型變化的影響[9]。Gorbunov等應(yīng)用自編程的色譜二維可視模型研究了不同吸附作用下嵌段共聚物在液相色譜中的分配系數(shù)[10]。Sarti等應(yīng)用數(shù)學(xué)模型描述了蛋白質(zhì)在膜親和色譜中的純化過程[11]。這些理論研究為色譜實驗提供了豐富的基礎(chǔ)理論，然而，由于計算量的限制，這些理論研究主要集中在色譜熱力學(xué)上，如吸附模型、兩相分配等。

色譜在分離科學(xué)中占據(jù)不可或缺的地位也和它的動力學(xué)內(nèi)容，如速率理論、流出曲線、待分離粒子的遷移過程等密不可分。色譜動力學(xué)不僅能促進提升色譜熱力學(xué)研究的理論層次，而且能夠?qū)ιV分離的本質(zhì)，操作條件的優(yōu)化等起到關(guān)鍵性的指導(dǎo)作用。在理論層面上研究色譜動力學(xué)，就需要對整個色譜體系構(gòu)建模型，繁復(fù)且計算量大。而隨機行走特別適合用來模擬這種大規(guī)模體系。隨機行走是Pearson在1905年引入的一項數(shù)學(xué)模型，用于描述在給定時間內(nèi)連續(xù)且隨機的運動軌跡[12]。隨機行走最顯著的優(yōu)勢在于它提取了運動行為中的關(guān)鍵性因素，并對其采取了適當且合理的近似，因而簡化了計算量，使得隨機行走能夠運用在一些大規(guī)模的模擬體系中，如土壤和生物組織中的擴散問題[13,14]、藥物輸運[15,16]、神經(jīng)傳導(dǎo)[17]和晶體生長[18,19]等。

在氣相色譜中，待分離粒子是限制在有限的空間環(huán)境內(nèi)的，粒子與流動相間沒有相互作用，它的運動行為跟隨機擴散密切相關(guān)?；谶@些相似性，我們應(yīng)用隨機行走模擬了氣相填充柱色譜中待分離粒子的運動行為。同時，還開展了相對應(yīng)的實際氣相色譜實驗。模擬與實驗得出了相一致的結(jié)果，驗證了本隨機行走模型在氣相色譜模擬上的可應(yīng)用性。本模擬方法對于色譜操作條件的優(yōu)化，及發(fā)展高效能色譜等具有參考價值。

1 計算方法與實驗

1.1 計算方法

為了研究外部因素對保留因子，柱效等與色譜動力學(xué)相關(guān)參量的影響，我們采用了一種微觀的隨機模型用來描述填充柱色譜中分子運動的過程。所用軟件為自主開發(fā)的隨機行走程序，可在Window系統(tǒng)上運行[20 - 23]。

1.1.1分離柱分離柱定義為一個二維平面空間。根據(jù)實際實驗中的色譜柱長，定義長度為2 m(X方向)，寬度為2 mm(Y方向)，見圖1。在整個空間內(nèi)，周期性排布著單元格，如圖1中虛框所示。單元格中多邊形為勢壘，用來模擬色譜中的固定相。待分離粒子可以在這些空腔內(nèi)(無勢壘區(qū)域)自由運動，但不能進入這些勢壘中。

1.1.2待分離粒子我們以忽略質(zhì)量和體積的粒子(圖1中的灰點)來模擬待分離粒子。粒子的位置用其所在單元格中的相對于單元格中心點的相對位置來描述：

Parrel=Parabs-Cellabs

(1)

其中，Parrel是粒子的相對位置對應(yīng)的位矢，Parabs是粒子的絕對位置對應(yīng)的位矢，Cellabs是勢壘中心點的絕對位置對應(yīng)的位矢。

圖1 算法模型示意圖。虛框為周期性排布的單元格，格內(nèi)的綠色多邊形代表著固定相?；疑珗A點為待分離粒子，紅色曲線為粒子的運動軌跡。桔色箭頭和藍色箭頭分別代表著橫向速度的方向和擴散速度的方向Fig.1 Schematic diagram of the simulated chromatographic system.The dashed box is the repeated arranged subcells,in which the green polygon represents the stationary phase.The grey dot is the point-like molecule,and the red curve is its’ mobile trajectory.The orange and the blue arrows represent the direction of transverse velocity and diffusion velocity,respectively

粒子的速度由兩部分組成。一個是由驅(qū)動力施加的橫向速度，定義為X方向的橫向速度(Vx，單位為mm·s-1)。該速度可以用來代表實際實驗中載氣的流速。另一個速度是擴散速度，該速度是與溫度相關(guān)的一個變量。在本文中該擴散速度的初始賦值是遵照氣體布朗運動估算得出的。粒子擴散的運動方向隨機生成,在以粒子為中心的0≤α<2π范圍內(nèi)等概率分布，α為方向角。粒子達到勢壘邊界上時，將按照彈性碰撞模型處理，即粒子的擴散速度大小不變，方向按照反射定律進行反射后繼續(xù)在色譜柱中運動。

1.1.3運動過程粒子和勢壘的相互作用可以用多種模型來描述，本模擬中采用彈性碰撞模型簡化了粒子與固定相的作用。假定粒子在某段時間之內(nèi)速度的方向和大小不發(fā)生變化,即速度和時間步(即粒子在程序中的行走時間單位與實際物理單位之間的比值,程序中采用0.005 s·step-1)決定了粒子一步隨機行走的步長。在不受到勢壘阻礙的情況下,根據(jù)速度、每步時間和粒子的當前位置即可確定粒子完成一步行走之后所處的位置。若粒子與勢壘發(fā)生碰撞,它的橫向速度不變，擴散速度大小不變，方向按照反射定律進行反射。時間也劃分為兩部分，分別為粒子運動到碰撞點前、后所需要的時間。碰撞前時間作用于碰撞前速度矢量，碰撞后時間作用與碰撞后速度矢量，據(jù)此可以確定碰撞之后的粒子運動軌跡。粒子若干個時間段內(nèi)位置參數(shù)的時間積分,即可以確定整個時間段內(nèi)粒子運動的軌跡：

(2)

1.2 色譜實驗

高純N2、高純H2、空氣和標準CH4(體積含量為3.2×10-6，剩余為Ar)，均購買于南京紅健氣體有限公司。實驗所使用的氣相色譜儀為Agilent GC-6890N系統(tǒng)，配備氫火焰離子化檢測器(FID)。填充柱為不銹鋼材質(zhì)，尺寸為2 m(X方向)×2 mm(Y方向)，內(nèi)填充固定相13XMS(載體粒徑60～80目)。爐溫和檢測器溫度分別設(shè)置為60 ℃和200 ℃。H2流速為30 mL·min-1，空氣流速為400 mL·min-1。載氣為N2，流速設(shè)為30 mL·min-1，在此載氣流速下樣品CH4氣體的保留時間為2.05 min。每次實驗以不分流的方式進樣，進樣體積均為10 mL。為了考察溫度和壓力的影響，實驗中的爐溫及載氣流速會進行相應(yīng)調(diào)整。

2 結(jié)果與討論

2.1 粒子的擴散行為

為了評價模擬中粒子的擴散行為，分析了粒子的均方位移(MSD)。均方位移是指粒子位移平方的平均值，用來表征粒子運動性質(zhì)的微觀物理量，同時它也是一個受溫度影響的物理量。在模擬程序中，均方位移定義為[24]：

(3)

其中，N代表粒子數(shù)目，t為時間，ri(t)-ri(0)為在t時間內(nèi)的位移矢量。

如果以時間為橫坐標，均方位移的斜率即為擴散系數(shù)D：

〈r2(t)〉=2dDt

(4)

其中，d是多維體系的維數(shù)。據(jù)此可以推導(dǎo)出在二維體系中擴散系數(shù)的表達式：

(5)

圖2 不同溫度條件下粒子的均方位移Fig.2 Mean square displacement under different temperature conditions

圖3 不同橫向速度條件下保留時間的柱狀統(tǒng)計分布圖Fig.3 Statistical distribution of the retention times with different transverse rates The particle size of stationary phase is 205 μm;temperature is 200 ℃.

圖2為在不同模擬溫度下，均方位移隨時間的變化曲線。由圖可見，當溫度比較低的時候，粒子的均方位移隨時間呈線性變化，說明粒子的自發(fā)運動是一個自由擴散的過程。從斜率上可以求出粒子的擴散系數(shù)。如，在溫度20、40、60 ℃(曲線g、h和i)時，擴散系數(shù)分別為2.72×10-3、3.12×10-3和3.52×10-3mm2·s-1。這些數(shù)值與文獻中報道的通過實驗方法求解出的氣體的擴散系數(shù)非常接近，均在同一個數(shù)量級內(nèi)[25,26]。在模擬的溫度區(qū)間內(nèi)，當溫度升高，均方位移隨時間呈現(xiàn)曲線變化，說明除了擴散行為外，粒子還有具有明顯的流動行為[27]。同時，曲線的曲率隨著溫度的升高而變大，表明擴散系數(shù)與溫度之間應(yīng)呈指數(shù)函數(shù)的變化關(guān)系。

2.2 壓力影響

氣相色譜對于壓力變化比液相色譜更加的敏銳。在實際色譜實驗中，如果壓力增大，死時間和CH4氣體的出峰時間會縮短。載氣的流速由柱長和死時間計算得出。相應(yīng)的，載氣的流速隨著壓力的增大而增大。因而在模擬中，采用橫向速度Vx作為壓力大小的度量。Vx的取值設(shè)定是在實際實驗結(jié)果的基礎(chǔ)上，略有擴展而來。

圖3為隨著壓力變化，粒子保留時間的柱狀統(tǒng)計分布圖。其中，溫度固定為200 ℃，勢壘粒徑為205 μm。勢壘粒徑的設(shè)置也是根據(jù)實際實驗中固定相載體的大小(60～80目)而定的。固定相填充率和形狀等對粒子運動行為的影響可以參看我們以前的報道[22]。圖3中，每個Vx條件下的柱狀圖都分別為1 000個粒子運動的統(tǒng)計結(jié)果。由圖可見，這些柱狀分布峰符合高斯正態(tài)分布，峰形光滑，無缺陷，且對稱性好。壓力增大勢必引起橫向速度的增加，因而在保留時間的統(tǒng)計圖上，保留時間減小與預(yù)期結(jié)果相符，出峰時間從319 s縮短至126 s。值得注意的是，壓力增大還引起了統(tǒng)計峰峰形的變化，峰形逐漸變窄變尖銳，說明柱效也相應(yīng)增加。半峰寬(FWHM)是色譜流出曲線的一項重要參數(shù)，涉及到色譜動力學(xué)和組分分離方程。在其它條件不變的情況下，半峰寬的大小與柱效的高低成反比[28]。

圖4 平均保留時間(a)和半峰寬(b)隨橫向速度的變化曲線；插圖分別為各自相對應(yīng)的實驗結(jié)果Fig.4 The plots of the average retention time (a) and the full width (b) at half maximum versus the transverse rates;Insets are corresponding experimental results,respectively

圖4為平均保留時間和半峰寬隨Vx變化的曲線圖，插圖分別為相對應(yīng)的色譜實驗的結(jié)果。在實際色譜實驗中，N2的壓力只能在一定范圍內(nèi)調(diào)節(jié)，流速過大或過小都會使得檢測器熄火。圖4顯示平均保留時間隨著Vx的增大而指數(shù)衰減。壓力增大，粒子受到的驅(qū)動力增加，通過色譜柱的時間也就越短。實驗中當Vx從1.8×104增大到2.6×104μm·s-1時，出峰時間從175 s縮短至114 s。在相同Vx變化區(qū)間內(nèi)，模擬結(jié)果從141 s縮短至127 s。模擬結(jié)果與實驗結(jié)果基本一致。模擬中隨著壓力的增大，半峰寬先有個陡降的過程，再緩慢減小。而緩慢減小的區(qū)間與實際實驗結(jié)果也是相吻合的。由于我們所采用的簡化模型，忽略了兩相間的分配平衡，所以模擬得出的半峰寬比實際結(jié)果要窄。根據(jù)速率理論，半峰寬減小的變化趨勢說明無論是在模擬中還是在實驗中，Vx都是設(shè)定在低流速區(qū)。半峰寬的減小對應(yīng)于板高的減小，也就是柱效的增加。從模擬結(jié)合實驗結(jié)果可知，在低流速區(qū)，我們可以通過施加一個適當?shù)母邏簛砜s短分析時間，提高柱效。

圖5 不同溫度條件下保留時間的柱狀統(tǒng)計分布圖Fig.5 Statistical distribution graph of the retention times with different temperature The particle size of stationary phase is 205 μm.The flow rate is 1.67×104 μm·s-1.

2.3 溫度影響

在模擬程序中，溫度的大小是由擴散速度反映，也就是分子的步速。溫度與步速之間為線性關(guān)系，因而通過步速的調(diào)控來實現(xiàn)柱溫的改變。圖5為當固定相粒徑和壓力固定不變時，不同溫度下粒子保留間的統(tǒng)計柱狀圖。峰形依然符合高斯正態(tài)分布。當溫度升高，統(tǒng)計峰逐漸左移，出峰時間縮短。我們在圖中觀察到，隨著溫度的升高，統(tǒng)計峰先變窄變尖銳，溫度繼續(xù)升高，峰形又會展寬。

圖6為平均保留時間和半峰寬隨溫度變化的曲線圖。插圖分別為相對應(yīng)的色譜實驗的結(jié)果。在實際色譜實驗中，為了保證體系的穩(wěn)定性和避免固定液的流失，溫度只能在小范圍內(nèi)調(diào)節(jié)(小于140 ℃)。由圖6可見，在實驗的溫度變化區(qū)間60～140 ℃內(nèi)，保留時間從123 s 縮短至 82 s。相同區(qū)間內(nèi)，模擬結(jié)果從173 s縮短至161 s。無論是模擬還是實驗結(jié)果，保留時間均隨著溫度的升高而指數(shù)級衰減的。對于半峰寬而言，實驗結(jié)果也是單調(diào)遞減的，但模擬結(jié)果給出先減小后增大的變化趨勢。這可能是由于溫度不僅影響著擴散速度，也影響著擴散系數(shù)，擴散系數(shù)對于板高的影響本身就是一個復(fù)雜的變量。在與實際實驗相同的溫度區(qū)間，半峰寬的變化是一致的。半峰寬在模擬結(jié)果中給出了一個最小值，也就是在427 ℃左右可以獲得最高柱效。

圖6 平均保留時間(a)和半峰寬(b)隨溫度的變化曲線;插圖分別為各自相對應(yīng)的實驗結(jié)果Fig.6 The plots of the average retention time (a) and the full width (b) at half maximum versus the temperature;Insets are corresponding experimental results,respectively

在氣相色譜中，溫度是一個相對復(fù)雜的影響因素，它在選擇性、分子擴散及柱效之間構(gòu)建了一個平衡。從模擬結(jié)果可以看出柱溫對于柱效的影響并不是單調(diào)變化的。如果在實際實驗中想提升柱效，可以在低流速區(qū)施加一個適當?shù)母邷亍Ｍ瑫r，模擬結(jié)果給出了最佳柱效的溫度條件。盡管由于實際操作條件的限制，我們不能此溫度下運行色譜柱，但是這對于其它色譜體系操作條件的選擇還是具有一定的參考意義。

3 結(jié)論

待分離粒子的動力學(xué)行為在氣相色譜的研究中至關(guān)重要，它能夠影響色譜動力學(xué)的一些基本內(nèi)容，如van Deemter方程、Giddings方程、分離過程等。本文中運用自主開發(fā)的隨機行走軟件對氣相填充柱色譜中的運動行為展開模擬，并考察了溫度和壓力因素的影響。同時，我們也在實際氣相色譜實驗中考察了這兩項影響因素。從模擬與實際實驗中我們得出了相一致的結(jié)果。由于實際實驗操作條件的限制，模擬結(jié)果給出了一個更寬泛的應(yīng)用區(qū)間。通過隨機行走模擬，我們可以從分子層面對氣相色譜中分子的運動行為給出合理解釋，并能夠?qū)崿F(xiàn)操作條件的優(yōu)化。