基于微流控-SERS的基底銀溶膠團聚動力學分析

丁家堅，楊蘭英，陳桂蘭，趙璐瑤，趙婧琳

（成都理工大學核技術與自動化工程學院，四川 成都 610059）

基于微流控-SERS的基底銀溶膠團聚動力學分析

丁家堅，楊蘭英，陳桂蘭，趙璐瑤，趙婧琳

（成都理工大學核技術與自動化工程學院，四川 成都 610059）

為探討團聚形成的二聚體與其相關參數的關系，該文以2，4-DNT作為仿真目標分子，運用COMSOL Multiphysics對微流道中基底銀溶膠團聚進行數值模擬，將銀溶膠濃度與流體流速參數化，以此研究銀溶膠濃度與流體流速對二聚體最大濃度、二聚體最大濃度位置的影響，得到二聚體最大濃度與銀溶膠濃度、二聚體最大濃度位置與流體流速的關系。該關系表明銀溶膠濃度決定二聚體的最大濃度值，流體流速是決定二聚體最大濃度位置的主要因素，可為基于微流控-SERS的痕量爆炸物探測技術提高拉曼散射光譜強度提供一定參考。

表面增強拉曼光譜；COMSOL；銀溶膠濃度；流體流速；二聚體濃度

0 引言

近年來，表面增強拉曼光譜因為具有分辨率高、溶液干擾性小、穩定性好及適合界面研究等特點，已經廣泛應用于分析化學、材料化學檢測[1]，如麻醉劑檢測及一些生物分子探測[2]。對表面增強拉曼散射（SERS）技術而言，其商業化應用于痕量爆炸物探測的技術瓶頸在于：1）分子的拉曼散射光譜強度遠遠低于其熒光光譜，目前以金、銀溶膠為代表的比較成熟、成本適宜于商業化推廣的增強技術也還存在溶膠增強能力的穩定性問題[3]；2）連續探測時SERS基底的更新問題。

仿生學[4]的出現為科技創新提供了新思路。江雷[5]團隊利用“荷葉”效應研發了大接觸角的疏水表面材料。利用“蛾眼”結構，科學家研發了寬帶抗反射表面[6]。而對于爆炸物探測來說，嗅爆犬一直是軍隊、警方搜索隱藏爆炸物的得力助手，其高靈敏度的探測引起了科研人員的注意。醫學研究表明，嗅爆犬具有特殊的鼻腔構造，使得它們具有特別敏銳的嗅覺，隨即，美國的科研人員提出通過微納加工技術MEMS[7]構建嗅爆犬的鼻腔近似結構，結合SERS技術應用于痕量爆炸物的探測。類“嗅爆犬鼻腔”結構的主體為自由表面微流道，該微流道工作原理為含有痕量目標分子空氣流流經微流道，在微流道自由表面進行氣-液（相）傳質交換。微流道自由表面液體區域具有濃縮器的功用，其最大能提高目標分子濃度達106倍[8]，并結合SERS技術的增強作用，最終能實現爆炸物分子的痕量探測，達到嗅探犬的靈敏度。

自由表面微流道結合SERS應用于痕量爆炸物分子的探測中，研究銀溶膠在微流道中的團聚行為有助于確定作為SERS主要載體的銀溶膠二聚體的特性，即二聚體最大濃度值及位置與銀溶膠濃度、流體流速的關系。根據該關系可實現于拉曼光收集位置出現最大濃度的二聚體，極大提高拉曼散射信號，從而提高痕量探測精度。

1 自由表面微流道結構及其工作原理

自由表面微流道是基于MEMS技術在硅芯片上刻蝕而成，如圖1所示。該微流道長為0.6mm，寬15μm，兩端流道深約為40μm，中間流道較淺約為4μm。銀溶膠溶液從左端入口流入，待測目標分子在c流道上表面經氣液（相）傳質交換進入流道，待測分子進入流道后，在b、c流道中開始團聚反應，最終于b流道中出現最大濃度的二聚體。該自由表面微流道具有很大表面體積比，具有預濃縮的功用，在一定條件下最大提高目標分子濃度達106倍。由于大多數爆炸物分子在常壓中濃度都很低，自由表面微流道的預濃縮效應很大程度的提升了探測精度。

圖1 自由表面微流道模型工作原理圖

2 模型建立

2.1 控制方程、邊界條件及初始條件

常見管道流體主要分為兩種類型：層流、紊流。一般認為，雷諾數Re≤2000，流體保持層流狀態。雷諾數逐漸增大會出現擾動，Re＞4 000時，形成紊流。對于流體流動，采用任意拉格朗日-歐拉（ALE）描述比較方便，某一時刻，在參考系中x位置流動的流體質點速度u的時間導數可以表示為

V·▽——遷移導數，在物理上表示由于流體微團從流場中的一點運動到另一點，流場的空間不均勻性引起的時間變化率。

流體域Ω的不可壓縮流動動量守恒定律以及連續方程可以表示為

式中：u——流體速度；

ρ、μ——流體的密度、動力粘度；

P——壓強；

F——體積力向量；

[0，T]——所考慮的時間間隔。

在微流道內流體邊界Γ處，邊界條件如下：

式（4）為流體壁面的無滑移條件，式（5）為流體入口邊界條件，方向為法向流入。

對于銀溶膠團聚，根據DLVO理論[9]，銀溶膠納米粒子因其表面電層引起的靜電斥力穩定懸浮于流體中。當目標分子進入流體后，將吸附于銀溶膠表面并減少其表面電荷，從而消除靜電排斥引起銀溶膠聚合，銀溶膠團聚反應方程[10]：

式中：k0——最大團聚反應常數；

V0——電勢，恒定；

β——隨溫度變化常數；

Kb——波爾茲曼常數；

mc——流體中目標分子濃度。

目標分子與銀溶膠在流體中對流擴散控制方程為

式中i=a表示目標分子濃度，i=1表示銀溶膠濃度，i=2表示二聚體濃度，i=3表示三聚體濃度。Di表示相應物質的擴散系數，擴散系數求解方程：D=KbT/6πμRs，Rs指相應物質半徑，得到Da=10-10m2/s，Di=1，2，3=1.13×10-11m2/s。反應變量Ri：

各個場的初始值都設為零，但為了計算易于收斂，將一個很小的初始值輸入到模型中。溶膠團聚時間：Ta=1/（k·C1）。團聚時間Ta主要由反應常數k、銀溶膠濃度C1決定。其中k=106m3/（mol·s）、C1= 5×10-6mol/m3，得到Ta=5 s；擴散時間Tdiff=L2/Di，L為流道的有效特征長度；分別計算得到2，4-DNT、銀溶膠顆粒擴散時間為103～104s。

2.2 仿真模型及材料屬性

為了易于收斂，選擇基于微流道的2D界面進行仿真，幾何模型如圖2所示。設置1為銀溶膠溶液入口，2為目標分子入口，3為流體出口。流體介質為等離子水。仿真材料屬性參數如表1所示。

利用自由剖分三角形網格對微流道進行劃分，單元數為15036個。通過參數化掃描設置如表2所示，研究各相關參數對二聚體濃度的影響。研究時間T=2000s，時間步長為10s，相對容差為0.01，求解時間約1h。

表1 與銀溶膠團聚相關的材料屬性參數

圖2 自由表面微流體2D截面示意圖

表2 相關參數的參數化設置

3 結果與分析

3.1 銀溶膠濃度與二聚體濃度的關系

通過對銀溶膠濃度參數化掃描（1.0×10-5，1.0×10-5，4.0×10-5），求解得到其中3種不同的銀溶膠濃度下二聚體濃度分布如圖3所示。

圖4為二聚體的最大濃度、平均濃度與銀溶膠濃度關系圖，圖中二聚體最大濃度值與其平均濃度值都隨銀溶膠濃度的增大而增大。同時，可以看到銀溶膠濃度與二聚體最大濃度存在線性增長關系。通過提取銀溶膠濃度與二聚體最大濃度關系數據，得出兩者的線性關系式：

式中：C2——二聚體最大濃度；

C1——銀溶膠濃度。

顯然，銀溶膠濃度是決定二聚體達到最大濃度的關鍵參數，且通過該結論能夠較為準確地確定二聚體最大濃度值，以實現最優探測。

3.2 流體流速與二聚體濃度的關系

通過對流速參數化掃描（4.0×10-6，1.0×10-6，1.0×10-5），求解得到其中3種不同的流速下二聚體濃度分布如圖5所示。

圖3 3種不同銀溶膠濃度下二聚體分布圖

圖4 二聚體最大濃度、平均濃度曲線

圖5 3種不同流速下二聚體濃度分布圖

圖6 所示為二聚體最大濃度、平均濃度與流速關系圖，圖中二聚體最大濃度隨流速的增大增長緩慢，且增長幅度較小；而二聚體平均濃度隨著流速的增大，增長幅度較大。由此，得出流速對銀溶膠的團聚速率影響很小，但其引起了銀溶膠的強制對流從而使銀溶膠平均濃度在增大的同時逐漸趨于平穩。進一步對不同的流速下最大濃度出現位置數據提取，得到如圖7所示的關系圖，從圖中得到流速與二聚體最大濃度位置同樣也存在著線性關系：x=8.357×103u+ 0.0004。顯然，流速是決定二聚體最大濃度于流道位置的關鍵參數。根據以上關系，可以有效控制二聚體濃度出現在激光探頭照射位置，實現高效探測。

圖6 二聚體最大濃度、平均濃度曲線

圖7 流速與二聚體最大濃度位置關系

4 結束語

本文通過研究銀溶膠濃度、流體流速對銀溶膠團聚運動的影響，得到以下結論：

1）團聚形成的二聚體最大濃度值隨著銀溶膠濃度的增大而增大，且兩者有線性關系：C2=0.1905×C1。

2）流體流速對團聚形成的二聚體最大濃度影響很小，對二聚體平均濃度影響較大。由此，我們確定流體流速不改變銀溶膠的團聚速率，僅引起強制對流使二聚體加速擴散，最終使得微流道中的二聚體濃度趨于平衡。此外，改變流體流速能夠改變二聚體最大濃度于微流道中的位置，且兩者形成線性關系：x=8.357×103u+0.0004。

3）本文通過研究銀溶膠濃度、流體流速對團聚形成的二聚體濃度的影響。得出了銀溶膠濃度決定二聚體的最大濃度值；流體流速決定二聚體最大濃度位置。由此，通過調節兩項關鍵參數，能夠實現二聚體最大濃度較為精準地出現在拉曼信號收集位置，提高拉曼光譜強度，進而實現最優探測。

[1]夏天紅.微流控技術結合先進校正模型用于SERS定量分析[D].長沙：湖南大學，2014.

[2]ANDREOU C，HOONEJANIM R，BRAMIM R，et al. Meinhart rapid detection of drugs of abuse in saliva using surface enhanced raman spectroscopy and microfluidics[J].ACS Nano，2013，7（8）：7157-7164.

[3]孟衛.高SERS活性的納米銀溶膠的制備及其增強機理研究[D].南京：南京理工大學，2014.

[4]孫毅.仿生學研究的若干重要進展[J].科技情報開發與經濟，2010，20（3）：143-144.

[5]江雷.從自然到仿生的超疏水納米界面材料[J].科技導報，2005，23（2）：4-7.

[6]董曉軒，申溯，陳林森.銀鏡反應制備納米蛾眼減反結構法[J].光子學報，2014，43（7）：178-183.

[7]方震華，黃慧鋒.微電子機械系統（MEMS）技術在軍用設備中的應用現狀[J].電子機械工程，2010，26（4）：1-4.

[8]PIOREK B D，LEE S J，MOSKOVITS M，et al.Freesurface microfluidics/SERS for real-time trace vapor detection of explosives anal[J].Chem，2012，84（22）：9700-9705.

[9]DERJAGUIN B V，LANDAU L.The theory of stability of highly charged lyophobic sols and coalescence of highly charged particles in electrolyte solutions[J].Acta Physicochim URSS，1941（14）：633-652.

[10]MOSKOVITS M，VICKOVA B.Adsorbate-Induced Silver Nanoparticle Aggregati on Kinetics[J].Journal of Physical Chemistry B，2005，109（31）：14755-14758.

（編輯：莫婕）

Aggregation kinetics of silver sol based on m icro fluidic-SERS

DING Jiajian，YANG Lanying，CHEN Guilan，ZHAO Luyao，ZHAO Jinglin
（School of Nuclear Technology and Automation Engineering，Chengdu University of Technology，Chengdu 610059，China）

It is of great significance to discuss the influence of the related parameters on the aggregation of silver sol.In this paper,based on 2,4-DNT as simulation target molecules, simulating the aggregation of silver sol in the micro-channel through COMSOL Multiphysics,the influence of fluid velocity and silver sol concentration on the maximum concentration of dimer was studied by parameterizing the flow rate and the concentration of silver sol.The simulation results show that the relationship between the dimmers of maximum concentration and silver sol,the position of dimmers and fluid velocity.The relationship shows the concentration of silver sol determines the maximum concentration of dimmers,and the fluid flow rate is the main factor of determining the maximum concentration of dimmers.The result provides a basis for improving the surface-enhanced Raman signal.

surface-enhancement Raman spectroscopy；COMSOL；silver sol concentration；fluid velocity；dimmers concentration

A

1674-5124（2017）05-0036-04

10.11857/j.issn.1674-5124.2017.05.008

2016-08-05；

2016-11-12

四川省科技計劃項目（2015HH0063）

丁家堅（1991-），男，廣西南寧市人，碩士研究生，專業方向為機械電子工程。