楊敬一， 何 蕭， 周秀歡， 袁佩青， 徐心茹

(華東理工大學(xué) 化工學(xué)院， 上海 200237)

咔唑類化合物在MoS2晶體表面吸附行為的第一性原理

楊敬一， 何 蕭， 周秀歡， 袁佩青， 徐心茹

(華東理工大學(xué) 化工學(xué)院， 上海 200237)

基于密度泛函理論，采用VASP方法研究了咔唑類化合物(咔唑、3-甲基咔唑和1,5-二甲基咔唑)在MoS2(0 0 -3)晶體表面的吸附行為，考察引起其吸附能差異的原因。結(jié)果表明，咔唑類化合物在MoS2晶體表面的吸附過程中為電子受體，其吸附能隨MoS2片層數(shù)的增多而減弱；由于空間位阻和電子效應(yīng)的共同作用，使得咔唑類化合物在MoS2晶體表面的吸附能隨著苯環(huán)上甲基取代數(shù)目的增加而變小；Bader電荷及高低占據(jù)軌道的分析表明，咔唑類化合物分子上N原子得電子多少的順序與它們在MoS2晶體表面的吸附能大小順序一致，吡咯環(huán)N原子上的電子流向催化劑Mo原子，同時咔唑類化合物分子的LUMO軌道接受S原子HOMO軌道的電子形成反饋鍵，使得咔唑類化合物在MoS2晶體表面穩(wěn)定吸附。

咔唑類化合物； MoS2晶體； 吸附； 密度泛函理論； Bader電荷

原油日趨重質(zhì)化劣質(zhì)化，加氫工藝成為生產(chǎn)清潔車用燃料的主要方法。在深度加氫脫硫過程中，原料油中的有機含氮化合物使加氫催化劑中毒，抑制加氫脫硫反應(yīng)活性，嚴(yán)重影響產(chǎn)品質(zhì)量，特別在生產(chǎn)國Ⅴ和國Ⅵ汽柴油時，影響尤為顯著[1-2]。

大部分原油和餾分油中的氮以有機雜環(huán)含氮化合物的形式存在，包括以喹啉、吡啶為代表的堿性含氮化合物和以咔唑、吲哚為代表的非堿性含氮化合物[3]。國內(nèi)外眾多學(xué)者基于密度泛函理論(DFT)研究了含氮化合物在加氫脫硫催化劑表面的吸附作用。Sun等[4-5]利用密度泛函理論計算了吡啶、喹啉、吖啶、吡咯、吲哚和咔唑在NiMoS晶體表面的吸附能，結(jié)果表明，堿性氮化物在晶體表面以垂直吸附為主，而非堿性氮化物只能平行吸附在晶體表面，且堿性含氮化合物的吸附能大于非堿性含氮化合物。孫煒等[6]同樣針對NiMoS晶體模型，選取廣義梯度近似(GGA)中的PW91交換關(guān)聯(lián)泛函計算了吡啶、喹啉、咔唑和吲哚的吸附構(gòu)型和吸附能，結(jié)果表明，吡啶和喹啉是通過Ni-N鍵垂直吸附在晶體表面，而吲哚通過吡咯環(huán)的β-C鍵吸附在晶體表面，咔唑則主要通過苯環(huán)吸附在晶體表面，吲哚和咔唑均以平行吸附為主。Abdallah等[7]采用DFT理論研究了吡咯在Mo(110)表面的平行、垂直和傾斜吸附，結(jié)果表明，吡咯以平行吸附構(gòu)型為主，不可能垂直吸附在Mo表面；多個吡咯吸附在Mo表面時，以傾斜吸附為主。Temel等[8]利用DFT理論結(jié)合掃描隧道顯微鏡(STM)，研究了吡啶分子在MoS2團(tuán)簇表面的吸附作用，結(jié)果表明，吡啶分子本身與MoS2團(tuán)簇的相互作用很弱，但在有H2存在的情況下，吡啶分子在MoS2團(tuán)簇邊緣穩(wěn)定吸附。Ren等[9]為了深入研究加氫脫氮過程，基于DFT理論考察了吡啶和吡咯在α-Mo2C(0001)表面的吸附，結(jié)果顯示吡啶和吡咯以大π鍵的形式吸附在3個Mo原子形成的三角型的位置。

筆者所在課題組前期研究表明，各種石化及煤焦油柴油餾分中咔唑類衍生物含量較高[10]，且在咔唑、喹啉、吲哚、吡啶、苯胺等5種含氮化合物中，咔唑在MoS2表面的吸附能最大[11]。因此筆者選擇咔唑及含量較高的2種咔唑類化合物3-甲基咔唑、1,5-二甲基咔唑，用Materials Studio 7.0構(gòu)建晶體結(jié)構(gòu)后，基于DFT理論采用Vienna ab initio Simulation Package(VASP)軟件包進(jìn)行了咔唑類化合物在MoS2表面吸附過程的優(yōu)化計算，再通過GAUSSIAN 03軟件對吸附分子全局參量(高低占據(jù)軌道能量、電負(fù)性和親電指數(shù))的計算揭示了引起其吸附能力差異的原因，以期為研究含氮化合物對加氫脫硫反應(yīng)影響的作用機理提供理論依據(jù)。

1 MoS2晶體模型和密度泛函計算方法

MoS2晶體模型取自Materials Studio 7.0的結(jié)構(gòu)數(shù)據(jù)庫，由Build模塊建立MoS2(0 0 -3)表面，采用(4×4)超晶胞模擬，晶胞間真空層厚度大于1 nm，圖1為優(yōu)化前后咔唑在3層MoS2表面吸附示意圖。采用DFT理論的VASP程序包進(jìn)行了表面及整個吸附過程的結(jié)構(gòu)優(yōu)化，計算了咔唑類化合物分子在不同層數(shù)的MoS2晶體表面的吸附能。周期性邊界條件的計算使用Projector-augmented Wave(PAW)結(jié)合PW91泛函，其中平面波基組的能量截斷為450 eV。咔唑類化合物分子在晶體表面的吸附能△E通過式(1)計算。

△E=Et-Es-Em

(1)

式(1)中，Et為分子吸附晶體穩(wěn)定體系的總能量，eV；Es為晶體晶簇結(jié)構(gòu)的能量，eV；Em為自由吸附分子的能量，eV。△E為正，表明吸附作用不能自發(fā)進(jìn)行，需要外加能量；△E為負(fù)，則表明吸附為放熱過程，吸附過程可以自發(fā)進(jìn)行，且△E絕對值越大，吸附越穩(wěn)定。

圖1 咔唑在3層MoS2晶體表面吸附模型Fig.1 Structure model of carbazole on 3 slabs of MoS2 crystalBlue—Molybdenum atom; Yellow—Sulfur atomA, B and C are x, y and z axises.(a) Before optimization； (b) After optimization

采用GAUSSIAN 03軟件包，用DFT理論中的雜化B3LYP方法，結(jié)合6-31G+(d,p)基組，計算了咔唑類化合物分子的最高占據(jù)軌道能量EHOMO和最低空軌道能量ELUMO，根據(jù)得到的EHOMO和ELUMO，通過公式(2)和(3)計算了咔唑類化合物分子的電負(fù)性χ和親電指數(shù)ω[12-13]。

(2)

ω=(EHOMO+ELUMO)2/2(ELUMO-EHOMO)

(3)

2 結(jié)果與討論

2.1 MoS2片層數(shù)對咔唑類化合物吸附能的影響

由文獻(xiàn)[4]和筆者所在課題組前期所做工作可知，咔唑是以平行吸附形式吸附在單片層MoS2晶體表面的[11]，筆者在此基礎(chǔ)上考察了咔唑類化合物在不同片層數(shù)的MoS2催化劑上的平行吸附行為。咔唑類化合物均是含有1個N原子、甲基取代苯環(huán)不同位置的H原子形成的衍生物，具有大的π-共軛剛性平面結(jié)構(gòu)。由圖1可見，優(yōu)化后咔唑分子整體從MoS2晶體表面向上移動，2個苯環(huán)距離催化劑表面較近，而咔唑的中心吡咯環(huán)距離催化劑表面較遠(yuǎn)，即苯環(huán)與MoS2表面的作用明顯強于吡咯環(huán)，說明咔唑分子主要通過苯環(huán)吸附在MoS2晶體表面，這一結(jié)果與文獻(xiàn)[4,6]所得的結(jié)論一致。

表1列出了咔唑分子在不同層數(shù)MoS2晶體表面吸附前后鍵長的變化以及吸附能。咔唑分子各原子編號如圖2所示。由表1可見，固定1層MoS2時，咔唑分子在其表面吸附后的C—C、N—C鍵比其在真空中的鍵長略??；當(dāng)MoS2層數(shù)增加到2時，由于底層原子固定，上層原子弛豫，使吸附后咔唑分子的鍵長顯著減小，同時咔唑在2層MoS2表面的吸附能絕對值為1.73 eV，也顯著小于其在1層MoS2表面的吸附能絕對值5.18 eV；當(dāng)MoS2層數(shù)增加至3時，底層原子固定，上面兩層原子弛豫，吸附后咔唑的鍵長略小于其在2層MoS2表面吸附后的鍵長，而咔唑在其表面的吸附能絕對值也相應(yīng)略小于2層MoS2時的吸附能絕對值。這表明催化劑弛豫原子層數(shù)對吸附分子的結(jié)構(gòu)影響較大，而固定層數(shù)的影響較小。隨著弛豫原子層數(shù)增多，吸附后吸附分子的鍵長變短，鍵能變大，其與催化劑表面相互作用趨于穩(wěn)定，從而使其在MoS2表面的吸附能也趨于恒定?？紤]計算效率，筆者選擇3層MoS2模型研究咔唑類化合物在其表面的吸附。趙巍等[14]計算了水分子在3層和4層Fe表面的吸附能，發(fā)現(xiàn)相對誤差不超過0.08%，印證了計算的合理性。

咔唑、3-甲基咔唑和1,5-二甲基咔唑在3層MoS2晶體表面的吸附能分別為-1.30 eV、-1.24 eV、-0.93 eV，由文獻(xiàn)[5]可知，單環(huán)含氮化合物吡啶和吡咯在NiMoS催化劑上不同吸附構(gòu)型的吸附能分別在-1.29～-0.10 eV和-0.71～-0.68 eV范圍內(nèi)，比較可知筆者計算所得的吸附能結(jié)果較為合理。

圖2 咔唑分子各原子編號示意圖Fig.2 Schematic of the atomic number of carbazole molecule

表1 咔唑在不同層數(shù)MoS2上吸附前后的鍵長和吸附能Table 1 Bond distances in free carbazole and adsorbed carbazole on different slabs of MoS2

2.2 甲基取代基對咔唑類化合物吸附能的影響

隨著咔唑類化合物苯環(huán)上甲基取代數(shù)目增多，其在MoS2晶體表面的吸附能由大到小的順序為咔唑、3-甲基咔唑、1,5-二甲基咔唑。當(dāng)咔唑分子處于自由狀態(tài)時，由于N的吸電子作用，離N原子較遠(yuǎn)的C2—C3和C3—C4鍵長較長，鍵能較小。吸附后，咔唑類化合物的2個苯環(huán)優(yōu)先吸附在MoS2表面，吡咯環(huán)凸起于催化劑表面，但吡咯環(huán)是1個剛性平面結(jié)構(gòu)的五元環(huán)，由于分子張力的作用，C11—C12鍵長最長，鍵能最小。因此咔唑分子在束縛狀態(tài)下，C11—C12鍵很可能先發(fā)生斷裂。當(dāng)苯環(huán)上有甲基取代基后，由于空間位阻效應(yīng)阻礙了吸附分子上苯環(huán)在MoS2晶體表面的吸附，使得吸附能變小。文獻(xiàn)[15]研究發(fā)現(xiàn)，二苯并噻吩(DBT)系列在NiMoS催化劑上平行吸附時，隨著DBT分子上氫原子數(shù)目增多，減少了吸附分子與催化劑表面的相互作用點，使得吸附能隨之減小。這一變化規(guī)律與本文研究結(jié)果類似。

計算了咔唑類化合物分子的最高占據(jù)軌道能EHOMO和最低非占據(jù)軌道能ELUMO，并由此得到咔唑類化合物分子的電負(fù)性χ和親電指數(shù)ω，電負(fù)性和親電指數(shù)是反映分子在反應(yīng)中吸引電子能力強弱的參數(shù)。表2為咔唑類化合物的全局參量。由表2可知，作為供電子基團(tuán)，隨著咔唑類化合物苯環(huán)上甲基取代數(shù)目增多，分子的電負(fù)性和親電指數(shù)有所減小，這可能是由于甲基的斥電子作用使得咔唑類化合物分子吸電子能力減弱，相應(yīng)的，其在MoS2表面的吸附能隨之變小，吸附越不穩(wěn)定。這表明咔唑類化合物在MoS2晶體表面的吸附過程中為電子受體，同時也從另一角度闡述了甲基取代基對吸附能的影響。綜上所述，隨著咔唑類化合物分子中苯環(huán)上甲基取代數(shù)的增加，由于空間位阻和電子效應(yīng)的共同作用，使其在MoS2晶體表面的吸附能變小。

表2 咔唑類化合物的全局參量Table 2 Part of global parameters of carbazole compounds

2.3 電子結(jié)構(gòu)對咔唑類化合物吸附能的影響

為了深入考察咔唑類化合物分子在MoS2表面的吸附能與原子間電荷轉(zhuǎn)移的關(guān)系，采用Bader電荷分析方法，通過在一個區(qū)域中整合電子密度從而量化了原子間電荷的轉(zhuǎn)移情況。表3列出了咔唑類化合物分子中C原子、N原子的Bader電荷數(shù)和轉(zhuǎn)移電荷數(shù)(凈電荷)，凈電荷為正表明該原子失電子，凈電荷為負(fù)表明該原子得電子。

表3 咔唑類化合物的Bader電荷Table 3 Bader charges of carbazole compounds

從表3可以看出，咔唑類化合物分子中電荷轉(zhuǎn)移比較顯著的位置是N原子和與其相鄰的2個C原子(10號C原子和13號C原子)。3種咔唑類化合物分子上C、N原子的電子轉(zhuǎn)移情況基本一致：N原子得電子，與其相鄰的10位和13位的C原子失電子，其他位的C原子電子基本不發(fā)生轉(zhuǎn)移。Sun等[16]認(rèn)為，對于咔唑分子，其親電反應(yīng)的活性中心在N原子上，親核反應(yīng)的活性中心分布在苯環(huán)上，這與筆者的研究結(jié)果相吻合。3種分子的N原子得電子由多到少的順序為咔唑、3-甲基咔唑、1,5-二甲基咔唑，與它們在MoS2晶體表面的吸附能大小順序一致，這一結(jié)果與Shiraishi等[17]研究的甲基對咔唑中N原子電子密度的變化規(guī)律類似。這進(jìn)一步證明了咔唑類化合物在與MoS2催化劑吸附過程中為電子受體。

圖3為咔唑類化合物分子的最高占有軌道(HOMO)和最低未占軌道(LUMO)的0.02等值面圖。由圖3可以看出，3種分子的 HOMO 軌道有大致相同的分布中心，對應(yīng)的吡咯環(huán)內(nèi)的π成鍵軌道分布一致，并且有N原子的p軌道貢獻(xiàn)。3種分子的LUMO軌道主要離域在吡咯環(huán)上，是典型的π反鍵軌道分布。這種分布使得咔唑類化合物分子N原子上的孤對電子轉(zhuǎn)移至催化劑Mo原子的缺電子中心形成配位鍵，同時吸附分子利用其LUMO軌道接受來自催化劑S原子HOMO軌道的電子形成反饋鍵，有利于咔唑類化合物分子在MoS2晶體表面穩(wěn)定吸附。

圖3 咔唑類化合物分子的最高占有軌道(HOMO)和最低未占軌道(LUMO)的0.02等值面圖Fig.3 The 0.02 isosurface figure of carbazole compounds molecular HOMO and LUMO The above three pictures are HOMO orbits and the following are LUMO orbits.(a) Carbazole; (b) 3-Methylcarbazole; (c) 1,5-Dimethylcarbazole

3 結(jié) 論

咔唑類化合物(咔唑、3-甲基咔唑和1,5-二甲基咔唑)在MoS2(0 0 -3)晶體表面的平行吸附能隨著MoS2晶體片層數(shù)增多而減弱。由于空間位阻和電子效應(yīng)的共同作用，使得咔唑類化合物在MoS2晶體表面的吸附能隨著苯環(huán)上甲基取代數(shù)目的增加而變小。咔唑類化合物在MoS2晶體表面的吸附過程中為電子受體，吡咯環(huán)N上的電子流向催化劑Mo原子，同時咔唑類化合物分子的LUMO軌道接受S原子HOMO軌道的電子形成反饋鍵，使得咔唑類化合物在MoS2晶體表面穩(wěn)定吸附。

[1] 邵志才， 高曉冬， 李皓光， 等. 氮化物對柴油深度和超深度加氫脫硫的影響I氮化物含量的影響[J].石油學(xué)報(石油加工)， 2006， 22(4)： 12-17. (SHAO Zhicai, GAO Xiaodong, LI Haoguang, et al. Effect of nitrogen-containing compounds on deep and ultradeep HDS for diesel oilⅠEffect of the content of nitrogen-containing compounds[J].Acta Petrolei Sinica (Petroleum Processing Section), 2006, 22(4): 12-17.)

[2] 張生， 李翔， 王安杰， 等. 有機含氮化合物對深度加氫脫硫反應(yīng)的影響[J].工業(yè)催化， 2005， 13(11)： 12-17. (ZHANG Sheng, LI Xiang, WANG Anjie, et al. Influence of organic nitrogen compounds on deep hydrodesulfurization(HDS)[J].Industrial Catalysis, 2005, 13(11): 12-17.)

[3] SHIRAISHI Y, TACHIBANA K, HIRAI T, et al. A novel desulfurization process for fuel oils based on the formation and subsequent precipitation of S-alkylsulfonium salts. 5. Denitrogenation reactivity of basic and neutral nitrogen compounds[J].Industrial & Engineering Chemistry Research, 2001, 40(22): 4919-4924.

[4] SUN M, NELSON A E, ADJAYE J. First principles study of heavy oil organonitrogen adsorption on NiMoS hydrotreating catalysts[J].Catalysis Today, 2005, 109(1): 49-53.

[5] SUN M, NELSON A E, ADJAYE J. Adsorption and hydrogenation of pyridine and pyrrole on NiMoS: An ab initio density-functional theory study[J].Journal of Catalysis, 2005, 231(1): 223-231.

[6] 孫煒， 方元， 王鵬， 等. 有機氮化物在鎳鉬硫催化劑的表面吸附[J].武漢工程大學(xué)學(xué)報， 2013， 35(12)： 35-39. (SUN Wei, FANG Yuan, WANG Peng, et al. Adsorption of organonitrogen on surface of nickel-molybdenum sulfide catalyst[J].Journal of Wuhan Institute of Technology, 2013, 35(12): 35-39.)

[7] ABDALLAH W A, NELSON A E. Density functional theory study of pyrole adsorption on Mo(110)[J].The Journal of Physical Chemistry B, 2005, 109(21): 10863-10870.

[8] TEMEL B, TUXEN A K, KIBSGAARD J, et al. Atomic-scale insight into the origin of pyridine inhibition of MoS2-based hydrotreating catalysts[J].Journal of Catalysis, 2010, 271(2): 280-289.

[9] REN J, WANG J, HUO C F, et al. Adsorption of NO, NO2, pyridine and pyrrole onα-Mo2C(0001): A DFT study[J].Surface Science, 2007, 601(6): 1599-1607.

[10] 楊敬一， 周秀歡， 蔡海軍， 等. 煤焦油和石油基柴油餾分中含氮化合物的分離鑒定[J].石油煉制與化工， 2015， 46(7)： 107-112. (YANG Jingyi, ZHOU Xiuhuan, CAI Haijun, et al. Separation and identification of nitrogen compounds in diesel fraction of coal tar and petroleum[J].Petroleum Processing and Petrochemicals, 2015, 46(7): 107-112.)

[11] 楊敬一， 周秀歡， 袁佩青， 等. 幾種含氮化合物在MoS2晶面的吸附研究[J].計算機與應(yīng)用化學(xué)， 2015， 32(7)： 769-773. (YANG Jingyi， ZHOU Xiuhuan， YUAN Peiqing， et al. The adsorption of nitrogen compounds on MoS2[J].Computers and Applied Chemistry， 2015， 32(7)： 769-773.)

[12] GEERLINGS P, DE PROFT F, LANGENAEKER W. Conceptual density functional theory[J].Chemical Reviews, 2003, 103(5): 1793-1874.

[13] CHERMETTE H. Chemical reactivity indexes in density functional theory[J].Journal of Computational Chemistry, 1999, 20(1): 129-154.

[14] 趙巍， 汪家道， 劉峰斌， 等. H2O分子在Fe(100)， Fe(110)， Fe(111)表面吸附的第一性原理研究[J].物理學(xué)報， 2009， 58(5)： 3352-3358. (ZHAO Wei, WANG Jiadao, LIU Fengbin, et al. First principles study of H2O molecule adsorption on Fe(100), Fe(110) and Fe(111) surfaces[J].Acta Physica Sinica, 2009, 58(5): 3352-3358.)

[15] 孫煒， 王鵬， 楊犁， 等. 二苯并噻吩加氫脫硫過程在鎳鉬硫催化劑表面的吸附[J].武漢工程大學(xué)學(xué)報， 2013， 35(2)： 51-56. (SUN Wei, WANG Peng, YANG Li, et al. Adsorption of dibenzothiophene hydrodesulfurization process over NiMoS catalyst[J].Journal of Wuhan Institute of Technology, 2013, 35(2): 51-56.)

[16] SUN M, NELSON A E, ADJAYE J. Correlating the electronic properties and HDN reactivities of organonitrogen compounds: An ab initio DFT study[J].Journal of Molecular Catalysis A: Chemical, 2004, 222(1): 243-251.

[17] RAYBAUD P, HAFNER J, KRESSE G, et al. Structure, energetics, and electronic properties of the surface of a promoted MoS2catalyst: An ab initio local density functional study[J].Journal of Catalysis, 2000, 190(1): 128-143.

First Principles of the Adsorption of Carbazole Compounds on MoS2

YANG Jingyi， HE Xiao， ZHOU Xiuhuan， YUAN Peiqing， XU Xinru

(ChemicalEngineeringCollege,EastChinaUniversityofScienceandTechnology,Shanghai200237,China)

The adsorption behavior and adsorption energy of three carbazole compounds (carbazole, 3-methylcarbazole, 1,5-dimethylcarbazole) on MoS2cluster were investigated by the VASP package based on the Density functional theory (DFT). The results showed that carbazole compounds were electron acceptors in the process of adsorption on the surface of MoS2and the adsorption energy of carbazole compounds on MoS2reduced along with growing amount of MoS2slabs. As the amount of methyl substitution on benzene ring increased, the adsorption energy of carbazole compounds on MoS2reduced due to the steric hindrance of methyl and electronic effect. According to the Bader charge, HOMO and LUMO analysis, the adsorption energy of carbazole compounds on MoS2increased along with increased amount of electrons on N atom. It was observed that the lone pair electrons on N atom of adsorbed pyrrole ring are transferred to the electron-deficient center on Mo atoms, and the LUMO orbit of carbazole compounds accepts electrons in the HOMO orbit of S atoms, leading to the stable adsorption of carbazole compounds on the surface of MoS2.

carbazole compounds; MoS2cluster; adsorption; density functional theory; Bader charge

2016-07-25

國家自然科學(xué)基金項目(21376075)資助

楊敬一，男，副教授，博士，從事石油與能源化工方面的研究；E-mail：jyyang@ecust.edu.cn

1001-8719(2017)03-0515-06

TQ015.9

A

10.3969/j.issn.1001-8719.2017.03.016