W(1-x)-TM(x)合金的结构、稳定性、力学及热物理性质——分子动力学(Molecular Dynamics, MD)研究

《Results in Materials》：Structural, Stability, Mechanical, and Thermophysical properties of W(1-x)-TM(x) - A Molecular Dynamics study

　　

研究人员采用分子动力学(Molecular Dynamics, MD)模拟研究了W1-x-TMx合金(TM＝Cr、Fe、Mo、Ni、V，x＝0～10 at.%)的力学、热力学及稳定性特征。与多数只关注单一合金体系的研究不同，本研究在同一计算框架与浓度范围内对多种W–TM合金进行了统一系统的对比，可直接对不同合金化效应进行排序。研究人员通过弹性常数(Cij)、体模量(Bulk Modulus, B)、剪切模量(Shear Modulus, G)、B/G比(Pugh's Ratio)、柯西压力(Cauchy Pressure, C12–C44)、泊松比(Poisson's Ratio, ν)评估钨基合金的力学性能，并通过德拜温度(Debye Temperature, ΘD)、最小热导率(minimum thermal conductivity, kmin)及熔点(melting point, Tm)分析其热学特征。结果表明二元W–TM合金的力学强度低于纯W，尤以W–Cr、W–Ni及W–Fe合金下降明显；但引入过渡金属可提高B/G比与泊松比，说明合金化能显著改善纯W的延展性(Ductility)，其中W–Ni合金对金属键(Metallic Bonding)的强化作用尤为突出。W–V掺杂对德拜温度影响最大，表现出较好的耐高温性。此外，仅W–Mo合金可在添加后提高纯W的热导率(Thermal Conductivity)，而过渡金属的加入均使纯W的熔点有所降低。

论文解读：W_(1-x)-TM_(x)合金的结构、稳定性、力学及热物理性质——分子动力学研究

钨(Tungsten, W)作为难熔金属(Refractory Metal)，因具备极高熔点(~3683 K)、优良热导率、低溅射产额及高热强韧性，被广泛视为聚变堆等离子面材料(Plasma-Facing Materials, PFMs)的首选候选。然而纯W存在低温辐照脆化、再结晶温度低、高韧脆转变温度(Ductile-to-Brittle Transition Temperature, DBTT)导致低温延展性差、断裂韧性不足等固有缺陷，限制了其工程化应用。传统解决方案是对纯W进行过渡金属(Transition Metal, TM)合金化以改善性能，但以往研究多聚焦单一W–TM体系，且所用计算方法或实验条件各异，缺乏统一框架下对不同TM掺杂效应的直接对比与排序，尤其W–Ni合金在分子动力学层面的力学与热物理性质研究尚不充分。为此，研究人员在同一分子动力学模拟框架下，系统考察了Cr、Fe、Mo、Ni、V五种TM以0～10 at.%置换掺杂W_1-xTM_x体心立方(body-centered cubic, bcc)结构时的力学、热力学及稳定性变化，为PFMs用钨基合金的合理成分设计提供统一对比数据库。该论文发表于《Results in Materials》。

主要关键技术方法：研究人员采用LAMMPS(Large-scale Atomic/Molecular Massively Parallel Simulator)软件包开展分子动力学模拟，借助Atomsk构建20a×20a×20a含16000原子的bcc超胞，按置换式随机分布引入TM原子(TM＝Cr、Fe、Mo、Ni、V，x＝0、2、4、6、8、10 at.%)。W–Cr与W–Fe体系选用Bonny等人Embedded Atom Method(EAM)势，W–V与W–Mo体系选用Chen等人Finnis-Sinclair形式势，W–Ni体系选用Shim等人Modified Embedded Atom Method(MEAM)势。先经共轭梯度法应力弛豫，NVT(300 K, 10 ps)与NPT(20 ps)系综平衡后施加微小应变求取弹性常数，进而计算体模量(B)、剪切模量(G)、杨氏模量(Young's Modulus, E)、Pugh比(B/G)、泊松比(ν)、柯西压力(C₁₂–C₄₄)、维氏硬度(Vickers Hardness)、声速(纵波v_l、横波v_t、平均声速v_m)、德拜温度(Θ_D)、最小热导率(k_min)及由C₁₁经验关联估算的熔点(T_m)。所有结果相对纯W归一化作图。

3.1. Validation of the pure W potentials

研究人员先计算纯W的晶格常数a₀、密度ρ、弹性常数C₁₁/C₁₂/C₄₄及模量B/G/E/ν并与实验值和前人理论值比对，三种势函数所得结果与文献吻合良好，验证了所选相互作用势及模拟参数的可靠性。

3.2. Mechanical property calculations

• •
  3.2.1. Elastic constants：随TM浓度升高，C₁₁整体线性下降，W–Mo降幅最缓且C₁₁最高(W–Mo＞W–V＞W–Ni＞W–Cr＞W–Fe)；C₁₂在W–Cr/Fe/Ni/V中上升而W–Mo基本不变；C₄₄变化幅度小，W–Cr与W–V近无变化。表明TM主要影响C₁₁，Mo对轴向刚度影响最小。
• •
  3.2.2. Mechanical stability：所有W–TM合金在0～10 at.%内满足Born–Huang立方晶系稳定性判据(C₁₁＞0, C₄₄＞0, C₁₁－C₁₂＞0)，bcc结构力学稳定，但稳定性随TM浓度微降。
• •
  3.2.3. Elastic moduli：体模量B、杨氏模量E、剪切模量G均随掺杂浓度增大而降低，降幅W–Fe＞W–Cr＞W–Ni＞W–V≈W–Mo，W–Mo与W–V的B和G下降极微。表明所有二元W–TM合金力学强度低于纯W，Fe弱化最强，Mo与V影响最小。B＞G说明剪切模量是结构稳定性主导因素。
• •
  3.2.4. Pugh's ratio(B/G)：所有W–TM合金B/G＞1.75(纯W≈1.917)且在阈值以上属延性区，随浓度升高B/G上升，W–Ni提升最显著，W–V最弱(10 at.%时B/G≈1.941)。说明TM掺杂可有效改善纯W延展性，效果排序：Ni＞Fe＞Mo≈Cr＞V。
• •
  3.2.5. Poisson's ratio(ν)：ν随TM浓度升高而增加，W–Ni最大，W–V最小，趋势与B/G一致，进一步证实合金化提升塑性。
• •
  3.2.6. Cauchy pressure(C₁₂－C₄₄)：五种W–TM合金柯西压力均为正且大于纯W，随浓度增大而增大，W–Ni增幅最大，W–Cr与W–V最小。正值表明金属键占主导，TM掺杂强化了W的金属键特征从而增强延性，与B/G及ν结论吻合。
• •
  3.2.7. Hardness：基于Chen模型及ν–E关系算得硬度随TM掺杂微降，W–V保留最高硬度接近纯W，W–Ni降幅最大。硬度降低对应延展性提升，符合强度–延性权衡规律。

3.3. Thermal property calculation

• •
  3.3.1. Sound velocity：平均声速排序W–V＞W–Mo＞W–Cr＞W–Fe＞W–Ni；W–V与W–Mo的纵/横波声速随浓度略升或基本不变，其余下降。高声速反映较强原子键合与较硬晶格。
• •
  3.3.2. Debye temperature(Θ_D)：Θ_D变化趋势同声速，W–V最高，W–Mo次高，W–Ni最低。说明V掺杂最显著提升晶格热稳定性与高温耐受性，W–V和W–Mo共价成分较强。
• •
  3.3.3. Minimum thermal conductivity(k_min)：仅W–Mo合金的最小热导率不低于纯W，其余TM掺杂使其下降或持平；各体系k_min在浓度范围内变化不显著。
• •
  3.3.4. Melting point(T_m)：纯W熔点最高，TM掺杂使T_m降低，W–Mo降幅极小，W–Fe降幅最大，W–Ni与W–V居中且相近。熔点降低源于TM打乱原始bcc点阵周期性与键强削弱。

结论(Conclusions)总结翻译：

研究人员通过统一分子动力学框架系统考察了五种过渡金属(Cr、Fe、Mo、Ni、V)在低掺杂量(0～10 at.%)下对钨力学、热力学及稳定性的影响，主要结论如下：(1)TM掺杂改变W弹性响应——C₁₁总体下降、C₁₂多上升，Mo对C₁₁影响甚微且C₁₂略降，C₄₄在W–Cr与W–V中基本不变；(2)所有W–TM合金在考察浓度内满足力学稳定性判据，低浓度合金化不破坏bcc W基体完整性；(3)TM杂质降低纯W力学强度，Fe/Cr/Ni影响最显著，以牺牲强度换取延性；(4)合金化显著提高Pugh比(B/G)与泊松比，延展性改善效果Ni＞Fe＞Mo＞Cr＞V，与金属键强化一致；(5)TM添加增强W金属键特征、削弱共价键组分，W–Ni金属键强化最明显；(6)Ni对声速与硬度影响最突出，V影响最小；(7)V对德拜温度提升最大，W–V具较好高温稳定性潜力；(8)仅W–Mo可维持或微增W热导率，其余TM无此效果且k_min变化不显著；(9)纯W熔点最高，TM降低T_m但W–Mo降幅很小，适合高温应用。