3. CALPHAD热力学建模基础

CALPHAD方法的核心，是用一套可计算的热力学模型描述多组元体系中各相的自由能随温度、压力和成分的变化，再通过相平衡计算得到相图与热力学性质。理解这一方法，关键要抓住四个层次：单个相的Gibbs自由能如何表达、溶液相如何描述混合行为、有序相与化合物如何用亚点阵模型处理、以及这些模型参数如何被评估并组织进数据库。本节围绕这四个层次展开，并给出一个可操作的参数评估与数据库组织示例。

一、Gibbs自由能模型的基本结构

对任一相φ，其摩尔Gibbs自由能通常写成三部分之和：

Gφ = Gφ_ref + Gφ_id + Gφ_ex

其中Gφ_ref是纯组元或端元相的参考自由能，Gφ_id是理想混合熵贡献，Gφ_ex是过剩自由能，用来描述实际溶液对理想行为的偏离。参考项一般写成各端元自由能的加权和，权重为端元在相中的摩尔分数；理想混合项为RTΣx_i ln x_i；过剩项则用Redlich-Kister多项式或其它展开形式表达。以二元A-B溶液为例，过剩自由能常写为：

G_ex = x_A x_B Σ_k L_k (x_A - x_B)^k

其中L_k是第k阶相互作用参数，通常设为温度的函数L_k = a_k + b_k T + c_k T ln T。这一形式的意义在于：L_0描述规则溶液行为，L_1、L_2等高阶项用于拟合非对称的混合焓与活度曲线。参数数量并非越多越好，需与实验数据量和物理合理性匹配，否则容易过拟合。

二、溶液模型：从理想溶液到亚正规溶液

溶液模型的选择取决于体系偏离理想行为的程度。理想溶液模型只保留参考项和理想混合项，适用于原子间相互作用近似相等的体系；规则溶液模型加入L_0，可描述对称的混合焓；亚正规溶液模型进一步引入L_1，能描述非对称的活度系数和偏摩尔性质。对于强相互作用体系，如含强化合物形成倾向的液态合金，还需考虑缔合模型或离子液体模型。

判断模型是否合适，不能只看拟合残差，还要检查外推行为。例如在二元体系评估得到的参数，外推到三元或更高组元时，若出现不合理的液相分层或负的熵，说明模型形式或参数约束需要调整。实际评估中常采用“先二元、后三元、再多元”的策略，并用实验相图、活度、焓和化学势数据交叉验证。

三、化合物能量模型与亚点阵描述

对于金属间化合物、有序相和陶瓷相，简单的置换溶液模型无法描述其晶体结构中的占位偏好。化合物能量模型（Compound Energy Formalism, CEF）用亚点阵描述这类相：每个亚点阵有固定的格点数，组元在各自亚点阵内混合。以B2结构的有序相为例，可写成(Fe,Ni):(Al,Fe)两亚点阵模型，其Gibbs自由能由端元化合物能量、理想混合熵和相互作用项组成：

G = Σ y_i^s y_j^t G_ij + RT Σ_s a_s Σ_i y_i^s ln y_i^s + G_ex

其中y_i^s是组元i在亚点阵s中的位点分数，a_s是亚点阵格点数，G_ij是端元化合物i:j的自由能。CEF的优点是能自然描述有序-无序转变、反位缺陷和成分范围，但参数数量随亚点阵数和组元数迅速增加，评估时必须依赖晶体结构信息和有序度实验数据。

四、参数评估与数据库结构

参数评估的本质是优化问题：以实验数据为目标，调整模型参数使计算值与实验值的加权残差最小。常用工具包括Thermo-Calc的PARROT、Pandat的PanOptimizer等。评估流程通常包括：收集并筛选实验数据、确定模型与初始参数、设置权重、迭代优化、检查残差与物理合理性。权重设置需反映数据精度，例如相边界温度权重可设为±2 K，活度数据权重可设为±5%。

评估完成后，参数需按数据库格式组织。以TDB（Thermo-Calc Database）为例，一个相的定义包括相名、亚点阵模型、端元自由能、相互作用参数和温度范围。下面给出一个简化的TDB片段示例，用于说明结构：

PHASE LIQUID % 1 1.0 !
CONSTITUENT LIQUID : A,B : !
PARAMETER G(LIQUID,A;0) 298.15 +5000 -2.5*T ; 6000 N !
PARAMETER G(LIQUID,B;0) 298.15 +3000 -1.8*T ; 6000 N !
PARAMETER L(LIQUID,A,B;0) 298.15 -12000 +4.0*T ; 6000 N !
PARAMETER L(LIQUID,A,B;1) 298.15 -3000 +1.0*T ; 6000 N !

该片段中，PHASE定义液相为单亚点阵置换溶液，CONSTITUENT列出组元A和B，PARAMETER给出端元自由能和相互作用参数，末尾的298.15–6000 K为有效温度范围。数据库的层次结构通常为：数据库→相→亚点阵→端元/相互作用参数。评估得到的参数必须与数据库中的相模型、组元定义和温度范围一致，否则计算时会出现参数缺失或外推异常。

五、评估中的常见问题与判断依据

参数评估不是纯数学拟合，必须结合物理约束。常见问题包括：参数外推后出现液相分层、端元自由能交叉导致负熵、亚点阵模型与晶体结构不符、以及权重设置不当使某一类数据主导优化。判断依据可以归纳为三条：第一，计算相图与实验相边界在误差范围内一致；第二，活度、焓、化学势等热力学量趋势合理；第三，参数在温度外推时保持连续且不产生非物理的相分离。若出现异常，应优先检查模型形式是否匹配晶体结构，再调整参数约束或权重。

六、数据库组织与使用要点

一个可用的CALPHAD数据库，除了参数本身，还需包含相模型定义、组元列表、温度范围、参考态和版本信息。使用数据库时，应先确认待计算体系的组元、相和温度区间是否被数据库覆盖；若需外推，应检查外推路径上是否存在参数缺失或相模型不兼容。对于工程应用，建议在数据库基础上建立验证用例，例如用已知相图关键点、活度数据或焓数据做回归测试，确保数据库更新后计算结果仍可复现。