英文原題:Identifying Photocatalytic Active Sites of C2H6 C?H Bond Activation on TiO2 via Combining First-Principles Ground-State and Excited-State Electronic Structure Calculations

通訊作者:郭慶,南方科技大學;胡偉,楊金龍,中國科學技術大學
作者:王曉寧,萬凌云,王子健,劉小峰,高蘊智,王蕾,劉杰,郭慶,胡偉,楊金龍
研究背景
TiO2表面的光催化反應對于太陽能驅動技術至關重要,而活性位點和過程識別仍然存在爭議。只有確定這些基本的光催化過程,才能更好地理解TiO2光催化作用,這對于開發新的光催化劑和光催化過程的表征至關重要。實驗中對光催化反應中吸附位點的影響以及光參與反應的步驟進行了大量研究,但始終沒有得出統一結論。之前的計算工作對TiO2表面光催化分子鍵的活化提供了許多有價值的見解,但仍然不足,因為這些計算大多數是在單粒子水平進行的,而光生電子和空穴之間的相互作用,即激子效應,很少被覆蓋到。本文結合密度泛函理論(DFT)和線性響應含時密度泛函理論(LR-TDDFT),詳細研究了C2H6在金紅石TiO2(110)上的兩個典型位點的吸附行為。
文章亮點
DFT 計算和LR-TDDFT計算揭示了吸附在金紅石TiO2(110)表面上C2H6分子光活性的吸附位點差異性。C2H6更容易吸附在Ti5c位點上并且基態的脫氫能壘較低,這與單粒子水平的計算一致。在考慮激子效應后的激發態情況則相反,吸附在Obr位點的C2H6具有更高的光活性,Obr原子作為光氧化過程中提供光生空穴的活性位點。在C2H6脫氫過程中,光生電子始終位于TiO2襯底中的Ti4+位點,而光生空穴可以被C2H6捕獲以激活C-H鍵。這一工作為構建新型高效的烷烴脫氫催化劑提供了重要的思想和理論基礎。
內容介紹
本項工作首先計算了C2H6吸附在金紅石TiO2(110)上的兩個典型位點Ti5c和Obr上的基態電子結構特征。隨后,在本組自主開發的KSSOLV代碼中實現的激發態LR-TDDFT方法計算了體系的光生載流子分布(圖1c,1d)。結合VBM的軌道分布(圖1a,1b右列),O原子的Py軌道和Pz軌道對S1激發態有主要貢獻。由于載流子在襯底中產生,需要轉移到吸附分子以催化氧化還原反應,因此C2H6吸附在Obr位點時光催化活性高于在Ti5c位點吸附。

圖1(a)Ti5c位點和(b)Obr位點吸附C2H6體系的VBM、CBM的基態投影態密度圖和電子空間分布(黃色)和空穴空間分布(白色)。C2H6(g)表示C2H6與TiO2襯底間的距離為大于10 ?。(c)Ti5c位點和(d)Obr位點吸附C2H6體系的第一、第二、第三和第四激發態的電子(藍色)和空穴(綠色)在實空間中的分布情況。等值面大小為0.0002 e/bohr3。括號中的數字表示相應軌道激發的貢獻值。

圖2 C2H6吸附在(a)Ti5c和(b)Obr位點的基態脫氫過程。C2H6在(c)Ti5c位點兩步脫氫過程和(d)Obr位點第一步脫氫的基態(黑線)和激發態(紅線)能量演化。等值面大小為0.0002 e/bohr3。
在Ti5c位點上C2H6脫氫轉化為C2H4(圖2a)的活化能壘出現在第二步C-H鍵的伸長過程。對于Obr位點吸附的C2H6,兩步C-H鍵活化的能壘分別為1.17 eV和0.03 eV。較高的能壘表明室溫下很難進行第一步C-H鍵的活化,需要光的參與才能完成,而第二步脫氫能壘可以在無光參與的室溫情況下自發進行。
從基態脫氫反應來看,C2H6吸附在Ti5c位點的脫氫能壘相比Obr位點時更小。激發態計算顯示,Ti5c位點吸附的C2H6脫氫過程中由于激發態能壘TS1的出現,使第一步C-H鍵活化傾向于發生在基態。當第一個H解離后,大量空穴被吸附分子C2H5捕獲,隨后的C-H鍵活化過程在激發態完成,能壘比基態反應時減小了0.24 eV。對于Obr位點吸附的C2H6,光生空穴在初始吸附時即被C2H6捕獲,整個C-H鍵活化過程傾向于在激發態進行。C-H鍵的活化能壘顯著降低,從基態的1.17 eV降低到激發態的0.4 eV。

圖3 當C2H6吸附在(a)Ti5c位點和(b)Obr位點時的前線軌道分布。(c、d)C2H6/TiO2沿垂直表面方向的平均靜電勢分布和電荷差分密度圖。黃色等值面表示電子增加,白色等值面表示電子損耗。綠色虛線與縱坐標交點之間表示C2H6所在的電位范圍(綠色陰影區域)。
機理部分圍繞兩個吸附位點的差異性進行解釋:從基態的表面軌道分析, Ti5c位點吸附的C2H6與TiO2(110)表面的軌道重合更多,電荷轉移更多,從而吸附越穩定。由于兩個吸附位點的C2H6所處的表面電勢不同,吸附態的C2H6能級出現了上下移動進而影響了其對襯底光生電子-空穴的捕獲能力。此外,TiO2(110)表面的光生空穴主要通過Obr原子遷移到吸附分子上。
綜上所述,本文對金紅石TiO2(110)表面光催化C2H6 C-H鍵活化的吸附位點差異性進行了探究,采用KSSOLV代碼實現的LR-TDDFT方法計算了吸附體系中C2H6 分子C-H鍵活化過程的光生電子-空穴的分布,并對光參與C-H鍵活化的步驟進行了識別。該研究為構建新型高效烷烴脫氫催化劑提供了重要思路和理論基礎。
KSSOLV(Kohn-Sham Solver)是一款基于平面波基組的求解Kohn-Sham密度泛函理論的代碼,利用MATLAB編程和計算,使用解釋性語言,無需額外編譯安裝,支持Windows、Mac OS X和Linux系統,同時支持英偉達GPU加速,能夠在個人PC端達到主流計算軟件的計算速度。目前KSSOLV中已經實現了分子和固體的結構優化,電子能帶結構計算(包括自旋軌道耦合效應SOC),局域密度近似(LDA)、廣義梯度近似(GGA)和雜化泛函(HSE06)的交換關聯泛函以及激發態電子結構(LR-TDDFT和GW)計算和動力學模擬(RT-TDDFT)。KSSOLV還可以用于驗證公式和算法的正確性,其中已經實現的三種加速算法,包括自適應壓縮交換算子(ACE)方法、投影的C-DIIS(PC-DIIS)方法和插值可分密度近似(ISDF)方法,加速了雜化泛函、激發態電子結構計算和動力學模擬,有希望實現針對數百原子中等尺度的分子固體材料體系激發態光電性質的計算模擬。
原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2024/03/19/b712e81a5d/