2.2 氣/液界面微觀行為的表征能力

2.2.1 界面結(jié)構(gòu)與組分分布：數(shù)密度分布分析

圖4為不同復(fù)配體系下z軸方向的數(shù)密度分布。由圖4可知，沿z軸的數(shù)密度分布可清晰界定氣相、界面區(qū)與液相的位置及界面厚度。通過該分析能直觀地揭示氟碳鏈、氫碳鏈、水分子及各復(fù)配組分在界面區(qū)域的富集程度、分布輪廓及其與純體系的差異，不同原子/基團(tuán)（如極性頭基與疏水尾鏈）在z軸上的位置偏移間接反映分子的取向偏好。

圖4   不同復(fù)配體系下z軸方向的數(shù)密度分布

（a） PFHxA/SDS/Li+ （b） PFHxA/SDS/Na+ （c） PFHxA/SDS/K+

2.2.2 分子取向與有序度：動(dòng)力學(xué)分析

為定量表征界面分子（特別是氟碳鏈與氫碳鏈）的排列有序性與取向，鍵序參數(shù)（P?）和分子軸序參數(shù)（S）被廣泛應(yīng)用。這些參數(shù)能精確量化鏈的排列有序度、平均傾斜角度（直立、平躺或傾斜）以及頭基（如羧基、磺酸基）插入水相的深度。通過直接對(duì)比發(fā)現(xiàn)，氟碳鏈通常具有比氫碳鏈更高的剛性與排列緊密性。

2.2.3 分子間相互作用空間關(guān)聯(lián)：徑向分布函數(shù)（RDF）分析

RDF可定量揭示界面區(qū)域分子間的特定相互作用。圖5為復(fù)配體系的徑向分布函數(shù)。其中，g（r）表示距離參考粒子r處找到目標(biāo)粒子的相對(duì)概率密度。通過分析特定原子對(duì)（如F…F、F…H（水）、C（氟碳）…C（氫碳）等）的RDF，可探究氟碳鏈間潛在的氟?氟作用、氟碳/氫碳鏈與水分子間的特征水合結(jié)構(gòu)（水合殼層）以及不同組分間的聚集傾向，為闡明復(fù)配體系中鏈間是傾向于分相聚集還是混合分布、是否存在特定分子間協(xié)同作用提供關(guān)鍵證據(jù)。

圖5   復(fù)配體系的徑向分布函數(shù)

2.2.4 微觀構(gòu)象與聚集狀態(tài)：可視化與輔助分析

對(duì)界面微觀結(jié)構(gòu)與形貌的可視化分析（包括軌跡快照、二維密度俯視圖）可提供分子構(gòu)象、排列及聚集狀態(tài)的直觀圖像，輔以界面漲落分析及團(tuán)簇分析，可有效探測(cè)界面在橫向（x?y平面）的均質(zhì)性、微相分離現(xiàn)象（如氟碳/氫碳富集區(qū)）、微小聚集體（如半膠束）的形成以及界面粗糙度。圖6為x?y平面內(nèi)表面活性劑分子的二維數(shù)密度分布圖。

圖6   x?y平面內(nèi)表面活性劑分子的二維數(shù)密度分布圖

（a） PFHxA/SDS/Li+ （b） PFHxA/SDS/Na+ （c） PFHxA/SDS/K+

2.2.5 分子動(dòng)態(tài)行為：動(dòng)力學(xué)分析

通過計(jì)算均方位移（MSD）、轉(zhuǎn)動(dòng)相關(guān)函數(shù)、氫鍵動(dòng)力學(xué)等指標(biāo)，可動(dòng)態(tài)評(píng)估界面分子的流動(dòng)性（平動(dòng)與轉(zhuǎn)動(dòng)擴(kuò)散）、鏈的柔性（擺動(dòng)、扭曲幅度與頻率）、吸附穩(wěn)定性（翻轉(zhuǎn)、脫附難易度）以及界面水分子的結(jié)構(gòu)動(dòng)態(tài)特征。圖7為不同復(fù)配體系中表面活性劑分子的關(guān)鍵動(dòng)態(tài)性質(zhì)隨時(shí)間的演變曲線。

圖7   不同復(fù)配體系中表面活性劑分子的關(guān)鍵動(dòng)態(tài)性質(zhì)隨時(shí)間的演變曲線

2.2.6 宏觀性質(zhì)關(guān)聯(lián)：界面張力計(jì)算

基于Kirkwood?Buff公式或壓力張量計(jì)算模擬體系的界面張力，將微觀模擬結(jié)果與關(guān)鍵宏觀性質(zhì)直接關(guān)聯(lián)。通過對(duì)比純組分模擬值、理論加和值與復(fù)配體系的模擬值，可定量驗(yàn)證復(fù)配體系降低界面張力的協(xié)同效應(yīng)及其程度。

2.2.7 吸附熱力學(xué)與協(xié)同機(jī)制：自由能計(jì)算

利用傘形采樣（Umbrella Sampling）等自由能計(jì)算方法，可定量評(píng)估不同組分（氟碳/氫碳鏈分子或基團(tuán)）在界面上的吸附自由能，揭示其自發(fā)吸附界面的熱力學(xué)驅(qū)動(dòng)力強(qiáng)弱。該方法能深入探究復(fù)配體系中一種組分的吸附如何影響另一種組分的吸附自由能，從而闡明協(xié)同作用（或拮抗作用）產(chǎn)生的熱力學(xué)微觀機(jī)制。

2.3 多組分體系協(xié)同機(jī)制的解析

MD模擬在原子尺度解析短鏈氟碳復(fù)配體系（如陰/陽離子型、氟碳/碳?xì)浠旌闲停┑亩嘟M分協(xié)同機(jī)制方面具有獨(dú)特優(yōu)勢(shì)，能定量揭示界面分子排布、相互作用及微觀?宏觀關(guān)聯(lián)。圖8為不同n（PFB?MC）/n（1?OA）復(fù)配體系的最終MD模擬模型。

圖8   不同n（PFB?MC）/n（1?OA）復(fù)配體系的最終MD模擬模型

（a） 1∶0 （b） 2∶1 （c） 1∶1 （d） 1∶2 （e） 0∶1

2.3.1 全原子模型揭示陰/陽離子靜電互補(bǔ)驅(qū)動(dòng)機(jī)制

X.L.QUAN等利用全原子模型（CHARMM力場(chǎng)）結(jié)合多尺度分析（密度分布、RDF、SASA），定量揭示了陰/陽離子表面活性劑在氣?液界面的協(xié)同機(jī)制。結(jié)果表明，靜電互補(bǔ)（-SO/-NH+）驅(qū)動(dòng)n（PFB?MC）/n（1?OA）=1∶1的復(fù)配體系，可形成致密單分子膜，顯著降低表面張力至28.0 mN/m（降幅>30%）；n（PFB?MC）/n（1?OA）=1∶1的復(fù)配體系中抗衡離子實(shí)現(xiàn)了均勻分散，且1?OA具有強(qiáng)氫鍵能力；通過動(dòng)態(tài)可視化技術(shù)可揭示非等物質(zhì)的量比下抗衡離子的競爭吸附行為，且可闡明比例調(diào)控機(jī)制。該研究可為環(huán)保短氟碳滅火劑設(shè)計(jì)提供原子尺度依據(jù)，指導(dǎo)復(fù)配比例優(yōu)化。

A.SHARMA等采用全原子分子動(dòng)力學(xué)模擬（OPLS?AA力場(chǎng)），系統(tǒng)研究了1～4頭陽離子表面活性劑（吡啶鹽類）在水溶液中的自組裝行為，揭示了多頭結(jié)構(gòu)對(duì)聚集行為的調(diào)控機(jī)制。結(jié)果表明，隨頭基數(shù)目增加，靜電排斥增強(qiáng)，導(dǎo)致聚集數(shù)顯著降低，且聚集體形狀偏離球形；4頭表面活性劑因頭基靈活性較高，可通過構(gòu)象調(diào)整減弱靜電排斥，形成較大聚集體（最大聚集數(shù)達(dá)15）。該研究從原子層次明確了陰/陽離子靜電互補(bǔ)是驅(qū)動(dòng)多頭表面活性劑自組裝及界面行為的主導(dǎo)因素，為理性設(shè)計(jì)離子型表面活性劑體系提供了關(guān)鍵依據(jù)。