基于密度泛函理論的活性炭孔徑分布改進算法

朱子文，鄭青榕，陳武，王澤浩，唐政

?

基于密度泛函理論的活性炭孔徑分布改進算法

朱子文，鄭青榕，陳武，王澤浩，唐政

（福建省船舶與海洋工程重點實驗室，集美大學輪機工程學院，福建廈門 361021）

為提高非定域密度泛函理論（NLDFT）預測活性炭孔徑分布（PSD）的精度，考慮了活性炭孔壁面晶體粗糙度對結果的影響。在傳統NLDFT基礎上，結合吸附壁面碳原子的密度分布，推導出改進NLDFT，預測了氬在光滑及具粗糙碳晶體表面的吸附平衡，并根據87.3 K、氬在活性炭上的吸附平衡數據，在由兩種NLDFT確定了不同孔徑的理論等溫線核后，由尋優函數確定活性炭在0.35～12 nm區間的PSD。結果表明，以改進的NLDFT預測活性炭的PSD時，確定的活性炭孔徑呈連續分布，預測平衡數據的相對誤差小于10%；傳統NLDFT確定的孔徑在1 nm處出現斷點，最大的預測相對誤差范圍達45%。改進NLDFT能較準確預測氬在具有粗糙晶體碳表面活性炭的PSD。

氬；活性炭；吸附；分布；密度泛函理論

引 言

由于活性炭具有較大的微孔容積和比表面積，其在物理吸附儲氫、吸附式天然氣存儲（ANG）以及二氧化碳的吸附捕集等領域具有光明的應用前景[1]。然而，由于吸附質分子直徑及極性的差異，適合于不同吸附質分子吸附的活性炭具有不同的孔徑分布（PSD），如何準確獲取用于描述吸附平衡的活性炭PSD極為關鍵。

當前，表征活性炭PSD的方法主要有BJH、t-plot、H-K、DR理論和非定域密度泛函理論（NLDFT），NLDFT具有較高的預測精度[2-3]。吸附過程中，近吸附劑壁面的吸附質分子受吸附勢作用，密度分布出現振蕩。在早期的研究中，Tarazona等[4-5]提出了光滑密度近似（SDA）的概念，Lastoskie等[6]將SDA應用于確定活性炭的PSD，Neimark[7]則引入不定Lagrange因子法（ILM）簡化求解過程，而Ustinov等[8]在求解過程時進一步考慮了一維能量分布不均勻項（EH）。隨后2D-NLDFT[9]、基于基本度量理論（FMT）[10-11]及統計締合流體理論（SAFT）[12-13]的NLDFT模型被逐漸應用。然而，由上述方式確定的流體在第一層內的局部密度過于密集，偏離了假設的硬球模型，PSD計算結果在1 nm處出現不連續的情況[14]。為此，Ustinov將非石墨化碳黑以及活性炭狹縫孔視為能量分布非均勻表面，通過對實驗數據進行非線性擬合確定了表面碳原子密度分布誤差函數，并將具有粗糙碳晶體表面近似為吸附質的淬火狀態（quenched），而靠近壁面區域的流體分子則近似為退火狀態（annealed），通過改進Carnahan-Starling(CS)方程求解混合狀態過剩自由能，最終由吸附數據擬合確定碳原子表面對流體分子施加的作用勢[15-19]。Ravikovitch等[20]則將具有粗糙碳晶體表面視為淬火相態的碳原子，提出了驟冷固體密度函數理論（QSDFT），利用FMT求解該二元硬球流體排斥項，預測了氮、氬在硅膠上的吸脫附平衡，并通過X射線衍射儀（XRD）衍射結果確定表面碳原子密度分布后，計算了活性炭的PSD[21-22]。鄭青榕等[23-24]引入碳原子表面密度分布誤差函數，結合氫在活性炭以及石墨化碳黑上的吸附數據，確定了活性炭的PSD。上述文獻引入碳原子表面密度分布的誤差函數后，均保證了在1nm處PSD的連續性，提高了計算精度。而對于處于亞臨界區域的吸附質，在壓力變化較大的區間，吸附相為液相或過冷液態，粗糙的吸附壁面使流體出現毛細凝聚和起泡等現象。為準確描述上述現象，改進FMT以及SAFT等理論被拓展應用于確定受限流體密度分布，但相關理論的泛函方程求解過程復雜[12-13,25]。

基于此，為精確預測活性炭的PSD，基于傳統NLDFT，參考Ustinov以及Ravikovitch等考慮活性炭孔壁面粗糙度對于預測PSD的研究，通過引入表面碳原子密度分布函數改進NLDFT，預測氬在具有光滑與粗糙晶體表面的碳黑上的吸附平衡后，依據氬在活性炭SAC-01上吸附平衡數據，比較了兩種NLDFT確定的活性炭PSD。

1 試 驗

活性炭（SAC-01）試樣由福建寧德鑫森碳業有限公司提供。為確定SAC-01活性炭的PSD，本文選用Micromeritics 3Flex固體表面分析儀于液氬溫度下測得吸脫附等溫線，結果如圖1。更多實驗信息參閱文獻[26]。

從圖1看出，氬在SAC-01活性炭上的吸脫附等溫線在/0>0.4區域具有“滯回”現象，表明SAC-01活性炭具有介孔和大孔。通過BET法以及脫附數據利用BJH法確定的結構參數如表1所示。

表1 SAC-01活性炭結構參數

2 基本理論

2.1 傳統理論

在給定化學勢和溫度的條件下，自由能是決定系統平衡特性最重要的熱力參數。對于非均質系統，一維分布的受限流體的自由能可表示為[8]

式中，()為流體沿壁面處垂直方向的密度分布曲線。吸附平衡時，吸附自由能為最小值，密度分布曲線滿足。

本征自由能[()]由理想氣體自由能id和流體分子的過剩自由能ex構成。id可由式(2)求得

其中，B為Boltzmann常數，為de-Broglie波長。而ex可區分為流體分子間吸引項att和根據硬球分子間排斥力HS[()]。根據平均場理論，流體分子間吸引項為[27]

(3)

其中，()與(′)為位于三維空間內與′局部位置的密度分布，對于一維模型，局部密度分布函數只與垂直壁面方向相關；相吸自由能att由WCA擾動理論近似[28]。

(5)

(7)

式中，w為加權系數，=0,1,2。對于一維平均密度求解過程可參閱文獻[9]。

傳統NLDFT理論認為碳材料表面可視為光滑壁面。壁面對于流體分子施加的作用勢ext可由10-4-3勢能方程確定

式（8）相關參數如表2，物理意義參閱文獻[29]。

表2 氬在碳表面上吸附的物性參數[15]

化學勢可由CS狀態方程求得[8]

其中，g為本體氣相密度，則；末項為由WCA擾動理論計算的本體氣相流體分子間相互作用，。

聯立式（1）～式（9）可得受限流體的自由能方程；基于平衡條件，即可對相應的巨勢[()]取最小值而獲得密度分布()。本文利用ILM法求解[()]最小值，詳細過程參見文獻[7]。

2.2 改進NLDFT

Ustinov等[15]認為活性炭孔墻為表面能量不均勻且為粗糙晶體的碳原子表面，并將其近似為非石墨化碳黑表面。對于非石墨化碳黑BP280，可附加關于壁面碳原子表面的密度分布誤差函數

其中，為碳原子表面分布修正系數，與Ustinov不同，本文依據Ravikovitch等對于BP280的XRD衍射結果，取0.026nm[16]。對于近壁面的硬球流體，受到的壁面作用勢呈現不均勻分布，該區域流體密度分布需要引入附加項加以修正

(11)

其中，e()為等效密度分布；在靠近壁面區域，可將壁面視為粗糙晶體的碳表面，吸附質分子與表面碳原子晶體呈相互作用的混合狀態，由于該區域表面碳原子的密度大于吸附質分子密度卻遠小于石墨層內碳原子密度，故可將表面碳原子密度等效為吸附質的淬火相態密度m，。同樣利用SDA求解其光滑密度，則吸附壁面對流體施加的排斥項為[16]

與Ustinov等利用最小二乘法擬合確定壁面對于流體分子施加的作用勢不同，本文利用WCA擾動理論確定碳原子與流體分子間的相吸自由能[16]

(13)

考慮吸引與排斥項的共同作用，通過積分求得流體受到壁面施加的作用勢ext

其中，s為碳黑石墨層內碳原子密度。

由式（1）～式（13）確定粗糙晶體碳表面吸附平衡時的自由能，并對相應的巨勢[()]取最小值即可得到密度分布()。將傳統NLDFT求解的密度分布設為初值，再由數值迭代求解粗糙表面的密度分布()。在此，針對一維分布的受限空間，以Δ0.005ff構建空間內相應位置分布的矩陣，通過Picard迭代利用Matlab的矩陣運算離散求解[25]。

3 吸附平衡分析

3.1 氬在碳黑上的吸附平衡

根據文獻[30]，87.3 K時，氬在石墨化碳黑CarbopackF（BET=6.2 m2·g-1）以及非石墨化碳黑BP280（BET=36 m2·g-1）上的吸附數據如圖2所示[圖2（b）為對數坐標]。在此，石墨化碳黑CarbopackF近似為表面光滑的碳表面，而非石墨化碳黑BP280則具有粗糙晶體碳表面。利用傳統NLDFT及改進NLDFT的預測結果如圖2所示，氬在不同碳表面上的密度分布的計算結果如圖3所示。

由圖2（a）可發現，傳統NLDFT的預測結果與氬在石墨化碳黑CarbopackF上吸附等溫線較為接近，而改進NLDFT的預測結果則與氬在BP280上的吸附等溫線更為吻合。從圖2（b）還可發現，在低壓區域，改進NLDFT的預測結果與CarbopackF的吸附數據間差距較大；在/0=0.01～0.1范圍內，改進NLDFT預測值與氬在BP280吸附數據有顯著誤差，但在其他區域具有較高的預測精度。如圖2（b）所示，/0=0.1時氬在碳黑表面上為單層吸附且接近最大值，圖2中改進NLDFT預測值小于實驗值，這說明僅引入一維誤差函數，仍無法準確描述碳黑表面的晶體分布，需進一步建立二維模型[9]。

由圖3更為直觀地發現兩者的區別：坐標原點為第1層碳原子的中心位置，由傳統NLDFT確定密度分布起始點距壁面約0.5ff，而改進NLDFT的密度分布變化則從壁面開始。這表明在粗糙晶體的碳表面，由于表面能量分布不均勻，氬將優先聚集在吸附作用勢較強的吸附位上。此外，傳統NLDFT確定的氬在碳黑表面的密度分布函數()明顯比改進NLDFT確定的振蕩幅度大，如第1層流體分子的局部密度最大值()3ff=7，而改進NLDFT確定的密度分布則較為平緩，其值在液氬（L3ff=0.76[15]）相態附近。顯然，改進的NLDFT更適合于描述氬在具有粗糙晶體表面碳基材料上的吸附平衡。

3.2 PSD的尋優確定

活性炭吸附空間可視為兩側由無限大碳黑平板構成的狹縫孔，其孔徑為兩側壁面間的距離[3]。在實際求解壁面對流體分子的作用勢時，計算距離為對應壁面碳原子間的中心距離*，其中*=+[31]。為石墨層厚度，=0.34 nm。

狹縫孔內的流體分子受到的壁面作用勢為[26]

(16)

求解過程中由于需要對密度分布函數作離散化處理，首先由梯形數值積分求得流體在不同孔徑的理論吸附等溫線核(,)[23]

其中，(,)是與平衡壓力和孔寬對應的吸附量。傳統NLDFT求解的、孔徑為0.35～12 nm的等溫線核()如圖4（a）所示，圖中的壓力與圖1中測試點的平衡壓力對應。由改進NLDFT計算結果如圖4（b）所示。

對比圖4（a）和圖4（b）可發現，改進NLDFT求解的等溫線隨壓力的變化較為平緩；孔徑較小、相對壓力/0>0.01時，氬的吸附相密度將逐漸增大并超過液氬密度；孔徑增大時，不同孔內吸附相發生相變的臨界壓力也隨之增大，尤其在介孔范圍，當/0>0.4時，吸附相密度呈直線急劇上升。

積分吸附等溫方程（IAE）為

其中，max=12 nm；()為孔徑分布函數。根據之前的研究，需要考慮氬在大于max的狹縫孔內的吸附量[23-24]

(19)

式中，為圖1等溫吸附線中吸附平衡壓力的點數，本文=77；附加項ext為氬在大于max的狹縫孔內的吸附量，可由圖2中氬在碳黑CarbopackF/BP280上吸附面密度插值計算獲得，ext為類似碳黑的外表面積。由此可以構建一個關于不同孔的理論吸附等溫線核(,)的×(-1)矩陣，將式（19）附加項加入該矩陣，變為維矩陣|(,)|。同樣地，增量體積Δv與外面積ext也可組成一個長度為的向量|Δ|，式（19）可變為=的矩陣運算。由于增量體積|Δ|為非負向量，直接求解該線性方程無法得到具有物理意義的解。

因此，需要確定關于增量體積|Δ|向量的尋優函數以進行數值求解，設誤差函數為，且滿足約束條件|err|<0.1，則尋優函數為

利用MATLAB中Optimtool工具箱中函數fmincon求出最優解|Δ|，兩種NLDFT預測的吸附等溫線如圖5所示。

圖5 87.3 K氬在活性炭SAC-01上的吸附等溫線測試值與預測值

Fig.5 Isotherms of argon adsorption on activated carbon SAC-01 at 87.3 K from experiments and calculations

由兩種NLDFT的預測值與實驗數據間的相對誤差|err|/exc×100%如圖6所示。

如圖6所示，由改進NLDFT預測值與實驗數據間的相對誤差范圍在±10%，當/0>10-5時，相對誤差約在±1%，而傳統NLDFT的相對誤差范圍在-20%～45%。這說明，將活性炭狹縫孔墻視為粗糙的碳原子表面，并由改進NLDFT預測其孔徑分布更為合理。

孔徑分布函數()與增量體積Δv的關系為[31]

其中，Δ為該孔孔徑的增量；ΔNLDFT為孔內增量面積。利用式(21)即可求出活性炭的孔徑分布PSD如圖7所示，累積孔容與比表面積如表3所示。

圖7 由不同算法確定的活性炭SAC-01的PSD

Fig 7 PSD of activated carbon SAC-01 determined by different approaches

表3 由PSD確定的SAC-01活性炭結構參數

從表3可看出，改進NLDFT的結果較為接近BET法確定的比表面積。從圖7也可看出，兩種方法的預測結果的變化趨勢相似，即具有較大的微孔孔容，介孔變化趨勢一致，但對于孔寬小于1 nm微孔的預測差異明顯，傳統NLDFT預測時，在1 nm處()為0，而改進的NLDFT則是連續的。這表明將活性炭孔墻視為具有粗糙晶體的碳表面是合理的。

4 結 論

根據活性炭具有表面能量分布不均勻的特點，附加碳原子粗糙表面密度分布改進傳統NLDFT，對比了兩種NLDFT對于氬在不同碳黑表面上的吸附行為的預測結果，得到如下結論。

（1）改進的NLDFT可更準確描述氬在具有粗糙晶體表面碳基材料上的吸附平衡。由傳統NLDFT確定的氬在碳黑表面上的密度分布函數()振蕩幅度大，如第1層流體分子的局部密度最大值()3ff7，而改進NLDFT確定的密度分布則較為平緩，其值在液氬相態附近波動。這說明由于在粗糙晶體的碳表面上，氬將優先聚集在吸附作用勢較強的吸附位上，呈現不均勻分布。

（2）改進NLDFT在預測活性炭的PSD時具有較高的精度。利用改進NLDFT確定的活性炭孔徑呈現連續分布，預測值與實驗數據相對誤差分布范圍在±10%，傳統NLDFT確定的活性炭孔徑分布在1 nm處出現斷點，相對誤差分布范圍在-20%～45%，說明將活性炭狹縫孔孔墻視為表面能量分布非均勻、碳表面粗糙的碳黑平板更為合理。

基于SDA理論的CS狀態方程預測二元混合硬球流體的精度有待提高，后續研究需要考慮建立二維模型或利用FMT以及SAFT等理論進一步提高改進NLDFT對于活性炭PSD的預測精度。

符 號 說 明

dHS——硬球分子直徑，nm F,Fid,Fex——分別為本征Helmholtz自由能、理想氣體自由能和來自流體分子的過剩自由能，J fex——體系內吸附質分子的過剩自由能，J H,H*——狹縫孔兩側壁面間的距離及壁面碳原子間中心距離，nm kB——Boltzmann常數，J·K-1 N——吸附量，mmol·g-1 r,r′——分別為流體分子位于三維空間內r與r′的位置，nm S——比表面積，m2·g-1 u——流體分子間相互作用勢，J V——壁面對流體分子施加的作用勢，J v,Δv——比孔容積以及容積增量，cm3·g-1 z——流體分子與壁面碳原子中心的距離，nm Γ——吸附質的表面密度，mmol·m-2 Δ——碳原子層間距，nm δ——碳原子表面密度分布參數，nm ε——勢阱，J ——堆積因子 μ——化學勢，J r,——吸附質密度分布及光滑密度，nm-3 σ——分子直徑，nm Ω——巨勢，J 上角標 AC——活性炭 rep——排斥項 下角標 att——吸引項 BET——比表面積 e——等效 ex——過剩 ext——外部附加 ff——流體分子間 i——序號 m——淬火相態 max——最大值 mic——微孔 min——最小值 s——固體壁面 sf——流體相對于固體壁面

References

[1] POLICICCHIO A, FILOSA R, ABATE S,. Activated carbon and metal organic framework as adsorbent for low-pressure methane storage applications: an overview[J]. Journal of Porous Materials, 2016, 23(6): 1-18.

[2] 楊侃, 陸現彩, 徐金覃, 等. 氣體吸附等溫線法表征頁巖孔隙結構的模型適用性初探[J]. 煤炭學報, 2013, 38(5): 817-821. YANG K, LU X C, XU J T,. Preliminary verification of common calculation methods of pore size distribution of shale based on gas adsorption isotherm[J]. JCSE Journal, 2013, 38(5): 817-821.

[3] LANDERS J, GOR G Y, NEIMARK A V. Density functional theory methods for characterization of porous materials[J]. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2013, 437(6): 3-32.

[4] TARAZONA P. Free energy density functional theory for hard spheres[J]. Physical Review A, 1985, 31(4): 2672-2679.

[5] TARAZONA P, MARINI U, EVANS R. Phase equilibria of fluid interfaces and confined fluids[J]. Molecular Physics, 1987, 60(3): 573-595.

[6] LASTOSKIE C, GUBBINS K E, QUIRKE N. Pore size distribution analysis of microporous carbons: a density functional theory approach[J]. Journal of Physical Chemistry, 1993, 97(18): 1012-1016.

[7] NEIMARK A V. The method of indeterminate Lagrange multipliers in nonlocal density functional theory[J]. Langmuir, 1995, 11(10): 4183-4184.

[8] USTINOV E A, DO D D. Application of density functional theory to analysis of energetic heterogeneity and pore size distribution of activated carbons[J]. Langmuir, 2004, 20(9): 3791-3797.

[9] JAGIELLO J, OLIVIER J P. 2D-NLDFT adsorption models for carbon slit-shaped pores with surface energetical heterogeneity and geometrical corrugation[J]. Carbon, 2013, 55(2): 70-80.

[10] ZENG M, TANG Y, MI J,. Improved direct correlation function for density functional theory analysis of pore size distributions[J]. Journal of Physical Chemistry C, 2009, 113(40): 17428-17436.

[11] BALZER C, CIMINO R T, GOR G Y,. Deformation of microporous carbons during N2, Ar, and CO2adsorption: insight from the density functional theory[J]. Langmuir, 2016, 32 (32): 8265-8274.

[12] SCHINDLER B J, MITCHELL L A, MCCABE C,. Adsorption of chain molecules in slit-shaped pores: development of a SAFT-FMT-DFT approach[J]. Journal of Physical Chemistry C, 2013, 117(41): 21337-21350.

[13] MALHEIRO C, MENDIBOURE B, PLANTIER F,. Density functional theory for the description of spherical non-associating monomers in confined media using the SAFT-VR equation of state and weighted density approximations[J]. Journal of Chemical Physics, 2014, 140(13): 108-112.

[14] USTINOV E A. Comparative features of Ar adsorption on smooth and amorphous surfaces examined by density functional theory[J]. Adsorption, 2008, 14(2): 171-179.

[15] USTINOV E A, DO D D, FENELONOV V B. Pore size distribution analysis of activated Carbopacks: application of density functional theory using nongraphitized Carbopack black as a reference system[J]. Carbon, 2006, 44(4): 653-663.

[16] USTINOV E A, DO D D. Modeling of adsorption in finite cylindrical pores by means of density functional theory[J]. Adsorption, 2005, 11(5): 455-477.

[17] USTINOV E A, DO D D, JARONIEC M. Application of density functional theory to equilibrium adsorption of argon and nitrogen on amorphous silica surface[J]. Applied Surface Science, 2005, 252(3): 548-561.

[18] USTINOV E A, DO D D, JARONIEC M. Adsorption of argon and nitrogen in cylindrical pores of MCM-41 materials: application of density functional theory[J]. Applied Surface Science, 2005, 252(4): 1013-1028.

[19] USTINOV E A, DO D D, JARONIEC M. Equilibrium adsorption in cylindrical mesopores: a modified Broekhoff and de Boer theorydensity functional theory[J]. Journal of Physical Chemistry B, 2005, 109(5): 1947-1958.

[20] RAVIKOVITCH P I, NEIMARK A V. Density functional theory model of adsorption on amorphous and microporous silica materials[J]. Langmuir, 2006, 22(26): 11171-11179.

[21] NEIMARK A V, LIN Y, RAVIKOVITCH P I,. Quenched solid density functional theory and pore size analysis of micro-mesoporous carbons[J]. Carbon, 2009, 47(7): 1617-1628.

[22] GOR G Y, THOMMES M, CYCHOSZ K A,. Quenched solid density functional theory method for characterization of mesoporous carbons by nitrogen adsorption[J]. Carbon, 2012, 50(4): 1583-1590.

[23] 鄭青榕, 蔡振雄, 陳武, 等. 非局域密度泛函理論表征活性炭孔徑分布的改進算法[J]. 低溫與超導, 2010, 38(1): 80-84. ZHENG Q R, CAI Z X, CHEN W,. An improved method to determine PSD of activated carbon by using non-local density functional theory [J]. Cryogenics and Superconductivity, 2010, 38(1): 80-84.

[24] 鄭青榕, DO D D, 陳武, 等. 應用超臨界溫度氫吸附數據表征活性炭結構[J]. 離子交換與吸附, 2010, 26(6): 551-558. ZHENG Q R, DO D D, CHEN W,. Determination of the pore size distribution of the activated carbon by adsorption data of supercritical hydrogen [J]. Ion Exchange and Adsorption, 2010, 26(6): 551-558.

[25] ROTH R. Fundamental measure theory for hard-sphere mixtures: a review[J]. Journal of Physics Condensed Matter, 2010, 22(6): 503-511.

[26] 朱子文, 馮玉龍, 鄭青榕. 甲烷在石墨烯和活性炭上的吸附[J]. 化工學報, 2015, 66(S2): 244-249. ZHU Z W, FENG Y L, ZHENG Q R. Methane adsorption on graphene sheets and activated carbon[J]. CIESC Journal, 2015, 66(S2): 244-249.

[27] BALL P C, EVANS R. Structure and adsorption at gas-solid interfaces: layering transitions from a continuum theory[J]. Journal of Chemical Physics, 1988, 89(7): 4412-4423.

[28] WEEKS J D, CHANDLER D, ANDERSEN H C. Role of repulsive forces in determining the equilibrium structure of simple liquids[J]. Journal of Chemical Physics, 1971, 54(12): 5237-5247.

[29] DO D D, DO H D, TRAN K N. Analysis of adsorption of gases and vapors on nonporous graphitized thermal carbon black[J]. Carbon, 2003, 19(14): 5656-5668.

[30] GARDNER L, MICHAL K A, JARONIEC M. Reference data for argon adsorption on graphitized and nongraphitized carbon blacks[J]. Journal of Physical Chemistry B, 2001, 105(50): 12516-12523.

[31] BIRKETT G, DO D D. New method to determine PSD using supercritical adsorption: applied to methane adsorption in activated carbon[J]. Langmuir, 2006, 22(18): 7622-7630.

Modified approach of determining pore size distribution of activated carbon by using density functional theory

ZHU Ziwen, ZHENG Qingrong, CHEN Wu, WANG Zehao, TANG Zheng

(Provincial Key Laboratory of Naval Architecture & Ocean Engineering, Institute of Marine Engineering, Jimei University, Xiamen 361021, Fujian, China)

In observation of accurately determining the pore size distribution (PSD) of activated carbons by employing thenon-local density functional theory (NLDFT), the effect of surface roughness resulted from the distribution of carbon atoms on the accuracy was evaluated in terms of comparison of the adsorption equilibrium of argon. The classical NLDFT, which was developed based on the approximation to the smooth distribution of the density near the surface, was modified by taking into account the influence from the uneven distribution of the carbon atoms on the surface. The accuracies of the results from the modified model were firstly evaluated by the adsorption data of argon on adsorbent surfaces which respectively have even and uneven distribution of carbon atoms. Experimental data of argon adsorption on the activated carbon at 87.3 K was then taken as a reference, and comparisons were accordingly made between the adsorption isotherms calculated by the classical and modified NLDFT. The PSD, which covers the range of 0.35—12 nm, was finally determined by employing optimization function. Results show that the modified model can bring about the continuous distribution of PSD, and the relative error between the results determined by calculationthe PSD and those from experiment can be less than 10%. However, the calculation from the classical model has an interrupt at a pore width about 1nm, and the resulted maximum relative error amounts to 45%. It suggests that the modified model can more accurately calculate the PSD of the activated carbon having uneven distribution of carbon atoms around the surface.

argon; activated carbon; adsorption; distributions; DFT

10.11949/j.issn.0438-1157.20170013

O 647.3; TQ 424.1

A

0438—1157（2017）09—3328—08

2017-01-06收到初稿，2017-06-12收到修改稿。

鄭青榕。

朱子文（1991—），男，博士研究生。

國家自然科學基金項目（51679107）；福建省自然科學基金項目（2015J01216）。

2017-01-06.

Prof.ZHENG Qingrong, qrzheng816@sina.com

supported by the National Natural Science Foundation of China (51679107) and the Department of Sciences & Technology of Fujian Province(2015J01216).