第三章 化合物相模型：从化学计量比到CEF
做CALPHAD这么多年，我最大的体会就是——相模型选对了，数据库就成功了一半。今天咱们聊聊化合物相模型的几个核心概念。说实话，刚入行那会儿我也被这些模型搞得头大，但后来发现，它们其实是一步步演化过来的，理解了这个脉络，一切就顺了。
3.1 化学计量比化合物：最朴素的模型
先说说最简单的。化学计量比化合物，说白了就是成分固定的化合物。比如Al₂O₃，铝和氧的比例就是2:3，变不了。你想想看，这种模型多简单——一个相就是一个点，在成分空间里动不了。
我记得刚做第一个数据库时，处理MgO就是这么干的。模型写成：
MgO: (Mg)1(O)1吉布斯自由能表达式也很直接：
G = G⁰ + ideal_term + excess_term这里G⁰是标准生成吉布斯自由能，通常从实验或第一性原理计算得到。ideal_term是理想混合熵——等等，化学计量比化合物没有混合熵，所以这项为零。excess_term也是零。
关键点：化学计量比化合物的自由度为零。成分固定，结构固定，能变的只有温度和压力。
实际项目中，纯物质、稳定化合物、中间相，很多都可以用这个模型。但问题来了——真实材料哪有这么完美？
3.2 线性化合物：成分可以滑动了
后来我发现，有些化合物成分不是完全固定的。比如某些碳化物，碳含量可以在一定范围内变化。这时候化学计量比模型就不够用了。
线性化合物模型应运而生。它允许成分沿着一条线变化，但限制在两端点之间。我习惯把它想象成一根绳子，两端固定，中间可以滑动。
举个例子，某二元体系中的化合物相，成分范围是A₀.₆B₀.₄到A₀.₄B₀.₆。模型写成：
(A,B)1(B,A)1吉布斯自由能表达式变成：
G = y_A * G_A + y_B * G_B + RT(y_A ln y_A + y_B ln y_B) + excess_term这里y_A和y_B是位点分数。注意，多了混合熵项！这就是线性化合物和化学计量比化合物的本质区别。
我的经验：线性化合物模型适合处理那些成分范围窄但确实存在的相。我曾经处理一个Ni基高温合金的数据库，γ'相就用这个模型，效果不错。但要注意，成分范围太宽就不行了，得用更复杂的模型。
3.3 Sublattice模型：把晶体结构拆开看
好了，接下来是重头戏。Sublattice模型，也叫亚点阵模型。这个想法其实很朴素——把晶体结构拆成几个子晶格，每个子晶格上可以放不同的原子。
你想想看，一个复杂的化合物，比如尖晶石结构，有四面体位和八面体位。不同原子占据不同位置，这怎么描述？Sublattice模型就是干这个的。
模型的一般形式：
(A,B,...)p(C,D,...)qp和q是位点比例。每个子晶格上的原子可以独立变化。
吉布斯自由能表达式：
G = Σ Σ y_i^s * y_j^t * G_ij + RT Σ (p^s * Σ y_i^s ln y_i^s) + excess_term看着复杂，其实逻辑很清楚：
- 第一项：端元化合物的自由能加权和
- 第二项：每个子晶格上的混合熵
- 第三项：相互作用参数，描述不同原子之间的"脾气"
注意：子晶格的数量和位点比例必须与晶体结构一致。我曾经见过有人为了拟合方便，随便加子晶格，结果预测出来的相稳定性完全不对。晶体学数据是基础，不能乱来。
3.4 CEF模型：CALPHAD的终极武器
CEF，全称是Compound Energy Formalism，化合物能量模型。说白了，它是Sublattice模型的升级版，也是目前CALPHAD领域最常用的模型。
CEF的核心思想：每个子晶格上的原子排列组合，形成一系列"端元化合物"。每个端元化合物都有其能量，整个相的自由能就是这些端元能量的加权和，再加上混合熵和过剩项。
我画个图帮你理解：
CEF模型的优势在哪？我总结三点：
- 灵活性高：可以描述任意复杂的晶体结构
- 物理意义明确：每个参数都有晶体学基础
- 可扩展性强：从二元到多元体系，框架不变
实际应用中，CEF模型最常见的写法：
// 一个两子晶格的CEF模型示例
// 子晶格1: (Al,Cr) 位点比例1
// 子晶格2: (C,N) 位点比例1
$ PARAMETER G(AL1CR1C1N1) 298.15 +30000; 6000 N !
$ PARAMETER G(AL1CR1C1V1) 298.15 +50000; 6000 N !
// ... 其他端元参数
// 相互作用参数
$ PARAMETER L(AL,CR:C,N) 298.15 -12000; 6000 N !
避坑指南：我曾经在构建一个多元数据库时，端元参数太多导致优化困难。后来学乖了——先用第一性原理计算端元能量，再拟合相互作用参数。这样既保证了物理合理性，又减少了拟合工作量。
3.5 四种模型的对比
说了这么多，咱们用个表格总结一下：
| 模型类型 | 成分自由度 | 混合熵 | 适用场景 | 复杂度 |
|---|---|---|---|---|
| 化学计量比化合物 | 0 | 无 | 纯物质、稳定化合物 | ★☆☆☆☆ |
| 线性化合物 | 1 | 有 | 窄成分范围相 | ★★☆☆☆ |
| Sublattice模型 | 多个 | 有 | 复杂晶体结构 | ★★★☆☆ |
| CEF模型 | 多个 | 有 | 通用，多组元体系 | ★★★★☆ |
嗯，这里要注意——不是越复杂的模型越好。我见过有人用CEF模型处理纯物质，这不是杀鸡用牛刀吗？选模型的原则是：够用就好，能简单就别复杂。
3.6 实际应用中的选择策略
说了这么多理论，到底怎么选？我个人的经验是这样的：
- 第一步：查晶体学数据，搞清楚结构有几个子晶格，位点比例是多少
- 第二步：看实验数据，成分范围有多宽？如果很窄，线性模型就够了
- 第三步：考虑组元数量，二元体系用Sublattice，多元体系直接上CEF
- 第四步：检查是否有有序-无序转变，有的话需要特殊处理
重要提醒：模型选择直接影响数据库的预测能力。我曾经接手一个项目，前人用了错误的子晶格比例，导致整个数据库的相图预测偏差很大。花了三个月才把问题找出来。所以，晶体学数据一定要核实，别偷懒。
好了，关于化合物相模型就聊到这儿。这些概念是CALPHAD数据库构建的基石，理解了它们，后面的内容就好办了。记住，模型是工具，不是目的。选对工具，事半功倍。
公众号:蓝海资料掘金营，微信deep3321