安 然,孔 明,胡萬(wàn)權(quán),倪義東,沈永鵬,李淑君,郭元茹,*,潘清江,*

1.東北林業(yè)大學(xué) 材料科學(xué)與工程學(xué)院, 黑龍江 哈爾濱 150040;2.黑龍江大學(xué) 化學(xué)化工與材料學(xué)院, 黑龍江 哈爾濱 150080

單原子催化劑由于其優(yōu)異的催化性能而引起了科學(xué)工作者的廣泛關(guān)注[1]。目前,已采用了多種方法,諸如構(gòu)筑缺陷、負(fù)載金屬和摻雜原子等[2-4],改性和調(diào)整催化劑的微觀結(jié)構(gòu)和電子性質(zhì)。單原子催化劑可以獲得最大的原子利用效率,提高催化活性和減少消耗;高度孤立的單個(gè)金屬原子可以提供更多的活性中心來(lái)增強(qiáng)催化活性,從而促進(jìn)多相催化的發(fā)展。

利用空間限域效應(yīng)將金屬原子或離子固定在具有配位作用的空穴中是制備單原子催化劑的有效策略之一。由于易于制備、成本低、電子結(jié)構(gòu)匹配和配位空穴大小適當(dāng),石墨型氮化碳(g-C3N4)可以作為單原子催化劑的優(yōu)良基質(zhì)材料,并表現(xiàn)出優(yōu)異的半導(dǎo)體光電催化活性[5-6]。最近,Gao等[7]發(fā)現(xiàn)g-C3N4負(fù)載Pd和Pt可以制備穩(wěn)定的光催化劑,并對(duì)CO2還原具有很好的性能。Huang等[8]開(kāi)發(fā)了g-C3N4負(fù)載Co單原子催化劑,實(shí)現(xiàn)了CO2的選擇性還原。文獻(xiàn)[9]合成了氮化碳納米管支撐稀土Er單原子的復(fù)合材料,具有高效的光催化CO2性能。

基于這些研究,選擇適當(dāng)大小和電子結(jié)構(gòu)匹配的金屬或離子錨定在g-C3N4空穴中,對(duì)增強(qiáng)其催化活性具有重要的意義。在本工作中,理論設(shè)計(jì)g-C3N4(CN)負(fù)載系列U、Sc族金屬離子,形成穩(wěn)定的配合物[M(CN)]z+,其中金屬M(fèi)分別為Sc、Y、La、Ac和U,z分別為2和3。通過(guò)探索配合物幾何和電子結(jié)構(gòu)、配位空穴對(duì)金屬離子的穩(wěn)定化作用(成鍵)和單電子還原性質(zhì),以期對(duì)相關(guān)單原子催化劑合成與應(yīng)用提供可靠理論支持。

1 計(jì)算細(xì)節(jié)

計(jì)算采用Priroda程序在氣態(tài)和不使用任何對(duì)稱性條件下優(yōu)化配合物結(jié)構(gòu)。運(yùn)用標(biāo)量相對(duì)論全電子(AE)哈密頓,并同時(shí)使用廣義梯度近似的GGA-PBE(Perdew-Burke-Ernzerhof)[10]泛函和全電子Gaussian基組。以優(yōu)化的結(jié)構(gòu)為基礎(chǔ)進(jìn)行頻率計(jì)算,結(jié)果未發(fā)現(xiàn)虛頻,說(shuō)明這些結(jié)構(gòu)為勢(shì)能面上的極小值點(diǎn)。計(jì)算還得到電子自旋密度、Mulliken電荷和Mayer鍵級(jí)。

2 結(jié)果與討論

2.1 幾何結(jié)構(gòu)

采用密度泛函(DFT)理論方法優(yōu)化了g-C3N4配體結(jié)構(gòu)單元,標(biāo)記為CN。將三價(jià)和二價(jià)金屬離子放置于配體空穴中,計(jì)算得到配合物[M(CN)]z+(M=Sc、Y、La、Ac和U;z=3和2)穩(wěn)定結(jié)構(gòu)。配體和代表性配合物示于圖1,而優(yōu)化的結(jié)構(gòu)參數(shù)(鍵長(zhǎng)、鍵角)、鍵級(jí)等列于表1。盡管有共軛存在,優(yōu)化顯示g-C3N4的各原子并不完全處于一個(gè)平面內(nèi),其三個(gè)3-s-三嗪結(jié)構(gòu)單元相互交錯(cuò)。然而引入金屬后,配體的平面性增強(qiáng),其六個(gè)N供體原子近似在一個(gè)平面內(nèi),標(biāo)記為N6。測(cè)量的M到N6平面鍵長(zhǎng)(DM-N6)顯示,[Sc(CN)]z+(z=3和2)的中心金屬幾乎位于N6平面內(nèi),對(duì)應(yīng)的鍵長(zhǎng)分別為0.003 ?(1 ?=0.1 nm)和0.004 ?(表1);其它金屬均不在N6平面內(nèi),其中金屬Ac偏離N6平面最遠(yuǎn),DM-N6分別是0.959 ?和0.907 ?。相比較而言,金屬相同時(shí),[M(CN)]3+的中心金屬比[M(CN)]2+的偏離N6平面更遠(yuǎn)些。另外,優(yōu)化的N-M-N角(α)也支持這一結(jié)論:即,當(dāng)α→60°時(shí),金屬位于平面內(nèi);偏離越大,金屬偏離平面越遠(yuǎn)。因此,計(jì)算的[Ac(CN)]z+(z=3和2)的α角分別為55.5°和56.0°,偏離平面最遠(yuǎn)。當(dāng)金屬相同時(shí),從z=3→2變化,即發(fā)生單電子還原時(shí),Sc配合物的DM-N6基本保持不變(變化僅為0.001 ?),而其它四個(gè)金屬時(shí)DM-N6則變短0.010～0.084 ?,這表明還原作用使得金屬中心更趨近于N6平面。

表1 計(jì)算得到的[M(CN)]z+的鍵長(zhǎng)(D)、鍵角、鍵級(jí)(括號(hào)內(nèi)數(shù)值)和FSRM-NTable 1 Optimized bond lengths, angles and bond orders(in parenthesis) of [M(CN)]z+, along with calculated formal shortness ratio(FSRM-N)

圖1 研究的含六個(gè)N空穴的g-C3N4配體單元(CN)(a)和配合物[M(CN)]z+(M=Sc、Y、La、Ac和U;z=2和3) (b)的結(jié)構(gòu)圖Fig.1 Investigated structures of molecular unit of g-C3N4 containing six nitrogen donors(marked as CN)(a), and [M(CN)]z+ (M=Sc, Y, La, Ac and U; z=2 and 3)(b)

2.2 配體與金屬相互作用

為了更好對(duì)比,將M-N鍵長(zhǎng)相關(guān)的系列參數(shù)繪制于圖2中。對(duì)于配合物[M(CN)]z+,從圖2(a)可以看出,當(dāng)帶有三個(gè)正電荷(z=3)時(shí),優(yōu)化的M-N鍵長(zhǎng)(平均值)先逐漸增長(zhǎng)后縮短:即中心金屬由Sc變到Ac,M-N鍵長(zhǎng)是逐漸增長(zhǎng)的,而從Ac到U的鍵長(zhǎng)縮短;在Ac原子處出現(xiàn)拐點(diǎn)。同樣的,當(dāng)z=2時(shí),M-N的鍵長(zhǎng)變化趨勢(shì)與z=3時(shí)一致,均為先增長(zhǎng)后縮短,在Ac處出現(xiàn)拐點(diǎn)。[M(CN)]3+從Sc→Y→La→Ac鍵長(zhǎng)增長(zhǎng)與同族元素隨著周期數(shù)增加原子半徑增大相符合,而Ac→U的M-N鍵長(zhǎng)縮短可歸結(jié)為錒系收縮概念。

圖2 計(jì)算的配合物[M(CN)]z+的M-N平均鍵長(zhǎng)(a)、鍵級(jí)(b)、FSRM-N(c)和還原前后配合物M-N鍵長(zhǎng)差(ΔDM-N)(d)Fig.2 Calculated distances of M-N(a), bond order(b) and formal shortness ratio FSRM-N(c) for [M(CN)]z+, along with difference of M-N distances between complexes before and after one-electron reduction(d)

由圖2(a)還發(fā)現(xiàn),二價(jià)配合物陽(yáng)離子(z=2)的M-N(M=Sc和Y)的鍵長(zhǎng)要稍長(zhǎng)于三價(jià)配合物(z=3)的。而有趣的是,隨著中心金屬原子序數(shù)增加,二價(jià)陽(yáng)離子配合物的M-N鍵長(zhǎng)逐漸與三價(jià)的相接近,最后二價(jià)的M-N鍵長(zhǎng)變得比三價(jià)的短。出現(xiàn)拐點(diǎn)的位置同樣是在Ac處。為了進(jìn)一步考查金屬和配體作用,繪制M-N鍵級(jí)示于圖2(b)。由圖2(b)可以看出,當(dāng)z=3時(shí),隨著中心金屬原子序數(shù)的增大,鍵級(jí)的趨勢(shì)是先減小后增大的,總體變化趨勢(shì)是上升的;當(dāng)z=2時(shí),鍵級(jí)的變化趨勢(shì)為單調(diào)遞增的。比較z=3和2時(shí)的鍵級(jí)發(fā)現(xiàn),在M為Sc、Y和La的時(shí)候,三價(jià)配合物的M-N鍵級(jí)高于二價(jià)的;當(dāng)中心金屬為Ac和U時(shí),則是三價(jià)配合物的M-N鍵級(jí)低于二價(jià)的。與圖2(a)中M-N鍵長(zhǎng)出現(xiàn)拐點(diǎn)的位置相同。

一般來(lái)說(shuō),鍵級(jí)大時(shí),表明形成的化學(xué)鍵作用越強(qiáng),而且相應(yīng)的鍵長(zhǎng)則是越短的。對(duì)比圖2(a)和圖2(b)可以看出,在M=Sc、Y、La時(shí),三價(jià)配合物的M-N鍵級(jí)大于二價(jià)的,相應(yīng)三價(jià)配合物的鍵長(zhǎng)應(yīng)該是短于二價(jià)的;而M=Ac、U時(shí)則規(guī)律相反。

在應(yīng)用M-N鍵長(zhǎng)比較化學(xué)鍵強(qiáng)度時(shí),會(huì)因?yàn)樵影霃讲煌斐杀容^困難。因此,引入了FSRM-N扣除原子半徑不同的影響。從圖2(c)可以看出,三價(jià)配合物的FSRM-N高于二價(jià)的。這與鍵長(zhǎng)相反,和鍵級(jí)的情形相一致。同時(shí)發(fā)現(xiàn),在La之后,FSRM-N趨于升高,與圖2(a)和圖2(b)中出現(xiàn)拐點(diǎn)的位置略有不同。在中心金屬原子相同時(shí),考慮[M(CN)]3+→[M(CN)]2+單電子還原對(duì)鍵長(zhǎng)的影響,對(duì)二價(jià)和三價(jià)陽(yáng)離子配合物的M-N鍵長(zhǎng)進(jìn)行了做差處理,得到還原前后配合物M-N鍵長(zhǎng)差(ΔDM-N)示于圖2(d)。由圖2(d)可以看出,在這五種配合物的還原過(guò)程中,M-N鍵長(zhǎng)的變化值較小。當(dāng)M=Sc、Y、La時(shí),還原作用使得M-N鍵長(zhǎng)增長(zhǎng),即還原作用使得金屬和配體作用減小;而當(dāng)M=Ac、U時(shí),還原作用使得M-N鍵長(zhǎng)變短。M(Ⅲ)的離子半徑小于M(Ⅱ),所以結(jié)合ΔDM-N結(jié)果,可以粗略判斷,Sc、Y和La配合物的單電子還原有可能發(fā)生在金屬中心,而Ac和U配合物的單電子還原更可能發(fā)生在配體。

表2 配合物[M(CN)]z+在M-N鍵臨界點(diǎn)的QTAIM參數(shù)和鍵相互作用能EintTable 2 Calculated QTAIM data at M-N bond critical points(BCPs) of [M(CN)]z+, together with bond interaction energy Eint

2.3 單電子還原與電子結(jié)構(gòu)

計(jì)算的配合物[M(CN)]z+(M=Sc、Y、La、Ac和U;z=3和2)各片段的電子自旋密度(S)與電荷(Q)列于表3。通過(guò)這些計(jì)算結(jié)果,一方面可以判斷中心金屬在單電子還原前后的氧化態(tài),另一方面,幫助分析各部分之間的相互作用,并與鍵長(zhǎng)、鍵級(jí)和QTAIM參數(shù)相印證。當(dāng)中心金屬M(fèi)為二價(jià)時(shí),Sc、Y、La和Ac離子的電子自旋密度理論期望值為1.0,而U(Ⅱ)的為4.0。當(dāng)中心金屬M(fèi)為三價(jià)時(shí),U(Ⅲ)離子的電子自旋密度理論期望值為3.0,其它M(Ⅲ)的均為0.0。根據(jù)表3數(shù)據(jù),可以看出[M(CN)]2+中心金屬如果為二價(jià),計(jì)算的中心金屬的電子自旋密度SM均小于期望值,而配合物中配體的SCN大于0。這說(shuō)明有部分電子從金屬轉(zhuǎn)移到了配體上,且在Sc族中隨著原子序數(shù)的增加,電子自旋逐漸偏離期望值。當(dāng)中心金屬為Ac原子時(shí),發(fā)生的電子流向最為明顯。這說(shuō)明Ac應(yīng)該為Ⅲ價(jià)。對(duì)于[M(CN)]3+系列配合物,中心金屬SM的計(jì)算值與M(Ⅲ)期望值相符,幾乎沒(méi)有發(fā)生電子轉(zhuǎn)移。

表3 [M(CN)]z+的金屬和配體片段的電子自旋密度和電荷Table 3 Electron-spin density and charge of each fragment for [M(CN)]z+

進(jìn)一步考察[M(CN)]3+單電子還原生成[M(CN)]2+過(guò)程中各組分的自旋密度變化ΔS,結(jié)果示于圖3。如圖3所示,當(dāng)M為Sc時(shí),發(fā)生單電子還原時(shí),金屬捕獲的電子相對(duì)較多,而配體部分捕獲的電子相對(duì)較少;這說(shuō)明在還原時(shí),金屬與配體均參與了反應(yīng),且金屬Sc起主導(dǎo)作用,可以認(rèn)為三價(jià)Sc還原為二價(jià)。當(dāng)中心金屬為Y時(shí),ΔSM為0.485,而ΔSCN為0.515,金屬和配體同時(shí)被還原。在中心金屬為L(zhǎng)a和U時(shí),單電子還原主要發(fā)生在配體上,對(duì)應(yīng)的ΔSCN值分別為0.725和0.777。而當(dāng)中心金屬為Ac時(shí),是發(fā)生配體還原,金屬部分僅捕獲0.035的電子。

圖3 配合物[M(CN)]3+單電子還原為[M(CN)]2+過(guò)程中金屬和配體片段的自旋密度變化ΔSFig.3 ΔS value of each fragment while reducing [M(CN)]3+ to [M(CN)]2+ by one electron

為了解析配合物電子結(jié)構(gòu),繪制了[U(CN)]z+(z=3和2)的α自旋軌道態(tài)密度(DOS)和電子云圖,結(jié)果示于圖4和圖5,并將部分前線軌道成分列于表4和表5。對(duì)于[U(CN)]3+,由圖4和表4可以看出,其HOMO、H-1和H-2為金屬U的軌道性質(zhì),5f分別約貢獻(xiàn)86.7%、86.5%和93.7%,U(7s)和U(6d)貢獻(xiàn)很少;這些軌道中沒(méi)有CN的貢獻(xiàn)。[U(CN)]3+中有三個(gè)自旋平行5f電子,U為Ⅲ價(jià),與其計(jì)算的SM≈2.99相符合。該配合物的H-4為配體貢獻(xiàn)軌道,為π(CN)性質(zhì)。其低能的非占據(jù)軌道LUMO、L+1和L+2均由U(5f)組成,貢獻(xiàn)超過(guò)92.6%。當(dāng)發(fā)生單電子還原生成[U(CN)]2+時(shí),由圖5和表5可以看到,出現(xiàn)了四個(gè)U(5f)為主的軌道,即HOMO—H-3,其中HOMO和H-1的5f成分超過(guò)74.6%,而在H-2和H-3中,金屬貢獻(xiàn)約為50%～60%。值得指出的是,在H-1、H-2和H-3中配體CN的貢獻(xiàn)分別約為2.3%、28.9%和17.9%;如果把這三個(gè)軌道的金屬和配體成分進(jìn)行歸一化處理,發(fā)現(xiàn)CN成分總和約為60%。這說(shuō)明[U(CN)]3+→[U(CN)]2+的金屬和配體均發(fā)生了還原,且主要發(fā)生在配體上。這與上面的電子自旋密度分析結(jié)論相一致。

表4 計(jì)算的配合物[U(CN)]3+的軌道成分分布(P)Table 4 Contribution of orbitals for [U(CN)]3+

表5 計(jì)算的配合物[U(CN)]2+的軌道成分分布Table 5 Contribution of orbitals for [U(CN)]2+

圖4 配合物[U(CN)]3+的α自旋軌道態(tài)密度和電子云圖Fig.4 Density of states of α-spin orbitals of [U(CN)]3+ and electron-density diagrams

圖5 配合物[U(CN)]2+的α自旋軌道態(tài)密度和電子云圖Fig.5 Density of states of α-spin orbitals of [U(CN)]2+ and electron-density diagrams

3 結(jié) 論

采用相對(duì)論密度泛函理論探索了氮化碳對(duì)系列金屬離子負(fù)載和穩(wěn)定化作用,系統(tǒng)研究了生成的配合物[M(CN)]z+(M=Sc、Y、La、Ac和U;z=2和3)的結(jié)構(gòu)、成鍵和還原性質(zhì),得到如下結(jié)論:

(1) M-N鍵長(zhǎng)隨著原子序數(shù)增加的變化趨勢(shì)為先增加后減小,并且發(fā)現(xiàn)對(duì)應(yīng)的鍵級(jí)也增大;值得注意的是,這一研究結(jié)果不同于一般的常規(guī)認(rèn)識(shí),即鍵長(zhǎng)和鍵級(jí)具有反比關(guān)聯(lián)性;

(2) 不同中心金屬配合物的被還原片段有很大差異:在M為Sc時(shí),還原作用主要發(fā)生在金屬部分,配體僅有極少的參與;而Y配合物則是金屬和配體均參與,且占比相近;La和U配合物是金屬修飾的配體還原機(jī)制;不同的是,Ac配合物的還原反應(yīng)完全發(fā)生在配體上,金屬未參與;

(3) 結(jié)合電子結(jié)構(gòu)計(jì)算,指認(rèn)[U(CN)]3+的中心U為三價(jià)氧化態(tài),而[U(CN)]2+的U則介于三價(jià)和二價(jià)之間;計(jì)算發(fā)現(xiàn)配合物為M-N6配位模式,N→M配位鍵強(qiáng)度值約在-0.47～-0.79 eV范圍。

總之,g-C3N4可以通過(guò)其N(xiāo)6配位空穴負(fù)載金屬原子,從而形成單原子催化劑。計(jì)算得到的幾何/電子結(jié)構(gòu)信息、化學(xué)鍵和氧化還原性質(zhì)等將為氮化碳基金屬單原子催化劑的合成與應(yīng)用提供準(zhǔn)確的理論依據(jù)。