【摘 要】本文是在 电子 、原子层次上进行多元合金设计的一个初步探索， 研究 表明，由原子间相互作用势 计算 确定合适的基体成分，选择尝试的合金总成分，通过经验公式计算碳 化物体积分数，取得与合适的相结构对应的合金总成分, 由试验检验设计的合理性。这样的 一种思路是可行的。
关键词: 第一原理; 相互作用势; 多元合金; 合金设计
Abstract
This is a preliminary investigation of alloy design for multi-element alloys at electronic and atomic
scale. The research shows that this idea of alloy design was feasible that first calculating the properties
of matrix with different compositions by interatomic potential, secondly selecting the attempt total composition of the alloy, then calculating the carbide volume fraction by empirical formulae, until obtaining the appropriate total composition of the alloy corresponding to the desired phase structure, finally testing the design by experiment.
Keywords: first principles, Interatomic potentials, Multi-element alloy, alloy design
1．引言
目前 从电子、原子层次上进行材料设计是材料 科学 领域的学者们广泛关注的热点 问题 ， 主要研究 方法 有第一性原理方法，第一性原理赝势方法，原子间相互作用势方法，分子力学 方法，分子动力学方法及蒙特卡罗方法等。其中前两种方法是在电子层次上进行材料设计的 方法，其方法的物理基础可靠，但由于计算工作量很大，因而所计算的体系受到一定的限制。 后几种方法是在原子及分子层次上的设计方法，这几种方法不考虑电子结构的 影响 ，虽然会 损失一些精度，但大体上反映出由相互作用势所决定的晶体结构，以及由晶体结构所决定的 材料性质，且计算速度明显提高[1]。本文在多元合金的电子、原子层次的 理论 计算上联合使 用了第一性原理方法和原子间相互作用势方法，根据 3 种系列合金的关键问题进行理论计算 并结合其他理论计算和经验计算进行了合金设计探讨。
2．Fe-Cr-Mn-C 系亚稳奥氏体基铸造合金
Fe-Cr-Mn-C 系亚稳奥氏体基铸造合金有优异的耐磨性和高的抗冲刷腐蚀能力。其合金 设计的关键理论问题是 B 对该多元合金奥氏体体系的影响，以及对含 B 多元合金奥氏体电 子、原子层次的计算研究。通过对奥氏体合金大体系的能量计算，既可以解释 B 元素在奥 氏体中占位、分布、固溶度、与 C 的替代作用、与其他合金元素的配合对奥氏体的影响， 又可进而解释 B 对摩擦诱发马氏体相变的作用[2]。Fe-Cr-Mn-C-B 系铸造合金一般为基体（奥 氏体或马氏体）和碳硼化物组成的双相系统，其中 C 元素和 B 元素在奥氏体基体中的作用 对合金材料的性能有重要影响。C 在奥氏体中的固溶度和占位已很清楚，而 B 的固溶度和 占位还不很清楚，采用量子化学从头计算方法，通过对含 C、B 的奥氏体小团簇电子结构计 算，来研究 B 在奥氏体中的固溶度和占位情况。
团簇的选取以从奥氏体的实测晶格结构出发，从中选取奥氏体中的八面体和四面体小团 簇。并在团簇中心分别加入一个 B 或 C 原子，计算团簇的结合能，见表 1。由表 1 可见奥 氏体中 Fe4 四面体团簇的原子平均结合能略大于 Fe6 八面体的，因为八面体中存在距离较远 的三对原子，使结合能降低。B 和 C 均使四面体体积增大很多（棱长增加 34.2%），使结合 能降低，表明它们在四面体间隙存在的可能性很小，尤其是 B 更小。B 使八面体的体积增 大（棱长增加 16.2%），结合能略有升高；在体积不变时，C 使八面体的结合能略有降低， 表明 C 在奥氏体八面体间隙中的溶解度要远大于 B。
表 1 量子化学从头算合金小团簇的电子结构计算结果
采用量子化学从头计算精确处理含 B 小团簇，以半经验原子间相互作用对势处理大团
簇，研究含微量 B 元素的合金奥氏体大体系。对含微量元素小团簇进行局部精确计算，对 大团簇采用低精度的计算方法，既能反映微量元素的作用，也使电子、原子层次的计算处理 多元合金大体系成为现实。
应用 量子化学从头计算方法对含 B、C 奥氏体大体系进行局部精确计算，计算结果：B 在八面体间隙中的原子平均结合能为 1.6978eV，最近键距为 0.29967nm；C 在八面体间隙中 的原子平均结合能为 1.3520eV，最近键距为 0.25780nm；每个八面体间隙 B 原子使奥氏体
Fe 团簇 总结 合能降低 116.91443eV，每个 C 原子使小团簇总结合能降低 0.459142eV，B 原
子对奥氏体能量的影响是 C 原子的 254.6 倍。利用量子化学从头计算方法计算了
Fe-Cr-Mn-C-B 系双原子团簇的势能，提出按势能曲线最低点与从头计算所得的势能最低值 相重合的拟合原则，得出半经验原子间相互作用对势的参数，并计算了奥氏体团簇的平均结 合能、平衡原子间距，结果与实验符合。将量子化学从头计算方法精确计算小团簇和半经验 原子间相互作用对势处理大团簇相结合，计算 分析 B、C 元素在奥氏体中的间隙固溶度，得 出了 B 原子处于奥氏体的八面体间隙中时的固溶度为 0.097wt%，此时对奥氏体大体系能量 的影响贡献最大，B 在晶界和缺陷中存在对体系能量影响很小。表 2 是随含量增加合金奥 氏体的原子平均结合能的变化。
表 2奥氏体中其他元素近似不变时 B 含量变化引起的团簇原子平均能的变化
根据半经验原子间相互作用对势计算结果，B 在晶界上与 Fe 较易结合，B 在固溶体中
晶界上的存在几率很大，而且 B 在硼碳化物中的含量比率也很高，由此，可估算出奥氏体 基铁合金中 B 的加入量范围约为 0.05~1.00wt%。计算了奥氏体大团簇中 Cr、Mn、C 元素含 量固定时，B 的加入对团簇原子平均结合能的影响。随奥氏体中 B 量的增加，原子平均结 合能降低，当 B 含量增至 0.0427wt%时，原子平均结合能与纯铁奥氏体相比降低 10%，将此 时的 B 含量定义为 B 在该团簇中的极限含量，当其他元素含量改变时，B 在奥氏体中极限 含量将有所变化。随着结合能的降低，奥氏体更易摩擦诱发马氏体相变。图 1 分别是无 B 和含 B 合金磨损表面的 XRD 衍射谱。所研究开发的 Fe-Cr-Mn-C-B 系亚稳奥氏体基耐磨铸 造合金为新型耐磨材料（图 2）。 
（a）合金 1 （b）合金 4
图 1 Fe-Cr-Mn-C-B 系合金摩擦表面 XRD 图
图 2 Fe-Cr-Mn-C-B 奥氏体基合金的组织(1000×),
(a)合金 2, (b)合金 3, (c)合金 4