热力学相图计算的原理，底层只有一句：在给定温度、压力和总成分下，体系会自发向吉布斯自由能最低的状态演化，平衡态就是自由能最小值对应的相组成和相成分分布。
这句话在本科物理化学教材里就已经写清楚了，但从原理到能实操的计算程序，中间隔了几层数学建模和算法实现。把这几个层次拆开来看，才能真正理解CALPHAD软件输出的是什么——以及它没有输出的是什么。
吉布斯自由能：你能测量的量构成了你能计算的边界
热力学相图计算的核心变量是吉布斯自由能G。对于一个单相，G是温度T、压力P和成分x的函数：
G = Σ x_i * G_i^0 + RT * Σ x_i * ln(x_i) + G^ex
三项分别对应：纯组元的参考态自由能(Reference)、理想混合熵(Ideal)、过剩自由能(Excess)。
第一项G_i^0来自纯组元的热化学数据——比热(Cp)、焓(H)、熵(S)对温度的函数。这些数据的主要来源是SGTE(Scientific Group Thermodata Europe)数据库，Larry Kaufman在1970年代开始系统收集和评估纯元素的热化学数据，至今仍在持续更新。SGTE纯元素数据库是热力学相图计算的”元素周期表”——它定义了每个元素在每种晶体结构下的吉布斯自由能基准线。
第二项理想混合熵是统计热力学的直接结果，不需要拟合参数。
第三项过剩自由能G^ex，是整个CALPHAD方法的”可调空间”——它吸收了所有非理想混合导致的吉布斯自由能偏差。G^ex本身不可直接测量，但它间接决定了相边界的位置和形状，而相边界是可以通过实验(XRD、DSC/DTA、扩散偶)测量的。所以，CALPHAD的本质是用相边界的实验数据反向约束G^ex中的参数。
热力学相图计算的原理能否在具体体系落地，取决于G^ex的模型选择是否匹配该相的物理特征。模型用错了，参数拟合可以收敛，但外推能力为零——因为你的模型在物理上就是错的。
三种相模型：正则溶体、亚晶格、化合物能量
热力学相图计算的原理在实现时，需要为不同类型的相配置不同的吉布斯自由能模型。常见的三种：
**正则/亚正规溶体模型**适用于置换固溶体(FCC、BCC、HCP、液相)。它假设所有原子在晶格点上随机分布，过剩自由能来源于不同原子对之间的交互作用能差异。参数少(L0、L1、L2)，外推能力强——这是CALPHAD能外推到多元体系的基础。代价是不能描述有序化现象，因为”随机分布”假设本身就是对有序化的否定。
**亚晶格模型**(Compound Energy Formalism，CEF)适用于有序固溶体和线性化合物。典型的如碳化物MC(M=Ti, Nb, V)：金属原子占据一个亚点阵，碳原子和空位共同占据间隙亚点阵。化学式写成(M)₁(C,Va)₁，空位Va的存在意味着这个相有一个有限的成分范围——碳亚点阵可以是全满(计量比碳化物)或部分空缺(亚计量比碳化物)。亚晶格模型比正则溶体模型多了一个维度——除了成分变量，还有每个亚点阵内部的原子/空位分布变量。
**化学计量化合物模型**适用于成分范围极窄的线化合物。比如Al₂O₃——它在相图上确实有线宽(line compound)，不严格是点化合物(point compound)，但在实际可用精度范围内，把它当作’成分固定’处理不会引入明显误差。这种模型最简单：一个固定成分对应一个固定的吉布斯自由能值，没有成分变量，没有混合熵，唯一的温度依赖来自纯组元的G_i^0(T)。
热力学相图计算的一个关键判断是：这个相该用哪种模型?而判断的依据来自于对晶体结构、缺陷化学和相宽度实验数据的理解。这不是建模技巧，是物理洞见。
吉布斯自由能最小化的算法实现
给定温度T、压力P、总成分x_bulk，热力学相图计算的数学问题是：
在满足质量守恒(各相成分加权平均等于总成分)的前提下，最小化体系的吉布斯总自由能：
min G_total = Σ f^φ * G^φ(x^φ)
其中f^φ是各相的摩尔分数，G^φ是各相的摩尔吉布斯自由能，x^φ是各相的成分向量。
这是一个带约束的多元非线性优化问题。早期算法用的是”两步法”：第一步枚举所有可能的相组合(单相、两相、三相……)，第二步对每个组合做吉布斯自由能最小化，最后选全局最小的那个组合。这个方法对二元、三元体系可以，对四元以上体系，相组合数量爆炸式增长，枚举法耗时过长。
现代CALPHAD软件(Thermo-Calc、Pandat)用的是全局优化算法——在成分-相分数空间中寻找全局自由能最小值。具体做法因软件而异：Thermo-Calc用的是基于Newton-Raphson的局部搜索配合多起点策略，Pandat用的是基于网格搜索(Grid Minimizer)的全局扫描法。两种方法对同一个输入数据库应该给出基本一致的结果——如果结果不一致，说明自由能曲面存在多个局部极小值，需要检查相模型的合理性。
有一个实际踩过的坑：含有序-无序转变的体系(如Ni-Al的γ/γ’)，自由能曲面上存在多个距离很近的极小值——一个对应有序态，一个对应无序态。优化算法有时会卡在无序态(更浅的极小值)，输出”没有γ’相”，而实际物理上γ’在低温区稳定存在。这时候需要人工干预——把γ’相的初始分数设为一个非零值，让优化器从有序态附近出发。这是热力学相图计算的原理和算法实现之间的一个缝隙：物理上存在的解，算法可能找不到。
从原理到实践的桥：数据库评估
热力学相图计算的原理清晰，实现稳定，但当计算结果与实验数据出现系统性偏差时，问题几乎总是出在数据库的模型参数评估质量上——而不是原理本身。
一个好的数据库评估(Assessment)长什么样?至少满足三条：
第一条，二元参数同时拟合了相边界数据(来自扩散偶、XRD、DSC)和热化学数据(形成焓、活度)。只拟合相边界、不顾热化学数据，会导致计算出的形成焓偏离量热实验值——外推三元时必定不准。反过来，只拟合热化学数据、不检查相边界，可能给出形状正确的自由能曲线但在关键温度窗口内相边界偏了十几度。
第二条，三元参数经过了至少两组独立实验数据的验证——比如扩散偶的相边界+等温截面的XRD相鉴定。单组数据拟合三元参数容易过拟合——参数多了，一条实验曲线可以拟合出多组不同的参数组合，但外推时差别显著。
第三条，相模型的物理基础是正确的。比如，不应该用正则溶体模型去描述有超晶格有序化的B2相(即便拟合能通过)，因为正则溶体模型在这种成分区域的函数形式就不对——它的结果在拟合区间内可以很准，出了拟合区间完全失真。
热力学相图计算的原理提供了方向，数据库评估提供了精度。方向对了，精度不够，出来的结果依然是误导。
计算结果的”不确定性标注”
热力学相图计算的一个常见错误是把计算结果当作精确值。相边界温度±10°C、相分数偏差±2%——这些是CALPHAD数据库评估质量的”典型公差”，不是软件菜单里的参数。
不同体系的评估质量差异悬殊：Fe-Cr-Ni体系经过40年以上的迭代评估，相边界预测偏差通常在±5°C以内，可靠性极高。Co-Cr-Fe-Mn-Ni(Cantor合金)体系，三元参数覆盖率不足，相边界偏差在±30°C甚至更高。
计算结果的输出应该像实验数据一样带上”误差估计”——至少标注出温度和成分窗口内哪些区域的预测经过了充分的实验验证、哪些区域靠外推。把这个意识贯彻到每一张输出的相图中，才是热力学相图计算的原理在工程实践中的完整落地。
理论教你怎么算。数据库评估告诉你算多准。两者之间的差距，就是实际项目中需要花时间判断的灰色地带。
