3. 相图计算的热力学模型与自由能函数

在多元合金相图计算（CALPHAD）中，相图本质上是各相吉布斯自由能随成分、温度变化关系的几何投影。要计算一个多元体系的相平衡，第一步就是为每一个可能出现的相建立自由能表达式。本节围绕理想溶液、正规溶液、亚正规溶液、Redlich-Kister多项式以及化合物能量模型（CEF）展开，重点说明这些模型的数学形式、参数物理意义，以及如何从参数值判断体系中的相互作用性质。

一、自由能函数的统一框架

对于任意一个相 φ，其摩尔吉布斯自由能通常写成三部分之和：

G_m^φ = G_ref^φ + G_ideal^φ + G_excess^φ

其中 G_ref 是纯组元机械混合项，由各纯组元自由能按摩尔分数加权得到；G_ideal 是理想混合熵贡献，形式为 RTΣx_i ln x_i，恒为负值，代表混合使自由能下降；G_excess 是过剩自由能，承载了所有偏离理想行为的相互作用信息。相图计算的核心工作，就是确定 G_excess 的函数形式并拟合其参数。

需要特别强调符号约定：本文所有自由能均以纯组元在相同温度、相同晶体结构下的自由能为参考态，过剩自由能的正负号表示对理想混合的偏离方向——正值表示同类原子偏好聚集（趋向分层），负值表示异类原子偏好结合（趋向有序或化合物形成）。

二、从理想溶液到正规溶液

理想溶液模型假设异类原子与同类原子间的键能完全相同，因此 G_excess = 0。二元理想溶液的 G_m 曲线是一条关于成分对称的、开口向上的光滑曲线，任意温度下只有一个极小值，对应完全互溶。理想溶液模型在真实合金中极少严格成立，但它是判断其他模型偏离程度的基准线。

正规溶液模型在理想溶液基础上引入一个与成分无关的相互作用参数 L：

G_excess = x_A x_B L

L 的单位是 J/mol，其物理意义是 1 mol 原子由纯 A 与纯 B 混合形成固溶体时的焓变。当 L > 0 时，混合焓为正，自由能曲线在低温下可能出现双谷，对应溶解度间隙与两相分离；当 L < 0 时，异类原子吸引，曲线更趋单谷，有利于形成连续固溶体或有序相。正规溶液模型只用一个参数就能描述对称的偏离行为，是理解相互作用参数含义的最佳入口。

三、亚正规溶液与Redlich-Kister多项式

真实体系的相互作用往往随成分变化，正规溶液的常数 L 不够用。亚正规溶液模型允许 L 随成分线性变化，即 L = L0 + L1(x_A − x_B)。更一般地，Redlich-Kister多项式将过剩自由能展开为：

G_excess = x_A x_B Σ_v L_v (x_A − x_B)^v

其中 v 为展开阶数，L_v 为第 v 阶相互作用参数。v = 0 时退化为正规溶液；v = 1 时即亚正规溶液；v 取到 2 或 3 可描述非对称的、甚至带拐点的自由能曲线。Redlich-Kister 形式之所以被 CALPHAD 广泛采用，是因为它在 x_A = 0 和 x_A = 1 两端自动为零，不破坏纯组元参考态，且各阶参数之间近似正交，便于逐阶拟合。

参数判读要点：L0 主导整体偏离方向与幅度；L1 控制不对称性，其符号决定曲线极小值偏向哪一侧；L2 及以上用于修正中段曲率。若拟合后发现 L1 的绝对值接近甚至超过 L0，说明体系非对称性很强，应检查是否遗漏了中间化合物或有序相。

四、化合物能量模型（CEF）

对于金属间化合物、有序相、固溶体中的亚点阵结构，需要采用化合物能量模型。CEF 把相划分为若干亚点阵，每个亚点阵可被一种或多种组元占据，自由能写成：

G_m = Σ_I Π_I y_i^I G_I^end + RT Σ_I a_I Σ_i y_i^I ln y_i^I + G_excess

第一项是端元化合物自由能加权和，Π_I 为各亚点阵位点分数乘积，G_I^end 是端元化合物能量；第二项是每个亚点阵内部的理想混合熵；第三项为亚点阵间的相互作用过剩项，同样可用 Redlich-Kister 形式展开。CEF 的关键参数包括亚点阵位点比、端元能量和相互作用参数，它们共同决定有序—无序转变温度、溶解度范围和化合物稳定性。

五、参数拟合与判读的实际操作

在实际计算中，建立自由能函数通常按以下顺序推进。首先确定相的结构与亚点阵划分，例如 FCC 相用单亚点阵 (A,B)，B2 有序相用双亚点阵 (A,B):(A,B)。然后收集该相的实验或第一性原理数据：混合焓、活度、相边界、有序化温度等。接着用最小二乘或贝叶斯优化拟合 L_v 参数，目标函数同时包含相平衡条件与热力学数据残差。

拟合完成后，判读参数是否合理可从三方面入手。其一，检查 L0 的符号与量级是否与已知混合焓一致，例如 Ni-Al 体系中 NiAl 形成焓约为 −60 kJ/mol 量级，若拟合得到正值则明显异常。其二，观察 L1/L0 的比值，若超过 1 应警惕模型阶数不足或数据覆盖成分范围过窄。其三，用拟合参数重新计算相图，与实验相边界对比，重点检查溶解度间隙、共晶点和有序转变温度是否落在合理区间。

六、模型选择建议

选择模型时，应遵循“从简到繁、按需升阶”的原则。若正规溶液已能复现主要相边界，就不必引入高阶项；若体系中存在明显有序化或化合物，则必须使用 CEF 并合理设置亚点阵。模型复杂度应与可用数据量匹配，避免参数过多导致外推行为失控。最终，自由能函数的建立不是终点，而是后续相平衡计算、相图外推和工艺设计的基础。