2. 热力学基础与相图计算

平衡凝固与Scheil计算的所有结果，最终都来自同一套热力学判据：体系在给定温度、压力和成分下，各相如何分配元素、各相的量是多少、成分如何随温度演化。要理解CALPHAD数据库为什么能同时支撑这两类计算，必须先回到吉布斯自由能、化学势和杠杆定律这三个基本工具，再看相图热力学优化如何把它们组织成可计算的数据。

吉布斯自由能与相稳定性

对于恒温恒压下的多元体系，判断某一相是否稳定、两相是否共存，依据是体系总吉布斯自由能是否达到最小。以二元A-B体系为例，某一相φ的摩尔吉布斯自由能可写成：

Gφ = x_A·G_A^φ + x_B·G_B^φ + RT(x_A ln x_A + x_B ln x_B) + G_ex^φ

其中前两项是纯组元参考项，中间是理想混合熵贡献，最后一项是过剩自由能，用来描述原子间相互作用偏离理想混合的程度。CALPHAD方法的核心，就是把G_ex^φ写成Redlich-Kister多项式等形式，用一组可调参数描述不同成分、不同温度下的相互作用。例如二元规则溶体近似下：

G_ex^φ = x_A x_B · Lφ
Lφ = ^0L + ^1L(x_A - x_B) + ^2L(x_A - x_B)^2 + …

这些L参数就是数据库中的热力学参数。它们不是从第一性原理直接算出的常数，而是通过拟合相图实验数据、活度数据、焓数据等优化得到的。参数一旦确定，任意温度、任意成分下的Gφ都可以计算，相平衡问题就转化为“在给定总成分下，寻找使总自由能最小的相组合与各相成分”。

化学势与相平衡条件

多相平衡的严格条件是：每个组元在各相中的化学势相等。对A-B二元系中α和β两相平衡，有：

μ_A^α = μ_A^β
μ_B^α = μ_B^β

化学势μ_i定义为在恒温恒压下，向体系中加入微量i组元所引起的自由能变化，即μ_i = ∂G/∂n_i。在摩尔自由能-成分图上，化学势对应自由能曲线在某成分处的切线截距。两相平衡时，两相自由能曲线存在一条公切线，切点成分就是平衡相成分。公切线斜率对应化学势差，截距对应各组元化学势。这个几何图像非常重要：它把“化学势相等”翻译成“公切线”，也把“相图上的相界线”翻译成“自由能曲线公切点的轨迹”。

对凝固计算而言，液相线和固相线就是液、固两相自由能公切点随温度移动形成的。温度降低时，公切点位置变化，液相成分和固相成分随之改变，这正是平衡凝固路径和Scheil路径中成分演化的热力学根源。

杠杆定律：相分数与成分分配

在两相区内，杠杆定律给出各相摩尔分数。设合金总成分为x_0，α相成分为x_α，β相成分为x_β，则：

f_α = (x_β - x_0) / (x_β - x_α)
f_β = (x_0 - x_α) / (x_β - x_α)
f_α + f_β = 1

杠杆定律本质是质量守恒，不依赖具体热力学模型。但它必须与相图或自由能计算配合：x_α和x_β由公切线条件确定，f_α和f_β再由杠杆定律求出。平衡凝固中，每一温度都允许固相和液相充分扩散、成分均匀，因此每一步都可用杠杆定律；Scheil计算中，固相无扩散、液相完全混合，固相成分沿固相线变化，液相成分沿液相线变化，但每一微小温度步内仍用杠杆定律分配新增固相与剩余液相的量。区别在于：平衡凝固是全局平衡，Scheil是局部界面平衡加固相无扩散。

相图热力学优化：从实验数据到可计算数据库

相图热力学优化（assessment）是CALPHAD方法的核心工作。它不是简单地把相图数字化，而是为每个相建立自由能模型，用实验数据反过来确定模型参数。典型流程如下：

1. 选择模型：液相通常用替代溶体模型；固溶体用亚点阵模型，例如FCC用(Fe,Cr,Ni)₁(C,N,Va)₁两个亚点阵；金属间化合物用有序-无序模型或化合物能量模型。
2. 收集数据：包括相界线、三相温度、活度、混合焓、热容、形成焓、第一性原理计算形成能等。数据按可信度加权。
3. 参数拟合：用最小二乘或贝叶斯优化，调整L参数、亚点阵相互作用参数，使计算相图与实验相图偏差最小。常用软件如Thermo-Calc、Pandat、FactSage。
4. 自洽性检查：检查所有相在全部成分温度范围内是否稳定、是否出现虚假相、是否满足热力学一致性（如Gibbs-Duhem关系）。
5. 数据库封装：把优化后的参数写成TDB文件，供平衡计算和Scheil模块调用。

以Al-Si二元系为例，液相和固相（FCC-Al、金刚石Si）的自由能参数经优化后，可同时复现共晶温度577°C、共晶成分12.6 at.% Si以及两端固溶度。用同一数据库做平衡凝固，得到的是平衡固相分数随温度的变化；做Scheil计算，则得到固相无扩散条件下的偏析曲线。两者共用同一套G函数，只是求解时施加的扩散假设不同。

数据库质量如何影响计算结果

CALPHAD数据库的误差会直接传递到凝固计算结果。常见问题包括：

• 相界线偏差：液相线温度偏差5–10°C，会导致Scheil终凝温度明显偏移，因为终凝温度对剩余液相成分非常敏感。
• 相互作用参数外推不当：在多元系中，二元参数向三元、四元外推时若缺少三元交互参数，可能高估或低估液相稳定区。
• 缺失相或错误相：若数据库遗漏了实际会析出的金属间化合物，Scheil计算会低估偏析或错误预测终凝相。
• 磁性贡献处理：Fe基合金中磁性转变对自由能贡献不可忽略，忽略G_mag会使相界线在居里温度附近出现明显偏差。

从热力学到凝固计算：一条完整的逻辑链

把上述内容串起来：CALPHAD数据库提供每个相的G(T, x)函数；平衡计算通过最小化总G或求解化学势相等，得到各温度下的平衡相组合与相成分；杠杆定律把相成分转换为相分数；平衡凝固假设全局平衡，Scheil计算假设固相无扩散、界面局部平衡。两类计算共享同一热力学基础，区别只在动力学假设。理解这一点，就能明白为什么同一个TDB文件既能算平衡凝固，也能算Scheil偏析，也能算后续的凝固路径、热裂敏感性和均匀化工艺窗口。

下一阶段进入具体计算时，需要重点关注：数据库版本、相的选择与抑制、计算温度步长、Scheil模块中固相扩散是否开启、以及终凝判据。这些设置会直接影响结果的可解释性。