氢在镁中的双重作用：从孪生促塑到晶界开裂
DOI：10.1016/j.ijplas.2024.104084
推荐理由
氢不总让镁更脆。Ji等在特定原子氢环境下测得纯Mg延伸率增加，同时强度下降。论文把表面入氢、晶内孪生和晶界开裂联系起来，说明判断氢的作用，必须看它在哪里、怎样参与变形。
作者研究思路
作者结合拉伸/EBSD、第一性原理（DFT）、原子模拟和晶体塑性（CPFEM）。贡献是用晶面能量变化约束含氢参数，并借助裂尖与多晶模型检验机制，再连接宏观响应。
CPFEM依据EBSD建立587晶粒、500×500μm二维模型；空气参数由曲线标定，含氢参数结合DFT缩放。它并未显式求解氢扩散与裂纹演化。
核心公式
式（7）：相对位移引起的能量变化除以界面面积，得到广义层错能；路径峰值用于比较变形势垒。
式（13—15）：采用已有晶体塑性框架，分别描述母相滑移、孪生转变与孪晶内再滑移。
式（16—18）：g为临界分解剪应力，孪生只取正向驱动；0.80是模型上限，不是实测孪晶含量。
核心结果图
图 2(a)：10⁻⁶s⁻¹、pH14、−1.6V相对SCE，各条件3次。延伸率19.46%→22.06%，增加2.60个百分点，抗拉强度降低。
图 3(b)(c)(e)：依次为空气、含氢断后组织，白色标尺均为100μm；孪晶面积分数1.06%→2.84%，支持孪生参与增塑的解释。
图 6：柱面I入氢能垒0.156eV/H，低于所算基面及锥面I；不能仅凭能垒判断整体耐蚀性。
图 9(a)(b)：基面势垒降18.5%，稳定层错能为−7.13mJ/m²；柱面I势垒反增4.5%，作用有晶面选择性。
图 10(b)(c)：近0K裂尖模型中，含氢后裂纹偏折并发射位错。原σxy色条单位为eV/ų。
图 11(a)(b)(d)(e)：2H沿路径与6H环绕裂尖的设置同时改变数量、位置，后者层错发射反而延后。
图 12(a)(f)：近0K将氢保留在晶内，含氢模型至50%模拟应变未见孔洞；这不是实测延伸率。
图 13(a)(c)：10%应变下，含氢模型孪生更多、应力较低；曲线吻合仍需局部实验进一步检验。
图 14(b)(c)：全局孪生体积分数低于1.8%；增长速率差异主要集中在前4%应变，勿与局部场混同。
图 16(c)：Mg-H仍取约200ppm。300K下氢向晶界富集，24%模拟应变出现晶界孔洞，与近0K情景相反。
作者核心结论
晶内氢可降低部分变形势垒、促进孪生；晶界富集又可能诱发损伤。氢效应取决于位置、晶面和变形条件，延伸率提高不代表强度与耐蚀性同步改善。
对我们的启示
可借鉴能量约束本构、局部组织检验机制的思路。原子模拟约10⁸s⁻¹，与慢拉伸相差巨大；应用时还需可测氢分布与晶界证据，避免把特定情景直接当成工程加氢方案。
来源：Ji等，International Journal of Plasticity 181（2024），104084。公式及结果图为原文节选、重排，版权归原权利人。
