实验测定三元相图——一个等温截面就要几十个合金、上百小时的热处理和表征。但热力学相图计算把这件事搬到了理论框架下：用少量实验数据拟合Gibbs自由能模型，然后通过自由能最小化外推到整个成分-温度空间。这个项目在分析Al-Mg-Si合金（6000系铝合金）的时效析出路径时，全程依赖CALPHAD方法做相图计算，过程中对几个核心原理的体会值得整理。
热力学相图计算的基本逻辑是：给定温度T、压强P和总成分xᵢ，系统平衡态是总Gibbs自由能G最小的那个状态。多相共存时，每相的Gibbs自由能都是成分和温度的函数——Gᵐ(T,x)——系统通过调整各相的成分和比例，找到G的最小值。
以Al-Mg-Si三元系为例：在500 K下，总成分为Al-0.8Mg-0.6Si（wt%）的合金，系统需要决策是形成α-Al单相，还是α-Al+Mg₂Si两相，还是α-Al+Mg₂Si+Si三相。这个决策靠的是对比三种相组合下的总自由能——相图计算软件（如Thermo-Calc、Pandat）做的就是这件事：在每个温度-成分格点上求解约束条件下的自由能最小化问题。
自由能最小化的前提是每相的Gᵐ(T,x)有解析表达式。对固溶体相（如α-Al），最简单的模型是理想溶液模型——原子完全随机混合，混合熵来自配置熵，混合焓为零。但实际合金中原子间相互作用不是零，需要加入超额项。
亚规则溶液模型（subregular solution model）是最常用的：它把Gᵐ分解为参考项、理想混合项和超额项。超额项用Redlich-Kister多项式展开到一阶或二阶：
G_ex = x_A x_B [L₀ + L₁(x_A−x_B)]
L₀和L₁是交互作用参数——来自实验数据的拟合。对Al-Mg-Si三元系，Al-Mg、Al-Si、Mg-Si三个二元系的L₀和L₁需要通过二元相图实验来评估。三元交互参数L_Al,Mg,Si如果二元描述不够，还要引入三元项——但这部分数据往往不充分，是CALPHAD数据库不确定性的主要来源之一。
回到项目实际需求：Al-0.8Mg-0.6Si合金在180℃时效时，Mg₂Si（β”相的前驱体）的析出驱动力有多大？用TCAL6数据库计算了该合金在200℃下的相平衡：平衡析出相为α-Al + β-Mg₂Si + 少量Si，Mg₂Si的摩尔分数驱动力（在形核阶段可用体积自由能变化近似）约-8.5 kJ/mol，对应临界形核半径约1.2 nm。
做了成分微调方案对比：Mg从0.8降到0.6 wt%，析出驱动力从-8.5 kJ/mol降到-5.1 kJ/mol，降幅40%。这个差异和实验观察到的时效硬化响应减弱是一致的——在热力学层面，时效强化潜力的天花板是Mg₂Si的析出驱动力决定的，工艺参数（时效温度和时间）决定的是能否到达这个天花板。
CALPHAD准不准，完全看数据库。TCAL6库对Al-Mg-Si三元系的评估基于2000年代中期的实验数据，对亚稳析出相（β’、β”等）的描述用的是稳定相模型的近似。这意味着你能准确算平衡相分数，但亚稳相的形核序列（GP zone → β” → β’ → β）不是纯热力学能回答的——需要结合形核理论和动力学模型。
另一个局限是晶粒尺寸效应。CALPHAD假设的是无限大块体相平衡，界面能可以忽略。但对纳米析出相（<5 nm），界面能对总自由能的贡献可能超过体积自由能项的5-10%，这会导致CALPHAD预测的析出相分数在低温下偏高。项目在析出相分数计算中加入Gibbs-Thomson修正（ΔG_GT=2γV_m/r），把纳米尺寸效应补了回去。
热力学相图计算的原理说起来就是”求自由能最小值”这一句话，但把这个原理从二元延伸到三元、从平衡相延伸到亚稳相、从块体延伸到纳米尺度，每一步都要知道模型的适用范围在哪里。
参考文献：Saunders & Miodownik, CALPHAD: Calculation of Phase Diagrams, Pergamon, 1998；Andersson et al., Calphad, 2002, 26, 273-312.
