📚 课程目录第2章热力学基础与数据库关联

2. 热力学基础与数据库关联

CALPHAD 方法的全部计算能力，最终都建立在同一个热力学框架之上：体系在给定温度、压力和成分下，通过最小化总吉布斯自由能来确定平衡状态。数据库则是这个框架的“数值载体”——它把每一种相的吉布斯自由能写成温度、压力和成分的函数，并把这些函数的系数以模型参数的形式存储下来。理解数据库，本质上就是理解“自由能函数如何被参数化”，以及“这些参数在物理上意味着什么”。

一、吉布斯自由能：相稳定性的统一判据

对恒温恒压下的封闭体系，热力学第二定律给出的平衡判据是总吉布斯自由能最小。对含 c 个组元、φ 个相的多元多相体系，总自由能可写为各相自由能的加权和：

G_总 = Σ_φ n^φ · G_m^φ(x^φ, T, P)

其中 n^φ 为第 φ 相的摩尔数，G_m^φ 为该相的摩尔吉布斯自由能。CALPHAD 计算的核心任务，就是在满足各相成分归一化约束和全局质量守恒约束的前提下，寻找使 G_总最小的相组合与各相成分。相图、相分数、相成分、活度、驱动力等所有输出，都是这一最小化问题的副产品。

对单个相，摩尔吉布斯自由能通常分解为三部分：

G_m = G_ref + G_ideal + G_excess
G_ref    = Σ_i x_i · °G_i        （纯组元参考项）
G_ideal  = R·T · Σ_i x_i · ln(x_i) （理想混合熵项）
G_excess = Σ_i Σ_j>i x_i·x_j·L_ij(x) （过剩项，含相互作用参数）

参考项 °G_i 描述纯组元 i 在该相结构下的自由能，通常以 SGTE 推荐的晶格稳定性多项式给出；理想项永远使自由能下降，是熵驱动的混合倾向；过剩项则刻画实际体系偏离理想混合的程度，是数据库参数最集中、也最需要谨慎处理的部分。

二、化学势与相平衡条件

化学势 μ_i 定义为在恒温恒压、其他组元不变时，向体系中加入 1 摩尔组元 i 所引起的自由能变化。对理想溶液，μ_i = °G_i + RT ln x_i；对实际溶液，则用活度 a_i 替代浓度：

μ_i = °G_i + RT ln a_i，a_i = γ_i · x_i

其中 γ_i 为活度系数，直接由过剩自由能对成分的偏导数给出。多元多相平衡的严格条件是：每一组元在所有共存相中的化学势相等，即

μ_i^α = μ_i^β = μ_i^γ = … （i = 1, 2, …, c）

这一组等式与质量守恒方程共同构成非线性方程组，是 Thermo-Calc、Pandat、OpenCalphad 等软件求解平衡的数学基础。对二元系，它退化为“公切线法则”：两相平衡时，两相自由能曲线被同一条公切线相切，切点横坐标即平衡成分，切线斜率即公共化学势。理解公切线，是理解多元平衡的几何入口。

三、多元多相体系的热力学描述

当组元数增加到三元、四元甚至更高时，自由能曲面无法再用平面曲线直观展示，但数学结构不变：每个相对应一个 c 维成分空间上的自由能超曲面，平衡态由所有共存相超曲面的“公共切超平面”确定。此时需要求解的未知量包括各相摩尔分数、各相成分，方程数随组元数和相数迅速增长，因此实际计算依赖数值优化而非解析求解。

多元体系的一个关键难点是相的选择：候选相可能有几十个，但平衡时通常只有少数几个稳定存在。CALPHAD 软件采用“全局最小化 + 相剔除”策略，先假设所有相都可能存在，再通过迭代把不稳定相的摩尔数压到零。判断某相是否稳定，可用其驱动力（driving force）——即该相在当前化学势下相对于平衡态的自由能盈余。驱动力为负说明该相可以稳定析出，为零说明处于临界。

四、数据库中的模型参数与物理意义

一个典型的 TDB（Thermo-Calc Database）文件中，参数以关键字和多项式形式书写。例如一个简单置换固溶体的相互作用参数可能写成：

PARAMETER L(LIQUID,A,B;0) 298.15 -25000 +8.5*T; 6000 N !
PARAMETER L(LIQUID,A,B;1) 298.15 -12000; 6000 N !
PARAMETER L(FCC_A1,A,B;0) 298.15 -18000 +6.0*T; 6000 N !

这里 L(...,;0) 是零阶相互作用参数，;1 是一阶参数，分别对应过剩项中 x_Ax_B·L₀ 与 x_Ax_B(x_A−x_B)·L₁ 的形式。参数值中的常数项反映焓贡献，乘以 T 的项反映熵贡献。负的 L₀ 通常意味着 A–B 之间吸引作用较强，倾向于形成有序或化合物；正的 L₀ 则意味着排斥，可能引起液相分层或溶解度间隙。

需要特别注意的是，数据库参数不是纯物理常数，而是拟合参数。它们由实验相图、热化学数据（如混合焓、活度）、第一性原理计算等多源信息联合优化得到。因此同一元素在不同数据库中的 °G_i 可能略有差异，跨数据库计算时必须保证参考态一致，否则会出现虚假的相稳定性偏移。

实用提示：拿到一个数据库后，先检查三件事——纯组元参考态是否与你的体系匹配（如 FCC_A1 还是 BCC_A2）、相互作用参数的阶数是否足够描述你的成分范围、以及是否有磁性或有序化模型被激活。这三项直接决定计算结果的可靠性。

五、从参数到相图：一个可复现的检查流程

假设我们要验证某 Al–Cu 数据库在 500 °C 附近的相平衡是否合理，可以按以下流程操作：

在 Thermo-Calc 中定义体系：GO DICTRA 或 GO TC，设置数据库 SW ALCO，定义组元 Al、Cu。
设置条件：S-C T=773.15 P=101325 N=1，成分 X(CU)=0.02。
执行平衡计算 COMPUTE-EQUILIBRIUM，观察输出相组合与相分数。
检查结果：若出现 α-Al + θ-Al₂Cu 两相区，且 θ 相分数随 Cu 含量单调增加，说明数据库在该区域行为合理。
若出现异常相（如液相在 500 °C 稳定），则需检查该相的自由能参数是否被错误外推，或参考态是否不一致。

这个流程的价值在于：它把“数据库参数—自由能函数—相平衡结果”这条链路变成可观察、可验证的操作。参数是否正确，最终要通过相分数、相成分、转变温度等可观测输出反推判断，而不是只看参数数值本身。

常见误区：把数据库参数当作绝对真理。实际上，参数具有适用范围（温度区间、成分区间、压力范围），超出范围外推时误差可能急剧增大。工程计算中应始终记录所用数据库版本、参数来源和适用区间，以便追溯和复现。

总结来看，吉布斯自由能是相稳定性的统一判据，化学势相等是多元多相平衡的数学条件，而数据库参数则是把这两者连接起来的数值桥梁。掌握“参数—自由能—相平衡”的因果链，才能在后续的相图计算、凝固模拟和扩散模拟中，正确判断结果是否可信。