高熵合金（HEA）的卖点是等原子比或近等原子比的多主元混合——五种或更多元素，没有哪一个元素明确占”溶剂”的位置。传统的CALPHAD方法建立在一元和二元基础之上，三元和四元通过外推得到，五元再加一层外推。外推层数越多，累积误差越大——到HEA这里，外推已经叠了至少三层。
HfNbTaTiZr是一个经典的难熔HEA体系。五个元素全是体心立方（BCC）结构，混合焓都是小幅正值或负值（Hf-Ta: +3 kJ/mol, Nb-Zr: +4 kJ/mol），预期应该是单相BCC固溶体。但铸态样品在SEM里看到了第二相——富Hf-Zr的析出物。CALPHAD能不能预测到这个第二相？
Thermo-Calc的高熵合金数据库（TCHEA）是目前对HEA覆盖最全的商业数据库。输入HfNbTaTiZr在1200°C的平衡相组成，CALPHAD给出：单相BCC，100%。没看到第二相。
问题在哪？Hf-Nb、Hf-Ta、Nb-Ta、Nb-Ti这些二元系的热力学描述（交互作用参数）都是完备的。但Hf-Ti、Ta-Ti、Zr-Ta这些二元的CALPHAD描述就稀疏了——有的二元系的液相参数是完备的，但BCC相的L12交互参数可能只是估值或从类似体系类推出的。三元系中，Hf-Ti-Zr、Ta-Ti-Zr的热力学描述几乎完全依赖Köhler外推——即假设三元交互作用可忽略，用三个边界二元的加权平均描述三元。
对一个完整三元系做Köhler外推的误差通常是可接受的（相边界偏差在50-100K以内），但如果四个三元系有描述，一个三元系纯外推——这第五个三元系在五元体系里的”权重”被稀释了，单看偏差不大。但这个被稀释的误差叠加在Hf-Zr的对等排斥效应上（Hf和Zr在BCC相中交互参数正且大），刚好在五元空间的Hf-Zr富集区域催生了混溶间隙。
混溶间隙是BCC固溶体内部的不互溶区——成分在这个区域内的BCC相会分解成两个BCC相（成分不同）。CALPHAD漏掉了这个混溶间隙，因为Hf-Ti、Hf-Zr、Ti-Zr这三对在二元层面都是完全互溶的——互溶间隙只出现在三元及以上的多组元空间（即混合焓正值的协同效应使互溶区在三元层面打开）。
追到根上，需要对Hf-Ti-Zr三元系的BCC相做不完备热力学描述的补充。最快的做法是用DFT算Hf-Ti-Zr三元空间几个关键成分点（Hf25Ti50Zr25, Hf50Ti25Zr25等，用SQS做化学无序模型）的形成焓，将这些DFT值作为三元L12交互参数的拟合目标。
在这个项目里，我们跑了12个成分点的SQS-DFT计算（每个64原子超胞），用VASP/PBE+PAW，ENMAX=350eV，k点4×4×4。12个点的形成焓全为正值（+2到+5 kJ/mol）——三元系的BCC形成焓为正，三元交互参数L12也取正值——这就意味着三元系确实存在混溶倾向。
把DFT拟合的三元参数回填到Thermo-Calc数据库中（作为用户自扩展数据库），重新跑HfNbTaTiZr五元相图，1200°C下BCC相区出现了一个狭窄的混溶间隙——富Hf-Zr侧析出一个BCC#2相（富Hf、Zr，贫Ta、Nb），基体BCC#1相富Ta、Nb。析出相的摩尔分数约8%，与SEM中第二相的面积分数（约7-10%）一致。
有DFT补充的三元参数后，跟踪从1400°C到800°C的平衡相分数演化。在1400°C以上是单相BCC；1350°C左右混溶间隙打开，BCC#2开始析出；1200°C时BCC#2的摩尔分数达到峰值（~8%）；降到900°C以下时，BCC#2开始回溶——因为温度降低后熵的稳定化效应减弱，混合焓正的组分倾向于形成有序化合物而非两个BCC相。
实验上，退火温度从1200°C降到900°C，SEM中第二相的面积分数确实从~8%减少到~3%。CALPHAD+DFT给出的趋势和实验一致。
高熵合金相图计算的核心经验：五元以上体系不能纯靠数据库的外推。当一个三元系的热力学描述不完整时，在多元空间里它可能刚好是那个打开混溶间隙的薄弱环节。DFT补充三元形成焓数据——12个成分点、几天的计算——是把这个薄弱环节补上的最直接手段。
