納米顆粒間相互作用對油-水界面張力的影響機制

文章編號：1671-3559（2024）02-0219-08DOI：10.13349/j.cnki.jdxbn.20240009.001

摘要： 為了研究油-水界面上納米顆粒的動態吸附過程及其對界面張力的影響，采用耗散粒子動力學模擬方法，建立納米顆粒在油-水界面物理模型，研究單顆粒的吸附動力學過程及多顆粒相互作用對界面張力的影響機制。結果表明：單顆粒在油-水界面的吸附分為自由擴散、 界面吸附、 動態平衡3個階段；單顆粒吸附過程自由能變化遠大于顆粒的動能，顆粒吸附可自發、 快速進行，且吸附后能穩定在界面上；多顆粒間的相互作用力隨顆粒間距離的增大而振蕩衰減，這是由顆粒間的溶劑粒子所產生的溶劑化效應所致；當顆粒間相互作用為引力時，界面張力增大，當顆粒間相互作用為斥力時，界面張力減小。

關鍵詞： 油-水界面； 納米顆粒； 界面張力； 吸附； 自由能

中圖分類號： O793

文獻標志碼： A

開放科學識別碼（OSID碼）：

Effect Mechanism of Nanoparticle Interactions on

Oil-Water Interfacial Tension

LIAN Yanfei， LI Yuxiu， CHEN Ying， ZHENG Jiajie， ZHENG Danjing

（School of Materials and Energy， Guangdong University of Technology， Guangzhou 510006， Guangdong， China）

Abstract： To explore the dynamic adsorption process of nanoparticles at oil-water interface on the interfacial tension， a physical model of nanoparticle adsorption at oil-water interface was established by using the dissipative particle dynamics simulation method， and the adsorption kinetic process of a single particle and the mechanism of multi-particle interaction on the interfacial tension were investigated. The results show that the adsorption of a single particle is divided into three stages， namely the free diffusion stage， the interfacial adsorption stage， and the dynamic equilibrium stage. The free energy change of particle adsorption is calculated by using umbrella sampling， which reveals that particle adsorption is spontaneous and rapid， and the particle can be stabilized at the interface after adsorption. The interaction force between the particles decreases with the increase of the distance between particles， which is due to the solvation effect of the solvent particles between the particles. The interfacial tension increases when the interaction of particles is attractive， and the interfacial tension conversely decreases when the interaction of particles is" repulsive.

Keywords： oil-water interface; nanoparticle; interfacial tension; adsorption; free energy

收稿日期： 2022-12-31""""""""" 網絡首發時間：2024-01-09T18：13：26

基金項目： 國家自然科學基金廣東省聯合基金項目（U20A20299）

第一作者簡介： 練燕菲（1997—），女，廣東清遠人。碩士研究生，研究方向為納米顆粒油-水界面自組裝。E-mail： 1123640828@qq.com。

通信作者簡介： 陳穎（1969—），女，重慶人。教授，博士生導師，研究方向為微尺度流動及傳熱。E-mail： chengying@gdut.edu.cn。

網絡首發地址： https：//link.cnki.net/urlid/37.1378.N.20240109.1405.002

作為典型的液-液分散體系，乳液在重油運輸［1］、 藥物輸送［2］、 納米材料制造［3］等方面都具有重要的作用。因為存在較大的界面面積，所以乳液是一種熱力學不穩定的體系。為了形成穩定的乳液，通常將表面活性劑添加到分散體系中，目的是減小2種互不相容液體間的界面張力。研究［4］發現，與傳統的表面活性劑相比，納米顆粒在穩定乳液方面具有很多優勢，例如成本低、 穩定性好、 毒性小等，因此越來越多研究者利用納米顆粒替代傳統表面活性劑分子來穩定乳液。

關于界面上的納米顆粒是否一定能減小界面張力的話題一直存在爭議［5-6］。大多數研究者認為，界面張力會隨著吸附在油-水界面的納米顆粒數量增多而減小。Zhou等［7］研究了功能性二氧化硅納米顆粒在油-水界面的吸附特性和乳化液的穩定性，實驗結果表明，隨著吸附在油-水界面的納米顆粒數量增加，界面張力逐漸減小。Sofia等［8］采用垂滴法研究了納米二氧化硅對油-鹽水界面張力的影響，發現納米二氧化硅能夠減小油-鹽水界面張力，原因在于納米粒子在油-水界面處的自組裝。隨著納米顆粒數目的增多，也可能使界面張力呈現先減小再增大、 甚至大于本征界面張力的變化規律。Dong等［9］研究了高pH條件下電荷穩定二氧化鈦分散體的界面張力，結果發現，隨著系統中顆粒數量的增加，界面張力先減小后增大，直至大于本征界面張力，可能是由界面處納米顆粒之間的毛細相互作用所致。Yang等［10］利用高速數碼相機研究了二氧化硅納米粒子在電場作用下對液滴界面聚集的影響，結果表明，隨著顆粒數量的增加，界面張力先減小后增大。雖然界面張力的減小促進了液滴下落過程的變形；但是，由于顆粒數量繼續增加，形成的顆粒膜抑制了液滴的變形，因此導致液滴聚結，起到靜電脫水的作用。還有一些研究者認為，納米顆粒在油-水界面的吸附不會對界面張力產生很大的影響，原因在于：一方面，納米顆粒表面活性有限，即使吸附在界面上也不會改變油-水間的相互作用；另一方面，納米顆粒之間的間距比納米顆粒的粒徑還大，顆粒間的相互作用對界面張力的影響很小。Moghadam等［11］研究氧化鋅納米顆粒在有、 無表面活性劑作用下對界面張力的影響，結果表明，當沒有表面活性劑時，納米顆粒可以吸附在界面上，但顆粒表面活性不足以改變界面上的油與水的相互作用，因此對界面張力沒有影響。Luu等［12］利用耗散粒子動力學模擬方法研究了納米顆粒在油-水界面上的動力學特性，結果表明，只有當界面上的納米顆粒覆蓋率足夠大且排斥力足夠大時，界面張力才會顯著減小。

納米顆粒能夠自發吸附在界面的條件是其在界面處的自由能最小，即顆粒從體相移動到界面處能量最低的位置，處于熱力學穩定狀態。吸附在界面上的顆粒對界面張力的影響主要有2個方面：一方面，界面處顆粒-油和顆粒-水的相互作用取代了原來的油-水相互作用，使得界面張力減小；另一方面，界面處顆粒間的相互作用也會對界面張力產生影響。本文中采用耗散粒子動力學模擬方法，建立納米顆粒在油-水界面物理模型研究納米顆粒的動態吸附的過程，以及界面處顆粒間的相互作用對界面張力的影響。

1" 模型建立與驗證

耗散粒子動力學（dissipative particle dynamic，DPD）［13］是采用球形粒子代表系統中的原子團， 以單個粒子作為基本計算單元求解系統中所有粒子運動軌跡的數值模擬技術， 計算的出發點是求解牛頓運動方程， 其中每個粒子所受到的合力fi的表達式為

fi=∑j≠i（FCij+FDij+FRij），（1）

FCij=aij（1-rij/rc）r^ij，（2）

FDij=-λ（1-rij/rc）2（r^ijvij）r^ij，（3）

FRij=σ（1-rij/rc）ξijr^ij ，（4）

式中： FCij、 FDij、 FRij分別為與粒子i、 j距離呈線性關系的保守力、 耗散力、 隨機力； aij為粒子i、 j之間的排斥參數； rij為粒子i、 j間的距離； rc為粒子相互作用的截斷半徑； r^ij為rij方向上的單位向量； λ為與流體黏度相關的模擬參數，取為常用值4.5； v→ij為粒子i、 j的相對速度； σ為噪聲振幅； ξij是均值為0、方差為1的高斯隨機變量。 σ和λ滿足漲落耗散關系σ2=2kBTλ， kB為玻爾茲曼常數， T為熱力學溫度，模擬溫度設置為298 K。DPD模擬中常采用約化單位， 模擬中選取kBT作為能量約化單位， rc作為長度約化單位， τ作為時間約化單位， 時間間隔為Δt=0.001τ， 將DPD模擬中水的擴散系數與實驗中的自擴散系數擬合， 得到τ=7.6 ps［14］。

以癸烷和水分別作為油相和水相建立初始模型。 根據癸烷和水的物性參數， 選取粗粒化程度Nm=5， 即一個DPD粒子代表5個水分子，每個粒子的體積Vb=0.15 nm3。

定義邊長為rc的立方體中的粒子個數為

粒子數密度ρ，取ρ=3，由rc=3ρVb得到rc=0.766 nm。一個癸烷分子的體積約為0.323 nm3， 由于一個系統所有DPD粒子的體積相同，因此一個癸烷分子需要用2個DPD粒子表示，粒子之間用彈性系數為350kBT/rc的簡諧鍵連接，平衡鍵長為r0=0.72rc。 模型中的納米顆粒為半徑為2rc的空心球，如圖1所示。 納米顆粒表面有162個小顆粒且中心有1個小顆粒。 圖2為初始的油-水界面物理模型，模擬域為底面邊長為20rc、 高度為40rc的正方體，其中紅色粒子為油相， 藍色粒子為水相，油相與水相的體積比為1∶1。

為了驗證所建模型的準確性，利用Irving-Kirkwood方法計算DPD單位下的油-水界面張力γ1，然后再轉換為標準單位（mN/m）的界面張力γ。

γ1=12［pzz-（pxx+pyy）/2］Lz，（5）

γ=（kBT/r2c）γ1，（6）

式中： pzz為法向壓力； pxx、 pyy為切向壓力； Lz為模擬域高度；系數1/2表示系統中存在2個油-水界面。

系統中水相、 油相、 納米顆粒相相同組分的相互作用參數aww=aoo=app=131.5kBT/rc、 油-水相互作用參數aow=198.5kBT/rc［15］，具體的初始相互作用參數如表1所示。x、 y、 z方向均采用周期性邊界條件。為了使計算體系達到平衡穩定狀態，模擬共運行106步，所有的DPD模擬均通過大規模原子分子并行模擬（LAMMPS）軟件包完成。模擬結束后采用式（5）、 （6）計算的界面張力為52.7 mN/m，與實驗值51.98 mN/m［16］的誤差僅為1.38%，因此所建立的油-水界面物理模型是可靠的。

2" 結果與討論

2.1" 單個納米顆粒的相界面吸附

體相中的納米顆粒總是時刻受到周圍溶劑分子的熱擾動， 這種熱擾動是顆粒搜索相空間的驅動力。 為了研究納米顆粒的吸附過程， 首先選取顆粒與油相的排斥參數apo=140kBT/rc， 顆粒與水相的排斥參數apw=130kBT/rc， 將單個納米顆粒置于水相并讓其自由擴散， 顆粒吸附過程中質心沿z方向的變化如圖3所示。 從圖中可以看出： 在模擬總時間為91.2 ns內， 顆粒由水相吸附至界面的過程可分為3個階段。 第Ⅰ階段是自由擴散階段，持續時間為54.18 ns， 此時顆粒在水分子的熱擾動下進行布朗運動，運動時間較長，布朗力是顆粒運動的主要動力，顆粒的質心因自由擴散而隨機振蕩；第Ⅱ階段是界面吸附階段，吸附過程持續的時間僅為213 ps，顆粒從體相運動到質心z坐標為16.7rc時開始進入吸附階段， 此時顆粒向界面運動的合力大于布朗力， 促使顆粒吸附至界面； 第Ⅲ階段是顆粒在界面上的動態平衡階段， 顆粒吸附至界面后， 質心z坐標在界面上、 下1.3rc的范圍內波動， 并且不會從界面解吸至體相。

當納米顆粒由體相吸附至界面動態平衡后，可采用接觸角表征顆粒吸附后的平衡狀態以及顆粒的整體疏水程度。接觸角的定義示意圖如圖4所示。由于納米顆粒表面疏水程度在一定程度上決定了顆粒在油-水界面上的吸附，因此可以通過調節顆粒與油相、 水相的排斥參數，進而調控顆粒的吸附行為。顆粒的吸附能也會影響接觸角大小，進而影響顆粒吸附在界面時的張力大小。接觸角θ計算公式為

θ=arccos（h/r），（7）

式中： r為顆粒半徑； h為顆粒質心到界面的距離。

θ—接觸角； h—顆粒質心到界面的距離。

油-水界面的位置可通過計算油、水粒子的密度分布曲線來獲得。首先將模擬域沿z方向平均分成4 000層，每層厚度為0.01rc。在油-水系統平衡后，計算每層區域內油粒子或水粒子的數密度，并將結果進行平均后得到如圖5中的油、 水粒子沿z方向的密度分布曲線，曲線交點即為油-水界面的位置。

rc—粒子相互作用的截斷半徑。

圖5" 油、 水粒子數密度ρ沿z方向的變化

固定納米顆粒與油相之間的排斥參數，調節顆粒與水相之間的排斥參數apw在（100～170）kBT/rc內變化，利用式（7）計算顆粒在界面處的接觸角隨排斥參數apw變化的規律，結果如圖6所示。由圖可見，隨著apw增大，即顆粒疏水性增強，顆粒的接觸角增大。當顆粒與水相以及顆粒與油相的排斥參數相等時即apw=apo=110kBT/rc，接觸角為90°，此時顆粒的親水性和疏水性相當；當apw＜apo時，此時接觸角小于90°，顆粒表現為親水性；當apw＞apo時，此時接觸角大于90°，顆粒表現為疏水性。

顆粒由體相吸附至界面的過程是熱力學自發過程，如果吸附后顆粒的自由能降低，顆粒就有能自發吸附并穩定在界面。通過計算顆粒越過油-水界面過程的自由能變化，可以全面了解吸附過程以及顆粒在界面的穩定性。本文中利用統計熱力學原理及傘形抽樣［17］方法計算納米顆粒由水相運動至油相的自由能變化。假設顆粒沿z方向坐標從14.7rc運動至25.2rc，將上述運動過程每隔0.3rc劃分一個窗口，并在每一個窗口上施加簡諧勢，使得顆粒在該窗口附近作簡諧運動并取得足夠多的采樣點。所添加的簡諧勢為W（x）=ks（x-xK）/2，其中ks為傘形抽樣彈性常數，xK為第K個反應坐標的位置， x為每次抽樣時顆粒的位置。在傘形抽樣結束后，采用加權直方圖分析（WHAM）算法對抽樣結果進行統計。自由能的計算過程中選取ks=500kBT/r2c，每個窗口的采樣次數為104，采樣時間間隔為0.5τ。

圖7所示為納米顆粒從水相運動到油相的自由能變化。 由圖可以看出： 顆粒由水相運動至界面的過程自由能減小了90.32kBT， 而顆粒的動能僅為1.5kBT，說明顆粒的吸附過程完全能自發進行；顆粒由界面運動至油相時，自由能增加了137.69kBT，說明顆粒從界面解吸到油相時需要克服該能量壁壘，該能量壁壘遠遠大于顆粒的動能，所以顆粒發生自發解吸可能性很小。自由能的計算結果也在一定程度上印證了圖3中納米顆粒的運動質心變化，原因在于：一方面，顆粒的吸附能遠大于顆粒本身的動能，因此吸附過程是自發且快速的；另一方面，顆粒具有遠大于動能的解吸能，因此顆粒可以長時間穩定在界面。

2.2" 顆粒間相互作用對界面張力的影響

在界面張力的模擬計算公式（5）中， 界面張力來源于界面附近處法向壓力和切向壓力的差值， pxx、 pyy和pzz是壓力張量p的3個對角分量， 壓力張量p又可以進一步分解為與原子速度和力有關的項［18-19］， 即p=（2EK+EV）/V， 其中EK為與原子速度有關的動能項， V為系統的體積，EV為與原子間相互作用力有關的維里項， EV=∑Ng=1Fgrg，其中Fg為系統中作用于原子g的力， rg為原子g的位置，N為系統中的原子總個數，表示張量積。由于Δp=pzz-0.5（pxx+pyy），也可以寫成Δp（2ΔEK+ΔEV）/V，其中動能的變化量ΔEK可以定義為，其中Δp為壓力的變化量，ΔEV為維里項的變化量，EKzz為z方向的動能項，EKxx為x方向的動能項，EKyy為y方向的動能項。為了觀察動能項對界面張力的貢獻，本文中計算了純油-水系統的動能沿z方向的分布，結果如圖8所示。從圖中可以看出，純油-水系統的動能在界面附近基本不變，即動能的貢獻很小，與文獻［19］中的研究結果一致，因此壓力張量p中與原子相互作用有關的維里項在界面張力中起主要作用。

因為壓力張量中的動能項對壓力影響很小，本文中主要關注與顆粒間相互作用有關的維里項對界面張力的影響。為了研究顆粒間相互作用對界面張力的影響，本文中建立了油-水界面的多顆粒吸附模型。界面張力與顆粒在界面上的覆蓋率有關，界面覆蓋率可定義為=nπd20/4A，其中n為界面吸附顆粒的個數，A為界面面積，d0為界面顆粒與顆粒之間的二維徑向分布函數第一個最高峰橫坐標，表示顆粒與顆粒間的最短距離［20］。圖9所示為3種不同覆蓋率時顆粒-顆粒的二維徑向分布函數。圖中曲線第一個最高峰值的橫坐標為4.8rc，表示界面上顆粒間的最短距離，也被視為顆粒的有效半徑，可通過該有效半徑獲得顆粒的真實覆蓋率。

為了研究界面顆粒覆蓋率對界面張力的影響，建立了底面邊長為30rc、 高度為40rc的正方體模擬區域來容納更多的顆粒，分別將不同數量的顆粒隨機放置于界面上，并采用2種方法計算不同顆粒覆蓋率下的界面張力。

第一種方法是不考慮顆粒間的相互作用，直接利用吸附能計算界面張力。單個球形顆粒從體相吸附至界面后，利用平衡接觸角計算的吸附能ΔE為

ΔE=-πr2γ0（1±cos θ）2 ，（8）

式中： γ0為裸界面的油-水界面張力； 括號里的正、 負號分別表示顆粒質心在油相、 水相。

顆粒吸附后的界面張力γ2可表示為

γ2=E/A=γ0+nΔE/A 。（9）

可以看出，在界面面積一定時，界面張力的變化Δγ2與顆粒數量呈正比，即Δγ2=γ2-γ0=nΔE/A。采用單顆粒的吸附能計算在不同顆粒數量時的界面張力的變化結果如圖10中的紅線所示。

kB—玻爾茲曼常數； T—熱力學溫度； rc—粒子

相互作用的截斷半徑； Δγ1、 Δγ2分別是利用

式（5）、 （9）計算的界面張力變化量。

第二種方法是考慮所有因素對界面張力的影響，利用式（5）計算界面張力的變化Δγ1，如圖10中的黑點所示，黑點中的誤差棒是獨立計算多次結果得到的誤差。隨著顆粒覆蓋率的增加，界面張力降幅增大。Δγ1和Δγ2的變化趨勢是一致的，都是隨著界面覆蓋率的增大，界面張力降幅增大。在覆蓋率較小時，由于顆粒數量較少，在界面分布較為稀疏，使得顆粒間相互作用不明顯，因此顆粒間相互作用對界面張力的計算沒有太大的影響，Δγ1和Δγ2變化趨勢基本一致。隨著覆蓋率的增大，Δγ2gt;Δγ1。由于直接利用吸附能和吸附顆粒數量計算界面張力存在明顯的局限性，因此接下來要進一步研究顆粒間相互作用對界面張力的影響。

本文中首先計算2個納米顆粒間的相互作用。 在油-水界面上模擬放置2個納米顆粒并用彈簧連接起來， 彈簧的彈性常數為kpp=5 000kBT/rc， 初始時刻2個顆粒的距離即為彈簧的平衡長度， 彈簧的平衡長度從4.5rc依次增加至7.5rc， 取樣間隔為0.1rc。在每個初始距離， 2個顆粒經過50 000步平衡后開始取樣，每次取樣共運行6×105步，每隔200步輸出一次，共有3 000個數據點。取樣結束后，利用胡克定律Fpp=kpp（rI-rI，0）計算2個顆粒在不同間隔距離下作用力Fpp，其中rI為第I時刻顆粒的位置， rI，0為初始時刻彈簧的平衡長度。 顆粒的間隔距離與作用力之間的關系如圖11所示。 作用力小于0時2個顆粒的相互作用表現為引力， 作用力大于0時2個顆粒的相互作用表現為斥力。2個顆粒在接近的過程中，需要不斷地一層層排擠出油、 水粒子，隨著顆粒間距增加，顆粒的受力呈現出一個幅值由大到小的振蕩衰減狀態。

隨顆粒距離的變化

以下基于油-水界面處2個納米顆粒之間油、 水粒子的分布說明顆粒間相互作用力波動變化的原因。 如圖12所示， 本文中采用界面上2個納米顆

粒之間局部動態圓柱區域來統計顆粒之間溶劑分粒子即油、 水粒子的數密度， 局部區域為2個納米顆粒的質心C1、 C2連線減去2個顆粒的半徑為母線L、 半徑Rd為顆粒直徑（3/4）D（D為顆粒直徑）即1.5rc的圓柱體積。

在顆粒間隔距離為（4.5～7.5）rc的區間內，每隔0.4rc統計一次局部區域中的油、水粒子數量分布，結果如圖13所示。為了更清晰地比較局部區域內油、水粒子密度，取界面上、 下2rc的矩形區域的粒子數密度作為基準，該矩形區域的粒子數密度為圖14中的紅線所示，約為2.85。當顆粒間隔距離為4.5rc時，顆粒間的局部區域粒子數密度大于顆粒外圍的基準粒子數密度，導致局部區域中的油、 水粒子對顆粒間的布朗撞擊大于顆粒外圍油、 水粒子的撞擊， 所以此時顆粒之間的相互作用力表現為斥力。當顆粒的間隔距離增加到4.9rc時，此時局部區域的粒子數密度小于基準粒子數密度，此時顆粒間的相互作用表現為引力。當間隔距離大于6rc時，顆粒間局部區域的粒子數密度在基準粒子數密度附近微弱波動，顆粒間的相互作用表現為引力與斥力的小幅度波動。

以油-水界面處的顆粒為研究對象，利用維里定理［21］研究顆粒間的相互作用對界面張力的影響。

MA=nkBT+〈12∑ne=1

∑nugt;1reu feu〉，（10）

式中： M為表面壓力， M=Δγ2=γ2-γ0; reu為顆粒e、 u間的距離； feu為顆粒e、 u間的相互作用力。

根據圖11中的顆粒間相互作用力的結果，利用式（10）計算不同間距的2個納米顆粒相互作用對界面張力的影響，結果如圖15所示。由圖可以看出，隨著顆粒間距減小，顆粒間相互作用導致界面張力的變化與顆粒間相互作用力具有相同的趨勢，顆粒間相互作用力表現為斥力時起到減小界面張力的作用，相互作用力為引力時起到增大界面張力的作用。圖15中還給出利用式（5）計算的2個納米顆粒間距不同時的界面張力變化量，該結果考慮了顆粒自身的吸附能和顆粒間相互作用對界面張力的影響。由圖可見： 2個顆粒在任何間距時均能起到減小界面張力的作用。當顆粒間隔距離大于6.5rc時，由顆粒間相互作用引起的界面張力變化量幾乎為0，界面張力降幅取決于顆粒自身的吸附能的影響；當顆粒間隔距離小于6.5rc時，界面張力的變化取決顆粒間相互作用力的大小。當顆粒間距使顆粒間表現為較強的斥力時，界面張力降幅明顯增大；當顆粒間距使顆粒間表現為較強的引力時，界面張力降幅較小；如果顆粒間的引力足夠大，則可能出現界面張力不降反增的現象。

3" 結論

本文中采用耗散粒子動力學模擬方法，研究了納米顆粒在油-水界面處的吸附動力學以及顆粒覆蓋率對界面張力的影響，得到以下主要結論：

1）采用傘形抽樣法計算了單個納米顆粒從水相運動到油相過程中的自由能變化。 顆粒從水相吸附到界面的自由能變化量為90.32kBT， 遠遠大于顆粒的動能， 使顆粒吸附能自發、 快速進行（約為10-10 s）， 如果顆粒要解吸回體相， 需要克服137.69kBT的能壘。 顆粒吸附到界面后的平衡接觸角隨顆粒與水相的排斥參數的增大而增大。

2）在界面隨機放置多個納米顆粒時的界面張力隨顆粒覆蓋率的增大而減小，2個納米顆粒間的相互作用力和二維表面壓力都隨著顆粒間距增加而呈現振蕩衰減的規律。當顆粒間相互作用力為斥力時，顆粒間維里作用與顆粒占據油-水界面面積效應一起促進界面張力的減小；但當顆粒間相互作用為引力時，顆粒間維里作用大于顆粒占據油-水界面面積效應，因此界面張力降幅減小，甚至發生界面張力增大的現象。界面張力的減小有利于改善乳液的穩定性，可應用于藥物運輸、 材料制造等領域，而界面張力的增大則可以應用在原油脫水、 廢水處理等去乳化的領域。

參考文獻：

［1］" 王子琛， 佘躍惠， 翁雪. 納米材料提高原油采收率的機理研究綜述［J］. 當代化工， 2018， 47（12）： 2612.

［2］" ZHOU S B， FAN J， DATTA S S， et al. Thermally switched release from nanoparticle colloidosomes［J］. Advanced Functional Materials， 2013， 23（47）： 5925.

［3］" GRZELCZAK M， VERMANT J， FURST E M， et al. Directed self-assembly of nanoparticles［J］. ACS Nano， 2010， 4（7）： 3591.

［4］" ZHU Y， JIANG J Z， LIU K H， et al. Switchable pickering emulsions stabilized by silica nanoparticles hydrophobized in situ with a conventional cationic surfactant［J］. Langmuir， 2015， 31（11）： 3301.

［5］" REZVANI H， KHALILNEZHAD A， GANJI P， et al. How ZrO2 nanoparticles improve the oil recovery by affecting the interfacial phenomena in the reservoir conditions？［J］. Journal of Molecular Liquids， 2018， 252： 158.

［6］" ZHANG H T， ZHOU M， GUO Z Y， et al. Effect of hydrophobicity on the interfacial rheological behaviors of nanoparticles at decane-water interface［J］. Journal of Molecular Liquids， 2019， 294： 111618.

［7］" ZHOU H D， DAI C L， ZHANG Q S， et al. Interfacial rheology of novel functional silica nanoparticles adsorbed layers at water-oil interface and correlation with Pickering emulsion stability［J］. Journal of Molecular Liquids， 2019， 293： 111500.

［8］" SOFIA S J D， JAMES L A， ZHANG Y H. Understanding the behavior of H+-protected silica nanoparticles at the oil-water interface for enhanced oil recovery（EOR） applications［J］. Journal of Molecular Liquids， 2019， 274： 98-114.

［9］" DONG L C， JOHNSON D. Surface tension of charge-stabilized colloidal suspensions at the water-air interface［J］. Langmuir， 2003， 19（24）： 10205.

［10］" YANG D H， SUN H Y， CHANG Q， et al. Study on the effect of nanoparticle used in nano-fluid flooding on droplet-interface electro-coalescence［J］. Nanomaterials（Basel）， 2021， 11（7）： 1764.

［11］" MOGHADAMTF，AZIZIANS.EffectofZnO nanoparticle and hexadecyltrimethylammoniumbromideonthe dynamic and equilibrium oil-water interfacial tension［J］. Journal of Physical Chemistry B， 2014， 118（6）： 1527.

［12］" LUU X C， YU J， STRIOLO A. Nanoparticles adsorbed at the water-oil interface： coverage and composition effects on structure and diffusion［J］. Langmuir， 2013， 29（24）： 7221.

［13］" 李紅霞， 強洪夫. 耗散粒子動力學模擬方法的發展和應用［J］. 力學進展， 2009， 39（2）： 165-175.

［14］" GROOT R D， WARREN P B. Dissipative particle dynamics： bridging the gap between atomistic and mesoscopic simulation［J］. Journal of Chemical Physics， 1997， 107： 4423-4435.

［15］" FAN H， STRIOLO A. Nanoparticle effects on the water-oil interfacial tension［J］. Physical Review E， 2012， 86（5）： 051610.

［16］" ZEPPIERI S， RODRGUEZ J， DE RAMOS A L L. Interfacial tensionofalkane-watersystems［J］.JournalofChemical amp; Engineering Data， 2001， 46（5）： 1086-1088.

［17］" SINGH H， SHARMA S. Free energy profiles of adsorption of surfactant micelles at metal-water interfaces［J］. Molecular Simulation， 2021， 47 （5）： 420-427.

［18］" WEN B Y， SUN C Z， BAI B F， et al. Ionic hydration-induced evolution of decane-water interfacial tension［J］. Physical Chemistry Chemical Physics， 2017， 19（22）： 14606-14614.

［19］" ZHAO L L， LIN S C， MENDENHALL J D， et al. Molecular dynamics investigation of the various atomic force contributions to the interfacial tension at the supercritical CO2-water interface［J］. Journal of Physical Chemistry B， 2011， 115 （19）： 6076-6087.

［20］" YONG X， QIN S Y， SINGLER T J.Nanoparticle-mediatedevaporationatliquid-vaporinterfaces［J］.ExtremeMechanicsLetters， 2016， 7： 90-103.

［21］" SUN J Z， STIRNER T. Molecular dynamics simulation of the surface pressure of colloidal monolayers［J］. Langmuir， 2001， 17（10）： 3103-3108.

（責任編輯：于海琴）