第一性原理表面参数与Cu-Ni功能梯度材料中的梯度依赖性塑性：一种分子动力学方法

《Materials & Design》：First-principles surface parameters and gradient-dependent plasticity in Cu-Ni FGMs: A molecular dynamics approach

　　

随着微纳尺度器件日益集成于微机电系统（MEMS）、纳机电系统（NEMS）及能量转换系统中，对具有空间可调力学性能材料的需求显著增长。功能梯度材料（FGMs）因其组成连续变化而成为这些应用的潜在候选材料。然而在纳米尺度下，经典连续介质模型无法捕捉尺寸依赖性效应，尤其是表面主导现象，这造成了一个关键空白：从原子视角理解组成梯度如何影响力学行为。本研究利用分子动力学（MD）模拟，研究了具有幂律组成梯度的Cu-Ni功能梯度材料的力学响应。特别关注了连续方法无法触及的表面效应、位错演化和相变机制。研究人员对具有不同Ni体积分数的Cu-Ni功能梯度材料进行了拉伸、压缩和纳米压痕模拟。拉伸模拟获得的不同组分的杨氏模量与连续介质理论结果接近，误差小于6%。从纯铜到纯镍，力学性能逐渐提升，且体心立方（BCC）结构的峰值比例与杨氏模量呈负相关。研究人员发现，当BCC作为结构前驱体时，可能出现位错密度增加等异常现象。通过压缩和纳米压痕模拟推导了表面参数，值得注意的是，这些参数首次在原子层面被确定。结果表明，铜较低的层错能有助于早期应力松弛，宏观表现为硬度降低9%且塑性区扩大2.4倍。这为连续介质力学理论下的微尺度表面效应提供了来自原子层面的可靠表面参数。三点微弯曲实验（Cu-Ni梯度膜）与模拟挠度曲线的误差小于6%。研究还探讨了变形过程中的位错扩展，并对大量位错堆垛层错活动进行了分析，为功能梯度材料的微观力学响应和变形机制提供了原子视角。

### 论文解读：Cu-Ni功能梯度材料的第一性原理表面参数与梯度依赖性塑性——分子动力学视角

#### 研究背景与问题

随着微机电系统（MEMS）、纳机电系统（NEMS）及能量转换器件的微型化发展，对具有空间可调力学性能的材料需求急剧增加。功能梯度材料（FGMs）因组成沿特定方向连续变化，能有效缓解界面应力集中、提高断裂韧性，在航空航天、生物医学植入体及能源电子等领域具有广阔前景。然而当器件尺寸降至纳米级时，经典连续介质理论无法捕捉表面/界面主导的尺寸效应，导致对FGM力学行为的预测失效。表面弹性理论（如Gurtin-Murdoch框架）虽能描述纳米尺度表面效应，但其应用高度依赖准确的、组成依赖的表面参数（如表面能、表面残余应力、表面弹性模量）。目前，这些参数极难通过实验测量，也难直接通过第一性原理（如密度泛函理论）为复杂合金体系获取，形成了“描述性空白”：纳米尺度连续组成梯度如何影响表面能量学，进而改变位错成核与传播等变形机制？为填补这一空白，研究人员以Cu-Ni合金为模型体系，利用分子动力学（MD）模拟深入揭示原子尺度机制，并将结果与连续介质理论及实验对比验证。该论文发表于《Materials》。

#### 关键技术与方法

本研究采用分子动力学（MD）方法，使用LAMMPS软件进行模拟。原子间相互作用由嵌入原子方法（EAM）势函数描述，该势函数基于第一性原理热力学/力学状态计算，而非经验拟合宏观数据。构建了沿厚度方向遵循幂律分布的Cu-Ni功能梯度纳米梁模型（原子总数176,000，晶向[100]、[010]、[001]），选取梯度指数n=0.1、0.5、1、5及纯Cu、纯Ni作为对照。在300K、1atm下经NPT系综弛豫100 ps后，切换至NVT系综300 ps稳定温度，然后施加应变率10⁹ s^-1的拉伸/压缩以及0.1 ?/ps的纳米压痕模拟。后处理使用OVITO软件进行共近邻分析（CNA）和位错分析（DXA）。表面参数通过局部能量微分法和Shenoy提出的病毒应力定义从原子层面提取，表面区域定义为最外三层原子（0.6 nm）。所有模拟无实验样本队列来源，均为纯计算模型。

#### 研究结果

**1. 表面偏析主导梯度依赖性行为**
通过计算不同梯度指数下上、下表面原子间距的能量与局部Ni原子分数，发现上表面（名义富Ni区）的表面能（0.325–0.342 J/m²）远低于纯Ni（1.88 J/m²），而下表面（富Cu区）表面能（0.181–0.211 J/m²）接近纯Cu（1.25 J/m²）。这一矛盾直接证明：Cu原子因其更低的表面能而自发迁移至自由表面，即使在上表面也发生明显富集。这种热力学偏析改变了真实表面成分，使得表面残余应力τ₀（上表面0.3343 N/m）低于纯Ni预期，直接影响位错成核应力。

**2. 通道?阻塞统一变形框架**
通过位错密度和位错类型演化分析，揭示了梯度依赖的塑性机制：Cu富集区具有较低层错能（SFE，约45 mJ/m²），促进位错解离为Shockley不全位错，形成“通道”供位错优先成核；Ni富集区具有较高SFE（约125 mJ/m²），阻碍位错滑移，形成“阻塞”边界。在中等梯度（n=0.5、1）模型中，位错在成分过渡区堆积产生背应力，导致持续较高的流动应力；而陡梯度（n=5）模型中，狭窄过渡区导致位错强烈局域化，早期应力松弛。

**3. BCC作为异常结构前驱体**
在n=5及纯Cu模型的拉伸过程中，研究人员在应变约10%时观察到亚稳态体心立方（BCC）相先于位错出现（图5e、5a）。这一应力诱发马氏体相变产生了剪切模量较低的“软点”，集中后续变形。定量分析表明，BCC峰值比例与屈服强度呈强负相关（R2=0.883，图13a），且BCC形成速率（dφ_BCC/dε）与总应力降高度线性相关（R2=0.865），证明结构转变动力学是宏观软化的主要因素。

**4. 表面参数提取与多尺度关联**
采用局部能量微分法从原子层面首次提取了四种梯度模型的表面能、表面残余应力τ₀、表面密度ρ₀和表面弹性模量E_s。纯Cu和纯Ni的表面能与文献值误差分别仅为0.89%和0.11%，验证了方法准确性。表面残余应力随局部Ni浓度线性变化：τ₀(X_Ni) = 0.183 + 0.151 X_Ni N/m。将MD所得表面参数嵌入1 μm尺度FGM梁有限元模型，三点微弯曲实验与模拟挠度曲线误差小于6%，证实了原子级表面效应在连续介质框架中的可传递性。

**5. 纳米压痕揭示梯度依赖性塑性**
对上、下表面分别进行纳米压痕模拟发现：上表面（富Ni）压痕产生紧凑塑性区，位错被高SFE基体钉扎，硬度较高；下表面（富Cu）压痕因铜的低SFE导致位错均匀成核并快速滑移，塑性区扩大2.4倍，硬度降低9%。原子应力分析显示富Cu区存在长程应力线传播，而富Ni区应力仅集中于压头下方。

#### 讨论与结论

研究讨论部分指出：表面偏析效应使得真实表面成分偏离名义幂律分布，是理解梯度FGM表面参数的关键；通道?阻塞机制统一解释了梯度指数对强度与塑性的调控；BCC作为结构前驱体虽增加位错密度却降低强度，其反常行为源于相变耦合塑性。此外，表面参数的成功提取为Gurtin-Murdoch表面弹性理论提供了第一性原理输入，实现了从原子尺度到宏观力学响应的跨尺度关联。研究同时承认局限性：EAM势在极端条件下可能失效；175,000原子模型不足以捕获长程应力场；未考虑表面重构与各向异性。

**研究结论**（原文翻译）：
本研究利用分子动力学（MD）模拟系统探究了具有幂律组成分布的Cu-Ni功能梯度材料（FGMs）的力学行为。进行了拉伸、压缩和纳米压痕模拟，获得了关键力学性能参数，包括杨氏模量和表面参数。直接从原子模拟中提取了连续尺度表面效应理论所需的表面参数，实现了原子模型与宏观模型的多尺度关联。主要结论总结如下：
（1）力学性能预测：MD模拟获得的杨氏模量与连续介质理论预测及文献值吻合良好，偏差小于6%。随着Ni体积分数从20%增加至80%，杨氏模量从68.4 GPa提升至82.7 GPa，增加约21%。
（2）表面效应：铜的较低层错能（SFE）促进了早期位错成核与滑移，导致应力提前松弛。宏观上表现为硬度降低9%且塑性区扩大2.4倍。
（3）微观变形机制：变形过程由位错成核、滑移和湮灭主导，伴生相变。Shockley不全位错引起的交错堆垛层错显著增强了材料强度。然而当BCC相作为结构前驱体时，虽位错密度增加，但未观察到相应的性能提升。
（4）梯度依赖性行为：通过调整组成梯度可有效控制位错和堆垛层错的分布，这为通过原子级调控微观结构定制宏观力学性能提供了途径。
（5）验证与应用：Cu-Ni梯度膜的三点微弯曲模拟挠度曲线与实验结果误差小于6%。此外，将MD模拟获得的表面参数成功应用于1 μm尺度的有限元模型，证明了原子级表面效应在宏观力学响应中的可传递性。
总之，本研究从宏观和微观两个角度全面揭示了Cu-Ni功能梯度材料的力学行为，所提出的多尺度方法为微纳尺度功能梯度材料的设计与优化提供了稳健框架。