《Computational Materials Science》：Influence of grain boundary energy on hydrogen segregation tendencies
编辑推荐：
为了解决氢脆性（氢致延性损失）这一结构合金耐久性的关键挑战，研究人员开发并应用了一个结合分子动力学与巨正则蒙特卡罗采样的原子模型，系统研究了氢在多晶镍中的偏析行为。该研究首次揭示了氢偏析与晶界能之间存在强烈的依赖性，晶界能越高，氢偏析越显著。这一发现为预测和减轻氢脆性提供了原子尺度的洞见与定量模型基础，对于能源存储、燃料电池和航空航天等关键领域的材料设计具有重要意义。
氢气在金属中的渗入会显著改变其机械性能，导致材料延展性下降，即所谓的“氢脆”。这种现象在储氢系统、燃料电池、核反应堆和航空航天部件等富氢环境中尤为突出，常常引发材料的过早失效。那么，氢究竟是如何导致脆化的呢？其根源在于氢与材料内部微观结构缺陷的相互作用，比如晶界和位错，这些相互作用会改变局部的变形机制。一个广为接受的机制是氢增强局部塑性理论（Hydrogen-Enhanced Localized Plasticity, HELP）。然而，要精确理解这些相互作用，一个关键的科学问题在于：当存在多种不同类型的晶界等复杂缺陷网络时，氢会偏析（即聚集）在哪些具体位置？以及氢的分布如何影响其后的变形过程？
直接通过实验手段探测氢在晶界等缺陷处的精确位置非常具有挑战性。例如，热脱附光谱可以定量氢的总含量，但缺乏空间分辨率来定位特定缺陷位点；电子背散射衍射和透射电子显微镜可以揭示位错、晶界等变形机制，但无法直接探测氢原子；原子探针断层扫描虽然能提供原子级的化学图谱，但其微小的样品体积、低温要求及背景信号限制了其对全局氢分布的代表性。因此，通过计算方法来补充实验研究，从原子尺度阐明氢的行为变得至关重要。
为了填补这一空白，美国卡耐基梅隆大学材料科学与工程系的研究团队在《Computational Materials Science》上发表了一项研究，他们开发并应用了一种结合分子动力学和巨正则蒙特卡罗采样的原子模型，来阐明多晶材料中的氢偏析现象。他们首先通过第一性原理计算验证了所选原子间势函数在预测晶界稳定性及其氢亲和性方面的性能。然后，他们模拟了具有27个随机取向晶粒的多晶样品，研究了氢-晶界相互作用随氢浓度和晶界取向差角的变化关系。分析结果揭示了一个核心发现：氢偏析行为强烈依赖于晶界能。这项工作为在结构合金中建模氢的掺入以及后续研究其变形机制以减轻氢脆奠定了基础。
为开展这项研究，研究人员主要运用了以下关键技术方法：首先，他们使用分子动力学（Molecular Dynamics, MD）和巨正则蒙特卡罗（Grand Canonical Monte Carlo, GCMC）相结合的混合模拟方法，高效处理了氢原子的长程扩散与局域平衡问题。其次，他们利用第一性原理密度泛函理论（Density Functional Theory, DFT）计算，对经验原子间势的准确性进行了严格验证，确保了模拟的可靠性。最后，他们构建了两种样品体系：一是利用aimsgb工具库生成的12种不同取向差（16.26°至155.38°）的双晶模型，用于系统研究晶界能与氢亲和性的关系；二是通过Voronoi镶嵌法生成的五个具有随机晶粒取向、包含约31万个镍原子的多晶模型，用于考察复杂微结构下的整体偏析趋势。
2.1. 样品生成
研究人员系统性地生成了两类模型。对于双晶样品，他们利用aimsgb库精确构建了12种不同晶界，涵盖从16.26°到155.38°的取向差范围，并通过移除原子重叠来优化结构。随后，他们将样品尺寸放大，以包含2500、10,000和40,000个镍原子，用于氢加载模拟，确保结果不受有限尺寸效应影响。对于多晶样品，他们采用Voronoi镶嵌法构建了五个三维多晶模型，每个模型包含27个晶粒，尺寸为15 nm × 15 nm × 15 nm，平均晶粒尺寸为5 nm，总计约31万个镍原子。这些模型的晶粒取向是随机分配的，从而在数据集中获得了低角晶界和高角晶界近乎均衡的分布，最大取向差角达到70.53°。这些模型为全面研究氢在不同晶界类型上的偏析行为提供了多样化的微结构平台。
2.2. 第一性原理计算
为了验证后续分子动力学模拟中所选原子间势的可靠性，研究人员使用Quantum ESPRESSO代码对所有双晶样品进行了第一性原理计算。他们采用投影级加波赝势和Perdew-Burke-Ernzerhof交换关联泛函，计算了晶界能和氢亲和能，为经验势函数提供了高精度的基准数据。
2.3. 分子动力学模拟
研究采用了一种专为镍和镍-氢系统开发的嵌入原子法势函数进行分子动力学模拟。在通过共轭梯度算法对样品进行能量最小化后，研究进行了氢加载模拟。首先，样品在等温等应力（NσT）系综中于300 K和零外应力下平衡5纳秒。随后，在μVT系综中，通过结合巨正则蒙特卡罗和分子动力学的混合方法进行平衡，初始氢化学势μ0设为-2.547 eV。之后，化学势被阶梯式增加至μ1= -2.447 eV，模拟了不同氢压下的氢掺入过程。该过程在双晶和多晶样品中均得以实施，从而系统分析了氢在晶界处的吸收行为。
2.4. 后处理
在双晶样品中，研究人员利用OVITO软件中的公共近邻分析工具，根据原子的局部环境（如面心立方、体心立方、密排六方或“其他”）来识别晶界原子。对于Σ3(1?11)[110]孪晶界，密排六方原子被识别为晶界原子。随后，基于晶界处的镍原子数量计算了局域氢浓度。在多晶样品中，同样使用公共近邻分析直观观察氢沿晶界的偏析。为了确定氢偏析偏好的取向差角，他们选取单个氢原子，并将其选择范围扩展到1纳米截止距离内的所有原子，从而捕捉氢原子的局部原子环境，识别相邻晶粒及相应的晶界取向差角。最后，利用OVITO中的“创建键”功能（设定配对截止距离为1纳米），识别出参与至少一个键的独特镍原子，作为晶界原子的代表，进而计算局域氢浓度。
3.1. 原子间势函数的验证
分子动力学模拟的准确性高度依赖于原子间势函数的选择。为了确保可靠性，研究人员将经验势函数计算得到的晶界能与氢亲和能与第一性原理结果进行了对比验证。计算了每个双晶模型的晶界能εGB和氢亲和能ΔHεGB。结果表明，所选用的嵌入原子法势函数能够合理地预测不同晶界结构的能量趋势，从而为后续大规模的氢偏析模拟提供了可信的物理基础。
3.2. 连接Sieverts定律与热力学原理
研究将Sieverts定律与基本热力学原理联系起来，建立了氢浓度与化学势之间的关联框架，为解释材料中的氢偏析行为奠定了理论基础。Sieverts定律描述了稀溶液中气体的溶解度与压力的平方根成正比。通过热力学推导，研究将实验测量的氢浓度映射到相应的化学势上，确保了模拟的氢加载条件与实际实验条件相符。
3.3. 多晶样品中的氢结合
对五个多晶样品的模拟分析提供了氢在代表性微结构中偏析趋势的概览。氢原子被清晰地观察到优先偏析在晶界区域。统计分析表明，氢偏析的趋势与晶界取向差角（进而与晶界能）存在相关性。高能晶界（通常对应较大的取向差角）显示出更强的氢捕获倾向。这一结果在复杂的多晶环境中初步验证了晶界能作为氢偏析关键驱动因素的假设。
3.4. 双晶样品的详细研究
为了更细致地揭示晶界取向差和晶界能如何控制偏析行为，研究对选定的双晶样品进行了深入分析。通过系统改变样品的尺寸（从约2500到40,000个镍原子），研究证实了观察到的氢偏析行为对系统尺寸不敏感，排除了有限尺寸效应的影响。核心发现是：氢偏析程度与晶界能呈正相关。高能晶界（如某些高角晶界）为氢提供了能量上更有利的偏析位点，导致该区域氢浓度显著高于低能晶界（如某些特殊重合点阵晶界）和晶粒内部。这种关系通过计算局域氢浓度与晶界能的关联得以量化，为理解氢脆的微观起源提供了清晰的物理图像。
结论与讨论
本研究通过结合分子动力学与巨正则蒙特卡罗采样的多尺度原子模拟，系统揭示了金属镍中氢偏析行为与晶界能量之间的内在联系。核心结论是：氢偏析的倾向强烈依赖于晶界能，能量越高的晶界对氢的亲和性越强，从而导致更显著的氢局域富集。 这一发现将晶界的微观结构特征（通过取向差角等参数体现）与其宏观的氢脆敏感性直接关联起来。
研究的意义在于多个层面：首先，在科学认知上，它定量地建立了晶界能量作为预测氢偏析关键描述符的地位，深化了对氢-缺陷相互作用机制的理解，特别是对氢增强局部塑性等理论的微观解释提供了支持。其次，在方法论上，研究展示了结合第一性原理验证、分子动力学和巨正则蒙特卡罗的混合模拟框架，是研究复杂多晶体系中溶质偏析行为的有力工具，能够克服单纯实验或单一计算方法的局限性。最后，在工程应用上，这项工作为设计抗氢脆合金提供了新的指导原则。通过调控材料的晶界特征分布（例如，通过晶界工程引入更多低能晶界），可能有效降低氢在关键区域的偏析浓度，从而提升材料在氢环境下的耐久性和可靠性。这对于面向清洁能源（如氢能储运、燃料电池）和极端环境（如航空航天）的高性能结构材料开发具有重要价值。
