淺談氫氣在共價有機骨架材料中的吸附機理論文

淺談氫氣在共價有機骨架材料中的吸附機理論文

　　引 言

　　化石燃料的大規模利用導致了日益嚴重的環境和能源枯竭問題。氫能作為新型綠色能源，具有資源豐富、清潔環保、發熱值高等優勢，為解決能源危機提供了一條可行途徑。作為氫能工程的關鍵一環，如何安全、經濟、高效地儲存氫氣成為研究的熱點，各國科學家致力于尋找有效的儲氫材料以達到美國能源部的技術指標。

　　共 價有機骨架(covalent organic frameworks,COFs)材料是近年來合成的一種新型的骨架結構材料，是由有機配體和B—O 團簇通過強共價鍵組裝而成的配位聚合物，僅含有B、O、C 和H 等輕元素。COFs 具有較大的孔隙率和表面積，晶體密度小，且可根據調控有機配體改變孔道結構和化學特性，熱穩定性好，是非常有潛力的儲氫材料。COF-102 和COF-103 的比表面積分別達到了3472、4210 m2·g1，遠遠超過沸石、活性炭等傳統材料。此外，COF 材料密度小，如COF-108密度為0.17 g·cm3，是目前報道晶體中密度最小的。COF-6、COF-8、COF-10 等具有類似石墨一樣的層狀二維結構[4]。Furukawa 等研究了二維COF-1、COF-5、COF-6、COF-8 和三維COF-102、COF-103 對H2 的吸附性能，結果表明三維COFs 的儲氫能力比二維COFs 好，且COF-102 和COF-103的吸附量接近目前儲氫量最高的MOF-177。Cao 等研究了COF-105 和COF-108 對氫氣的吸附行為，發現在77 K 和10 MPa 的情況下，儲氫量可達18.05%和17.80%。Assfour 等采用分子動力學模擬研究了多種COFs 的主要吸附位，結果表明COFs能夠穩定吸附氫氣。Yang 等采用巨正則MonteCarlo 方法研究了低溫下H2 在二維COFs 吸附，發現H2 在二維COFs 的吸附等溫線不存在階梯現象。Klontzas 等采用多尺度理論方法研究H2 在三維COFs 的吸附能，發現H2 在三維COFs 團簇中的吸附能與具有良好儲氫性能的IRMOFs 接近。

　　現有研究表明COFs 具有良好的儲氫能力，但針對H2 在COFs 上機理的研究仍缺乏系統性。計算化學作為新興的研究手段，可以突破傳統研究的局限性，從微觀角度對機理做出解釋，可為最佳吸附材料的設計提供理論依據。本文主要計算了H2 分子在典型COFs 材料(COF-6 、COF-8 、COF-10 、COF-102、COF-103、COF-105、COF-108)主要團簇的吸附能，摻雜改性后COFs 的吸附性能，揭示 H2在COFs 中的吸附機理。本文研究對于COFs 吸附和儲氫性能提高具有一定的指導意義。

　　1 計算模型和方法

　　1.1 密度泛函理論方法的選擇與驗證

　　計算主要采用密度泛函理論(density functionaltheory)方法。計算過程使用Material Studio 軟件包中DMol3 模塊。計算使用All Electron 核處理方法，采用DNP (double numeric polarization basis set)基組和優等精度收斂。

　　為了選取描述范德華弱作用力最優的泛函方法，選擇了苯環與H2 的吸附計算進行驗證，如表1所示。Lee 等計算得到氫氣在苯環的吸附能分布于1.55～3.97 kJ·mol1，通過比較發現GGA/HTCH過高估計H2 與苯環之間的吸附距離，BLY3P、GGA/BLYP 過高估計H2 與苯環的吸附能，GGA/PBE 和GGA/PW91 是描述范德華弱作用力的良好泛函。PBE 泛函已經成功地應用在氫氣和MOF的弱作用體系計算中，本文采用GGA/PBE 的泛函方法。

　　1.2 COFs 團簇模型的建立

　　COFs 模型中的各個原子位置采用XRD 技術測定的實驗數據[3-4]確定并構建模型。COF-6、COF-8、COF-10、COF-102、COF-103、COF-105、COF-108分別由兩個團簇組成。由于元胞含有數百個原子，體系過大，難以計算每個可能的吸附構型，因此本文截取各團簇作為計算模型，根據截斷部位周圍環境使用氫原子進行飽和，得到如下7 種結構，其中存在兩種COFs 共有的團簇，進行結構優化得到穩定構型。前期使用了與本文相同的方法截取團簇模型研究金屬有機骨架材料對CO2 的吸附機理。上述 7 類團簇模型基本包括了COFs 的主要骨架部分。其中，COF-6 由團簇a 和團簇e 兩部分構成，COF-8 由團簇a 和團簇g 兩部分構成，COF-10由團簇a 和團簇f 構成，COF-102 由團簇b 和團簇d 構成，COF-103 由團簇b 和團簇c 構成，COF-105由團簇a 和團簇c 構成，COF-108 由團簇a 和團簇d 構成。

　　1.3 H2 在COFs 中的吸附

　　本文以7 種團簇為計算模型，并考慮H2 分子在團簇中的不同吸附位。以團簇a 為例，主要吸附位為中心苯環、側苯環、BOC 環、BOC 環氧位、BOC環硼位。H2 分子在吸附位考慮垂直(V)和平行(P)兩種吸附取向。對每一種吸附構型進行幾何優化和能量計算，并計算吸附能大小。吸附能定義為吸附前后各物質的總能量的變化，其符號和大小可表示發生吸附的可能性和吸附的程度。吸附能的計算公式如下ΔE=ΔEcluster + H2 (ΔEcluster + ΔEH2) (1)式中，Ecluster 表示團簇的能量，EH2 表示H2 分子的能量，Ecluster+H2 表示吸附H2 分子后團簇與H2的總能量，ΔE 表示吸附能的大小。當吸附能ΔE 為負值的時候表示吸附可以發生，吸附能的絕對值小于40 kJ·mol1 時為物理吸附。吸附能的絕對值越大，則吸附越穩定。

　　1.4 兩個 H2 分子的競爭協同吸附

　　為了研究兩個H2 分子在二維COFs 和三維COFs 的吸附，本節選取二維的COF-10 和三維的COF-105 元胞的完整團簇作為研究對象，團簇模型如圖3 所示。選取主要計算位置有中心苯環、側苯環、BOC 環、有機鏈，將單個H2 分子用不同取向放置于吸附位，計算過程與1.2 節相同。在單個分子吸附位確定的基礎上，將第2 個H2 分子放置在另一吸附位，計算吸附能。

　　成對能量(pairing energy)的計算為Epairing=Eij(Ei+Ej) (2)式中，Ei 和Ej 表示單個分子在吸附位i、j 上的吸附能，Eij 表示兩個分子分別在吸附位i、j 時的吸附能，Epairing 表示成對能量。當Epairing 為正值時表示兩個分子同時的吸附能比單個分子值的和要小，意味著多個氫氣吸附具有競爭效應。當Epairing 為負值時表示兩個分子同時的吸附能比單個分子值的和要大，意味著多個氫氣吸附具有協同效應。

　　1.5 COFs 改性之后的吸附

　　Snurr 等發現吸附能是影響材料儲氫性能的主要因素之一，在低壓時，H2 的吸附量與吸附能呈正比。Cao 等實驗研究表明，Li 摻雜的COF-105和COF-108 在298 K、10 MPa 下氫氣的吸附量分別為6.84%和6.73%，比未改性COFs 在常溫下吸附量有明顯增加。因此對COFs 進行改性是提高儲氫量、增加H2 與COFs 的吸附能的重要途徑之一，然而目前對改性COFs 的具體吸附機理還不清楚。本文研究了COFs 改性后，如金屬摻雜、官能團替代對H2 吸附的影響。主要對團簇a 和團簇b 進行改性，對團簇a 摻雜堿金屬Li、Na、K，確定穩定摻雜位，分析吸附能的變化。對團簇b 進行官能團替代，用甲基—CH3、巰基—SH、氨基—NH2、氯基—Cl 等替換單個苯環的氫位，分析H2 吸附能的變化。

　　2 結果與討論

　　2.1 不同 COFs 團簇對H2 的吸附機理

　　H2分子在7 種團簇的吸附構型以及密度泛函理論計算的吸附能結果。結果表明H2 的吸附能大小分布于0.52～7.65 kJ·mol1 之間，H2和團簇的相對位置也大于范德華半徑，說明H2 在COFs 團簇的吸附主要為物理吸附。

　　比較 H2 和7 種團簇的32 種不同吸附構型下吸附能大小，可以得出如下結論：當H2 垂直接近時，苯環的中心位置為最優先吸附位，存在中心苯環的團簇a 和團簇g 對H2 的吸附能分別為5.84 kJ·mol1和5.88 kJ·mol1，側苯環對于H2 的吸附能小于中心苯環。作為COFs 的典型團簇，BO 環、BOC 環以及位于BOC 環上的氧位對H2 的吸附能也超過了5 kJ·mol1，為H2 垂直吸附時的次優先吸附位。BOC 環上的硼位和氧位作為單原子吸附位，由于接觸面積相較于苯環、BO 環等多原子吸附位小，因此 H2 分子的吸附角度會出現一定偏斜。當H2 平行接近時，BO 環的中心以及氧原子(O)為優先的吸附位，吸附能分別為7.65 kJ·mol1 和7.06 kJ·mol1，BO 環硼位和BOC 環的硼位的吸附能分別為6.31kJ·mol1 和6.53 kJ·mol1，為次優先吸附位，團簇c 的中心Si 位對H2 的吸附能最小。

　　此外，比較只含C、H 元素的有機鏈團簇(團簇e 苯環、團簇f 苯環、團簇d 側苯環、團簇g 中心苯環) 吸附H2 的吸附能數據，發現隨有機鏈團簇C原子數目增加，吸附能增加，這與Snurr 等[17]在研究IRMOF 系列吸附甲烷的模擬研究結果一致。說明有機鏈團簇也是COFs 的主要吸附位，有機鏈團簇C 原子數目增加有助于吸附能的增加，可能與有機鏈改變了孔尺寸從而增加了H2 與有機鏈的有效接觸面積有關。因此，在合成COFs 時，有機配體的選擇是提高比表面積和H2 吸附的主要途徑之一。

　　2.2 兩個 H2 分子在COFs 中的吸附機理

　　在密度泛函理論研究中，計算對象多為單個分子。本文以COF-10 和COF-105 為例，計算兩個H2 分子在團簇不同部位同時吸附的情況，進一步揭示H2 在COFs 中的吸附機理。Lee 等曾采用與本節相同的方法研究H2分子在MOF-5 的協同吸附作用。

　　兩個H2 分子分別在COF-10 和COF-105 的吸附能和成對能量。如表3 所示，單個H2 分子在COF-10 首先吸附于中心苯環、側苯環、BOC 環、有機鏈，再增加一個H2 分子后，吸附能有明顯增加，在7.59～14.39 kJ·mol1 范圍。如表4 所示，單個H2 分子在COF-105 首先吸附于有機鏈Si 位、苯環和BO 環時，增加一個H2 分子后，吸附能分布于3.29～13.23 kJ·mol1。而成對能量Epairing 的值基本小于0，說明當1 個H2 分子發生吸附時，它有利于第2 個H2 分子的吸附，兩者存在協同吸附的關系。

　　2.3 通過改性提高COFs 材料吸附性能的機理

　　本文分別考察了堿金屬摻雜、官能團替代改性方法，采用堿金屬元素Li、Na、K 在COFs 主要團簇a 摻雜，采用甲基—CH3、巰基—SH、氨基—NH2、氯基—Cl 等替換團簇b 單個苯環的H 原子，優化得到穩定構型，計算改性后COFs 對H2分子的吸附能。結果表明，堿金屬在團簇a 中心苯環正上方時，此時堿金屬和團簇的結合能最低，構型最為穩定，其中K、Na、Li 摻雜改性后團簇a 對于H2 的吸附能分別為13.79、12.39、11.04 kJ·mol1，吸附能有明顯的增加，吸附效果為K>Na>Li。原因是堿金屬摻雜產生孤立的金屬原子，強烈的電荷力增加了吸附能。

　　采用官能團替代改性后，采用—CH3、—SH、—NH2、—Cl 改性的團簇b 對H2 的吸附能分別為12.42、7.96、7.34、4.84 kJ·mol1，官能團替代后吸附效果為：甲基(—CH3)>巰基(—SH)>氨基(—NH2)>氯基(—Cl)，氯基對于吸附能提高無明顯變化。H2 分子具較弱的四極矩，巰基和氨基對于吸附能的提高可能是因為其對H2 分子極性產生改變，而甲基的存在增加了苯環的碳原子數目，導致相互作用力增強，說明對COFs 團簇進行化學改性，引入具有較大偶極矩的化學官能團或者含碳基團可以增加COFs 材料對于H2 的吸附能力。

　　為了進一步深入了解改性的影響，本文計算了團簇改性前后的軌道，分析其最高占據軌道(highestoccupied molecular orbital，HOMO)分布。根據前線軌道理論分析，團簇b 改性前后的HOMO 和H2 的最低空軌道( lowest occupiedmolecular orbital，LUMO)接近時，彼此對稱性不匹配，不能發生反應，因此屬于物理吸附，這與2.1節中結論是一致的。通過比較改性前后HOMOs 軌道變化，發現經過改性后的團簇b 電荷密度離域化明顯地非對稱性增強，HOMOs 分布于改性官能團附近，有助于為骨架提供更多的極性電荷，從而增強與H2 分子軌道的相互作用。

　　3 結論

　　本文采用密度泛函理論計算了H2 分子在COFs材料的7 種主要團簇中的吸附位置和吸附能。計算結果表明，H2 在COFs 中的相互作用力是色散力，屬于物理吸附。當H2 垂直吸附時，團簇中心苯環的正上方為最優先吸附位;當H2 平行吸附時，BO 環的中心以及氧原子上為優先的吸附位。有機鏈碳原子數目的增加有利于吸附能增加。兩個氫氣吸附時，吸附能總體增加，屬于協同吸附。摻雜了Li、Na、K 等堿金屬團簇，由于金屬原子的強電荷作用，吸附能明顯增加。采用甲基—CH3、巰基—SH、氨基—NH2 官能團改性后的團簇，電荷密度離域化非對稱性增強，極性增強，對H2 的吸附能增加，說明金屬摻雜和官能團替代是COFs 增加H2 吸附能的有效方法。本文研究可為設計新型儲氫材料提供理論依據。

【淺談氫氣在共價有機骨架材料中的吸附機理論文】相關文章：

淺談養心氏片的臨床運用機理論文08-02

有機化學復習策略淺談07-01

淺談有機生姜豐產栽培技術10-09

淺談工程造價中的甲供材05-27

動機理論下的音樂教學論文10-06

有機化學與生活論文08-19

淺談壯醫藥浴在預防疾病中的作用機理及應用08-04

淺談教育論文寫作09-29

淺談生命教育的論文06-27

淺談藝術的意義論文05-28