十二烷基苯磺酸鈉氣/液界面聚集行為的分子動力學(xué)模擬*

韓 穎，曲廣淼，薛春龍，梁 爽，丁 偉

(東北石油大學(xué)化學(xué)化工學(xué)院 石油與天然氣化工省重點(diǎn)實(shí)驗(yàn)室，黑龍江 大慶 163318)

近年來，表面活性劑在工業(yè)中的應(yīng)用愈加廣泛[1-3]。例如在三次采油中降低油/水界面張力，提高石油的采收率，也可以改善農(nóng)藥使用性能，在紡織業(yè)中做洗凈劑等。采用不同的實(shí)驗(yàn)手段可對表面活性劑分子在界面處的聚集形態(tài)進(jìn)行研究，如熒光探針法、小角X射線衍射法和中子散射法等[4]。

分子模擬方法作為一種重要的手段受到了廣泛的關(guān)注，通過分子模擬研究表面活性劑分子在表面、界面的聚集形態(tài)，尤其適用于研究各種表面、界面的微觀現(xiàn)象，從而解釋表面活性劑微觀作用機(jī)理[5-7]。如Aksornnarong Ritwiset[8]等通過分子動力學(xué)模擬在分子水平上研究了山梨醇酐單硬脂酸酯在氣/液界面聚集的單分子膜的穩(wěn)定性，研究了界面張力、氫鍵動力學(xué)性質(zhì)等，并發(fā)現(xiàn)氫鍵作用單分子層的穩(wěn)定性有很大的影響。Li Zhang[9]等通過極性頭的分布、表面活性劑的有序性等性質(zhì)研究了非離子氟碳表面活性劑在氣/液界面的聚集，并將其與烷烴進(jìn)行了對比，發(fā)現(xiàn)烷烴有希望替代氟碳型表面活性劑。Li[10]等采用分子動力學(xué)模擬研究了非離子表面活性劑在氣/液界面的聚集，通過界面張力對其結(jié)構(gòu)性質(zhì)進(jìn)行研究，并分別研究了非離子表面活性劑親水部分和疏水部分的性質(zhì)。李曉峰[11]等通過分子動力學(xué)模擬的方法研究了十二烷基醇聚氧乙烯醚系列非離子表面活性劑單分子層在氣/液界面的微觀結(jié)構(gòu)，考察了表面活性劑分子結(jié)構(gòu)與性能的關(guān)系，并且通過表面活性劑親水基和憎水基的傾斜角研究表面活性劑的取向，發(fā)現(xiàn)親水基更傾向于氣/液界面平行排列。苑世領(lǐng)[12]等用分子動力學(xué)模擬的方法研究了陰離子表面活性劑十二烷基硫酸鈉(SDS)在氣液界面上單層膜的結(jié)構(gòu)，對比了不同界面濃度下單層膜以及疏水烷烴鏈的排列情況，通過均方位移根研究了不同界面濃度下的水分子動力學(xué)信息。李亞娉[13]等采用分子動力學(xué)模擬研究了十二烷基硫酸鈉和甜菜堿以及復(fù)配體系在氣/液界面和油/水界面的排布行為，分別考察了溫度、無機(jī)鹽以及復(fù)配比例對表面活性劑界面活性以及泡沫穩(wěn)定性的影響，并與實(shí)驗(yàn)進(jìn)行了對比。

上述研究均在分子動力學(xué)模擬方法研究表面活性劑界面聚集行為方面取得了一定的進(jìn)展。作者以十二烷基苯磺酸鈉(SDBS)為研究對象，采用分子動力學(xué)模擬的方法對其氣/液界面聚集行為進(jìn)行研究，通過體系中各組分的密度分布、徑向分布函數(shù)等參數(shù)對其在氣/液界面的聚集行為進(jìn)行表征。

1 模擬體系模型的建立及模擬方法

以十二烷基苯磺酸鈉(SDBS)為研究對象，采用GAMESS(US)[14]量子化學(xué)計算軟件對其進(jìn)行分子結(jié)構(gòu)優(yōu)化，并得到優(yōu)化后聯(lián)合原子結(jié)構(gòu)的電荷分布見圖1，分子模擬所需要的力場參數(shù)文件參照參考文獻(xiàn)[15]生成。采用Gromacs 4.6.7程序，選擇Gromos53a6聯(lián)合原子力場進(jìn)行分子動力學(xué)模擬。首先利用packmol[16-17]構(gòu)建模擬體系，構(gòu)建好的模擬體系見圖2。堆砌的盒子尺寸為5.0 nm×5.0 nm×15.0 nm，采用周期性盒子邊界條件。中間部分為水相，水分子采用SPC/E模型，水相兩側(cè)分別為相同數(shù)量的SDBS分子，SDBS分子的親水基朝向水相，研究中所用模擬體系見表1，模擬過程參照文獻(xiàn)[4]。

圖1 SDBS量子化學(xué)優(yōu)化后聯(lián)合原子分布

中間部分-水相；綠色-碳原子；紅色-氧原子；黃色-硫原子；藍(lán)色-Na+圖2 采用packmol軟件堆砌的模擬體系示意圖

模擬體系SDBS分子數(shù)H2O分子數(shù)離子種類及數(shù)量NA30305000Na+30NA45455000Na+45NA60605000Na+60NA60CA20605000Na+60Ca2+20Cl-40NA60M20605000Na+60Mg2+20Cl-40

2 結(jié)果與討論

2.1 單分子層結(jié)構(gòu)

首先對表1模擬體系分別模擬至10ns時各組分的密度分布進(jìn)行了研究，結(jié)果見圖3，可以看出水相密度在996.98～987.90 kg/m3，接近于密度真實(shí)值，說明模擬條件、力場參數(shù)以及水分子模型等是準(zhǔn)確的，可以用于SDBS分子在氣/液界面聚集行為的分子動力學(xué)模擬。

Z/nma NA30體系中平衡狀態(tài)下各組分的密度剖面圖

Z/nmb NA45體系中平衡狀態(tài)下各組分的密度剖面圖

Z/nmc NA60體系中平衡狀態(tài)下各組分的密度剖面圖圖3 各體系中平衡狀態(tài)下各組分的密度剖面圖

由圖3可見，不同SDBS分子數(shù)的密度剖面圖可知，經(jīng)過長時間的模擬SDBS分子已經(jīng)形成了良好穩(wěn)定的單分子層界面膜。單分子層界面膜厚度通過“10%～90%厚度原則”來定義[18]，NA30體系界面膜厚度為3.04 nm，NA45體系界面膜厚度為3.90 nm，NA60體系界面膜厚度為4.15 nm。上述結(jié)果說明，隨著SDBS分子數(shù)的增加，界面膜厚度明顯增加，這是由于SDBS分子的增多使疏水尾鏈之間的相互作用增強(qiáng)，呈現(xiàn)出更加舒展的形態(tài)，進(jìn)而增加了界面膜厚度。

SDBS分子數(shù)量為60時初始狀態(tài)、模擬至0.1 ns和模擬至10 ns時的界面排布形態(tài)示意圖見圖4,經(jīng)過長時間的模擬，SDBS分子已經(jīng)不像初始狀態(tài)那樣整齊的排列在氣/液界面，而是呈現(xiàn)了彎曲波動的形態(tài)，并且隨著模擬時間的增加，少量水分子滲透進(jìn)入SDBS單分子層中。此外，在模擬之初，Na+均勻分布在水相中，隨著模擬時間的增加，Na+逐漸向界面處遷移，這主要是由于SDBS分子中磺酸基帶負(fù)電，與帶正電的Na+存在靜電作用。

中間部分—水相；綠色—碳原子；紅色—氧原子；黃色—硫原子；藍(lán)色—Na+圖4 NA60體系中SDBS分子不同狀態(tài)下的界面排布示意圖

2.2 SDBS分子與水分子的相互作用

徑向分布函數(shù)g(r)表示一個粒子或基團(tuán)在指定半徑內(nèi)出現(xiàn)另一個粒子或基團(tuán)的概率[19]，通過徑向分布函數(shù)可以研究SDBS分子和水分子以及Na+之間的相互作用。SDBS分子磺酸基中氧原子與Na+和水中氫原子的徑向分布函數(shù)gOS-Na(r)和gOS-Hw(r) 見圖5。

r/nm圖5 磺酸基中氧原子與Na+及水中氫原子的徑向分布函數(shù)(NA60,10 ns)

由圖5可知，gOS-Hw(r)的第一個峰值出現(xiàn)在0.188 nm處，說明磺酸基中氧原子與水中氫原子存在較強(qiáng)的氫鍵作用形成了第一水化層，第二個峰值出現(xiàn)在 0.334 nm處出現(xiàn)，說明隨著距離增加，磺酸基中氧原子與水中氫原子隨著氫鍵作用的逐漸減弱而形成第二水化層，第三個峰值出現(xiàn)在0.418 nm處，此時磺酸基中氧原子與水中氫原子由于較弱的氫鍵作用而形成第三水化層。由此進(jìn)一步說明SDBS分子中親水基磺酸基與水分子存在著較強(qiáng)的相互作用從而使水分子更易于聚集在SDBS磺酸基周圍。在gOS-Na(r)中，第一個峰值出現(xiàn)在0.246 nm處，是由Na+和磺酸基的靜電作用產(chǎn)生的，第二個峰值出現(xiàn)在0.458 nm處，數(shù)值大小明顯大于第一個峰值，是由磺酸基中氧原子優(yōu)先和水中氫原子相互作用而形成氫鍵的緣故，也表明少部分Na+已經(jīng)進(jìn)入了水化層，第三個峰值出現(xiàn)在0.634 nm處，是由于磺酸基和Na+較弱的靜電作用而形成的。通過對比gOS-Na(r)和gOS-Hw(r)可以發(fā)現(xiàn)，gOS-Na(r)的峰值總和大于gOS-Hw(r)的峰值總和，說明磺酸基和Na+之間的靜電作用強(qiáng)于磺酸基中氧原子和水中氫原子的氫鍵作用。

2.3 SDBS分子與金屬離子間的相互作用

采用各組分密度分布和徑向分布函數(shù)分別考察了金屬離子Ca2+和Na+對SDBS分子在氣/液界面聚集行為的影響，見圖6。

Z/nm

Z/nm圖6 NA60CA20體系中SDBS分子與Ca2+和Na+的密度剖面圖(10 ns)

圖6為NA60CA20體系添加Ca2+后SDBS分子、Na+和Ca2+密度分布，體系中加入Ca2+后，Ca2+會優(yōu)先聚集到磺酸基周圍，從而迫使更多的Na+進(jìn)入水相，上述結(jié)果可以說明SDBS分子中磺酸基與Ca2+之間的作用優(yōu)于與Na+間的相互作用。

模擬體系中存在Na+、Ca2+和Mg2+的條件下，SDBS分子磺酸基中氧原子與水中氫原子的徑向分布函數(shù)對比見圖7，NA60體系的g(r)在0.188 nm處出現(xiàn)第一個峰值，表明磺酸基中氧原子與水中氫原子形成氫鍵，作用比較強(qiáng)，在0.330 nm處出現(xiàn)第二個峰值，表明磺酸基中氧原子與水中另外一個氫原子形成氫鍵。兩個峰值間距為0.142 nm，與水分子中兩個氫原子之間間距0.140 nm基本符合，說明磺酸基與水分子中的兩個氫原子分別作用形成氫鍵，在超過0.330 nm后出現(xiàn)的小的峰值是與體相中其它水分子作用減弱形成數(shù)量較少的氫鍵。

r/nm圖7 SDBS分子中磺酸基中氧原子與水中氫原子的徑向分布函數(shù)

三個體系峰值高低的排序?yàn)镹A60體系、NA60CA20體系、NA60M20體系，表明Ca2+和Mg2+的加入使磺酸基中氧原子和水中氫原子之間的作用減弱，并且Mg2+的作用大于Ca2+的作用。通過Gromacs中g(shù)_hbond功能可以計算模擬體系SDBS分子中磺酸基中氧原子與水分子間氫鍵的數(shù)量，結(jié)果見表2，加入Ca2+和Mg2+后的體系中氫鍵數(shù)量有所減少，上述結(jié)果均可說明表面活性劑SDBS耐二價離子能力較差。

表2 不同離子存在下氫鍵數(shù)量

3 結(jié) 論

采用分子動力學(xué)模擬對十二烷基苯磺酸鈉在氣/液界面的聚集行為進(jìn)行了研究。得到如下結(jié)論。

(1) SDBS分子可以在氣/液界面形成緊致的單分子層界面膜，降低表面張力；通過對比不同SDBS分子數(shù)的模擬體系，發(fā)現(xiàn)隨著SDBS分子數(shù)的增加，界面膜的厚度增加；

(2) 由于Na+與SDBS分子中的磺酸基之間存在靜電吸引力，體系中Na+主要分布在SDBS單分子層處；

(3) 通過對磺酸基中氧原子和水中氫原子的徑向分布函數(shù)的研究，發(fā)現(xiàn)磺酸基和Na+的靜電作用要強(qiáng)于磺酸基中氧原子和水中氫原子的氫鍵作用；

(4) 通過對不同離子存在下模擬體系中各組分密度分布研究可知，SDBS分子中磺酸基與Ca2+之間的作用優(yōu)于與Na+間的相互作用；結(jié)合不同離子存在下SDBS分子中磺酸基中氧原子與水中氫原子的徑向分布函數(shù)和氫鍵數(shù)量變化情況，發(fā)現(xiàn)Ca2+和Mg2+的加入使磺酸基中氧原子和水中氫原子之間的作用減弱，說明SDBS分子耐二價離子能力較差。

參 考 文 獻(xiàn)：

[1] 張昌輝,謝瑜,徐旋.表面活性劑在紡織工業(yè)中的應(yīng)用及發(fā)展[J].日用化學(xué)品科學(xué),2008(01):19-23.

[2] 趙成英,王利民.表面活性劑在制藥工業(yè)中的應(yīng)用[J].日用化學(xué)工業(yè),2007(06):389-392.

[3] MUHAMMAD,SHAHZAD,KAMAL.A review of gemini surfactants:potential application in enhanced oil recovery[J].Journal of Surfactants and Detergents,2016(2):223-236.

[4] 王峰,艾池,曲廣淼.十二烷基磺丙基甜菜堿表面活性劑的氣液界面聚集行為分子動力學(xué)模擬[J].計算機(jī)與應(yīng)用化學(xué),2014,31(1):83-89．

[5] 肖紅艷,甄珍,孫煥泉,等.陰離子表面活性劑在水/正烷烴界面的分子動力學(xué)模擬(英文)[J].物理化學(xué)學(xué)報,2010(02):422-428.

[6] 劉國宇,顧大明,丁偉,等.表面活性劑界面吸附行為的分子動力學(xué)模擬[J].石油學(xué)報(石油加工),2011(01):77-84.

[7] HUA Y F, LV W J,ZHAO S L,et al.Effect of surfactant SDS on DMSO transport across water/hexane interface by molecular dynamics simulation[J].Chemical Engineering Science,2015(134):813-822.

[8] AKSORNNARONG RITWISET,SRIPRAJAK KRONGSUK,JEFFREY JOHNS.Molecular structure and stability of the sorbitan monostearate (Span60)monolayers film at the water-air interface:a molecular dynamics simulation study[J].Journal of Molecular Liquids,2014(195):157-164.

[9] ZHANG L,LIU Z P,REN T,et al.Understanding the structure of hydrophobic surfactants at the air/water interface from molecular level[J].LANGMUIR,2014(46):13815-13822.

[10] LI,XF,LI,et al.Molecular dynamics simulations of nonionic surfactant at the air/water interface[J].ACTA CHIMICA SINICA,2011(19):2235-2240.

[11] 李曉峰,李應(yīng)成,吳志勇，等.非離子表面活性劑在氣液界面的分子動力學(xué)模擬[J].化學(xué)學(xué)報,2011,19:2235-2240.

[12] 苑世領(lǐng),崔鵬,徐桂英，等.氣液界面上陰離子表面活性劑單層膜的分子動力學(xué)模擬[J].化學(xué)學(xué)報,2006,64(16):1659-1664．

[13] 李亞娉,呂韋欽,曹緒龍,等.十二烷基硫酸鈉與甜菜堿在氣液和油水界面的復(fù)配協(xié)同作用研究[J].化學(xué)學(xué)報,2014,05:615-623.

[14] PERRI M J,WEBER S H.Web-based job submission interface for the GAMESS computational chemistry program[J].Journal of Chemical Education,2014,91(12):2206-2208.

[15] 丁偉,高翔.表面活性劑GROMAOS53a6 力場參數(shù)文件的構(gòu)建[J].精細(xì)石油化工進(jìn)展,2012,14(3):1-4．

[16] MARTíNEZ J M,ANDRADE R,BALDRIDGE E G,et al.PACKMOL:a package for building initial configurations for molecular dynamics simulations[J].Journal of Computational Chemistry,2009,30(13):2157-2164．

[17] MARTíNEZ L,ANDRADER R,BIRGINE E G.Packmol:a package for building initial configurations for molecular dynamics simulations[J].Journal of Computational Chemistry,2009,30(13):2157-2164.

[18] JANG S S,LIN S T,MAITI P K,et al.A molecular dynamics study of a surfactant-mediated decane-water interface:effect of molecular architecture of alky benzene sulfonate[J].J Phys Chem B,2004,108(32):12130-12140.

[19] 曲廣淼,薛春龍,程杰成,等.碳酸鈉對2,3-二甲基-5-(5'-十六烷基)苯磺酸鈉在油水界面聚集行為的影響[J].精細(xì)化工,2016,33(3):283-288.