黃建花，張偉

浙江理工大學(xué)化學(xué)系，杭州 310018

氣體在固體表面的吸附是物理化學(xué)課程中“界面現(xiàn)象”章節(jié)的重要教學(xué)內(nèi)容之一。和液體表面一樣，固體表面層的分子也處于不對(duì)稱力場(chǎng)中。但是，固體無法像液體一樣通過自動(dòng)收縮減小表面積；固體表面只能通過吸附外部氣體或液體分子減小表面分子受力不對(duì)稱程度，從而降低表面吉布斯函數(shù)。根據(jù)固體和氣體分子的相互作用，氣體分子在固體表面的吸附可分為物理吸附和化學(xué)吸附。國(guó)內(nèi)大部分物理化學(xué)教材均詳細(xì)介紹Langmuir單分子層吸附理論。該理論基于四點(diǎn)基本假設(shè)：(1) 單分子層吸附；(2) 固體表面均勻；(3) 被吸附在固體表面上的氣體分子之間無相互作用；(4) 吸附平衡是動(dòng)態(tài)平衡。Langmuir吸附等溫式能很好地描述化學(xué)吸附及低壓氣體的物理吸附，即物理化學(xué)教材中提到的第I類吸附等溫線。

但是大量實(shí)驗(yàn)發(fā)現(xiàn)隨著氣體壓強(qiáng)升高，氣體在固體表面的吸附常呈現(xiàn)多分子層吸附，吸附等溫線與Langmuir單分子層吸附理論不符。1938年，Brunauer、Emmett和Teller在Langmuir單分子層吸附模型的基礎(chǔ)上提出了多分子層吸附理論，即著名的BET理論[1]。BET理論不僅能描述第I類吸附等溫線，還能很好地描述第II和第III類吸附等溫線，是目前應(yīng)用最廣的吸附理論，BET法是目前測(cè)量固態(tài)物質(zhì)比表面積的標(biāo)準(zhǔn)方法。由于BET理論比較抽象，而且BET公式的推導(dǎo)過程涉及較多數(shù)學(xué)知識(shí)，因此大部分物理化學(xué)教材都沒有詳細(xì)介紹BET理論和BET公式的推導(dǎo)過程，而是直接給出BET公式[2–4]。學(xué)生只能被動(dòng)接受BET公式，很難透徹理解BET公式的物理意義。我校應(yīng)用化學(xué)專業(yè)和材料化學(xué)專業(yè)均開設(shè)了測(cè)量固體催化劑比表面積的綜合化學(xué)實(shí)驗(yàn)，面對(duì)不同材料，全自動(dòng)比表面分析儀輸出各類不同吸附脫附等溫線和比表面積，學(xué)生無法正確分析實(shí)驗(yàn)數(shù)據(jù)，只能生搬硬套。由于學(xué)生無法正確理解物質(zhì)比表面積與結(jié)構(gòu)、結(jié)構(gòu)與性能之間的關(guān)系，因此，很難達(dá)到預(yù)期的教學(xué)效果。

筆者參考BET多分子層吸附理論的原始文獻(xiàn)[1]及劉國(guó)杰和黑恩成教授編著的《物理化學(xué)導(dǎo)讀》[5]，在課堂教學(xué)過程中向?qū)W生詳細(xì)介紹BET多分子層吸附理論和BET公式的推導(dǎo)過程。BET理論保留Langmuir單分子層吸附理論的第(2)、(3)和(4)點(diǎn)假設(shè)，只對(duì)第(1)點(diǎn)假設(shè)作了修正。BET理論認(rèn)為已被固體吸附的氣體分子可通過氣體分子之間的范德華作用在垂直方向上吸附其他氣體分子，形成第二、第三等多分子層吸附。第一層即固體表面層吸附的吸附熱qa與氣體分子和固體表面的相互作用相關(guān)；第二及以上各層即氣體吸附層的吸附均源于氣體分子之間的范德華作用，吸附熱可視為氣體的凝聚熱qc。當(dāng)吸附達(dá)到平衡時(shí)，逐層的吸附速率和脫附速率相等，即由第i? 1層變成第i層的吸附速率等于由第i層脫附變成第i? 1層的脫附速率?；谏鲜黾僭O(shè)，得到著名的二參數(shù)BET公式：

BET公式中的Na是壓力為p時(shí)的氣體吸附量，p*和分別是飽和蒸氣壓和氣體在固體表面的單分子層飽和吸附量，參數(shù)c是與固體、氣體相關(guān)的特性參數(shù)，

其中k1和k?1分別是第一層的吸附和脫附常數(shù)，k2和k?2分別是第二及以上各層的吸附和脫附常數(shù)，kB和T分別是玻爾茲曼常數(shù)和以開爾文溫標(biāo)表示的吸附溫度。顯然c與吸附熱相關(guān)，吸附作用越強(qiáng)，吸附等溫線在低壓區(qū)上升越快，c值可通過分析BET實(shí)驗(yàn)數(shù)據(jù)得到[6]。BET公式(1)也可以寫成，

為了幫助學(xué)生正確理解BET吸附理論，筆者結(jié)合自身科研特長(zhǎng)，在課堂教學(xué)過程中結(jié)合計(jì)算機(jī)模擬氣體分子在固體表面的吸附脫附過程。通過改變固體表面層和氣體吸附層的吸附和脫附常數(shù)，分別模擬得到第I、II、III類吸附等溫線；結(jié)合可視化軟件使學(xué)生能直觀看到氣體分子在固體表面的吸附脫附動(dòng)力學(xué)過程，從而正確理解BET多分子層吸附理論和Langmuir單分子層吸附理論，為后續(xù)綜合化學(xué)實(shí)驗(yàn)提供知識(shí)儲(chǔ)備。此外，計(jì)算機(jī)模擬速度快、成本低、不受實(shí)驗(yàn)條件和課堂教學(xué)時(shí)間限制，學(xué)生隨時(shí)可以通過在計(jì)算機(jī)上調(diào)整各種參數(shù)模擬各類體系，從而深入理解各知識(shí)點(diǎn)，計(jì)算機(jī)模擬輔助化學(xué)教學(xué)能有效提高教學(xué)質(zhì)量[7–12]。例如，為了幫助學(xué)生正確理解碰撞理論和過渡態(tài)理論，趙力維博士等[7]在親核取代反應(yīng)機(jī)理的課堂教學(xué)中，引入商用Gaussian 09計(jì)算化學(xué)軟件計(jì)算氯離子與溴代甲烷的一步親核取代反應(yīng)中的能量變化和反應(yīng)途徑，通過Gaussian View軟件對(duì)反應(yīng)物、活化絡(luò)合物和生成物結(jié)構(gòu)實(shí)現(xiàn)可視化，從而幫助學(xué)生直觀理解親核取代反應(yīng)機(jī)理。PHREEQC軟件也被用于課堂教學(xué)中構(gòu)建多種化學(xué)過程的反應(yīng)模型，增進(jìn)學(xué)生對(duì)化學(xué)反應(yīng)動(dòng)力學(xué)的理解[8]。X射線衍射和拉曼光譜是分別表征物質(zhì)的晶體結(jié)構(gòu)、化學(xué)鍵及分子振動(dòng)等微觀信息的重要分析技術(shù)，這些微觀信息肉眼無法觀察，對(duì)學(xué)生而言過于抽象。溫大尉博士[10]將免費(fèi)的晶體學(xué)軟件和第一性原理計(jì)算軟件引入課堂教學(xué)中，把抽象難懂的微觀信息以可視化的形式在計(jì)算機(jī)上實(shí)時(shí)展現(xiàn)，幫助學(xué)生深入理解物質(zhì)的微觀結(jié)構(gòu)。計(jì)算機(jī)模擬還能將宏觀的色譜分離抽象成粒子在色譜柱中的微觀運(yùn)動(dòng)，直觀展現(xiàn)待分離粒子在色譜柱內(nèi)的擴(kuò)散行為，幫助學(xué)生理解色譜分離原理[12]。山東大學(xué)苑世領(lǐng)教授等[9]利用商用Materials Studio軟件開設(shè)了計(jì)算化學(xué)實(shí)驗(yàn)，旨在幫助學(xué)生從分子水平理解物質(zhì)的結(jié)構(gòu)–性質(zhì)關(guān)系，培養(yǎng)學(xué)生從微觀角度理解宏觀現(xiàn)象，并解決化學(xué)問題的能力。

本文將計(jì)算機(jī)模擬技術(shù)引入到抽象難懂的BET吸附理論的課堂教學(xué)中，通過自編程序模擬氣體分子在固體表面的吸附，結(jié)合實(shí)時(shí)可視化使學(xué)生直觀看到氣體分子在固體表面的吸附脫附動(dòng)力學(xué)。在物理化學(xué)課程之前學(xué)生已學(xué)習(xí)了“C語言設(shè)計(jì)”和“C語言設(shè)計(jì)實(shí)訓(xùn)”兩門必修課，部分學(xué)生還選修了“Python語言程序設(shè)計(jì)”等編程課程，因此學(xué)生已具備一定的編程能力。筆者的科研方向?yàn)楦叻肿芋w系的計(jì)算機(jī)模擬，擁有豐富的編程經(jīng)驗(yàn)，可以引導(dǎo)學(xué)生修改程序并合理調(diào)整參數(shù)模擬各類吸附等溫線，從而更好地理解BET理論及其應(yīng)用原理。

1 計(jì)算機(jī)模擬氣體分子在固體表面吸附

1.1 模型和模擬方法

模擬在立方格子上進(jìn)行，固體表面由L×L的格點(diǎn)組成，氣體分子落在格點(diǎn)上。氣體的壓強(qiáng)與氣體分子數(shù)密度n成正比，即p=nkBT。氣體以平均速度＜v＞向各個(gè)方向隨機(jī)運(yùn)動(dòng)，氣體分子的運(yùn)動(dòng)、吸附、脫附用Monte Carlo方法模擬，假定氣體分子在每個(gè)Monte Carlo模擬時(shí)間(1/＜v＞)運(yùn)動(dòng)一個(gè)格點(diǎn)長(zhǎng)度，那么每個(gè)Monte Carlo模擬時(shí)間有M=nL2/6個(gè)氣體分子向下運(yùn)動(dòng)到固體表面或吸附層表面。格點(diǎn)長(zhǎng)度作為長(zhǎng)度單位，L可隨模擬體系的尺寸而調(diào)整，本文選取L= 50。模擬系統(tǒng)的高度沒有限制。

氣體分子可吸附在固體表面或吸附層上方，吸附的氣體分子也可以沿垂直方向脫附但沒有水平方向的移動(dòng)。氣體分子吸附到固體表面層(稱為第一層)的吸附常數(shù)為k1，氣體分子從固體表面層脫附的脫附常數(shù)為k-1；同理，氣體分子吸附到其他氣體吸附層的吸附常數(shù)為k2，氣體分子從氣體吸附層脫附的脫附常數(shù)為k-2。這些吸附、脫附常數(shù)與固體、氣體的性質(zhì)相關(guān)，它們的取值沒有直接關(guān)聯(lián)。圖1是氣體分子在光滑固體表面的吸附脫附示意圖和程序設(shè)計(jì)圖。每一個(gè)格點(diǎn)的吸附氣體數(shù)目設(shè)為N(x,y)，N= 0代表空白表面；N(x,y)的平均值即為單位面積的氣體吸附量Na。當(dāng)取k2=k-2= 0時(shí)，我們可以模擬Langmuir單分子層吸附。

圖1 氣體分子在光滑固體表面吸附脫附示意圖(上)和程序設(shè)計(jì)圖(下)

如氣體分子吸附脫附過程的程序設(shè)計(jì)圖(圖1)所示，氣體分子以吸附常數(shù)k1或k2隨機(jī)地吸附到格點(diǎn)的最上層，當(dāng)吸附發(fā)生時(shí)N(x,y)加1；同樣地，格點(diǎn)上氣體分子的最上層分子以脫附常數(shù)k-1或k-2隨機(jī)地脫附，當(dāng)脫附發(fā)生時(shí)N(x,y)減1。當(dāng)吸附-脫附達(dá)平衡后，單位面積的氣體吸附量Na達(dá)到動(dòng)態(tài)平衡。吸附時(shí)，N(x,y) = 0表示表面層吸附，N(x,y) ＞ 0表示氣體層吸附；同理，脫附時(shí)，N(x,y) = 1表示從表面脫附，N(x,y) ＞ 1表示從氣體層脫附。

飽和蒸氣壓p*值取決于氣體分子凝聚的吸附常數(shù)k2和氣體分子的脫附常數(shù)k?2。模擬過程中我們不計(jì)算p*，而是計(jì)算約化壓強(qiáng)p/p*。根據(jù)BET理論，吸附脫附達(dá)到平衡時(shí)，第i層吸附分子數(shù)θi是第i? 1層吸附分子數(shù)θi-1的p/p*倍，即θi/θi-1=p/p*。因此模擬中，我們通過計(jì)算θ2/θ1和θ3/θ2的平均值得到p/p*= (＜θ2/θ1＞ + ＜θ3/θ2＞)/2。

計(jì)算機(jī)模擬程序流程如圖2所示。模擬程序采用Python語言編寫，方便學(xué)生學(xué)習(xí)和修改。吸附過程可以調(diào)用提前編譯好的Fortran可執(zhí)行程序，以便更快完成模擬。

圖2 計(jì)算機(jī)模擬流程示意圖

1.2 模擬結(jié)果和討論

吸附等溫線的形狀與吸附熱和凝聚熱相關(guān)，不同體系可通過調(diào)節(jié)吸附常數(shù)和脫附常數(shù)k1、k?1、k2和k?2實(shí)現(xiàn)。我們分別模擬了三類吸附等溫線：(1) 只存在單分子層化學(xué)吸附(qa＞＞qc)；(2) 低壓區(qū)以單分子層吸附為主(qa＞qc)；(3) 多分子層物理吸附為主(qa＜qc)。k1、k?1、k2和k?2四個(gè)參數(shù)與吸附熱qa、凝聚熱qc沒有直接關(guān)聯(lián)，為了明顯給出三類吸附等溫線的差別，本文選擇了不同的k1、k?1、k2和k?2進(jìn)行模擬。

當(dāng)qa＞＞qc時(shí)，體系近似于單分子層吸附(第I類吸附等溫線)。我們?cè)O(shè)置固體表面第一層的吸附和脫附常數(shù)分別為k1= 0.9和k?1= 0.05，氣體吸附層的吸附和脫附常數(shù)為k2=k?2= 0。因?yàn)闅怏w分子只有單層吸附，我們無法得到p/p*，所以采用氣體分子數(shù)密度n表示壓強(qiáng)p=n。模擬結(jié)果如圖3所示，在較低壓強(qiáng)p下，單位面積的平衡吸附量Na隨壓強(qiáng)升高而增加，當(dāng)壓強(qiáng)足夠大時(shí)，吸附達(dá)到飽和，不再隨壓強(qiáng)上升而增加。模擬結(jié)論能很好地?cái)M合成Langmuir單分子層吸附等溫線：，其中單位面積的飽和吸附量，吸附常數(shù)b= 2.9。飽和吸附量表示97%的固體表面被氣體分子覆蓋。對(duì)于不同性質(zhì)的固體表面和氣體分子，Nam和b值可通過改變k1和k?1進(jìn)行調(diào)節(jié)。模擬中選擇大的k1值和小的k?1值可以使系統(tǒng)更快地達(dá)到吸附飽和。

圖3 固體表面上單位面積的氣體吸附量Na隨氣體壓強(qiáng)p的變化關(guān)系

對(duì)于qa＞qc的體系，我們固定固體表面層的吸附和脫附常數(shù)k1= 0.9和k?1= 0.05，氣體吸附層的吸附和脫附常數(shù)分別取k2= 0.01和k?2= 0.05。模擬結(jié)果見圖4(a)，很明顯，我們的模擬結(jié)果符合第II類吸附等溫線，能很好地用二參數(shù)BET公式進(jìn)行擬合：

圖4 (a) 氣體吸附層的吸附常數(shù)k2 = 0.01時(shí)模擬得到的第II類吸附等溫線；(b) k2不同時(shí)，單位面積吸附量Na隨壓強(qiáng)p的變化

當(dāng)qa＜qc時(shí)，以多分子層物理吸附為主，固體表面層的吸附和脫附常數(shù)分別取k1= 0.1和k?1= 0.05，氣體吸附層的吸附和脫附常數(shù)分別取k2= 0.2和k?2= 0.05。模擬結(jié)果(圖5(a))能很好地再現(xiàn)第III類吸附等溫線，擬合參數(shù)和c= 0.5。我們進(jìn)一步模擬了氣體吸附層的脫附常數(shù)k?2對(duì)吸附的影響，脫附常數(shù)k?2增加導(dǎo)致脫附作用增強(qiáng)，在相同壓強(qiáng)下氣體分子的吸附量減少，結(jié)果如圖5(b)所示。小的k?2值代表弱脫附，因此小的k?2值有利于觀察到第III類吸附等溫線。

圖5 (a) 氣體吸附層的脫附常數(shù)k?2 = 0.05時(shí)模擬得到的吸附等溫線；(b) k?2不同時(shí)，單位面積吸附量Na隨壓強(qiáng)p的變化

我們還可以實(shí)時(shí)跟蹤氣體分子的吸附動(dòng)力學(xué)過程。以k1= 0.9，k2= 0.01，和k?1=k?2= 0.05的體系(第II類吸附等溫線)為例，如圖6所示，隨著吸附時(shí)間延長(zhǎng)，吸附量Na快速增加，并達(dá)到吸附平衡。平衡吸附量隨氣體壓強(qiáng)增大而增加；吸附達(dá)到平衡所需的時(shí)間也隨氣體壓強(qiáng)增大而增加。模擬過程中，我們還實(shí)現(xiàn)了氣體分子在固體表面的吸附過程以及吸附脫附動(dòng)力學(xué)等微觀信息的實(shí)時(shí)可視化，如圖7所示的不同壓強(qiáng)下的氣體實(shí)時(shí)吸附圖能讓學(xué)生直觀看到氣體分子在固體表面的吸附機(jī)制和吸附脫附動(dòng)力學(xué)，更好地理解各種吸附模型。將計(jì)算機(jī)模擬穿插到課堂教學(xué)中能充分活躍課堂氣氛，實(shí)現(xiàn)教與學(xué)互動(dòng)，提高教學(xué)質(zhì)量。

圖6 不同壓強(qiáng)下氣體吸附量隨吸附時(shí)間的變化過程

圖7 不同壓強(qiáng)下的氣體實(shí)時(shí)吸附圖

2 教學(xué)方案

計(jì)算機(jī)模擬輔助BET吸附理論教學(xué)包含課堂演示和課后自學(xué)兩部分。課堂教學(xué)安排1學(xué)時(shí)，教師按如下方案開展課堂教學(xué)：

(1) 介紹氣體分子在固體表面吸附的計(jì)算機(jī)模型、模擬方法、模擬程序的設(shè)計(jì)流程和參數(shù)調(diào)節(jié)；

(2) 課堂演示程序運(yùn)行方法和步驟，學(xué)生在教室的投影屏幕上能實(shí)時(shí)觀察氣體分子在固體表面的吸附過程以及吸附脫附動(dòng)力學(xué)；

(3) 程序運(yùn)行結(jié)果通過Origin軟件處理后確定吸附等溫線的類型。Origin軟件學(xué)生已在“化學(xué)信息學(xué)”課程中學(xué)過。

由于課堂教學(xué)時(shí)間有限，我們將模擬程序源代碼、程序說明和參數(shù)調(diào)節(jié)方法分享到班級(jí)QQ學(xué)習(xí)群，方便學(xué)生在課后自學(xué)練習(xí)。Python程序簡(jiǎn)單易學(xué)，非計(jì)算機(jī)專業(yè)的學(xué)生也能很快掌握和使用。Python程序可以實(shí)時(shí)輸出氣體吸附量Na隨吸附時(shí)間的變化過程、吸附等溫線，以及氣體實(shí)時(shí)吸附圖等內(nèi)容。學(xué)生在課后可以方便地調(diào)節(jié)參數(shù)和修改程序，激發(fā)學(xué)生的求知欲和學(xué)習(xí)積極性。如前所述，BET公式假設(shè)了固體表面均勻、表面各吸附位點(diǎn)的吸附熱相等、多層吸附時(shí)每層凝聚熱都相同等條件，因此，實(shí)驗(yàn)發(fā)現(xiàn)BET公式只在p/p*= 0.05–0.35范圍內(nèi)適用[13]。但模擬可以遵守這些BET假設(shè)，因此可以在更大的壓強(qiáng)范圍內(nèi)適用，如圖4(a)和5(a)所示。學(xué)生還可以通過修改程序，模擬更接近實(shí)際的體系，如粗糙固體表面、吸附位點(diǎn)不均勻的固體表面以及被吸附的同層氣體分子之間存在相互作用等，模擬結(jié)果在QQ學(xué)習(xí)群進(jìn)行交流，從而開闊學(xué)生視野，激發(fā)學(xué)生的科研興趣。

該教學(xué)方案率先在2020級(jí)應(yīng)用化學(xué)專業(yè)試行，學(xué)生反響積極，普遍反映有助于加深理解BET吸附理論，很多同學(xué)希望物理化學(xué)的其他抽象難懂的知識(shí)點(diǎn)以及涉及復(fù)雜數(shù)學(xué)推導(dǎo)的公式也能開展類似的計(jì)算機(jī)模擬輔助教學(xué)。

3 結(jié)語

物理化學(xué)課程教學(xué)中有很多看不見摸不著、很抽象但又具廣泛應(yīng)用背景的知識(shí)點(diǎn)，通過把計(jì)算機(jī)模擬技術(shù)引入到物理化學(xué)課堂教學(xué)中，教師能將抽象概念和理論表現(xiàn)為定量數(shù)據(jù)和直觀圖像，激發(fā)學(xué)生的求知欲；還能充分活躍課堂氛圍，提升教學(xué)趣味性，實(shí)現(xiàn)教與學(xué)互動(dòng)。這種將理論和計(jì)算機(jī)模擬相結(jié)合的教學(xué)方法豐富了現(xiàn)代教學(xué)方法，能激發(fā)學(xué)生的學(xué)習(xí)主動(dòng)性和科研興趣，提高教學(xué)質(zhì)量，為培養(yǎng)學(xué)生的創(chuàng)新能力和建立正確的科學(xué)研究方法奠定基礎(chǔ)。