第二章：吉布斯自由能模型
各位同学，今天我们来聊聊CALPHAD方法里最核心的东西——吉布斯自由能模型。说实话，我刚入行那会儿，觉得这些模型就是一堆数学公式，背下来就行了。后来做项目做多了才明白，选对模型比算对参数重要得多。
你想想看，我们做热力学计算，本质上就是在求解吉布斯自由能的最小值。但问题来了——不同体系的自由能长什么样？纯组元、溶液、化合物，它们的自由能表达式能一样吗？显然不能。所以这一章，我们就来拆解这些模型。
2.1 纯组元的吉布斯自由能
先说最简单的——纯组元。一个元素，没有混合，没有相互作用。它的吉布斯自由能怎么表达？
我习惯用SGTE（Scientific Group Thermodata Europe）数据库里的标准形式。说白了，就是一段分段多项式：
G(T) = a + bT + cT·ln(T) + Σ d_n·T^n这里面，a、b、c、d_n都是拟合参数。T是温度。为什么要有ln(T)项？因为热容C_p和温度的关系不是线性的，ln(T)能更好地描述低温行为。
举个例子，纯铝的吉布斯自由能，在298.15K到933.47K之间，表达式长这样：
G_Al(T) = -7976.15 + 137.093T - 24.3672T·ln(T) - 1.884662×10⁻³T² - 0.877664×10⁻⁶T³ + 74092T⁻¹嗯，看着挺吓人。但实际用的时候，这些参数都是现成的，你直接从数据库里调就行。我自己做项目时，最怕的就是手动输入这些参数——抄错一个数字，整个计算就废了。所以我建议，能用代码读数据库就别手敲。
关键点：纯组元的自由能是参考态。所有溶液模型、化合物模型，最终都要回归到纯组元的自由能上。没有这个基础，后面全是空中楼阁。
2.2 理想溶液模型
接下来，我们往体系里加第二种元素。最简单的混合情况是什么？就是两种原子谁也不搭理谁，完全随机排列。这就是理想溶液模型。
它的吉布斯自由能表达式：
G_mix = x_A·G_A + x_B·G_B + RT(x_A·ln x_A + x_B·ln x_B)前两项是纯组元贡献，最后一项是混合熵。注意，这里没有混合焓——因为理想溶液假设原子间相互作用和纯组元一样，混合时既不吸热也不放热。
我在项目中遇到过一种情况：用理想溶液模型去拟合液态合金，结果偏差特别大。后来发现，那套体系里A原子和B原子其实有很强的相互作用，根本不能用理想模型。所以，理想溶液模型只适用于：
- 同位素混合物（比如²³⁵U和²³⁸U）
- 化学性质非常相似的元素（比如Co和Ni在某些条件下）
- 稀溶液（溶质浓度极低时）
我的经验：如果你不确定能不能用理想模型，先算一下混合焓。如果实验测出来混合焓接近零，那可以试试。否则，老老实实用更复杂的模型。
2.3 正规溶液模型
理想溶液太理想了，现实世界哪有那么完美？原子之间总会有点"小情绪"。正规溶液模型就是在理想溶液的基础上，加了一个混合焓项。
表达式变成：
G_mix = x_A·G_A + x_B·G_B + RT(x_A·ln x_A + x_B·ln x_B) + x_A·x_B·L这个L就是相互作用参数。如果L是负的，说明A-B键比A-A和B-B键更稳定，体系倾向于形成化合物。如果L是正的，说明A和B互相嫌弃，容易发生相分离。
我记得有一次做钢中夹杂物的热力学计算，用正规溶液模型拟合CaO-Al₂O₃二元系，L参数调了好几天才收敛。后来发现，这个体系其实需要亚正规模型——因为相互作用参数不是常数，而是随成分变化的。
注意：正规溶液模型假设L与温度、成分无关。这个假设在很多体系中不成立。如果你发现拟合残差很大，先检查一下是不是L需要写成温度的函数。
2.4 亚正规溶液模型
亚正规溶液模型，说白了就是正规溶液的升级版。它允许相互作用参数随成分变化。最常见的做法是用Redlich-Kister多项式展开：
G_mix = x_A·G_A + x_B·G_B + RT(x_A·ln x_A + x_B·ln x_B) + x_A·x_B·Σ L_ν·(x_A - x_B)^ν这里ν从0开始取，ν=0就是正规溶液。ν=1、2、3...就是亚正规。每一项L_ν都可以是温度的函数。
我个人的习惯是：先试ν=0，如果不行再加ν=1。不要一上来就用高阶展开——参数越多，拟合越容易过拟合。你想想看，用5个参数去拟合10个实验数据点，结果看起来很美，但外推时可能完全跑偏。
| 模型 | 参数数量（二元系） | 适用场景 |
|---|---|---|
| 理想溶液 | 0 | 同位素、稀溶液 |
| 正规溶液 | 1 | 弱相互作用体系 |
| 亚正规（ν=1） | 2 | 不对称体系 |
| 亚正规（ν=2） | 3 | 复杂相互作用 |
避坑指南：我曾经在一个Fe-Cu二元系里用了ν=3的亚正规模型，拟合效果确实好。但后来做三元外推时，发现相图出现了不合理的"波浪形"。这就是过拟合的代价。所以，能用简单模型就别用复杂的。
2.5 化学计量相模型
最后说说化学计量相。这类相的特点是——成分固定，比如Al₂O₃、MgO、Fe₃C。你不能说Al₂O₃里铝多一点、氧少一点，它的成分就是严格2:3。
对于化学计量相，吉布斯自由能表达式很简单：
G_compound = G_pure_A + G_pure_B + ΔG_fΔG_f是生成自由能。它可以是常数，也可以是温度的函数。通常写成：
ΔG_f = a + bT这里a是生成焓，b是生成熵。注意，化学计量相没有混合熵项——因为成分固定，原子排列方式也固定，没有混合的"自由度"。
我做陶瓷材料热力学时，经常和化学计量相打交道。比如Al₂O₃，它的自由能数据非常成熟，直接拿SGTE数据库里的参数用就行。但有些复杂氧化物，比如Ca₃Al₂O₆，实验数据很少，这时候就需要用第一性原理计算来辅助拟合。
小技巧：如果你发现某个化学计量相的自由能参数拟合不出来，先检查一下它是不是真的化学计量相。有些相其实有固溶范围，只是实验没测出来。比如MgAl₂O₄尖晶石，名义上是化学计量，实际上Al和Mg可以互换位置，有一定的固溶度。
本章小结
好了，这一章的内容就这些。我们讲了五种模型：纯组元、理想溶液、正规溶液、亚正规溶液、化学计量相。它们的关系其实是一个递进的过程——从最简单到越来越复杂。
我个人觉得，学CALPHAD最重要的不是背公式，而是理解每个模型的假设和适用范围。你选对了模型，参数优化就是水到渠成的事。选错了模型，再好的优化算法也救不了你。
嗯，下一章我们会讲如何用这些模型去构建多元体系。到时候你会发现，二元系搞明白了，多元系其实就是"搭积木"。