《Materials & Design》：Influence of nitrogen flow on the mechanical properties of multilayers for the design of metal nitride coatings
编辑推荐：
自20世纪80年代末以来，多层涂层已被广泛用作高性能保护涂层，特别是在切削刀具领域。研究表明，在涂层中交替堆叠韧性材料与刚性材料可增强其力学性能[1], [2]。这些涂层的力学性能在很大程度上取决于堆叠周期：通常周期越小，力学性能越优。Misra等人[3]的研
自20世纪80年代末以来，多层涂层已被广泛用作高性能保护涂层，特别是在切削刀具领域。研究表明，在涂层中交替堆叠韧性材料与刚性材料可增强其力学性能[1], [2]。这些涂层的力学性能在很大程度上取决于堆叠周期：通常周期越小，力学性能越优。Misra等人[3]的研究表明，Cu/Nb多层膜在不同周期下表现出显著不同的硬度，同时保持相似的压痕模量。基于Hall-Petch效应[4]，研究人员建立了一个考虑层间界面引起的受限塑性的模型。这种硬度增强还可伴随韧性的提高[5], [6], [7], [8], [9], [10]。在某些情况下，压痕弹性模量也会受到影响[11]，表明层间界面不仅充当位错屏障。因此，界面不能再被视为二维表面边界，而应被视为具有自身厚度和性能的三维区域。当堆叠周期达到纳米尺度时，这种通常具有纳米厚度的界面层对整个多层膜性能的影响将占据主导地位。目前，对这种纳米尺度的界面层进行力学和结构表征变得至关重要[12]。
对于纳米级周期，界面层厚度过薄，无法通过实验手段进行力学表征。近期，Giljean等人[13]提出了一种通过标定不同堆叠周期多层膜压痕实验曲线的有限元模型来评估界面层厚度的方法。然而，该方法需要对界面层所采用的力学本构规律做出较强假设，并且由于微观结构的精细程度需要较大的计算资源。本研究采用了不同的方法。与Giljean等人类似，研究人员提出一种基于压痕曲线的模型标定方法，但使用的是由多层涂层两种组成相（TiAl和TiAlN）构成的等效均匀涂层模型，并调整两者的体积分数。
密度泛函理论（DFT）计算成为评估这些极薄界面层微观结构和力学性能的另一种选择[14], [15]。近期，Zhang等人[12]利用此类计算研究了面向多层结构的多种类型界面的力学响应。然而，其研究集中于不同层间具有突变界面的理想界面，未考虑两层之间可能形成第三种材料的情况，尽管此类计算方法完全可以处理这一问题。事实上，通过堆叠两层形成超结构后，利用分子动力学（MD）计算[16]进行弛豫以实现热力学稳定，可以重现涂层两层之间的界面。随后，可借助半经验方法对此类结构进行力学表征[17], [18], [19], [20], [21], [22], [23]。本研究中，DFT计算用于预测多层膜两种组成相（TiAl和TiAlN）的结构和力学性能，并以努氏硬度为优化目标对结构进行了优化。
本文的结构安排如下：第2节概述了研究中所使用的材料与方法，其中2.1节描述所研究样品的沉积方法及其特征；2.2节给出了X射线衍射（XRD）表征的细节；2.3节给出了电子能量损失谱（EELS）和透射电子显微镜（TEM）表征的细节；2.4节介绍了通过纳米压痕对样品进行力学表征的过程，以及用于确定其力学性能的压痕有限元模型（FEM）和有限元模型更新（FEMU）方法的描述；2.5节详细介绍了从头计算（ab-initio）方法，特别是通过Voigt、Reuss和Hill方法计算努氏硬度以及弹性张量得出材料性能的方法，以及采用Oganov等人模型计算维氏硬度的方法；最后，2.6节描述了界面的特征和分子动力学计算的条件。第3节描述了本工作中获得的不同结果及相关讨论。首先，3.1节展示了与微观和纳米结构表征相关的XRD、TEM和EELS结果；其次，3.2节介绍了纳米压痕结果，包括样品的压痕曲线；随后，3.3节展示了均匀化模型获得的结果；其后，3.4节展示了TiAl和TiAlN结构的DFT计算结果；最后，3.5节展示了分子动力学的结果，包括界面弛豫内能演化、不同类型原子随时间的位移，以及单独界面与TiAlN和TiAl等比例材料力学性能的对比。结论部分总结了所获得的所有结果。
一、研究背景与问题
自20世纪80年代末以来，多层涂层作为高性能保护涂层在工业界得到广泛应用，尤其在切削刀具领域。研究表明，在涂层中交替堆叠韧性材料与刚性材料可以显著增强其力学性能，且这种性能高度依赖于堆叠周期——周期越小，力学性能越优异。Misra等人基于Hall-Petch效应提出了层间界面导致的受限塑性模型，解释硬度随周期变化的规律。随着周期降低至纳米尺度，界面层不再仅仅作为位错屏障作用于二维表面，而是成为具有自身厚度与性能的三维区域，其对整个多层膜性能的影响趋于主导。然而，纳米尺度界面层过薄，无法通过实验手段直接进行力学表征。Giljean等人曾提出通过有限元模型标定压痕实验曲线来评估界面层厚度的方法，但该方法对界面层力学本构规律假设较强，且计算成本较高。为此，研究人员提出了一种基于等效均匀涂层模型的实验—数值联合标定方法，并结合密度泛函理论和分子动力学计算，从微观尺度考察界面结构与性能演变，以揭示多层涂层力学性能增强的深层机理。
二、研究内容与结论
本研究以射频磁控溅射结合反应气体脉冲工艺（RGPP）沉积的TiAl/TiAlN多层涂层为对象，系统研究了不同堆叠周期（Λ=4、10、16和50 nm）下多层涂层的微观结构、化学分布和力学性能。研究发现，随着周期减小，TiAlN层在多层结构中的体积分数显著增加（由XRD、TEM/EELS和两种模型标定方法共同确认），表明氮气脉冲关闭时间内未能完全清除腔室中的残留氮气，导致TiAl金属层被部分氮化。DFT计算获得了TiAl和TiAlN的晶体结构及力学参数，MD弛豫模拟则显示界面处存在氮原子从TiAlN向TiAl层的迁移，且迁移幅度不超过0.5 ?。界面力学性能计算表明，按[200]c方向生长的I200界面比按[111]c方向生长的I111界面具有更优的力学性能，且两者的性能均优于按等比例混合物规则预测的TiAlN/TiAl复合材料，说明界面结构本身带来了额外的力学增强效应。更重要的是，研究建立了TiAlN体积分数Vf与多层周期倒数1/Λ之间的线性关系，并给出了硬度与屈服强度的线性关系，使得仅通过周期值或硬度测试即可预测任意TiAl/TiAlN多层涂层的完整力学性能参数。
三、关键技术方法
本研究综合运用了多种实验表征与理论计算方法。沉积方面采用Alliance Concept AC450反应腔体通过射频磁控阴极溅射结合反应气体脉冲工艺（RGPP）制备单层和多层涂层，以(100)硅为衬底，TiAl烧结靶（Ti/Al=2）为溅射源；结构表征采用X射线衍射（XRD）确定晶体结构和择优生长方向，使用配备CuKα1源的Rigaku Smartlab衍射仪进行θ–2θ扫描；化学与微观结构分析采用扫描透射电子显微镜（STEM）结合电子能量损失谱（EELS）沿生长方向获取纳米尺度元素分布；力学表征采用Berkovich压头的超纳米压痕仪进行恒应变速率加载和连续刚度测量，并结合有限元模型更新（FEMU）方法标定塑性参数；理论计算方面利用USPEX软件耦合VASP进行基于DFT的结构预测和努氏硬度优化，通过Voigt-Reuss-Hill方法计算多晶弹性模量，采用Mazhnik/Oganov半经验模型计算维氏硬度和断裂韧性；界面研究采用PyMatgen构建TiAlN/TiAl界面，基于SIESTA程序进行分子动力学（MD）弛豫以获得热力学稳定界面结构。
四、主要研究结果
（一）微观和纳米结构表征（XRD、TEM和EELS）
XRD结果表明，单层TiAl呈现六方结构，主要衍射峰对应[002]h方向；单层TiAlN呈面心立方结构，具有[111]c和[200]c两个择优生长方向。10 nm和16 nm周期多层膜呈现典型的超晶格卫星峰特征。STEM图像清晰显示多层结构交替排列。EELS浓度分析表明，50 nm周期样品中金属层和氮化物层界面清晰、过渡区仅数纳米；而4 nm和10 nm周期样品中，TiAl金属层内检测到显著的氮元素存在（约为35%甚至60%以上），说明短周期条件下金属层被部分氮化。
（二）深度传感压痕的解析分析
纳米压痕实验表明，硬度与压痕模量均随周期减小而增加。研究人员采用受King模型启发的薄膜模型对压痕数据进行解析拟合，通过混合规则（rule of mixtures）确定TiAlN体积分数Vf。结果表明，所有多层涂层的Vf均大于50%，从50 nm周期的59%增加到4 nm周期的98%，反映了短周期下TiAlN占比的系统性增大。将多余TiAlN转化为等效界面附加厚度后，发现该厚度并非恒定值，暗示不同周期下纳米尺度的增厚机制存在差异。
（三）均匀化多层涂层模型得到的TiAlN体积分数
为克服逐层建模在有限元计算中的高昂成本，研究人员提出了单一均匀化涂层模型，通过TiAl和TiAlN串联混合规则描述等效弹塑性行为，并以FEMU方法标定各周期多层膜的压痕力—位移曲线，获得VfFEM值。FEMU识别的Vf值与解析模型结果吻合良好（差异范围为2%–22%），进一步验证了两种方法的可靠性。两种方法一致表明短周期涂层中TiAlN占主导，由此解释了超硬效应的来源。
（四）DFT计算TiAl和TiAlN的性能
DFT优化的TiAl和TiAlN晶体结构参数与XRD测定值接近，TiAl六方晶格参数与实验差异分别为c方向2%和a方向7%，TiAlN面心立方晶格参数与实验完全一致。弹性性能方面，DFT计算的TiAl/TiAlN弹性模量比值（0.48）与纳米压痕实验比值完全一致，验证了理论方法在相对比较中的有效性。然而硬度方面，DFT计算的TiAlN硬度（15.8 GPa）约为实验值（31 GPa）的一半，比值差异（理论0.47 vs 实验0.32）主要由半经验塑性理论模型的局限所致。
（五）分子动力学重排界面的力学与结构特征
分子动力学弛豫显示，I200界面比I111界面具有更低的内能绝对值，即I200在热力学上更稳定。原子位移分析表明，N和Al主要在Z方向迁移，而Ti在X和Z方向均有较大位移，N原子与Ti/Al原子在Z方向出现约0.5 ?的差分位移。氮浓度沿深度分布显示，无论何种界面取向，从金属到氮化物的氮浓度梯度均小于从氮化物到金属的方向，表明前者界面过渡更缓和。界面力学性能计算表明，I200界面（E=294.2 GPa，H=14.3 GPa）优于I111界面（E=260.7 GPa，H=12.7 GPa），两者的性能均明显高于等比例串联或并联混合规则预测的TiAlN/TiAl复合材料性能，体现了界面结构带来的额外增强效应。
（五）体积分数与周期及力学性能的关联
研究人员将四种周期多层膜的硬度H和屈服应力σy分别对1/Λ和1/√Λ作图，发现对1/Λ的线性相关性明显优于1/√Λ，支持了短周期下强化机制更符合界面主导而非纯粹Hall-Petch效应的观点。Vf与1/Λ呈现良好的线性关系，且MD计算的Vf范围与实验标定结果一致，为仅由周期预测多层膜的完整力学参数提供了经验基础。此外，硬度与屈服应力之间存在线性关系（斜率1.9），与经典的Tabor关系（斜率3.3）不同，反映了TiAl/TiAlN材料体系的特殊性，但该关系可用于从硬度测量直接推算屈服应力。
五、讨论总结
研究人员综合利用实验表征、解析模型、有限元标定和DFT/MD多尺度计算，系统揭示了TiAlN/TiAl多层涂层的力学性能增强机制。首先，自上而下的均匀化方法在宏观层面定量评估了TiAlN/TiAl之间的体积失衡，解析模型与FEMU两种方法在独立实验观测量下给出了一致的结果，验证了该方法的可靠性。其次，自下而上的DFT/MD方法证实了氮原子从TiAlN层向TiAl金属层的迁移，且这种迁移可能改善整个多层涂层的力学性能。然而MD结果表明氮原子迁移距离不超过0.5 ?，提示沉积过程中的氮污染是主要的强化机制。结合STEM观察确认多层结构、以及EELS浓度曲线显示短周期下各层中均存在氮元素，可以得出结论：在最短周期下，这些涂层不再是严格意义上的金属/陶瓷多层膜，而更像是具有TiAlNx/TiAlN（x<1）结构的准均匀材料。研究最后提供了完整的预测工具，仅通过多层膜的周期值或硬度测量即可获得弹性模量、初始屈服应力和硬化模量等全部力学参数，为TiAlN/TiAl多层涂层的设计与应用提供了便捷而可靠的工程途径。
六、结论翻译
本文通过实验、解析和数值方法相结合的多尺度策略，系统研究了射频磁控溅射沉积TiAlN/TiAl多层涂层的力学响应。主要结论如下：第一，基于多层膜均匀化方法的自上而下方法在宏观层面评估了TiAl与TiAlN之间的体积失衡，解析模型（结合硬度和压痕弹性模量随深度演变曲线的拟合）和FEMU有限元标定（对不同周期多层膜的压痕曲线）两种途径给出了非常相似的结果，且基于不同的实验观测量，相互验证了有效性。第二，通过DFT和分子动力学的自下而上方法证实了氮原子从TiAlN层向TiAl金属层的迁移，且这种迁移可能增强整个多层涂层的力学性能。然而，MD结果表明氮原子迁移不超过0.5 ?，提示主要的强化机制是沉积过程中发生的氮污染。因此，结合STEM对涂层多层结构的确认以及EELS浓度分布显示短周期下各层中均存在氮元素的事实，可以认为对于最短周期而言，研究对象已不再是严格意义上的金属/陶瓷多层膜，而是一种具有TiAlNx/TiAlN（x<1）结构的准均匀材料。最后，本研究提供了完整的预测工具，可仅由周期值或硬度测量值预测任意TiAlN/TiAl多层涂层的力学性能（弹性模量、初始屈服应力、硬化模量）。
