4. 亚点阵模型与有序相描述

在理想溶体模型里，我们假设所有原子随机混合、只用一个摩尔分数就能描述成分。但真实合金中大量存在有序相：不同元素倾向于占据晶格中不同的位置，例如面心立方（FCC）晶格中，Ni 倾向占据角位、Al 倾向占据面心位，形成 L1₂ 结构的 Ni₃Al。此时单一摩尔分数无法区分“Ni 在角位”还是“Ni 在面心位”，必须引入亚点阵模型（Sublattice Model）。它是 CALPHAD 描述有序相、金属间化合物、间隙相、离子熔体乃至陶瓷氧化物的核心工具。

一、亚点阵划分与位点分数

亚点阵模型的基本思想是：把一个晶体结构按晶体学位置等价性拆成若干亚点阵，每个亚点阵内部的位置对原子而言是等价的，不同亚点阵之间则不等价。例如 FCC 的 L1₂ 有序结构可拆成两个亚点阵：角位亚点阵（记为 A 位）和面心亚点阵（记为 B 位），其位置数之比为 1:3。用通式写作 (Ni,Al)₁(Ni,Al)₃，括号内是该亚点阵上允许存在的组元，下标是位置数。 描述成分不再用摩尔分数，而用位点分数（site fraction） y_i^s，表示组元 i 在亚点阵 s 上占据的位置比例，满足每个亚点阵内 Σy_i^s = 1。位点分数与摩尔分数 x 之间通过位置数换算，例如对 (Ni,Al)₁(Ni,Al)₃： x_Ni = (1·y_Ni^A + 3·y_Ni^B) / 4 当 y_Ni^A = y_Ni^B 时，两个亚点阵成分相同，即为无序态；当两者不等时，即为有序态。因此有序度可以用长程有序参数 η = y_Ni^A − y_Ni^B 来度量。这一“同一模型同时描述有序与无序”的能力，是亚点阵模型最优雅之处。

二、端元与化合物能量

每个亚点阵上各选一个组元，就构成一个端元（endmember）。对 (Ni,Al)₁(Ni,Al)₃，共有 2×2 = 4 个端元：Ni:Ni、Ni:Al、Al:Ni、Al:Al。每个端元对应一个化合物能量 G^end，即假设整个相就是该端元成分时的摩尔 Gibbs 能。这些能量通常以“化合物”形式给出，例如 G_Ni:Al 表示 Ni 占 A 位、Al 占 B 位的能量，它并不对应一个真实存在的稳定化合物，而是模型参数。 相的 Gibbs 能由三部分构成： G_m = G^ref + G^ideal + G^ex 其中参考项 G^ref = Σ_端元 (Π_s y_i^s) · G^end，即各端元能量按位点分数乘积加权；理想混合项 G^ideal = RT Σ_s a_s Σ_i y_i^s ln y_i^s，a_s 为亚点阵位置数；过剩项 G^ex 用 Redlich-Kister 型多项式描述同一亚点阵内组元间的相互作用，例如 L_Ni,Al:A 表示 A 位上 Ni-Al 相互作用参数。当相互作用参数足够负时，模型会自发产生有序化。

三、有序–无序转变与反位缺陷

有序–无序转变在亚点阵模型中表现为位点分数的连续变化。以 Ni-Al 的 L1₂/无序 FCC 为例，高温下 y_Ni^A ≈ y_Ni^B ≈ 0.75，两个亚点阵等价，相为无序 FCC；降温时 y_Ni^A 升高、y_Ni^B 降低，长程有序参数 η 从 0 连续增大，发生二级有序化。若把有序与无序写成同一亚点阵模型的两个不同解，则可在同一相内描述连续转变，避免人为设置两个相。 反位缺陷（antisite defect）指原子占据“错误”亚点阵位置，例如在化学计量 Ni₃Al 中，Al 占据角位（A 位）或 Ni 占据面心位（B 位）。在亚点阵模型中，反位缺陷直接由位点分数偏离理想值来体现：理想 Ni₃Al 为 y_Ni^A=0、y_Al^B=0；一旦 y_Ni^A > 0，就表示 Ni 反位到 A 位。反位缺陷浓度随温度升高而增大，是解释有序相成分偏离化学计量、影响扩散与力学性能的关键。

四、典型应用与建模要点

亚点阵模型在 CALPHAD 数据库中有几类典型应用。第一类是金属间化合物，如 Ni₃Al、TiAl、σ 相，用多亚点阵描述其成分范围与有序结构。第二类是间隙相，如钢中的渗碳体 Fe₃C 写作 (Fe)₃(C, Va)₁，其中 Va 表示空位，用于描述碳含量偏离化学计量。第三类是离子熔体与氧化物，如炉渣用 (Ca²⁺, Fe²⁺)(SiO₄⁴⁻, O²⁻) 等亚点阵描述，并需引入电荷平衡约束。第四类是有序–无序转变，如 Cu-Zn 的 B2 与 A2 转变、Ni-Al 的 L1₂ 与 A1 转变，用同一模型的不同位点分数解来描述。 建模时需注意几个要点。首先，亚点阵划分必须与晶体学结构一致，位置数之比要正确，否则端元能量失去物理意义。其次，端元能量通常由第一性原理计算或实验相图优化得到，不能随意设定。第三，相互作用参数 L 的符号和大小决定有序化倾向：L 越负，有序化温度越高。第四，若有序与无序用同一模型描述，需保证高温极限下位点分数趋于相等，否则会出现虚假的相分离。第五，反位缺陷的浓度可通过位点分数直接读出，是连接热力学与扩散、力学性能的桥梁。