機器學習原子間勢揭示α-Fe一般晶界氫脆起源

1  研究概述

氫脆是制約高強結構材料在綠氫制備、輸運與利用全鏈條中應用的關鍵瓶頸，尤其沿一般晶界的氫致開裂因其缺乏晶格對稱性而更難預測與防范。傳統經驗勢（如EAM）在描述含氫晶界的變形斷裂行為時存在顯著精度不足，而密度泛函理論（DFT）雖精度高卻受限于計算規模，無法直接模擬含氫晶界的大尺度力學響應。機器學習原子間勢（MLIP）為彌合這一精度-尺度的鴻溝提供了新途徑，但現有Fe-H MLIP的訓練集多局限于對稱傾轉晶界，難以準確描述一般晶界的氫偏聚與力學行為。

針對上述問題，Kazuma Ito、Takashi Otaki、Tatsuya Yokoi、Katsutoshi Hyodo等人在“Communications Materials"期刊上發表了題為“Machine Learning Interatomic Potential Reveals Hydrogen Embrittlement Origins at General Grain Boundaries in α-Iron"的研究論文。該研究通過并行學習（concurrent learning）策略，系統構建了覆蓋廣義層錯、表面、對稱傾轉晶界、位錯及空位團簇等所有與氫致脆化相關原子環境的Fe-H二元MLIP訓練集（圖1）?；?/span>Behler-Parrinello神經網絡勢（BNNP）框架，最終實現了對α-Fe中一般晶界氫偏聚、<111>/2全位錯發射、變形孿生及裂紋萌生擴展的DFT精度描述。大規模分子動力學模擬揭示了氫在不同變形機制晶界中的差異化影響，為高強合金抗氫脆設計提供了原子尺度的力學判據。

圖1  Fe-H二元體系MLIP訓練數據集構建流程：(a) DFT初始數據集；(b) 四個并行深度勢（DP）訓練；(c) 多條件分子動力學候選構型采集；(d) 候選構型DFT標記與迭代

2  MLIP構建方法與精度驗證

2.1  并行學習策略與訓練集

訓練集構建經歷169輪并行學習迭代，最終包含36,438個構型（約155萬局部原子環境）。訓練集系統涵蓋了：純α-Fe與應變Fe中的H溶解能、H-空位團簇、相鄰四面體間隙位的雙氫對、廣義層錯能、七種對稱傾轉晶界的氫偏聚能、以及納米多晶中一般晶界的結構與偏聚評估等。采用節點并行訓練與修正的n2p2超參數（截斷半徑從6.5 ? 優化至6.0 ?），在保證精度的前提下顯著提升計算效率。

2.2  基本物理性質與缺陷-氫相互作用精度

圖2展示了構建MLIP對Fe-H體系關鍵物理量DFT結果的復現能力。氫在四面體間隙位（TIS）的溶解能可準確匹配單軸應變與靜水壓條件的變化（圖2a,b）；螺位錯核心對氫的偏聚能預測與DFT高度一致（圖2c），盡管訓練集中未明確包含螺位錯-氫構型，驗證了MLIP優異的可遷移性；廣義層錯能曲線與氫-空位團簇相互作用能同樣達到DFT精度（圖2d,e）。全訓練集總能量與原子力的均方根誤差分別為0.13 meV/atom和35.8 meV/?。

圖2  構建MLIP的基本物性精度驗證：(a-b) TIS氫溶解能；(c) 螺位錯氫偏聚能；(d) (112)廣義層錯能（氫位于TIS）；(e) 空位-氫團簇相互作用能；(f) 折射率表面分離能；(g)能量的鍵長依賴性；(h) 四面體位置原子相互作用能；

2.3  對稱傾轉晶界氫偏聚精度

圖3展示了七種對稱傾轉晶界（Σ3-Σ11）的氫偏聚能預測結果。MLIP整體準確復現了DFT的偏聚能趨勢與相對次序，僅在高Σ值晶界的極低偏聚能區出現輕微系統偏差，但偏差幅度在可接受范圍內。這表明MLIP能夠精確描述氫在不同晶界能量/結構下的偏聚行為差異。

圖3  MLIP對七種對稱傾轉晶界氫偏聚能的預測精度：(a-g) 分別為Σ3(112)、Σ3(111)、Σ5(210)、Σ5(310)、Σ9(114)、Σ9(221)、Σ11(113)晶界；(h) DFT和原子間勢在七個晶界位點的偏析能對比。

2.4  一般晶界能量與氫偏聚精度

針對最難描述的一般晶界，研究從含8晶粒的納米多晶（28.8 nm）中提取計算單元，評估MLIP對一般晶界氫偏聚能的預測能力（圖4）。MLIP預測的平均晶界能與實驗外推值（圖4a）基本吻合，遠優于EAM勢。一般晶界氫偏聚能的RMSE僅為0.042–0.103 eV（圖4b,c），表明即使在訓練集中未出現相關構型，MLIP仍對一般晶界氫偏聚展現出優異的外推精度。這一結果為后續大規模MD模擬提供了必要的置信度。

圖4  MLIP對一般晶界氫偏聚的預測精度：(a) 平均晶界能與模型尺寸關系；(b) 納米多晶中一般晶界氫偏聚能的評估方法；(c-e) MLIP與DFT氫偏聚能對比

3  一般晶界氫偏聚行為

基于構建的MLIP，研究建立了六個隨機取向的雙晶一般晶界模型（約30萬原子），考察不同體氫濃度條件下的晶界氫偏聚狀態（圖5）。Voronoi分析表明，氫優先偏聚于晶界附近體積膨脹的八面體間隙位（圖5a,b），晶界氫占據數與體濃度呈近似正比關系（圖5c）。在中等與高體濃度條件下，一般晶界的氫面密度分別達到1.4和6.7 atoms/nm2，與Sato等人在高純Fe中陰極充氫實驗測得的2.3–2.9 atoms/nm2處于同一數量級，驗證了模擬結果的合理性。

圖5  一般晶界氫偏聚狀態：(a) 雙晶模型示例（GGB 1）；(b) 原子能量著色圖；(c) 體氫濃度與晶界氫濃度的關系；(d) GGB1晶界附近的氫濃度分布；(e-f) 氫在GGB1晶界中心的一個原子層（2.5 nm）區域的Voronoi體積和配位數的概率密度；(g) 體積氫濃度為0.0044 at.%時，GGB1晶界中心附近的原子結構

4  氫對一般晶界變形斷裂行為的影響

4.1  兩種變形模式

單軸拉伸MD模擬揭示了六個一般晶界模型的兩種本征變形機制（圖6a）：GGB 1–3由晶界處的<111>/2全位錯發射主導塑性；GGB 4–6由晶界附近連續發射<111>/6不全位錯形成變形孿晶主導塑性。前者的最大應力（11.2–13.0 GPa）通常略高于后者（9.3–9.8 GPa）。與Σ3(111)（15.4 GPa）和Σ3(112)（13.3 GPa）晶界相比，低Σ值晶界的最大應力顯著更高，體現了晶界有序度對應力響應的重要影響。

圖6  氫偏聚對一般晶界變形與斷裂行為的影響：(a) 清潔晶界應力-應變曲線；(b-c) GGB 1和GGB 3受氫影響的應力-應變曲線；(d-f) 氫對最大應力與斷裂應變的影響；(g) GGB 3晶界遷移與位錯發射；(h) GGB 4變形孿晶形成；(i) GGB 4-6斷裂應變隨氫濃度的變化；(j) GGB 4最大應力附近應變處的原子結構；(k) GGB4晶界變形孿晶的形核機制

4.2  氫對位錯主導型晶界的影響：抑制位錯發射、促進晶界斷裂

對于位錯主導塑性的GGB 1–3模型（圖6b–g），氫偏聚呈現一致趨勢：增大最大應力（抑制位錯發射的宏觀表現），并降低晶界斷裂應變。以GGB 1為例，中等氫濃度下最大應力顯著增大，斷裂應變顯著降低（圖6b）；位錯分析（DXA）表明，位錯發射位置不受氫影響，但晶界內位錯線伸入晶內即觸發發射，氫在位錯線附近偏聚提升了位錯發射的能壘（圖6e）。在GGB 2的中等氫濃度條件下，最大應力后并未立即開裂，而是繼續以位錯發射維持變形，約在14%應變達到第二應力峰后才發生斷裂。在GGB 3的低/中氫濃度條件下，晶界從約6%應變起發生顯著遷移，由遷移形成的突出部位發射位錯，但氫偏聚使應力單調下降并最終斷裂。高氫濃度時晶界遷移被抑制，彈性應變能無法釋放，導致在位錯發射之前即發生低應力脆斷。

4.3  氫對孿生主導型晶界的影響：弱化顯著降低

GGB 4–6模型由變形孿生負責塑性，其氫敏感性顯著弱于位錯主導型晶界（圖6h,i）。氫偏聚幾乎不改變最大應力與斷裂應變，即使晶界氫濃度與GGB 1–3相當（甚至更高）。變形孿生通過相鄰(112)晶面上連續發射<111>/6不全位錯形成（圖6j,k），該機制對氫的敏感度遠低于全位錯發射。這一發現與實驗報道的氫添加改善CoCrFeMnNi高熵合金力學性能的現象具有內在一致性——孿生主導變形的材料對氫脆的敏感性天然較低。

5  拉伸過程的外推精度驗證

為確認大變形條件下的預測可靠性，研究基于局部原子環境計算了拉伸過程中所有約30萬原子的外推等級（extrapolation grade）（圖7）。對位錯主導的GGB 1和孿生主導的GGB 4，無論氫濃度與應變條件，絕大部分原子始終處于內插（interpolation）區域（外推等級<1），僅有極少量原子接近精確外推邊界（圖7a–d）。這表明MLIP在拉伸過程中未遭遇超出訓練數據范圍的新原子環境，所有模擬結果均具有可靠的DFT精度。

圖7  拉伸過程MLIP外推精度驗證：(a-b) GGB 1和GGB 4全部原子的外推等級隨應變的變化；(c-d) 12%應變時GGB 1和GGB 4的原子結構與外推等級著色圖；(e-f) GGB 4單軸拉伸試驗中原子結構中所有原子平均和最大外推級的應變依賴性

6  討論與工程意義

6.1  氫抑制位錯發射的原子機制

模擬結果直接支持了實驗上關于氫偏聚于晶界可抑制位錯發射的推測。氫占據晶界附近的膨脹八面體間隙位后，提高了位錯線從晶界向晶內突出所需應力，這與HELP機制（氫降低位錯運動阻力促進滑移）看似矛盾，實則反映了氫在不同應力/濃度階段的差異化影響：低濃度促進位錯運動，高濃度偏聚反而抑制位錯發射，迫使晶界在低應變下發生脆性斷裂。這一發現為理解HEDE與HELP在高強合金晶界氫脆中的協同競爭提供了新視角。

6.2  對高強合金抗氫脆設計的啟示

本研究從原子尺度揭示了晶界本征變形機制是決定氫脆敏感性的核心參數：位錯主導變形的晶界為高敏感性，變形孿生主導的晶界為低敏感性。這為抗氫脆微結構設計指明了方向：（1）通過晶界工程增加特殊晶界（低Σ值）比例；（2）利用合金化調控層錯能，促進變形孿生而非位錯滑移；（3）引入納米析出相阻礙氫向一般晶界的擴散和偏聚。此外，團隊正基于本研究訓練數據構建精度-效率更優的新一代MLIP，未來結合加速分子動力學方法，有望在更接近工程應變率的條件下系統探索晶界特征與氫脆敏感性的構效關系。

7  總結與展望

本研究通過并行學習策略成功構建了Fe-H二元體系的BNNP型MLIP，實現了一般晶界氫偏聚、位錯發射、變形孿生及裂紋萌生擴展的DFT精度描述。主要成果可歸納為三點：（1）證明了系統覆蓋氫脆相關原子環境的訓練集設計策略可使MLIP對訓練集外構型（螺位錯-氫、一般晶界-氫）保持優異的外推精度；（2）從原子尺度揭示了一般晶界的兩種本征變形模式——全位錯發射與變形孿生，以及氫對兩者的差異化影響；（3）氫通過偏聚于晶界膨脹八面體位、提升位錯發射能壘的方式抑制全位錯發射并促進晶界斷裂，而對變形孿生的影響極其有限。

本研究為高強鐵基合金的抗氫脆晶界工程提供了原子尺度力學判據。未來工作將聚焦于：（1）研發精度-效率更優的MLIP以覆蓋更廣泛的應變率與晶界特征空間；（2）引入加速分子動力學方法實現更接近工程實際的低應變率模擬；（3）將MLIP構建方法學推廣至Fe-C-H等多元素體系，為碳中性背景下高強結構合金的理性設計提供高通量計算工具。

需注意的是，當前分析受限于高應變率條件（MD固有局限）和有限數量的晶界模型，且晶界遷移等因素在某些條件下（如GGB 3高氫濃度）亦發揮額外作用。實驗驗證方面，Sugiyama等人的觀測為本研究的位錯發射抑制結論提供了間接支持，但變形孿生對氫不敏感的預言有待進一步實驗檢驗。

原文閱讀：Machine learning interatomic potential reveals hydrogen embrittlement origins at general grain boundaries in α-iron | Communications Materials