一种计算高效的方法用于预测高温下过渡金属合金的输运性质

《Computational Materials Science》：A computationally efficient approach for predicting the transport properties of transition-metal alloys at elevated temperatures

　　

本研究提出了一种新颖的现象学框架，用于高效估计体心立方（BCC）过渡金属浓缩合金在室温及高温下的热电性质。该方法被用于预测BCC系统的电阻率，与实验数据表现出极好的相关性。该框架进一步扩展到预测BCC合金在300–1300 K温度范围内的电阻率ρ、热导率κ和定容热容C_v，结果通过实验观测进行了验证。研究人员通过使用该模型预测一种浓缩W₅₃Ta₄₂V₅合金的热电性质，展示了模型的能力，该合金在300–1300 K温度范围内表现出ρ的饱和。随后，该模型被用于预测另一种浓缩Nb₄₀Mo₄₀Ta₂₀合金在相同温度范围内的性质。

本文解读一篇发表在《Computational Materials Science》上的研究，题为《一种计算高效的方法用于预测高温下过渡金属合金的输运性质》。该研究聚焦于过渡金属合金，尤其是难熔合金（如W、Ta、Nb、Mo基合金），这些材料因高熔点和优异高温机械性能在高温应用中极具潜力。然而，传统实验合成与表征方法因成本高昂和原材料稀缺而难以规模化，亟需可靠的计算框架来辅助甚至替代实验。现有第一性原理方法（如基于密度泛函理论（DFT）求解玻尔兹曼输运方程（BTE））虽精确，但计算成本高昂，不适用于高通量筛选复杂合金成分空间。因此，研究人员提出了一种半经验现象学框架，旨在高效预测体心立方（BCC）过渡金属浓缩合金在室温至高温（300–1300 K）范围内的热电性质，包括电阻率ρ、热导率κ和定容热容C_v。该模型基于Mott-Jones理论，将过渡金属中s-d轨道散射（即传导电子从s轨道散射到未占据d轨道的过程）与电阻率温度依赖性关联，并引入构型熵（S_conf）描述合金无序度对电阻率的贡献。通过实验数据拟合参数，模型成功预测了Nb₉₅Mo₅、Nb₉₀Mo₁₀、Ta₉₀W₁₀等合金的电阻率，并与实验值高度吻合。进一步，该框架结合BoltzTrap2和ShengBTE软件包，预测了电子和晶格热导率及热容，验证了Wiedemann-Franz定律的适用性，并检查了电阻率饱和（Ioffe-Regel极限）现象。研究人员还应用模型预测了W₅₃Ta₄₂V₅合金（已知实验数据）和Nb₄₀Mo₄₀Ta₂₀合金（设计用途为加热元件）的热电性质，其中W₅₃Ta₄₂V₅表现出零温度系数（ZTC）行为，即电阻率随温度变化极小。该研究为高效筛选高温热电性能优异的复杂合金提供了新途径，具有重要的工程应用价值。

**主要关键技术方法**：研究人员采用基于DFT的VASP（Vienna Ab initio Simulation Package）软件进行结构弛豫和电子结构计算，使用GW-sv赝势和500 eV平面波截断能，构建72原子超胞。通过BoltzTrap2求解BTE（固定弛豫时间近似）获得电子输运性质（如κ_e），并利用ShengBTE数值求解声子BTE（基于Phonopy计算二阶力常数和Thirdorder.py计算三阶力常数）获取晶格热导率κ_l和热容C_v。电阻率预测基于经验公式ρ = ρ_o + ρ_r^alloy，其中ρ_o通过s-d散射参数A（从DFT态密度几何光滑化得到）计算，ρ_r^alloy通过构型熵S_conf拟合实验数据得到。所有合金超胞使用ATOMSK生成，计算在零压下进行。

**研究结果**：

**4.1 电阻率ρ**：通过DFT计算A参数并利用经验公式，研究人员预测了Nb₉₅Mo₅、Nb₉₀Mo₁₀和Ta₉₀W₁₀在300–1300 K范围内的电阻率，与实验数据相比，均方根误差（RMSE）分别为5.92、4.14和3.22 μΩ·cm，平均绝对百分比误差（MAPE）分别为11.9%、8.87%和6.29%，相关系数R²分别为0.71、0.89和0.92，表明模型预测与实验高度相关。

**4.2 热导率κ和热容C_v**：对Ta₉₀W₁₀，研究人员结合BoltzTrap2和ShengBTE计算了电子和晶格热导率，预测总热导率κ与实验值对比，RMSE为5.69 W/m·K，MAPE为11.93%，Pearson相关系数为0.994，验证了模型的有效性。

**4.3 电阻率饱和与Wiedemann-Franz关系有效性检查**：通过计算弛豫时间τ与温度T的幂律指数K（τ ∝ T^K），研究人员发现W₅₃Ta₄₂V₅的K值接近0.13（低于阈值0.2），表明该合金在300 K以下已进入Ioffe-Regel极限，电阻率饱和。而其他合金（如Nb₉₀Mo₁₀）K值均大于0.77，未饱和。同时，Wiedemann-Franz定律的验证显示所有合金的Lorenz数L在2–2.6×10^-8 W·Ω·K^-2范围内，符合要求。

**4.4 W₅₃Ta₄₂V₅系统的实验验证**：基于实验κ和计算κ_l，研究人员反推κ_e并利用Wiedemann-Franz定律计算ρ，发现ρ在300–1300 K内几乎恒定（约20.7–24.2 μΩ·cm），证实了电阻率饱和行为，类似于零温度系数（ZTC）合金，具有精密应用潜力。

**4.5 合金设计Nb₄₀Mo₄₀Ta₂₀**：研究人员预测该合金在300–1300 K内电阻率较高（112–134 μΩ·cm），热导率近乎恒定（19–25 W/m·K），且未出现电阻率饱和（K值未计算但假设高），适合作为高温加热元件。该合金的A参数为正（1.81×10^-7 K^-2），表明电阻率随温度上升缓慢，且s-d散射在1400 K以下未饱和。

**讨论与结论总结**：讨论部分指出，模型忽略了热膨胀效应、A参数的温度依赖性、德拜温度Θ_D的变量、短程有序（SRO）以及晶体缺陷（如位错、晶界）和磁性元素的影响，这些因素可能引入误差，但当前近似在温度范围内捕获了正确趋势。模型仅针对BCC结构开发，但物理基础可扩展至其他晶体结构，需进一步工作。结论部分翻译：本研究提出了一种物理启发的框架，用于高效估计BCC浓缩合金的热电性质，以极低计算成本实现了对BCC过渡金属及其合金的预测。通过少数单相合金的室温电阻率验证（见表1、表3），结果良好。由于缺乏300–1300 K范围内单相BCC固溶体的实验数据，仅对稀释二元BCC合金进行了温度依赖性验证，模型与这些系统高度相关，但仍需对更复杂浓缩合金进行验证以增强鲁棒性。W₅₃Ta₄₂V₅已合成并有κ和C_v数据，计划后续测量ρ。Nb₄₀Mo₄₀Ta₂₀被预测为高温应用（1300 K）的可行候选，但需实验验证电阻率饱和。若成功验证，该方法可加速开发具有特定热电性质的浓缩合金。