3. 溶液模型与过剩自由能
在 CALPHAD 数据库里，溶液相（液相、固溶体、有序/无序固溶体等）的 Gibbs 自由能几乎都写成“理想混合项 + 过剩自由能项”的形式。理解这两部分各自承担什么、相互作用参数如何进入模型，是后续读写 TDB、拟合实验数据、判断外推合理性的基础。本节从理想溶液出发，逐层引入规则溶液、亚规则溶液与 Redlich-Kister 多项式，并说明相互作用参数的物理含义与工程用法。
一、理想溶液：只考虑构型熵
理想溶液假设不同组元在晶格或液体中随机混合，原子间相互作用完全相同，混合焓为零。对二元 A-B 体系，摩尔 Gibbs 自由能为：
G_m^ideal = x_A·G_A° + x_B·G_B° + R·T·(x_A·ln x_A + x_B·ln x_B)其中 G_A°、G_B° 是纯组元在相同参考态（同一相、同一温度、同一晶体结构）下的摩尔 Gibbs 能，通常来自 SGTE 纯组元数据库；R 为气体常数 8.314 J·(mol·K)⁻¹；x_A、x_B 为摩尔分数，且 x_A + x_B = 1。第二项即理想混合项 G_m^id = R·T·Σx_i·ln x_i，它恒为负值，是熵驱动的“混合倾向”。
理想溶液模型能给出对称的、以 x=0.5 为最低点的自由能曲线，但它无法描述实际体系中常见的偏析、短程有序、液相分层等现象。原因在于它忽略了原子间结合能的差异，也就是忽略了混合焓。
二、过剩自由能：把“非理想”集中表达
CALPHAD 的做法是把实际自由能写成理想部分加一个修正项，即过剩 Gibbs 自由能 G_m^E：
G_m = G_m^ref + G_m^id + G_m^E
G_m^ref = x_A·G_A° + x_B·G_B°
G_m^id = R·T·(x_A·ln x_A + x_B·ln x_B)G_m^E 承载了所有非理想效应：混合焓、非构型熵、原子尺寸差异、有序化倾向等。它通常用组元摩尔分数的多项式展开表示，展开系数就是“相互作用参数”。在 TDB 文件中，这些参数以 L 符号出现，例如 L(LIQUID,A,B;0)、L(LIQUID,A,B;1)。
三、规则溶液：一个相互作用参数
规则溶液（Regular Solution）是最简单的非理想模型，只保留一个与成分无关的相互作用参数 L₀：
G_m^E = x_A·x_B·L₀对应的混合焓 ΔH_mix = x_A·x_B·L₀，混合熵仍取理想值。参数 L₀ 的符号决定体系行为：
L₀ < 0：A-B 结合强于 A-A、B-B，混合放热，自由能曲线下凹，倾向于形成化合物或有序相；L₀ = 0：退化为理想溶液；L₀ > 0：同类原子结合更强，混合吸热，当L₀ > 2RT时自由能曲线出现双谷，发生液相分层（如 Al-Zn 液相在低温下的不混溶间隙）。
规则溶液的局限在于它强制 G_m^E 关于 x=0.5 对称，而真实体系（如 Fe-Ni、Cu-Zn）的过剩自由能往往不对称，需要更多参数。
四、亚规则溶液与 Redlich-Kister 多项式
亚规则溶液（Sub-regular Solution）允许相互作用参数随成分变化，最通用的表达是 Redlich-Kister 多项式：
G_m^E = x_A·x_B · Σ_{v=0}^{n} L_v · (x_A − x_B)^v
= x_A·x_B · [L₀ + L₁(x_A−x_B) + L₂(x_A−x_B)² + …]其中 L_v 是第 v 阶相互作用参数，通常写成温度的函数：
L_v = a_v + b_v·T + c_v·T·ln T + …在 TDB 中，L(LIQUID,A,B;0) 对应 L₀，L(LIQUID,A,B;1) 对应 L₁，依此类推。各阶参数的物理意义可以这样理解：
L₀：平均相互作用强度，决定整体偏离理想的程度；L₁：一阶不对称项，使过剩自由能曲线向某一侧倾斜，描述 A 端与 B 端行为差异；L₂：二阶项，用于描述曲线中部的“鼓包”或“凹陷”，常见于存在短程有序或团簇的体系；- 更高阶项：一般不超过 3～4 阶，阶数过高会导致外推到实验数据范围之外时出现非物理振荡。
需要强调：L_v 是拟合参数，不是第一性原理量。它们通过拟合相图边界、活度、混合焓、化学势等实验或计算结果获得，因此同一体系在不同评估版本中参数可能不同。
五、从参数到相图：一个可复算的示例
设某二元液相在 1000 K 下取 L₀ = −20000 J/mol、L₁ = −5000 J/mol，其余阶数为零。则：
G_m^E = x_A·x_B·[−20000 − 5000·(x_A − x_B)]
x_A = 0.4, x_B = 0.6
x_A·x_B = 0.24
x_A − x_B = −0.2
G_m^E = 0.24 × [−20000 − 5000×(−0.2)]
= 0.24 × [−20000 + 1000]
= 0.24 × (−19000) = −4560 J/mol负的 G_m^E 说明该液相在 1000 K 下对理想混合有额外稳定化作用，A-B 亲和力较强。若把 L₀ 改为 +20000 J/mol，则 G_m^E = +4560 J/mol，体系趋于排斥混合，低温下可能出现不混溶间隙。实际评估时，还要把 G_m^E 与 G_m^id 相加，再对 x 求导得到化学势，进而通过公切线构造计算相平衡。
L 参数的单位是 J/mol，温度项系数 b_v 的单位是 J/(mol·K)。若把 L₀ 的符号写反，液相线形状会整体翻转，这是数据库调试中最常见的低级错误之一。六、参数关系与选择建议
下图把理想溶液、规则溶液、亚规则溶液与 Redlich-Kister 多项式的关系、参数构成和典型自由能曲线形态放在同一张图中，便于对照理解。
选择模型时，可遵循以下经验：若体系接近理想（如某些高温液相），用规则溶液即可；若两端行为明显不对称，至少保留 L₀ 与 L₁；若存在短程有序、团簇或明显的中部曲率变化，再引入 L₂。每增加一阶参数，都需要有相应的实验信息（活度、混合焓、相边界）支撑，否则容易过拟合。
G_m^E 在全成分范围内的单调性与曲率合理性。掌握本节内容后，读 TDB 时看到 L(LIQUID,A,B;0)、L(FCC_A1,A,B;1) 就不再是陌生符号，而能直接对应到“平均相互作用强度”“不对称修正”等物理图像，也能判断某个参数的正负号是否与已知相图特征一致。这是后续进行参数拟合、数据库合并与热力学计算的前提。