【純計算論文】J. Alloys Compd:富勒烯中封裝的雙金屬合金助力氧還原/析氧反應 2023年12月16日 下午10:23 ? 計算 ? 閱讀 17 研究背景 隨著工業的發展,化石燃料的快速消耗和相關的環境問題激發了人們對開發綠色和可持續能源系統和設備的研究興趣,如水分解、鋅空氣電池和燃料電池等。雖然這些能量體系和器件的研究取得了很大的進展,但氧還原反應(ORR)或析氧反應(OER)的緩慢反應動力學限制了它們的廣泛應用。鉑(Pt)和氧化釕(RuO2)通常分別被認為是ORR和OER的高效催化劑,但它們的高成本和稀缺性限制了這些材料的大規模應用。最近,一些研究表明,封裝在碳材料中的金屬團簇對ORR或OER具有良好的催化性能。對于這種類型的催化劑,金屬核與碳殼之間的相互作用使電子很容易從金屬核轉移到碳殼上,從而影響催化性能。 西南石油大學陳鑫等人基于密度泛函理論(DFT)方法,對富勒烯M1xM24?x@Cn中封裝的雙金屬合金的ORR和OER催化活性(M1xM24?x代表FexCo4?x,FexNi4?x,CoxNi4?x;x = 1、2、3;n = 40、50、60)進行系統的研究,揭示催化劑的構效關系。并通過計算反應中間體的結合能和各基本步驟的吉布斯自由能來確定其活性。 模型與計算方法 采用Materials studio軟件包的DMol3模塊對其電催化析氫性能及電子性質進行理論研究,計算基于廣義梯度近似(GGA)下的BLYP泛函計算電子交換相關性,通過DNP基組對內核電子進行處理,應用有效核近似法(ECP)來改善計算的準確性。將能量、最大力和最大位移的收斂公差分別設為2× 10-5 Ha、0.004Ha·?-1和0.005 ?。 結果與討論 圖1為M1xM24?x@Cn催化劑的結構圖,分別在C40,C50,C60上進行雙元合金的封裝。不同催化劑的封裝能如圖2所示,封裝能越小,結構越穩定。可以清楚地看到,除了Fe1Ni3@C60外,M1xM24?x@Cn的所有Eenc值都為正,這表明將雙金屬合金封裝在富勒烯中還需要額外的外界能量。此外也能發現碳殼尺寸越大,封裝能越低。 圖1 催化劑結構圖 圖2 催化劑的封裝能 對于封裝金屬團簇的碳材料,碳殼將是吸附反應中間體的位置,帶正電荷的碳原子更有利于成為活性位點。圖三為在C活性位點上進行ORR/OER的反應中間體的結合能。對于M1xM24?x@C40,*OOH的結合強度強于Pt(111)表面(3.66 eV),有利于ORR過程中*OOH的形成。但與Pt(111)表面的ΔE*OH值相比,所有M1xM24?x@C40的ΔE*OH值要小得多,說明所有M1xM24?x@C40與*OH之間的相互作用過于強烈。因此,所有M1xM24?x@C40的催化活性都受到*OH極強的結合作用的影響,對于OER,M1xM24?x@C40表面的ΔE*都比RuO2(110)表面的吸附能小得多,說明M1xM24?x@C40的OER催化活性可能較差。 M1xM24?x@C50上的ΔE*OOH值與Pt(111)表面上的值很接近,*OH在M1xM24?x@C50上的結合強度相對較強,因為ΔE*OH值較小。在所有M1xM24?x@C50中,Co3Ni1@C50上的ΔE*OOH和ΔE*OH值與Pt(111)表面的相對接近,表明此催化劑比其他M1xM24?x@C50具有更高的ORR活性。大多數M1xM24?x@C50上的ΔE*物種值大于M1xM24?x@C40上的相應值,因此將雙金屬合金封裝在更大尺寸的富勒烯中,可以通過調節中間體的結合,使催化體系具有更好的ORR和OER活性。 對于M1xM24?x@C60,Fe2Ni2@C60和Fe3Ni1@C60的ΔE*OOH值與Pt(111)相對接近,有利于后續的ORR過程??梢灶A測,CoxNi4?x@C60可能具有很好的OER催化活性。一般來說,隨著Cn尺寸的減小,反應中間體在M1xM24?x@Cn上的結合強度大致增強。這可能是由于較小的富勒烯導致合金芯和碳殼之間更大的電荷轉移所致。 圖3 催化劑在ORR/OER上的反應結合能 以往的研究表明,反應中間體的吸附自由能之間存在線性關系。因此,建立了M1xM24?x@Cn上ΔG*OOH和ΔG*O與ΔG*OH的比例關系,如圖4所示。ΔG*OOH與ΔG*OH之間的相關性可描述為ΔG*OOH = 0.82ΔG*OH + 3.32,決定系數(R2)為0.87,說明ΔG*OOH與ΔG*OH之間存在較強的線性關系。斜率和截距與其他材料非常相似。同時,ΔG*O和ΔG*OH顯示出ΔG*O = 1.32ΔG*OH + 0.96(R2 = 0.72)的標度關系。 圖4 吸附自由能線性關系 過電位與?G*OH之間的火山關系如圖5所示。速率決定步(RDS)定義為ΔG的最大值步長。對于ORR,與*OH結合較強的M1xM24-x@Cn的RDS (ΔG*OH <0.88 eV)為*OH→H2O的反應步驟。隨著*OH在M1xM24-x@Cn上的結合強度減弱,數據點落在紫色線上,即RDS從*OH還原步驟轉換為*OOH生成步驟。令人鼓舞的是,Co1Ni3@C50的ηORR值最小,為0.35 V,位于火山的頂峰(倒置)。此外,Co3Ni1@C50、Fe3Ni1@C60和Fe2Co2@C50的ηORR值分別為0.36、0.39和0.43 V,表現出優異的ORR活性。 同樣,對于OER,當ΔG*OH值為1.57 eV時,ηOER值為0.34 V,出現火山頂。結果表明,除Fe1Ni3@C60 (H2O→*OH)外,大部分M1xM24?x@Cn的RDS都是*O→*OOH的階。此外,在M1xM24?x@Cn中,Co2Ni2@C60距離火山頂最近,ηOER值為0.36 V。Co1Ni3@C60 的OER催化活性也與RuO2(110)表面(ηOER = 0.37 V)相當。因此,M1xM24?x@C50和M1xM24?x@C60大部分位于ORR和OER火山的近頂部,因此C50和C60比C40更適合包封雙金屬合金,更適合于高效ORR和OER催化劑。 圖5 自由能與過電位的關聯 為了深入了解催化活性,計算了M1xM24?x@Cn上每個反應步驟的自由能。Co1Ni3 @C50和Co2Ni2@C60的自由能圖如圖6所示,其余見支撐數據。對于ORR,可以清楚地發現除Co2Ni2@C40外,M1xM24 – x@Cn上每一步的自由能曲線都是下坡的。結果表明,這些催化劑上的ORR是自發發生的。Co2Ni2@C40上的*OH還原步驟能量向上,這是由于*OH的結合過于強烈造成的。Co1Ni3@C50對ORR的催化活性最高,優于先前報道的性能。對于OER, Co2Ni2@C60的ΔGOER值最小(1.59 eV),與RuO2(110)表面的ΔGOER值(1.60 eV)接近,催化活性顯著。 圖6 催化劑反應自由能 結論與展望 本文采用DFT方法研究了富勒烯包封雙金屬合金的ORR和OER催化活性。發現碳殼尺寸越大,封裝能越低,并且大部分催化劑存在優異的ORR/OER活性。本研究為雙金屬合金在富勒烯中封裝的新型催化劑設計提供了研究思路,并有助于指導高性能電催化劑的合理設計。 文獻信息 Chen, X., Huang, S., & Zhang, H. (2022). Bimetallic alloys encapsulated in fullerenes as efficient oxygen reduction or oxygen evolution reaction catalysts: A density functional theory study. Journal of Alloys and Compounds, 894, 162508. 原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2023/12/16/1b43c3bde0/ 贊 (0) 0 生成海報 相關推薦 【純計算】AEM: 利用靜電勢溶劑描述符實現鋰電池電解液的合理設計 2024年3月13日 郭再萍等EES:Cu單原子抑制鋅-碘電池的穿梭效應 2023年11月4日 【DFT+實驗】Nano Lett.:外延MnBi2Te4薄膜的柵極可調諧本征反常霍爾效應 2023年12月29日 【DFT+實驗】港中文Jimmy Yu課題組聯合港理工黃渤龍與勞子恒 JACS:電催化二氧化碳還原中的高分散鉍納米顆粒的尺寸效應 2024年3月5日 JPCL快訊 | 基于集合變量的增強采樣算法:從人類學習到機器學習 2024年3月26日 【DFT+實驗】中科大Angew.:構建CdS/Ti3C2 MXene肖特基結,用于高選擇性和活性光催化脫氫還原氨基化 2023年11月7日