于先超， 亓文帥， 鄧全花， 侯萬國，3

水-短鏈醇二元溶液的表面吸附

于先超1， 亓文帥1， 鄧全花2， 侯萬國1，3

（1. 山東大學膠體與界面化學教育部重點實驗室， 2. 化學與化工學院， 3. 國家膠體材料工程技術研究中心, 濟南 250100）

基于表面聚集吸附(SAA)模型和Gibbs吸附微分方程, 建立了二元液體混合物Gibbs表面過剩(2)和吸附層厚度()代數(shù)方程, 可在全濃度范圍內預測2和隨體相組成(如組分2的摩爾分數(shù)2,b)的變化. 采用所建立的模型, 研究了水(組分1)-短鏈醇(組分2)二元溶液的表面吸附行為, 其在水表面的吸附趨勢按甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇的順序依次增強, 與其疏水性增強趨勢一致, 且其2和在低2,b區(qū)均依次增大, 而在高2,b區(qū)均依次減小, 是短鏈醇間吸附趨勢和分子尺寸的差異所致. 對于給定的短鏈醇體系, 隨著2,b的增大,2先急劇增大然后降低, 呈現(xiàn)一個最大值, 而其持續(xù)降低.

液體混合物； 表面吸附； 表面聚集； 吸附層厚度； 熱力學模型

液體混合物的表面組成與其體相不同， 稱為表面吸附（或富集）， 這是一個基本的界面現(xiàn)象， 因其具有重要的理論意義和應用價值而受到廣泛關注［1～4］. 表面層組成（或表面過剩量）和吸附層厚度是表征表面吸附的兩個重要物理量. 雖采用多種技術［如橢圓偏振［5，6］、 中子反射［7，8］、 質譜（聯(lián)合表面層分子簇束產生裝置）［9］、 頻振動光譜［10］和亞穩(wěn)態(tài)誘導電子光譜［11］等］探索測定表面層組成和厚度， 但目前仍難以（獨立）獲得準確數(shù)據(jù)［7，8］. 因液體混合物的表面張力易于準確測定， 因而目前相關研究主要集中于建立表面層組成和厚度與表面張力的相關性模型， 依據(jù)表面張力隨體相組成的變化預測表面層組成和厚度的變化， 可以認識和描述表面吸附行為［8，12～15］.

基于熱力學分析及經驗或半經驗歸納等， 已提出了很多模型用于描述或預測表面性質（如組成和厚度）與體相組成的相關性［16～23］. 但由于表面現(xiàn)象的復雜性， 現(xiàn)有模型還缺乏普適性， 緣于模型構建時的一些假設與體系實際狀態(tài)不符［23～27］. 特別是， 對一些性質的預測甚至存在相反的結果. 如， 水-乙醇溶液作為模型體系得到廣泛研究， 關于其表面層厚度隨體相乙醇濃度增大的變化還存在爭議， Azizian等［14］和Hyde等［8］認為逐漸降低（呈線性或指數(shù)衰減）， 而Bagheri等［12，13］和Santos等［15］認為是先升高后降低， 這主要是所用模型及其假設不同所致.

Butler方程、 Gibbs吸附方程和Langmuir等溫式是表面熱力學的3個基本方程［22～28］， 其中Gibbs和Langmuir方程相互獨立， 但均可由Butler方程推出［28］. Butler方程將表面組成、 表面張力與體相組成相關聯(lián)， 理論上可獲得表面組成和厚度信息， 但目前尚未得到廣泛應用， 因其涉及各組分的活度系數(shù)和偏摩爾面積［28，29］， 而這兩個物理量也與體系組成有關， 目前還難以實驗測定和理論預測［30～32］. Gibbs方程可給出表面過剩（吸附）量隨體相組成的變化， 但不能直接給出表面組成信息［28］， 也存在溶質活度系數(shù)問題［25，26］. 另外， Gibbs方程是一個微分方程， 包含表面張力對溶質分子濃度的導數(shù)， 但其處理結果人為誤差大. Langmuir等溫式是一個關聯(lián)表面相和體相組成的簡單方程式［22，24，28］， 僅涉及一個物理意義明確的模型參數(shù)（即吸附平衡常數(shù)）， 但僅適用于組分性質相近的混合體系（或理想溶液， 因其假設活度系數(shù) 為1）. Pi?eiro等［22］考慮體系的非理想性（組分分子尺寸差異及特殊相互作用）， 提出了“擴展Langmuir方程”（Extended Langmuir equation）. 需說明的是， 經典的和改進的Langmuir方程均未涉及表面張力， 在依據(jù)表面張力數(shù)據(jù)推測其吸附平衡常數(shù)時， 需與一個關聯(lián)表面張力和表面組成的獨立方程聯(lián)用［22，23］. Eberhart［18］提出一個半經驗“簡單混合律”（Simple mixing rule）， 經Laaksonen和Kulmala［20］以及Pi?eiro等［22］改進后， 常與擴展Langmuir方程聯(lián)合描述二元液體混合物的表面張力等溫線［22，23］. 表面張力和表面組成等溫線（恒定溫度下， 表面張力和組成隨體相組成的變化線， 可統(tǒng)稱為表面吸附等溫線）可分為L型和S型， 擴展Langmuir方程難以準確描述S型等溫線［23］.

現(xiàn)有預測表面組成和厚度的熱力學模型或方法均基于上述3個基本方程（選取其中的兩個方程聯(lián)用）， 其差異主要是選取的方程組合以及假設（如活度系數(shù)和偏摩爾面積）不同［8，12～15］. 最近， 我們基于熱力學吸附平衡原理和表面膠束化模型［33］， 提出了“表面聚集吸附”模型（Surface aggregation adsorption model， SAA）［23，34］. 活性組分易于表面（界面）聚集， 這已被理論和實驗所證實［20，26，35］. SAA模型與改進的Eberhart混合律聯(lián)用， 可準確描述液體混合物的L型和S型表面吸附等溫線. SAA模型能否與Gibbs方程相關聯(lián)， 合理預測液體混合物的表面吸附層厚度是一個值得研究的問題. 另外， 表面吸附（吸附能、 吸附量和吸附層厚度）的影響因素目前還不清晰， 有必要探索研究.

本文基于SAA模型和Gibbs微分方程， 推導出了Gibbs表面過剩和吸附層厚度代數(shù)方程， 可在全濃度范圍內預測Gibbs表面過剩和吸附層厚度. 采用所建立的模型方程對水-短鏈醇（甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇）二元溶液的表面吸附進行了研究， 考察了其吸附行為的差異， 探究了其內在原因， 以期加深對液體表面吸附現(xiàn)象的認識.

1 理論分析

假設液體混合物的表面層為具一定厚度（）的均勻表面相［18， 23］（這是一個常用假設， 可簡化模型處理； 實際上， 表面層在平行表面方向上是均勻的， 而在垂直表面方向上是不均勻的）. 為了表述方便， 按慣例將表面張力高的組分記為組分1（常稱為溶劑）， 將表面張力低的組分記為組分2（常稱為溶質或表面活性組分， 因其在表面富集）.

1.1　SAA模型

對二元液體混合物體系， SAA模型假設［23］： （1） 表面活性組分在表面層吸附而形成聚集體； （2） 表面層由吸附位（或準晶格）組成， 每個吸附位由一個溶劑分子或一個活性組分聚集體占據(jù).

組分2的表面吸附（或與組分1的交換吸附）平衡可表示為

式中： Site表示吸附位（或準晶格）； C1和C2分別表示組分1和2的分子， aggr表示組分2的聚集體，為其平均聚集數(shù)（一個聚集體的分子數(shù)）， 下標“bulk”表示體相，為組分2與組分1的表面交換吸附平衡常數(shù).

設各物種在體相和表面相的活度系數(shù)均為1， 依據(jù)熱力學吸附平衡原理可知：

式中：1，s和aggr分別為組分1分子和組分2聚集體在表面相中的摩爾分數(shù)；1，b和2，b分別為組分1和組分2分子在體相中的摩爾分數(shù)， 下標“s”和“b”分別表示表面相和體相.

設聚集體的偏摩爾體積（aggr）與組分2的偏摩爾體積（2）關系為aggr=2， 則式（2）可寫為

式中：1，s和2，s分別是組分1和組分2在表面相中的體積分數(shù)；*≡2/1（1為組分1的分子尺寸，2為組分2的分子尺寸）. 因1，s+2，s=1和1，b+2，b=1， 故由式（3）可推得SAA模型方程為［23］

Connors和Wright提出“對比表面壓”（Reduced surface pressure，r）概念［19，27］， 定義為

式中：為液體混合物的表面張力，10和20分別為純組分1和純組分2的表面張力. 用r可方便對比 不同液體混合物體系間表面張力（或吸附）行為的差異（因其值在0～1間變化）.

對二元體系， 經Laaksonen和Kulmala改進的Eberhart混合律方程為［20］

由此式可得，r=2，s， 這是改進Eberhart方程的一個優(yōu)點， 由實驗測定的值可得出表面相組成.

由式（4）～式（6）可得［23］：

式（8）可寫成線性式：

此式表明， lg［r（1?2，b）/（1?r）］對lg2，b作圖為直線， 由其截距和斜率可求K和值. 兩個參數(shù)（K和）均可由線性擬合得到唯一值， 這是該模型的一個優(yōu)點［23］.

1.2　Gibbs表面過剩

表面過剩量（或表面吸附量）是一個相對概念， 表征表面相和體相組成的差異. 表面過剩有不同的定義［13，17］， 常用的是Gibbs表面過剩［12～14］. Guggenheim和Adam［17］揭示了Gibbs表面過剩概念的物理本質： 即單位面積表面層中溶質的量與其等量溶劑在體相時所含溶質的量之差（即基于等量溶劑定義的溶質表面過剩量）. Gibbs基于溶劑（組分1）過剩為零的分界面， 推出溶質（組分2）的表面過剩量［1］， 用2表示.

設溶質的活度系數(shù)為1， 則Gibbs吸附方程為［25，36］

這是一個微分方程， 在獲取d/d2，b值時人為誤差較大［25， 36］.

將式（7）與式（11）聯(lián)合可推導出：

此式即為Gibbs過剩量的代數(shù)方程， 利用SAA模型得到的和值可在全濃度（2，b=0～1）范圍內計算2值.

另外， 由式（12）可看出，2，b→0和2，b→1時的極限2值［分別用2（x→0）和2（x→1）表示］， 即

在2，b→0時，2=0， 而當2，b→1時，2>0， 這由2的物理本意所致［17］. 在全濃度范圍內，2≥0， 溶質（低表面張力組分）為正吸附.

1.3　吸附層厚度

現(xiàn)基于Gibbs方程和SAA模型建立吸附層厚度（Adsorption layer thickness）方程， 簡稱ALT方程. 根據(jù)Gibbs表面過剩的物理本意［17］，2可表示為［12］

由式（3）， 式（4）， 式（12）和式（15）可得ALT方程為

由式（16）可看出，2，b→0時的極限值（τ→0， 設n→0=1）為

在2，b→1時的極限值（τ→1）為

由式（16）～式（18）可以看出，隨（吸附趨勢）的增大而減小；τ→0（或低2，b區(qū)的）隨2的增大而增大（若>>1時，則與無關）， 而τ→1與2無關（與1有關）. 需說明的是， 當=1和=1時， 可視為理想體系， 不發(fā)生吸附， 其→∞（或不存在明確的吸附層； 吸附層不是表面層）.

1.4　模型擬合

采用Origin 2023b（學習版）軟件對實驗數(shù)據(jù)進行模型方程擬合， 采用平均相對偏差（Average relative deviation， ARD）和平均絕對偏差（Average absolute deviation， AAD）評判模型預測與實驗數(shù)據(jù)間的符合度. ARD和AAD計算式為［23，34］

式中：exp和cal分別是物理量的實驗值和模型計算值；為數(shù)據(jù)點數(shù).

2 結果與討論

采用所建立的模型方程研究了短鏈醇甲醇（MeOH）、 乙醇（EtOH）、 異丙醇（PrOH）、 正丙醇（PrOH）和叔丁醇（BuOH）與水（Water）的二元液體混合物的表面吸附行為. 表1列出了所涉及純液體在25 ℃下的摩爾體積（0）、 表面張力（0）以及溶度參數(shù)（）值［25，37～39］. 其中，參數(shù)包括Hildebrand參數(shù)（t）及其Hansen色散（d）、 極性（p）和氫鍵分量（h）， 用于計算組分間的Flory-Huggins參數(shù)（12）， 以考察其相互作用強度（或相似度）［40］. 另外， 在模型計算時， 各組分的偏摩爾體積用其0代替.

Table 1　Molar volumes, surface tensions, and solubility parameters of pure liquids

2.1　表面聚集吸附平衡

圖1（A）和（B）分別為水（1）-短鏈醇（2）混合物的和r隨2，b的變化， 其中數(shù)據(jù)引自文獻［25，37，38］. 可以看出， 各體系均為LI型等溫線［23］， 表明短鏈醇在混合溶液表面發(fā)生正吸附. 甲醇、 乙醇、 異丙醇和正丙醇體系的線依次降低，r線依次增高， 表明其界面活性（或吸附趨勢）依次增強. 叔丁醇和正丙醇體系的等溫線相近， 表明二者界面活性相近.r隨2，b增大先急劇升高后趨于平緩， 甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系其轉折點對應的2，b分別約為0.3， 0.2， 0.15， 0.1和0.1. Raina等［9］采用質譜測定了甲醇、 乙醇和正丙醇水溶液2，s隨2，b的變化， 其等溫線依次增高； 另外， 將2，s-2，b曲線的轉折點視為“吸附飽和”（相應2，s分別約為0.65， 0.65和0.85）， 其對應的2，b分別約為0.3， 0.2和0.1. 本文模型預測結果與之一致.

Fig.1　Change of surface tension(σ)(A) and reduced surface pressure(πr)(B) with bulk composition(x2,b) for water(1)?alcohol(2) mixtures

The symbols represent the experimental data， and the solid lines represent the model predictions.

采用式（9）對水（1）-短鏈醇（2）體系的π-2，b數(shù)據(jù)進行了線性擬合. 作為示例， 圖2（A）和（B）分別 給出了甲醇和叔丁醇體系的擬合結果（其它3個體系的結果見本文支持信息圖S1）， 可見， lg［r（1?2，b）/（1?r）］與lg2，b呈良好的線性關系， 與SAA模型預測一致. 需說明的是， 叔丁醇體系（還有乙醇和正丙醇體系）在高2，b區(qū)出現(xiàn)偏離， 可能是假設各物種活度系數(shù)為1所致［23，34］. 表2列出了各體系所得，K和值， 其=0.89～1.06（即≈ 1）， 表明短鏈醇在表面相無明顯聚集［23］.

Fig.2　Linear fitting plots of Eq.(9) for πr?x2,b data of water?methanol(A) and water?tert?butanol mixtures(B)

The symbols represent the experimental data， and the solid line represents the model fitting.

Table 2　Model parameter values of water-alcohol mixtures

采用式（7）和式（8）及所得和K值計算了各體系的-2，b和r-2，b曲線， 示于圖1. 可以看出， 模 型預測與實驗數(shù)據(jù)間具良好的符合度， 其中甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系計算值的ARD分別約為0.2%， 1.0%， 1.7%， 1.1%和2.6%， 其AAD值分別為0.09， 0.32， 0.60， 0.30和0.80 mN/m. 這表明用式（7）及所得和K（或）值可準確描述-2，b曲線， 因而其和值可用于 式（12）和式（16）計算2和值.

液體混合物表面吸附趨勢應由組分間的相互作用（或相似度）所決定. 二元液體混合物組分間的相互作用可用12表征［40］， 可由下式計算：

各體系的χ12值列于表2中， 可見甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇的χ12值有依次增大趨勢. χ12值越大， 組分間的親和力或相似度越低（即短 鏈醇的疏水性越強）， 越易發(fā)生表面 吸附.

2.2　Gibbs表面過剩

采用所建立的Gibbs代數(shù)方程（式12）計算了各體系2隨2，b的變化， 由圖4（A）可以看出， 各體系呈現(xiàn)相同的變化趨勢. 隨著2，b增大，2先急劇增大后逐漸 降低， 呈現(xiàn)一個最大值（記為2，max）， 與文獻［17，36］報道一致.2，max值在4.77～5.60 μmol/m2間（表3）， 也與文獻［7，36］報道值相近（若考慮活度系數(shù)，2，max值有所增大［7］）. 另外， 繪制了2-lg2，b圖， 以清晰觀察其2，max的位置［圖4（B）］， 可以看出甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系的2，max所對應的2，b值（記為2b，mΓ）依次降低（表3）， 分別約為0.15， 0.05， 0.02， 0.01和0.006， 這與其吸附趨勢依次增強一致. 對比圖1（B）和圖4（B）可看出，2，max時表面吸附未達飽和， 其相應表面相2，s值（記為2s，mΓ）也依次降低（表3）， 分別約為0.36， 0.26， 0.20， 0.20和0.17. Raina等［9］采用質譜測得甲醇在2，b=0.15時的2，s≈0.41， 乙醇在2，b=0.05時的2，s≈0.27， 正丙醇在2，b=0.01時的2，s≈0.2， 本文結果與其基本一致， 表明本模型預測數(shù)據(jù)是合理的. 下面分析影響2，max，2b，mΓ和2s，mΓ的主要因素.

Fig.4　Change of Gibbs excesses(Γ2) with bulk composition(x2,b)(A) and its semi?logarithmic coordinate plots(B) for water?alcohol mixture

Table 3　Adsorption parameter values of water-alcohol mixtures

為了簡便， 取=1， 由式（12）可得：

2為最大值時， d2/d2，b=0， 則

由式（2）和式（23）可得：

可以看出，和2越大，2b，mΓ和2s，mΓ均越小；2，max與2和（10-20）有關，2越大， （10-20）越小， 則2，max越大. 這可解釋短鏈醇間2，max，2b，mΓ和2s，mΓ的變化［式（23）～式（25）計算值與所觀察值可能略有差異， 因各體系的≠1］.

有趣的是， 從圖4（B）還可看出， 甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系的2值在低2，b區(qū)呈依次增大趨勢， 而在高2，b區(qū)與之相反； 當2，b→1時，2（x→1）也呈依次減小趨勢（表3）， 與文獻［36］報道 結果一致. 這些結果均由2值依次增大所致. 在高2，b區(qū)， 由式（14）可知，2增大，2減小. 另外， 由式（22）可知：

這表明在低2，b區(qū)，2值越大， 等溫線極限斜率（d2/d2，b）→0越大， 故2越大.

需要說明的是， 采用多項式=∑aexp（?bx2，b）擬合-2，b曲線， 其中a和b是擬合常數(shù)（令0=0）， 再采用Gibbs微分方程（式11）計算Г2值， 其結果與所建立的Gibbs代數(shù)方程（式12）結果基本一致（見本文支持信息圖S2、 圖S3及表S1）. 由以上討論可知， 相較Gibbs微分方程， Gibbs代數(shù)方程可以提供更多的信息.

2.3　吸附層厚度

采用所建立的ALT方程（式16）計算了各體系隨2，b的變化， 由圖5（A）可以看出， 各體系呈現(xiàn) 相同的變化趨勢. 對于給定的短鏈醇體系， 隨著2，b增大，降低（呈指數(shù)衰減趨勢）， 這與Li等［7］、 Azizian等［14］和Hyde等［8］（所研究水-乙醇體系）的結果一致. 有趣的是， 從-lg2，b圖可以清晰看出， 在低2，b范圍（lg2，b<0.001）內， 甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系的值依次增大， 而在高2，b區(qū)（lg2，b>0.02）依次減小［圖5（B）］， 這與其極限（τ→0和τ→1）值變化一致（表3）. Bagheri等［13］和 Azizian等［14］研究了甲醇、 乙醇、 異丙醇、 正丙醇和正丁醇水溶液在低2，b范圍內的值， 呈依次增大趨勢， 本文結果與之一致. 由式（16）～式（18）可以看出， 在低2，b范圍內（和τ→0）值的依次增大主要緣于醇分子尺寸（或2）的依次增大（因增大導致減?。?而在高2，b范圍內的依次減小可歸因于吸附趨勢的依次增強（值增大）.

Fig.5　Change of surface adsorption layer thickness(τ) with bulk composition(x2,b)(A) and its semi?logarithmic coordinate plots(B) for water?alcohol mixtures

另外， 所研究短鏈醇體系的τ→0≈1.19～1.88 nm，τ→1≈0.01～0.19 nm， 其中τ→0結果與文獻報道 值（甲醇： 1.2～1.7 nm［13，14］， 乙醇： 1.3～1.5 nm［8，13］， 異丙醇： 1.85 nm［13］）相近. Li等［7］采用中子反射 測得水-乙醇體系在2，b=0.022， 0.045和0.10時的值分別約為0.6， 0.52和（0.46±0.03） nm， 本文 結果（分別約為0.81， 0.47和0.42 nm）與之基本一致. 所得τ→1值很小， 緣于小的2（x→1）值［實際上τ→1≈22（x→1）］.

需說明的是， 模型預測的τ→0和τ→1值不是純組分1和純組分2的表面層厚度（純液體表面層一般由1～3個分子層形成［6］）， 而是Gibbs過剩層厚度的極限值， 可視為模型特征參數(shù)值， 反映吸附趨勢強弱（極限值越小， 吸附趨勢越強）. 純液體表面層單層厚度可用液體分子等當立方體邊長（ec）表征［6］， 其值列于表1中. 據(jù)極限值和溶質ec值可推算溶質的吸附分子層數(shù)（2=/ec）. 甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇體系的2（x→0）分別約為2.9， 2.9， 3.3， 3.2和3.5， 其2（x→1）值分別約為0.46， 0.19， 0.07， 0.02和0.04.2（x→0）值的相對變化率較?。ǜ黧w系可視為基本相近， 即2（x→0）=3.2±0.3）， 而2（x→1）值的相對變化較大， 表明前者可能主要由水（溶劑）的性質決定， 而后者可能與短鏈醇（溶質）性質有關（如表面分子排布形態(tài)等）， 這還有待進一步研究. 另外， 聚集體吸附不等同于多層吸附， 當2>1時可視為溶質發(fā)生多層吸附.

總之， 所建立的ALT模型方程可合理預測吸附層厚度隨體相組成的變化， 吸附層厚度由吸附趨勢和分子尺寸所決定.

3 結 論

基于SAA模型和Gibbs微分方程， 建立了Gibbs表面過剩（2）和表面過剩層厚度（）代數(shù)方程， 其優(yōu)點是采用SAA模型參數(shù)（和）可在全濃度范圍（2，b=0～1）內計算2和值. 甲醇、 乙醇、 異丙醇、 正丙醇和叔丁醇在水表面的吸附趨勢依次增強， 與其疏水性依次增強趨勢一致； 在低2，b區(qū)， 其2值和值均依次增大， 而在高2，b區(qū)均依次減小.值越大，值越小. 在低2，b區(qū)，2值的增大由和2依次增大所致， 而值的增大主要由2增大所致； 在高2，b區(qū)，2值的減小主要緣于2的依次增大， 而值的減小主要由值增大所致. 對于給定的短鏈醇體系， 隨著2，b的增大，2先急劇增大后降低， 呈現(xiàn)一個最大值， 而其持續(xù)降低. 本研究加深了對水-短鏈醇液體混合物表面吸附行為的認識.

支持信息見http://www.cjcu.jlu.edu.cn/CN/10.7503/cjcu20230316.

［1］ Marmur A.，， 2017，， 164—173

［2］ Kaptay G.，， 2020，， 102212

［3］ Stephan S.， Hasse H.，， 2020，， 319—349

［4］ Vegh A.， Korozs J.， Kaptay G.，， 2022，， 4906—4912

［5］ Privat M.， Bennes R.， Tronel?Peyroz E.，， 1988，， 198—207

［6］ Douillard J. M.，， 2009，， 307—310

［7］ Li Z. X.， Lu J. R.， Styrkas D. A.， Thomas R. K.， Rennie A. R.， Penfold J.，， 1993，， 925—939

［8］ Hyde A. E.， Ohshio M.， Nguyen C. V.， Yusa S. I.， Yamada N. L.， Phan C. M.，， 2019，， 111005

［9］ Raina G.， Kulkarni G. U.， Rao C. N. R.，， 2001，， 10204—10207

［10］ Chen H.， Gan W.， Lu R.， Guo Y.， Wang H. F.，， 2005，， 8064—8075

［11］ Kirmse K.， Morgner H.，， 2006，， 2193—2199

［12］ Rafati A. A.， Bagheri A.， Khanchi A. R.， Ghasemian E.， Najafi M.，， 2011，， 252—258

［13］ Bagheri A.， Fazli M.， Bakhshaei M.，， 2016，， 442—451

［14］ Azizian S.， Moghadam T. F.，， 2011，， 67—71

［15］ Santos M. S. C. S.， Reis J. C. R.，， 2018，， 419—428

［16］ Butler J. A. V.，， 1932，， 348—375

［17］ Guggenheim E. A.， Adam N. K.，， 1933，， 218—236

［18］ Eberhart J. G.，， 1966，， 1183—1186

［19］ Connors K. A.， Wright J. L.，， 1989，， 194—198

［20］ Laaksonen A.， Kulmala M.，， 1991，， 6745—6748

［21］ Wang Z. L.， Chen Z. X.，，1995，（2）， 191—195（王兆林， 陳鐘秀. 高等學?；W報， 1995，（2）， 191—195）

［22］ Pi?eiro á.， Brocos P.， Amigo A.， Gracia?Fadrique J.， Lemus M. G.，， 2001，， 4261—4266

［23］ Qi W.， Yu X.， Du N.， Hou W.，， 2023，， 507—518

［24］ Ross S.， Morrison I. D.，， 1983，， 121—134

［25］ Strey R.， Viisanen Y.， Aratono M.， Kratohvil J. P.， Yin Q.， Friberg S. E.，， 1999，， 9112—9116

［26］ Menger F. M.， Shi L.， Rizvi S. A. A.，， 2009，， 10380—10381

［27］ Santos M. S. C. S.， Reis J. C. R.，， 2021，， 21571—21578

［28］ Kaptay G.，， 2015，， 5796—5804

［29］ Kaptay G.，， 2019，， 10987—10992

［30］ Wilson G. M.，， 1964，， 127—130

［31］ Brocos P.， Pi?eiro á.， Amigo A.， Gracia?Fadrique J.，， 2007，， 343—353

［32］ Santos M. S. C. S.， Reis J. C. R.，， 2019，， 525—535

［33］ Zhu B.， Gu T.，， 1991，， 1—32

［34］ Qi W. S.， Hou W. G.，， 2023，， 1236—1246（亓文帥， 侯萬國. 中國科學， 2023，， 1236—1246）

［35］ Bermu?dez?Salguero C.， Gracia?Fadrique J.，， 2015，， 10304—10315

［36］ Yano Y. F.，， 2005，， 255—259

［37］ Vázquez G.， Alvarez E.， Navaza J. M.，， 1995，（3）， 611—614

［38］ Gliński J.， Chavepeyer G.， Platten J. K.，， 1995，， 2113—2117

［39］ Hansen C. M.， Hansen Solubility Parameters?a User’s Handbook， CRC Press， Boca Raton， FL， 2007

［40］ Qin J.， Wang X.， Jiang Q.， Cao M.，， 2019，， 1069—1097

Surface Adsorption of Water-alcohol Binary Solutions

YUXianchao1， QIWenshuai1， DENGQuanhua2， HOUWanguo1，3*

（，，，，，250100，）

Algebraic equations of the Gibbs surface excess（2） and adsorption layer thickness（） were developed for binary liquid mixtures based on the surface aggregation adsorption（SAA） model and the Gibbs differential equation， which can predict the change of2andwith the bulk composition（such as the molar fraction of component 2，2，b） in the whole concentration range. The model equations were used to investigate the surface adsorption of water（1）-alcohol（2） binary solutions. The alcohols involved include methanol（MeOH）， ethanol（EtOH），-propanol（PrOH），-propanol（PrOH）， and-butanol（BuOH）. It was found that their adsorption trend on the water surface increases in turn， consistent with the enhancement trend of their hydrophobicity. Their2andincrease in turn in the low2，bregion， while decrease in turn in the high2，bregion， which arise from the differences in adsorption trend and molecular size between alcohols. In addition， for a given alcohol system， with an increase in2，bfrom 0 to 1， its2initially sharply increases and then decreases， showing a maximum， but itsexhibits a continuous reduction. This work provides a better understanding of the surface adsorption behavior of water-alcohol liquid mixtures.

Liquid mixture； Surface adsorption； Surface aggregation； Adsorption layer thickness； Thermodynamic model

2023-07-05

侯萬國, 男, 博士, 教授, 主要從事膠體與界面化學研究. E-mail: wghou@sdu.edu.cn

國家自然科學基金(批準號: 22272088)資助.

O647

A

10.7503/cjcu20230316

2023-09-11.

Supported by the National Natural Science Foundation of China(No.22272088).

（Ed.： V， K， S）