【機器學習+DFT計算】JEC:基于DFT和機器學習探索穩定高效的OER電催化劑 2024年3月10日 下午12:21 ? 計算 ? 閱讀 173 研究背景 機器學習(ML)、數據科學和分析的革命性發展,以及它在材料科學中的應用,是過去十年科學界的一個重要里程碑?由于析氧反應(OER)在一系列電化學能量轉換過程中的重要性,研究活性高、穩定性好和經濟效益優異的OER催化劑至關重要。韓國順天鄉大學Hwanyeol Park教授和漢陽大學Ho Jun Kim教授等人設計了一種高效的高通量篩選方法,通過密度泛函理論(DFT)計算和機器學習算法,在酸性條件下為高性能析氧反應(OER)催化劑尋找穩定的過渡金屬氧化物? 研究亮點 1.?作者設計了一種高效的高通量篩選方法,以尋找具有優異OER性能的潛在催化劑,并介紹了一種利用Materials Project數據庫和DFT計算的方法來評估特定反應條件下催化劑的穩定性? 2. 此外,本研究利用梯度增強回歸(GBR)模型和ML方法,深入了解吸附性質與物理和化學描述符之間的關系,發現TM-O的鍵長和第一電離能是與吸附行為相關的兩個值得注意的描述符? 3. 最終,確定了十種具有良好前景的催化劑(IrO2、Fe(SbO3)2、Co(SbO3)2、Ni(SbO3)2、FeSbO4、Fe(SbO3)4、MoWO6、TiSnO4、CoSbO4和Ti(WO4)2),從而拓寬了OER電催化劑的選擇范圍? 計算方法 作者使用VASP計算軟件包進行密度泛函理論(DFT)計算,采用廣義梯度近似(GGA)下的PBE泛函計算電子交換相關性,同時通過投影增廣波(PAW)方法描述價電子與離子之間的相互作用? 考慮到電子之間的庫倫排斥作用,此文采用了Dudarev等人提出的Hubbard+U校正方法,并且在z方向上由8 ?的真空層以避免周期性系統之間的相互作用。同時在能量和力的收斂準則上分別設置為1×10-4eV和0.06eV/?。最后,DFT的自動化計算通過Python包Pymatgen中的各種模塊實現,機器學習方法是通過GBR模型實現,由Scikit-learn和PyTorch包提供。 圖文導讀 作者設計的高通量篩選過程的示意圖如圖1所示,篩選由Materials Project數據庫啟動,考察了含過渡金屬氧化物(TM-oxides)的材料,篩選出了2366個適合分析電化學穩定性的化合物,并對其性能進行了分析?最初的篩選步驟包括計算它們的相穩定性Ehull < 0.1 eV/原子和水穩定性(Pourbaix圖上穩定域的吉布斯自由能)ΔGpbx < 0.1 eV/原子的標準?其次,TM-oxides應該表現出高的熱力學穩定性(吉布斯形成能ΔGf ≤0.1 meV/原子)?第三步,進行電化學分析,對所有968種金屬組合的ΔGpbx值的分析表明,其中只有159種在酸性OER環境中保持穩定? 圖1 高通量選擇標準 對于這些穩定的組合,作者構建了其低指數面((100)、(001)、(110)、(101)和(111))的模型,如圖2所示?然后,作者對所構建的3786個結構的低指數面模型進行了ΔGO*、ΔGOH*和ΔGOOH*的DFT計算,得到OER過電位(η = 0.3 V時的ηTD和Gmax)?最后,作者指出了10個最有希望的候選結構? 圖2 計算流程示意圖 目前與OER相關的研究主要基于中間態(熱力學),而忽略了過渡態(動力學),使得速率決定步驟變得與實際不符?作者采用微動力學方法,結合熱力學和動力學闡明了ηOER = 400 mV時的決速步(圖3a),結果發現,作為O2生成催化劑的施加過電位ηOER與活性描述符ηTD不相對應,使得ηOER與ηTD的比較不恰當。 圖3 OER動力學相關關系 在圖3a中從火山圖的左側,在動力學方面考察了OOH與OH的尺度關系?當OOH的形成為決速步時,其ΔG3的降低導致生成OOH中間體的自由活化能降低,從而提高電催化效率?相反,當OOH的分解成為決速步時,會導致催化活性降低? 隨后作者使用描述符Gmax(η)來計算OER活性,其采用自由能跨度模型來估計自由能形貌中的RDS?圖4顯示了Gmax(η = 0.3 V)與OH分解時的ΔG2的關系? 圖4 Gmax(0.3 V) 與ΔG2的標度關系 表1總結了最有前景的催化劑材料IrO2、Fe(SbO3)2 (P42/mnm)、Co(SbO3)2 (Imma)、Ni(SbO3)2 (P42/mnm)、FeSbO4 (I41md)、Fe(SbO3)4 (P1)、MoWO6 (P21/c)、TiSnO4 (Cmma)、CoSbO4 (Imma)和Ti(WO4)2 (P2/c)的活性數據?在這些有希望的OER候選催化劑中,CoSbO4的最低Gmax(η = 0.3 V)估計為0.15 eV,帶隙約為0.14 eV?衡量材料穩定性的Ehull值為0.09 eV/原子,表明CoSbO4易于合成,帶隙約為0.14 eV,這意味著其導電性容易調節?另一個值得注意的候選材料—Co(SbO3)2,其(101)面在0.28 eV時具有較低的Gmax(η = 0.3 V),并且 Ehull值為0 eV/原子。IrO2和Ni(SbO3)2也表現出高活性和0 eV的帶隙,表明其金屬性質,使它們有很大的發展前途。 表1 過渡金屬氧化物OER的電催化性能 描述符φ和εd雖然對吸附能有一定的定量描述,但沒有充分建立明確的線性關系,吸附能與各種原子和結構性質之間的關系尚不明確?為了解決這個問題,作者采用ML模型來闡明描述符與吸附能之間的內在相關性,如圖5(a)所示。 圖5 機器學習方法的流程圖及部分特征比較 作者選取159種穩定的過渡金屬氧化物(3786種結構)在η = 0.3 V時的Gmax作為目標數據。選取與TM-oxides的鍵長、配位原子(dTM-O、dO-O)、d帶中心(εd)、電負性(Nm)、電子親和性(Am)、第一電離能(Im)、價電子數(Nv)和TM-d電子數(Ne)等8個與原子和結構性質相關的描述符作為輸入數據? 從數據生成開始,將364個輸入數據隨機分為292個訓練集和72個測試集。如圖S8所示,GBR模型在訓練過程中表現出優異的性能,R2為0.96,最小RMSE為0.06 eV。在測試中,模型保持了良好的性能,驗證了ML方法用于特征分析的可靠性。進一步對比了這些描述符的重要性,如圖6(c)所示,從而對吸附能力的起源有了更深入的了解。觀察發現TM-O的鍵長(dTM-O)和第一電離能(Im)與吸附能力的相關性較強,其特征重要性值分別為39.3%和33.5%,而其余6個描述符的特征重要性值相對較低? 在分析了催化活性、吸附能和電子結構之間的關系之后,仍然存在兩個未能解決的問題:(1)如何將雙官能團OER活性與過渡金屬氧化物的內在化學和物理性質聯系起來;(2)如何將此研究的方法擴展到不同反應的替代催化劑體系。通過在跨學科研究中的應用,該方法有可能解決這些挑戰。 文獻信息 Park, H., Kim, Y., Choi, S., & Kim, H. J. (2024). Data driven computational design of stable oxygen evolution catalysts by DFT and machine learning:Promising electrocatalysts. Journal of Energy Chemistry. https://doi.org/10.1016/j.jechem.2023.12.048 原創文章,作者:計算搬磚工程師,如若轉載,請注明來源華算科技,注明出處:http://www.zzhhcy.com/index.php/2024/03/10/96a9338b7f/ 贊 (0) 0 生成海報 相關推薦 加速材料屬性預測:結構感知圖神經網絡 2024年2月25日 北理工白瑩團隊EES綜述:自上而下方法論推動錳基正極在水系鋅離子電池中的性能優化——電子結構調控 2023年11月27日 【DFT+實驗】Nano Lett. | 湖南大學潘安練教授團隊: 單層WS?中氧吸附誘導的熒光增強及其擴展行為 2024年4月4日 崔屹/鮑哲南,最新PNAS!鋰硫電池再現“中庸之道”! 2024年2月6日 【DFT+實驗】廈大熊海峰/福大林森:羥基和Rh單原子之間的距離決定Rh?/CeO?催化劑在CO氧化中的活性 2024年1月2日 彭慧勝院士,Angew! 2025年2月19日