二元平衡相图计算是CALPHAD体系的细胞。三元相图、多元相图、Scheil模拟——一切高阶特性都叠加在二元体系之上。二元体系的参数评估错了,三元体系必然错:这不是简单的1+1+1=3.是误差在Gibbs-Duhem关系下非线性传递,三元相图的某些区域会被扭曲到面目全非。

公切线法则是唯一的数学基础
二元平衡相图计算的数学底层很简单:在给定温度和压力下,两相平衡的条件是组元在两相中的化学势分别相等。翻译成几何语言——在吉布斯自由能-成分曲线上,两相平衡对应的成分由公切线(Common Tangent)的两个切点确定。
看起来简单。但这里的陷阱是:吉布斯自由能曲线不是实验室里直接测出来的。它是模型拟合出来的。
常用的是亚正规溶体模型(Subregular Solution Model),过剩自由能用Redlich-Kister多项式展开:
G^ex = x_A * x_B * [L0 + L1*(x_A – x_B) + L2*(x_A – x_B)^2 + …]
L0、L1、L2是交互作用参数,需要从实验数据(相边界、活度、形成焓)反推拟合出来。正则溶体模型只留L0一项——假设混合焓随成分对称变化。这个假设对Cu-Ni这种连续固溶体还说得通,对Fe-C、Al-Mg这种形成中间化合物的体系,L0一项根本拟合不出正确的相图拓扑。
项目里有过一个案例:Al-Cu体系的二元平衡相图计算,用正则模型(L0 only)算出的Al₂Cu(θ相)的析出温度比实验值偏了超过60°C。原因很简单:Al-Cu在液相中的短程有序(SRO,Short Range Ordering)导致了混合焓的成分不对称性,必须在Redlich-Kister展开式中保留L1甚至L2项才能描述。加上L1之后,θ相的析出温度回到实验范围内(偏差<10°C)。少一个参数,60°C的差距——这不是细节,是定性错误。
亚稳相的取舍是计算判断力的一道考题
二元平衡相图计算给出的”平衡”相图,实际上包含了一个默认假设:所有可能存在的相都能被自由能模型所描述,并且允许所有相参与平衡竞争。
但实际问题中,某些相因为动力学原因在实际冷却中从未出现,而数据库中仍然给定了它的吉布斯自由能参数。这时候,二元平衡相图计算会输出一个在热力学上成立、在动力学上不成立的”假平衡”——石墨在铁碳相图中的析出就是一个典型例子。含硅较低的铁碳体系,渗碳体(Fe₃C)从动力学上抢在石墨之前析出,平衡相图上写的是石墨+奥氏体,实际金相里看到的是渗碳体+奥氏体。
这两个都”正确”——一个热力学正确,一个动力学正确。
在Thermo-Calc/Pandat里,可以通过”Suspend Phase”功能将动力学上不会出现的相排除在平衡计算之外。操作上就是一个右键菜单,但背后的判断——哪个相该排除、哪个该保留——需要理解体系的凝固动力学。这个判断力,是二元平衡相图计算从”绘图”升级为”建模”的分界线。
固溶体相的亚晶格模型:原子占位决定相稳定性
有序固溶体(如B2、L1₂)的稳定性和相边界,在二元平衡相图计算中需要用亚晶格模型(Sublattice Model)来描述,而不是简单的正则溶体模型。
以Ni-Al系的Ni₃Al(γ’相,L1₂结构)为例。这个相的结构是:面心位置被Ni占据,角顶位置可以同时容纳Ni和Al。亚晶格模型将这个相拆分成两个亚点阵:(Ni,Al)₀.₇₅(Ni,Al)₀.₂₅——面心亚点阵和角顶亚点阵,每个亚点阵都有自己的成分变量。
这种描述不是理论洁癖。不拆分成亚晶格模型,γ’相的有序-无序转变(order-disorder transition)根本无法计算——因为有序化能来自Ni原子偏聚到面心位置的能量降低,而这个效应在”均匀混合”的正则溶体框架里是不可见的。
对于含有序相的二元体系,二元平衡相图计算中的模型选择不是”精度偏好”,是”能不能”的问题——选错模型,γ’相要么不出现,要么占据了整个浓度范围。
实验数据的权重分配
二元平衡相图计算的参数拟合中,输入数据不止一种:
– 相边界数据(XRD、DSC/DTA、扩散偶) – 热化学数据(形成焓、活度、比热) – 第一性原理数据(0K形成焓、特殊准随机结构SQS的总能)
拟合器(如Thermo-Calc的PARROT模块)会同时拟合以上所有数据,每个数据点有一个权重。权重的分配是潜在的手艺活儿:热化学数据通常误差较小(几kJ/mol),相边界数据的实验误差在不同文献间可能相差几十°C。如果给一份来源不明、误差偏大的相边界数据配了高权重,整个二元体系的交互作用参数会被一根不可靠的杠杆撬偏。
经验是:优先用Calphad期刊和Landolt-Börnstein数据库中的评估数据集,其次用原始实验文献——并且至少用两个独立来源的相边界数据互相校验。单篇文献的相边界,不管来自哪本期刊,先存疑,再采用。
边界之外
二元平衡相图计算产出的结果是一张相图和一组经过优化的热力学参数。但这组参数的可适用区间是有限的——通常在经过实验数据覆盖的温度-成分窗口内可靠,超出窗口就是外推。
参数外推的后果可能很隐蔽。比如Ni-Ti二元体系的参数在500~1500K范围内经过评估,但在300K下外推时,B2-NiTi马氏体相变的自由能差可能从实验值约-2kJ/mol变成+5kJ/mol——符号反了,相稳定性判定从”可能发生”变成”不可能发生”。
二元平衡相图计算的最后一步,应该是标注出这张相图的”可信区域”和”外推区域”,而不是把计算结果当作在全温度-成分空间都成立的固件输出。材料热力学的自尊,来自对自身边界的清醒。