五角石墨烯吸附H2S的第一性原理計(jì)算

孟 浩， 李得生， 王向嶺， 王軍凱

(河南理工大學(xué) 材料科學(xué)與工程學(xué)院 河南省深地材料科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室， 焦作 454003)

1 引 言

隨著工業(yè)的發(fā)展，油氣開采和大量化石燃料的使用，造成了大量的有毒氣體被排入大氣，不僅嚴(yán)重危害人體健康而且還限制了工業(yè)的進(jìn)一步發(fā)展.其中硫化氫作為一種神經(jīng)毒物和腐蝕性氣體污染物，其廣泛存在于生活污水和工業(yè)生產(chǎn)中.一方面，硫化氫的存在會(huì)嚴(yán)重傷害人體機(jī)能，低濃度的硫化氫會(huì)損傷人的嗅覺，導(dǎo)致細(xì)胞缺氧，引起呼吸困難、惡心、頭暈和四肢乏力[1]，而高濃度的硫化氫則會(huì)直接導(dǎo)致人體窒息甚至死亡[2].另一方面，硫化氫在潮濕環(huán)境下還具有強(qiáng)烈的腐蝕性，會(huì)對(duì)金屬制品產(chǎn)生嚴(yán)重腐蝕，造成以金屬為材質(zhì)的設(shè)備和管道使用壽命降低[3, 4].而硫化氫的燃燒產(chǎn)物SO2、SO3等酸性氣體也會(huì)對(duì)空氣造成嚴(yán)重的污染[5]，因此設(shè)計(jì)一種有效的H2S吸附材料是非常重要的.

五角石墨烯是一種新型的石墨烯類似物，最早由張等在2015年發(fā)現(xiàn)，這是一個(gè)只有五元環(huán)的準(zhǔn)二維亞穩(wěn)碳結(jié)構(gòu)[6].五角石墨烯具有許多有趣的特性，包括機(jī)械穩(wěn)定性、動(dòng)態(tài)穩(wěn)定性、熱穩(wěn)定性、高理想強(qiáng)度、負(fù)泊松比和寬的準(zhǔn)直接帶隙[6-9]，因而不需要像石墨烯那樣通過化學(xué)或物理的修飾來打開帶隙，因此理論上更適合用來吸附氣體.目前已經(jīng)有部分學(xué)者研究了五角石墨烯對(duì)氣體分子的吸附性能.Lu等人[10]利用密度泛函理論(DFT)研究了五角石墨烯對(duì)燃燒后氣體中CO2的捕獲和分離，結(jié)果表明五角石墨烯具有高效的CO2捕獲和分離能力，是一種很有前途的CO2捕獲和分離材料.Enriquez等人[11]采用第一性原理的方法研究發(fā)現(xiàn)了五角石墨烯對(duì)H2吸附性能優(yōu)于石墨烯，H2在五角石墨烯上的吸附能(-0.092 eV)高于石墨烯上的吸附能(-0.067 eV)[12].這些研究都表明了在氣體吸附能力上五角石墨烯具有比石墨烯更好的效果.

研究表明，元素?fù)诫s是提高二維材料吸附性能的有效手段.Chen等人[13]采用密度泛函理論(DFT)研究了SF6分解產(chǎn)物在原始和B/N摻雜五角石墨烯上的吸附，結(jié)果表明原始五角石墨烯對(duì)SF6分解產(chǎn)物各分子的吸附能較小，而B和N原子摻雜后，吸附能均有不同程度的提高.Zhang等人[14]采用第一性原理計(jì)算的方法研究了CO和CO2在Fe摻雜五角石墨烯上的吸附行為，結(jié)果表明，相對(duì)于純五角石墨烯，F(xiàn)e的摻雜可以有效地增強(qiáng)CO和CO2氣體的吸附.Sathishkumar 等人[15]采用DFT計(jì)算的方法研究了氫分子在原始、B和N摻雜五角石墨烯上的吸附，結(jié)果表明將硼或氮加入到五角石墨烯薄片中，可以顯著提高H2分子和H原子在五角石墨烯上的吸附能.Villagracia等人[16]采用第一性原理計(jì)算的方法研究了3d過渡金屬(Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu和Zn)的摻雜對(duì)五角石墨烯氣敏性能的影響，結(jié)果表明，相對(duì)于原始的五角石墨烯，摻雜過渡金屬元素后，其帶隙由寬到窄，顯著增強(qiáng)了其對(duì)H2分子的吸附.Feng等人[17]利用DFT計(jì)算的方法研究了NO在鋁摻雜五角石墨烯上的吸附行為，結(jié)果表明NO氣體分子的取向?qū)ξ叫袨橛幸欢ǖ挠绊?，與原始五角石墨烯相比，鋁摻雜五角石墨烯體系具有更高的吸附能和更短的吸附距離.

然而目前尚缺乏五角石墨烯對(duì)H2S氣體吸附行為的研究，因此本文采用第一性原理計(jì)算的方法研究了五角石墨烯/摻雜五角石墨烯對(duì)H2S氣體的吸附行為，旨在為設(shè)計(jì)研發(fā)基于五角石墨烯的H2S吸附材料提供理論指導(dǎo).

2 計(jì)算方法

文中所有計(jì)算均采用Materials Studio軟件中基于密度泛函理論的DMol3軟件包[18,19]，電子和電子相互作用的交換關(guān)聯(lián)能函數(shù)采用的是廣義梯度近似(GGA)的PBE交換關(guān)聯(lián)泛函[20]，軌道占據(jù)參數(shù)采用0.005 Ha.布里淵區(qū)的特殊點(diǎn)積分使用Monkhorst-Pack方法[21].自洽場(chǎng)(SCF)計(jì)算收斂標(biāo)準(zhǔn)為電子密度變化小于1×10-6Ha，能量的收斂標(biāo)準(zhǔn)為1×10-5Ha，Hellmann-Feynman原子收斂判據(jù)為每個(gè)原子受力不大于2×10-3Ha/?，最大位移為5×10-3nm.

計(jì)算模型采用含96個(gè)C原子的4×4×1五角石墨烯超晶胞，真空層厚度設(shè)為30 ?.其結(jié)構(gòu)如圖1所示.

圖1 原始五角石墨烯晶胞幾何優(yōu)化結(jié)構(gòu)(金色框代表PG單元，綠色框代表4×4超晶胞，其中C1代表的是sp2C原子，C2代表的是sp3C原子，兩個(gè)sp2C原子之間的距離d2是1.341 ?，sp2C和sp3C原子之間的距離d2是1.551 ?，側(cè)視圖中相鄰兩個(gè)原子之間的距離d3為1.216 ?)Fig. 1 Geometrically optimized structures of pristine pentagonal graphene. (the golden frame represents the PG unit, the green frame represents the 4 × 4 supercell, C1 represents the sp2C atom, while C2 represents the sp3C atom, the distance d2 between two sp2C atoms is 1.341 ?, the distance d2 between sp2C and sp3C is 1.551 ? and the distance d3 between two adjacent atoms in the side view is 1.216 ?.)

為了定量描述氣體分子在純五角石墨烯(penta-graphene)表面以及元素?fù)诫s五角石墨烯(doped penta-graphene)表面的吸附強(qiáng)度，我們用吸附能Ead來定義：

Ead=Egas+penta-graphene-Epenta-graphene-Egas

(1)

Ead=Egas+doped penta-graphene-Edoped penta-graphene-Egas

(2)

其中，Epenta-graphene、Edoped penta-graphene、Egas、Egas+penta-graphene以及Egas+ doped penta-graphene分別代表純五角石墨烯、元素?fù)诫s五角石墨烯、單個(gè)小氣體分子、整個(gè)純五角石墨烯吸附體系以及元素?fù)诫s五角石墨烯吸附體系的能量.根據(jù)公式(1)、(2)中Ead的定義，Ead數(shù)值越負(fù)，吸附作用越強(qiáng).采用Mulliken方法[22]分析襯底與吸附物之間的電荷轉(zhuǎn)移.

3 結(jié)果與討論

3.1 硫化氫在純五角石墨烯表面的吸附

首先研究了純五角石墨烯對(duì)硫化氫氣體分子的吸附行為，考慮了硫化氫分子的不同吸附位點(diǎn)和取向.圖2為硫化氫分子的5種可能取向模型：氫原子與硫原子的連線垂直于五角石墨烯所在平面朝下(a)，氫原子與硫原子的連線垂直于五角石墨烯所在平面朝上(b)，硫化氫分子所形成的平面與五角石墨烯平行(c)，硫化氫分子所形成的平面垂直于五角石墨烯所在平面成“V”字型(d)，硫化氫分子所形成的平面垂直于五角石墨烯所在平面成倒“V”字型(e).可能的吸附位點(diǎn)分別是：sp2C原子上方的頂位點(diǎn)(T1)，sp3C原子上方的頂位點(diǎn)(T2)，sp2C-sp3C鍵中心位置上方的橋位點(diǎn)(B0)，sp2C-sp2C凸鍵中心位置上方的橋位點(diǎn)(B1), sp2C-sp2C凹鍵中心位置上方的橋位點(diǎn)(B2)，五邊形中心位置上方的空穴位點(diǎn)(H).表1為30種吸附構(gòu)型對(duì)應(yīng)的吸附能計(jì)算結(jié)果.可以看出，吸附能變化范圍為-0.292 eV到-0.079 eV.該結(jié)果表明，在純五角石墨烯表面，無論硫化氫是何種取向，吸附在何種位置，吸附都比較弱，因此純五角石墨烯不適合作為硫化氫氣體分子的吸附材料.

圖2 硫化氫分子的初始吸附構(gòu)型在五角石墨烯上的不同取向: 氫原子與硫原子的連線垂直于五角石墨烯所在平面朝下(a)，氫原子與硫原子的連線垂直于五角石墨烯所在平面朝上 (b)，硫化氫分子所形成的平面與五角石墨烯平行 (c)，硫化氫分子所形成的平面垂直于五角石墨烯所在平面成“V”字型(d), 硫化氫分子所形成的平面垂直于五角石墨烯所在平面成倒“V”字型(e)Fig. 2 Different orientations of initial adsorption configuration of hydrogen sulfide molecule on pentagonal graphene: the connecting line of hydrogen atom and sulfur atom is perpendicular to the plane of pentagonal graphene and faces downward (a) ，the connecting line of hydrogen atom and sulfur atom is perpendicular to the plane of pentagonal graphene and upward (b)，the plane formed by hydrogen sulfide molecule is parallel to pentagonal graphene (c)，the plane formed by hydrogen sulfide molecule is perpendicular to the plane of pentagonal graphene, forming a "V" shape (d)，the plane formed by hydrogen sulfide molecule is perpendicular to the plane of pentagonal graphene, and forms an inverted "V" shape (e)

3.2 非金屬摻雜對(duì)硫化氫在五角石墨烯表面吸附的影響

為了增強(qiáng)五角石墨烯對(duì)硫化氫分子的吸附能力，首先采用氮(N)、磷(P)和硫(S)三種非金屬元素對(duì)純五角石墨烯進(jìn)行摻雜(分別考慮了sp2C和sp3C兩種摻雜位置，如圖3所示)，研究了非金屬摻雜五角石墨烯對(duì)硫化氫氣體分子吸附的影響.其中，硫化氫氣體分子初始吸附位置為摻雜非金屬原子的上方.對(duì)應(yīng)的吸附能計(jì)算結(jié)果如表2所示.從表中可以看出，對(duì)所有非金屬摻雜而言，無論硫化氫是何種取向，吸附在何種位置，摻雜五角石墨烯對(duì)硫化氫氣體分子的吸附能都很低，電荷轉(zhuǎn)移量趨近于0.該結(jié)果表明：在非金屬摻雜的五角石墨烯表面，無論硫化氫是何種取向，吸附在何種位置，吸附都比較弱，即氮(N)、磷(P)、硫(S)摻雜的五角石墨烯都不適合作為硫化氫的吸附材料.

表1 優(yōu)化后不同構(gòu)型的硫化氫在純五角石墨烯表面上的吸附能和電荷轉(zhuǎn)移量

表2 硫化氫在非金屬摻雜五角石墨烯表面的吸附能和電荷轉(zhuǎn)移

圖3 (a)sp2C位置摻雜的模型代表圖(b) sp3C位置摻雜的模型代表圖 (X代表N、P及S)Fig. 3 (a) model representation of sp2c site doping diagram.(b) model representation of sp3c site doping diagram. (X represents N, P and S)

3.3 金屬摻雜對(duì)硫化氫在五角石墨烯表面吸附的影響

本部分采用鈷(Co)、鎳(Ni)、銅(Cu)及鈦(Ti) 四種金屬元素?fù)诫s五角石墨烯(分別考慮了sp2C和sp3C兩種摻雜位置)，研究了金屬摻雜五角石墨烯對(duì)硫化氫氣體分子吸附的影響.其中，硫化氫氣體分子初始吸附位置為摻雜金屬原子的正上方.對(duì)應(yīng)的吸附能和電荷轉(zhuǎn)移計(jì)算結(jié)果如表3所示.從表中可以看出，當(dāng)金屬原子摻雜在五角石墨烯的sp2C位置時(shí)，摻雜五角石墨烯對(duì)硫化氫氣體分子的吸附能更負(fù)，電荷轉(zhuǎn)移量更多.該結(jié)果表明金屬原子替換五角石墨烯的sp2C時(shí)更有利于硫化氫氣體分子的吸附.

對(duì)比分析了不同金屬元素?fù)诫s的五角石墨烯對(duì)硫化氫氣體分子的吸附行為，結(jié)果表明，四種金屬元素?fù)诫s五角石墨烯sp2C后，對(duì)硫化氫氣體分子的吸附能變化范圍為-1.185到-1.824 eV，為強(qiáng)的化學(xué)吸附.這表明，鈷(Co)、鎳(Ni)、銅(Cu)和鈦(Ti)摻雜的五角石墨烯可能都適合作為硫化氫的吸附材料.進(jìn)一步對(duì)比發(fā)現(xiàn)，不同金屬摻雜的五角石墨烯對(duì)硫化氫氣體分子的平均吸附能從小到大為Ead(Cu)

此外，以Cu摻雜(sp2C位置)五角石墨烯為例，分析了硫化氫分子不同吸附取向(圖4)對(duì)其在摻雜石墨烯表面吸附行為的影響.可以看出，無論硫化氫分子是何種取向進(jìn)行吸附，其吸附能數(shù)值均在-1.535和-1.603 eV之間微弱變化，都屬于化學(xué)吸附.這充分表明，Cu摻雜五角石墨烯是探測(cè)/吸附硫化氫氣體分子的理想材料.

4 結(jié) 論

采用第一原理計(jì)算研究了純五角石墨烯和元素?fù)诫s五角石墨烯對(duì)H2S氣體的吸附性能.得到的主要結(jié)論如下：

(1)純五角石墨烯對(duì)H2S氣體僅有微弱的物理吸附作用，無法直接用于吸附H2S氣體.

(2)采用非金屬N、P及S摻雜的五角石墨烯對(duì)H2S氣體的吸附也屬于較弱的物理吸附，仍然無法用于吸附H2S氣體.

(3)采用金屬Co、Ni、Cu及Ti摻雜的五角石墨烯對(duì)H2S氣體的吸附作用顯著增強(qiáng).當(dāng)金屬摻雜在sp2C位置時(shí)，摻雜五角石墨烯對(duì)H2S氣體的吸附效果較好，此時(shí)的吸附均為化學(xué)吸附.

(4)Co、Ni、Cu及Ti摻雜的五角石墨烯均可作為H2S氣體的傳感/捕獲材料.其中Cu摻雜五角石墨烯吸附H2S氣體的效果最好.H2S分子的不同吸附取向?qū)ζ湮叫袨橛绊懖淮?