第2章吉布斯自由能模型

第二章：吉布斯自由能模型——纯组元、理想溶液与正规溶液

各位好，欢迎来到第二章。上一章我们聊了CALPHAD的基本哲学，说白了就是把热力学数据塞进一个自洽的数学模型里。这一章我们得动点真格的了——聊聊吉布斯自由能到底怎么建模。

我个人习惯，讲模型之前先问一个问题：为什么是吉布斯自由能？ 答案其实很简单：在恒温恒压下，一个体系往哪个方向走，全看G值。G越小，体系越稳定。CALPHAD的核心工作，就是给每个相配一个靠谱的G(T, x)表达式。

核心思想： 整个CALPHAD数据库，本质上就是一堆吉布斯自由能模型的集合。每个相——不管是液相、固溶体还是金属间化合物——都得有一个数学表达式来描述它的G随温度、成分的变化。

先别急着搞混合体系。咱们得先把单个组元的行为摸清楚。纯组元的G只跟温度有关，不涉及成分变量。

我记得刚入行那会儿，导师跟我说：“纯组元的自由能是数据库的基石，基石不稳，后面全白搭。” 这话我到现在都记得。

最常用的表达式是SGTE（Scientific Group Thermodata Europe）推荐的格式：

G(T) = a + bT + cT·ln(T) + dT² + eT³ + f/T

这里面每一项都有物理意义：

你可能会问：这么多项，怎么确定系数？嗯，这就是CALPHAD评估的第一步——用纯组元的热容、相变焓、熵等实验数据去拟合这些系数。

避坑指南： 我曾经在评估一个高熔点金属时，忽略了高温段的电子比热容贡献，结果算出来的液相线偏差了50多度。后来加上f/T项，问题才解决。所以，高温数据一定要检查电子比热容项。

另外，纯组元通常有多个同素异构体。比如Fe有bcc、fcc、hcp三种结构。每种结构都得单独建一个G(T)表达式。它们之间的差值，就是相变驱动力。

好，纯组元搞定了。现在把两个组元混在一起，怎么描述混合相的G？

最简单的起点是理想溶液模型。它假设：

于是，理想溶液的吉布斯自由能写成：

G_ideal = ∑x_i·G_i⁰ + RT∑x_i·ln(x_i)

第一项是纯组元自由能的加权平均，第二项是理想混合熵的贡献。

说实话，真实体系里几乎没有严格意义上的理想溶液。但为什么还要学它？因为它是所有非理想模型的参考基准。你想想看，所有超额项（excess term）都是在这个基础上加的。

个人经验： 我在做Al-Cu二元体系时，一开始直接用理想溶液算液相线，结果跟实验差了200多度。这很正常——理想溶液只是起点，不是终点。但它给了你一个“零级近似”，让你知道后面需要加多少修正。

理想溶液太理想了。真实体系里，A原子和B原子之间是有“脾气”的——它们可能互相喜欢（负偏差），也可能互相嫌弃（正偏差）。

正规溶液模型就是在理想溶液的基础上，加了一个混合焓项。它假设：

其中Ω是相互作用参数。如果Ω < 0，表示A-B键比A-A和B-B键更强，体系倾向于混合（负偏差）。如果Ω > 0，则相反。

正规溶液的G表达式为：

G_reg = x_A·G_A⁰ + x_B·G_B⁰ + RT(x_A·ln(x_A) + x_B·ln(x_B)) + Ω·x_A·x_B

Ω可以是常数，也可以随温度变化：Ω = Ω₀ + Ω₁·T。我建议初学者先用常数Ω试试，拟合效果不好再考虑温度依赖。

注意： 正规溶液模型虽然简单，但它有一个致命缺陷——它假设混合熵是理想的。对于很多体系，特别是液态合金，这个假设并不成立。我曾经用正规溶液拟合一个液态体系，死活拟合不好，后来发现是短程有序导致的熵偏差。

正规溶液假设Ω是常数，但真实体系的相互作用往往随成分变化。于是就有了亚正规溶液模型（也叫Redlich-Kister模型）。

它的核心思想是：把相互作用参数Ω展开成成分的多项式：

Ω(x) = ∑ L_v · (x_A - x_B)^v,  v = 0, 1, 2, ...

其中L_v是Redlich-Kister系数。v=0对应正规溶液项，v=1及以上是非对称修正项。

完整的亚正规溶液模型写成：

G_sub = x_A·G_A⁰ + x_B·G_B⁰ + RT(x_A·ln(x_A) + x_B·ln(x_B)) 
        + x_A·x_B · ∑ L_v · (x_A - x_B)^v

这个模型是目前CALPHAD中最常用的固溶体模型。为什么？因为它足够灵活——你可以通过增加L_v的阶数来拟合任意复杂的实验数据。

实用建议： 我一般从v=0开始，也就是正规溶液。如果拟合残差太大，再加v=1项。最多用到v=3，再高就容易过拟合了。记住：模型越复杂，物理意义越模糊。

下面我用一张图来总结这几种模型的关系：

从这张图可以看得很清楚：从纯组元到亚正规溶液，每一步都是在增加模型的表达能力。但代价是什么？参数越来越多，物理意义越来越模糊。

我个人建议初学者这样选模型：

重要提醒： 不要为了拟合而拟合。我见过有人把亚正规溶液用到v=5，拟合R²确实漂亮，但外推到三元体系时完全崩溃。记住：CALPHAD的目标是预测，不是插值。

好了，这一章的内容就到这里。四种模型——纯组元、理想溶液、正规溶液、亚正规溶液——构成了CALPHAD数据库的骨架。下一章我们会聊更复杂的模型，比如缔合溶液和离子液体模型。但不管多复杂，底层逻辑都是今天讲的这些。