超越平均化：δH——多组分合金单相形成预测的化学异质性描述符

《ACTA MATERIALIA》：Beyond Averaging : δ H, a Chemical Heterogeneity Descriptor for Single-Phase Formation Prediction in Multicomponent Alloys

　　

摘要

研究人员提出了一种用于多组分合金（MCA）和高熵合金（HEA）单相微观结构预测的新型无量纲描述符：二元混合焓的归一化离散度δ_H。该描述符量化了组成元素对之间化学相互作用的内部异质性——这是经典描述符无法捕捉的贡献，因为经典描述符依赖于加权平均值，从而掩盖了强吸引对与强排斥对共存的事实。与原子尺寸错配参数δ相结合，δ_H–δ二元组合显著优于Yang和Zhang（2012）建立的Ω–δ参考判据，在来自三个独立来源、涵盖二元至九元成分的1675种合金实验数据集上，交叉验证准确率提升了+3.16个百分点（83.16%对80.00%），马修斯相关系数（MCC）提升了+0.063（0.663对0.600）（McNemar检验，p = 0.006）。尽管δ_H与δ在整体水平上呈中等程度相关（Spearman ρ = 0.27），部分原因在于两者均随成分复杂性增加而增大，但δ_H–δ相较于任何单变量分类器的显著改进证实，δ_H贡献了δ单独无法捕捉的预测信息。关键阈值δ_H ≤ 1.1具有直接的物理意义：当二元化学相互作用的异质性不超过凝固温度下可获得的平均热能（R_Tm）的1.1倍时，固溶体是稳定的。

论文解读：化学异质性描述符δ_H——多组分合金单相形成预测的新范式

**一、研究背景与问题提出**

高熵合金（HEA）及多组分合金（MCA）的设计核心挑战之一，在于准确预测给定成分体系将形成单相固溶体（SS）还是多相微观结构（Multi）。自Yeh和Cantor的开创性工作以来，这一问题始终占据着HEA研究的中心地位。由于成分空间极其庞大——仅涉及40种金属元素的等原子五元体系就有超过65万种组合——而实验成本高昂，因此建立简单可靠的预测准则具有毋庸置疑的实用价值。微观结构对合金的力学与功能行为具有决定性影响，这使得预测问题不仅具有学术意义，更直接关系到材料设计的效率与成败。

2008年，Zhang等人提出了三个经典描述符：原子半径错配参数δ、混合熵ΔS_mix和混合焓ΔH_mix，并分别给出了有利于单相微观结构的判据范围。2012年，Yang和Zhang引入无量纲参数Ω，将熵效应与焓效应统一为单一表达式，建立了Ω ≥ 1.1且δ ≤ 6.6%的联合判据，这一准则迄今仍是该领域的基准。然而，这些经典描述符在结构上存在一个根本性局限：它们均采用加权平均的方式处理所有元素对的二元混合焓，导致正负焓对在平均值中相互抵消。典型案例如CoCrCuFeNi体系，其整体ΔH_mix约-3 kJ/mol，看似利于互溶，但Cu-Fe对的二元混合焓高达+13 kJ/mol，Cu的偏析在该体系中系统性地发生。这种平均化处理掩盖了强吸引对与强排斥对共存的内部异质性，而后者恰恰是驱动相分离的关键因素。

**二、研究内容与核心发现**

针对上述局限，本研究提出了一种全新的无量纲描述符δ_H，定义为二元混合焓相对其加权均值的离散度，并以凝固温度下的平均热能RT_m进行归一化。δ_H量化了组成元素对之间化学相互作用的内部异质性——这是经典描述符无法捕捉的信息维度。研究表明，δ_H–δ二元组合在1675种合金的实验数据集上显著优于Yang-Zhang基准判据Ω–δ，交叉验证准确率提升+3.16个百分点（83.16%对80.00%），马修斯相关系数（MCC）提升+0.063（0.663对0.600），McNemar检验p = 0.006，差异具有统计学显著性。该发现发表于《ACTA MATERIALIA》。

**三、主要技术方法**

研究人员构建了包含1675种实验合金的综合数据集，整合了三个独立来源：Mendeley HEA Phase Dataset v1d（419种合金）、Kaggle High Entropy Alloy Properties数据集（608种合金）以及Pei等人的MOESM2_ESM补充资料（1135种合金）。经过去重和冲突标签清理后获得最终数据集，涵盖63种元素，成分范围从二元到九元。描述符计算采用统一的输入数据源：Takeuchi和Inoue的混合焓数据、Teatum等人的CN12原子半径、Brandes和Brook的熔点数据。评估采用5折分层交叉验证，通过网格搜索在每个训练折上优化阈值，并运用McNemar检验进行分类器间的统计比较。

**四、研究结果**

*4.1 个体判别能力*

Welch t检验显示，所有五个描述符在SS与Multi两组间的差异均高度显著（p < 0.001）。δ的t统计量绝对值最大（t = -25.6），确认其作为最佳单变量描述符的地位。δ_H达到t = -18.4，为第二大判别差异：SS合金的平均δ_H为0.43，而Multi合金为1.51，证实单相体系中的内部化学异质性显著更低。相比之下，VEC的组间差异最小（t = +4.5），统计显著但实际区分能力有限，这与VEC最初设计用于预测固溶体内部晶体结构（BCC/FCC）而非区分SS与Multi的初衷一致。

*4.2 交叉验证与比较分析*

5折交叉验证结果表明，δ_H–δ分类器的准确率达83.16% ± 1.75%，MCC为0.663，而参考判据Ω–δ为80.00% ± 1.18%，MCC为0.600。McNemar检验（配对1）确认δ_H–δ显著优于Ω–δ（b = 152，c = 205，χ2 = 7.57，p = 0.006）：在两类判据预测不一致的合金中，205种为δ_H–δ独有的正确预测，而反向独有正确预测为152种。在单变量描述符中，δ表现最佳（79.22%），非常接近Ω–δ二元参考判据（80.00%），表明Ω单独使用时对δ的增量贡献有限。δ_H单独准确率为69.97%，低于δ，但其与δ的组合产生了所有研究分类器中最大的性能提升（较δ单独提升+3.94个百分点），确认两者尽管存在中等程度相关，却捕捉了互补信息。涉及Λ的组合（Λ–δ和Λ–δ_H）均未超越δ_H–δ，且与参考判据Ω–δ接近，表明Λ的弹性修正对本分类任务无额外贡献。关于δ_H可能仅作为组分数n代理变量的担忧，按n分组的分析显示：对二元合金（n = 2），Ω–δ略优（88.4%对86.9%）；但对所有n ≥ 4的多组分体系，δ_H–δ均达到或超越Ω–δ，其中五元体系优势最大（+6.2个百分点），确认δ_H捕捉了超越组分数量的真实成分信息。全数据集上δ与δ_H的Spearman相关系数为ρ = 0.27（p < 0.0001），呈中等程度相关，部分归因于两者随成分复杂度共同增大的趋势。

*4.3 相图与象限分析*

在δ_H–δ相图中，四个象限展现出截然不同的行为特征。最具有决定性的象限是δ > 6.6%且δ_H ≤ 1.1，包含98%的Multi合金（N = 229）——当原子尺寸错配超过经典阈值时，几何不兼容性本身就足以排除单相形成。δ > 6.6%且δ_H > 1.1的象限产生89%的Multi（N = 79），纯度略低的原因在于高δ_H体系中强吸引对的存在可能局部补偿了不利因素。对δ_H贡献最具信息量的象限是δ ≤ 6.6%且δ_H > 1.1，包含82%的Multi（N = 528）——这些合金具有有利的几何错配，会被Ω–δ判据分类为SS，但实际主要为多相，内部化学异质性超过了可用热能（δ_H > 1.1），在几何错配适度的条件下驱动了相分离。该象限代表了δ_H–δ相对于Ω–δ改进的核心。最后，δ ≤ 6.6%且δ_H ≤ 1.1的象限包含79%的SS（N = 839），定义了可能形成单相的区域。在固定阈值下的完整数据集评估中，分类器准确率为83.0%，MCC为0.662，与交叉验证结果高度一致，证实δ_H = 1.1接近最优阈值。

*4.4 改进机制与物理诠释*

采用固定阈值（δ = 6.6%、δ_H = 1.1、Ω = 1.1）的分析表明，δ_H–δ正确分类了414种被Ω–δ误判为SS的Multi合金，同时损失了161种Ω–δ正确分类的合金，净增益为253种。对被挽回合金的成分家族分析揭示了两个主要类别：含Cu的Al-Co-Cr-Fe-Ni体系占总挽回Multi合金的36%，其中Cu-Ni对适度有利（ω_CuNi = -16 kJ/mol），但Cu-Cr和Cu-Fe对具有排斥性（分别为+12和+16 kJ/mol），产生了超过阈值的内部不对称性；Al-Nb-Ti-Zr及含Hf变体体系占6%，其中Al-Nb对强烈吸引（ω_AlNb = -150 kJ/mol），而Nb-Ti和Nb-Zr对几乎中性，尽管整体ΔH_mix为负且适中，仍产生了极高的内部化学离散度。其余58%的挽回案例分布在各种体系中，共同特征是至少存在一个强不对称二元对，使δ_H超过阈值，而全局ΔH_mix保持适中。

*4.5 局限性与未来方向*

数据集存在发表偏倚和合成方案偏倚两类已知局限性，但均等影响所有被比较的分类器，不改变相对比较的有效性。VEC对SS/Multi区分任务无额外判别力（单独准确率54.0%，接近随机），这与VEC的设计初衷一致，但不排除电子效应在单相形成中的更普遍作用。残留误差分析揭示了两种特征性错误模式：107个假阳性（6.4%）主要对应于富含Al且含Nb、Zr或Hf的高不对称成分合金，强吸引对抬高了δ_H，尽管其余元素高度兼容且体系实际形成固溶体；178个假阴性（10.6%）中53%为二元体系（δ_H恒为0），其余为所有二元交互均为均匀排斥的合金——此时围绕均值的离散度很小，尽管均值本身为正并指向不混溶。δ_H的输入数据来自基于Miedema模型的估算值，对稀土和锕系元素不确定性较大，但对数据集中常见元素验证良好。未来工作方向包括：与CALPHAD热力学计算耦合、探索δ_H对固溶体晶体结构的预测能力、以及在新合成合金上的前瞻性实验验证。

**五、讨论与结论**

讨论部分指出，δ_H与Ω并非严格独立——两者均源于相同的二元交互焓ω_ij——但它们捕捉该信息的不同方面：Ω关联总熵稳定化能T_mΔS_mix与平均混合焓量级|ΔH_mix|之比，而δ_H量化单个二元焓围绕其加权均值的离散度。研究结论明确指出：δ_H–δ组合的显著改进源于δ_H能够检测具有高内部化学异质性的多相体系——这类体系因有利对与不兼容对的贡献被平均化而被Ω–δ错误分类为SS。建议的操作阈值为δ ≤ 6.6%（与经典Yang-Zhang判据一致）和δ_H ≤ 1.1（与Ω ≥ 1.1阈值在数值上具有类比性），两者均具有直接的物理诠释：当内部化学异质性不超过凝固时可用热能RT_m的1.1倍时，固溶体是稳定的。两种物理机制——强负ω_ij对诱导的有序金属间化合物形成倾向和正ω_ij对诱导的偏析——正是δ_H能够检测而Ω在平均化过程中掩盖的效应。研究确认，在近等原子成分占主导的数据集中，焓异质性比Ω所嵌入的熵项贡献了更多的判别信息，但对于远离等原子比的成分，这一平衡可能发生转变。