【MS論文解讀】相鄰雙NiNx位點對提高燃料電池質子交換膜氧還原反應活性的協同作用 2023年11月5日 下午1:32 ? 計算 ? 閱讀 23 研究背景 隨著化石燃料燃燒引起的全球環境問題的發展,清潔可再生能源受到越來越多的關注。燃料電池被認為是未來汽車的主要動力源。在燃料電池中,陽極氧還原反應(ORR)的動力學緩慢的問題亟待解決。 重慶理工大學張志剛和宋潔等人通過第一性原理密度泛函理論(DFT)計算研究了石墨烯中嵌入N原子的雙Ni飽和配位中的反應機理和相應的熱力學能量演化。 計算方法 本文使用Materials Studio中的Dmol3模塊,采用廣義梯度近似(GGA)下的PEB泛函進行計算交換和相關能量。所有模擬均以摻雜Ni-N配位結構的5×5單層石墨烯為模型。為了避免周期單元之間可能的相互作用,將石墨烯片z方向的真空厚度設置為20 ?。原子軌道用雙數值加極化函數(DNP)表示。采用k點為4×4×1的Monkhorst-Pack網格方案。燃料電池陰極催化劑的水環境模擬采用Cosmo模型,介電常數設為78.54。考慮了范德華相互作用,采用OBS方法進行DFT-D校正。能量、力和最大位移的收斂標準分別為1.0×10-5 Ha、0.002 Ha/?和0.005?。 結果與討論 相鄰雙NiNx位點的優化幾何圖形如圖1所示。在溶劑環境下,壓力為0.035 bar,溫度為298.15 K時,整個反應的自由能為-4.92 eV,與實驗結果相吻合。在298.15 K的溶劑環境下,H2O2的還原自由能變化為-3.46 eV。 圖1. 相鄰雙NiNx位點優化示意圖 圖2.??(a) Mulliken電荷(b)靜電勢面(c) HOMO和(d) LUMO Mulliken電荷分布如圖2所示,Ni和N之間明顯的電荷差異表明與其他位置相比,Ni的頂部位置對含氧中間體的吸附活性最高。 所有結構中Ni的分波態密度(PODS)計算結果如圖3所示。所有的催化劑都是非磁性的,這些催化劑在費米能級以下的Ni的3d軌道PDOS有兩個明顯的峰,其中較高的峰更接近費米能級,說明其有較好的電子遷移性質。 圖3. Ni的3d軌道的分波態密度 作者研究了這含氧中間體O2*、OOH*、O*、OH*、OH*+OH*、O*+OH*、H2O*、H2O2*在催化劑的活性位點上的吸附,結果如圖4和表2所示。 表2. 吸附構型中含O中間體的鍵長( ?) 對于H2O和H2O2,Ni-O的鍵長接近3.0A,吸附能小于-0.3 eV。因此,兩者都是物理吸附。這意味著H2O分子的產物很容易從活性位點解吸,有利于另一個還原過程。這也意味著H2O2的吸附和隨后的還原是難以處理的。 幸運的是,Ni2N6-ST2/G和Ni2N6-ST3/G中的活性位點一旦與活性位點結合,就可以催化H2O2的還原過程。 圖4. 含氧中間體在催化劑上最穩定的吸附構型 含O中間體的吸附能排序如圖5所示。含O中間體在這些活性位點上的吸附能幾乎按照O*、OH*、OOH*、O2*、H2O*、H2O2*的順序遞減。此外,含氧中間體在Ni2N6-ST2/G和Ni2N6-ST3/G上的吸附能明顯大于其他催化劑上對應的中間體。 圖5. 含氧中間體在催化劑上的吸附能 這些催化劑上可能的基本反應步驟和相應的自由能見表3。在催化劑上吸附后,發生氧還原反應。 表3. 反應過程中的自由能(eV)變化 催化劑的氧還原途徑有不同的反應途徑,但它們的RDS都是第一步還原(O2*/OOH*)。Ni2N6-ST2/G和Ni2N6-ST3/G的氧還原路徑為O2*→OOH*→O*→OH*→H2O, 通過O*機制促進了四電子途徑。 Ni2N6-ST1/G和Ni2N7-ST4/G的氧還原路徑為O2*→OOH*→ 2OH*→OH*→H2O,通過2OH*機制促進了四電子途徑。NiN4/G的氧還原路徑為O2*→OOH*→H2O2,促進了雙電子途徑。圖6a顯示了計算得到的這些活性位點上氧還原途徑的自由能演化圖。這些催化劑上的所有反應步驟都是下坡的,表明在U = 0 V下反應很容易。 雖然NiN4/G對O2的吸附能為正,自由能為+0.15 eV,但其對ORR在NiN4/G上0.69 V的起始電勢仍然過大。在燃料電池運行狀態下,H2O2的生成會導致電池性能的下降。因此,減少孤立Ni活性位點,增加雙Ni活性位點有利于提高燃料電池的壽命。特別是Ni2N6-ST2/G和Ni2N6-ST3/G都能吸附H2O2并將其還原成H2O。 圖6. 酸性條件下施加電位(a) U = 0 V和(b) U = 0.68 V時氧還原反應過程的自由能變化 作者分別考察了外加電位和pH對對Ni2N6-ST3/G吸附含氧中間體穩定性的影響。如圖7a所示,這些含氧中間體的自由能隨外加電位的增加而線性減小,O2*的穩定性隨著外加電位的增加而迅速提高。pH值在0 ~ 14之間變化對含o中間體穩定性的影響,如圖7b所示。 這些含O中間體的臨界電位也隨pH值的增加而線性降低,而OH*和O*的電位差是恒定的,不隨pH值的變化而變化。 圖7. Ni2N6-ST3/G上ORR工藝含O中間體的穩定性 結論與展望 作者采用DFT計算研究了雙鎳位點催化劑吸附的含氧中間體上氧還原反應的ORR機理及其ORR性能。計算得到的形成能表明,與NiN4/G的活性位點相比,雙NiNx催化劑更有利于穩定性。這些催化劑表面Ni和N之間明顯的電荷差異也表明,與其他位置相比,Ni的頂部位置對含氧中間體的吸附活性最高。 Ni2N6-ST2/G和Ni2N6-ST3/G的3d軌道相鄰PDOS峰和對應的d帶中心明顯更接近費米能級,有利于O2和OH的吸附過程。Ni2N6-ST1/G、Ni2N6-ST2/G和Ni2N6-ST3/G上ORR的導通電位分別為0.36 V、0.65 V和0.68 V。此外,還發現氧還原途徑中最穩定的中間體取決于它們在給定pH值和外加電位下的最低自由能。這些含氧中間體的自由能隨外加電位的增加而線性減小。這些含氧中間體的臨界電位也隨pH的增加而線性降低。 文獻信息 Qi, W., Huang, J., Zhang, Z., Huang, Z., Song, J., Fu, L., … & Ran, J. (2023). The synergistic effect of the adjacent dual NiNx sites to boost the oxygen reduction reaction activity for proton exchange membrane fuel cell.?International Journal of Hydrogen Energy,?48(9), 3632-3642. https://doi.org/10.1016/j.ijhydene.2022.10.047 原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2023/11/05/961c8f85fc/ 贊 (0) 0 生成海報 相關推薦 【全文解讀+計算重現】Advanced Materials經典工作:DFT計算+實驗表征+電池! 2023年12月22日 復旦大學、中國科學院團隊綜述,化學機器學習:基礎知識和應用 2023年9月17日 【實驗+計算】Angew.:Bi-ZMOF高效電催化CO2還原為甲酸! 2023年10月2日 【計算+實驗】頂刊案例解讀:JACS、EES、Angew、Nature子刊、ACS Catalysis、ACS Nano等成果! 2023年12月27日 電池管理系統:物理模型與機器學習集成 2024年4月1日 Nano Res.│山西師大賈建峰課題組和中科大葛君杰課題組綜述:酸性介質中兩電子氧還原制備過氧化氫——單原子催化劑的優越性 2024年3月16日