热力学计算：CALPHAD相图外推与DFT+热力学在Al-Cu-Mg合金中的真实推演

接手这个项目时，合作单位已经做完了三轮Al-Cu-Mg合金热处理实验。第三轮的DSC曲线在275 ℃附近出现一个意料外的吸热峰——既不是文献里报过的θ’（Al₂Cu）相，也不是S’（Al₂CuMg）相。问题被甩到计算端：要回答”这个未知峰是什么相析出导致”。这正是热力学计算最擅长、也最容易翻车的领域——把实验观察翻译成相图的语言。

一、CALPHAD方法的起点：相图外推的边界条件

热力学计算的第一件事不是DFT，而是建一个稳健的CALPHAD（CALculation of PHAse Diagrams）数据库。Al-Cu-Mg三元系并非全新体系——TCGB、SGTE、PanAl2024等商业库都包含这一三元系的基础描述。但”基础描述”和”能预测未知相”之间有一段距离。

项目组从TCGB2023版库导入Al-Cu-Mg三元参数，使用Thermo-Calc 2023b的TCHEA2模块，设定外推压力为标准大气压，温度扫描300 K到900 K。在300到500 K区间，Al-3.5Cu-1.2Mg（质量百分比，下同）合金的平衡相图预测存在α-Al基体 + θ（Al₂Cu） + S（Al₂CuMg）三相共存。认定这一平衡相组态更可信的判断依据是：它与Monroe、Hornbogen等人在1980年代发表的多组实验相图一致。

二、亚稳相的DFS能垒与动力学约束

CALPHAD给的是热力学平衡相图，但它不能直接回答”为什么275 ℃的DSC峰不在平衡相图上”。这一类析出峰往往对应亚稳相——其吉布斯自由能比平衡相略高（通常高0.5-3 kJ/mol），但在特定温度和成分范围内因形核动力学优势而优先析出。

热力学计算在这里需要做两件事：先用CALPHAD算出各相在300 K下的形成自由能ΔGf，再与DFT点计算结果做交叉验证。S相的ΔGf = -22.4 kJ/mol（CALPHAD），DFT用GGA-PBE+U（U_d=4.5 eV处理Cu 3d电子）算出ΔGf = -23.1 kJ/mol。两者差距0.7 kJ/mol，差距不会说谎——这是CALPHAD经验参数与DFT第一性原理之间的”系统差异量级”，与Liu、Chen等人2020年关于Al-Cu-Mg体系的方法学对照工作的结论一致。

三、DFT+热力学：从端基组元到多元溶液模型

更系统的做法是DFT+热力学联合方法。第一步：DFT计算端基组元的形成焓。FCC-Al的晶格能、FCC-Cu的结合能、HCP-Mg的晶格能、Al₂Cu（θ）、Al₂CuMg（S）、Al₃Mg₂（β）、AlCuMg（τ）共6个端基相的ΔH全部用VASP计算。交换关联泛函用PBE，PAW势处理所有金属元素，ENCUT = 520 eV，k点密度按12×12×12 Monkhorst-Pack网格。

第二步：把DFT的ΔH数据用Compound Energy Formalism（CEF）拟合进CALPHAD数据库，参数包括L⁰_{Al,Cu:Mg}、L¹_{Al,Cu:Mg}等相互作用参数。这一步认定”采用CEF而非简单亚规则溶液模型”更适合本场景的判断依据是：Al-Cu-Mg三元系中Mg在θ相里实际存在可观的固溶度（0.5-1.2 at%），而亚规则模型无法处理这种”两相区中的第三组元溶解”。

四、亚稳相预测：被”卡住”的两周

把DFT数据回灌到CALPHAD后，Thermo-Calc重新计算的Al-3.5Cu-1.2Mg在300-500 K的亚稳相图（TOW=10⁻⁶，禁用S相作为平衡相）显示：在该亚稳约束下，θ’相和S’相分别在约320 K和约360 K附近为优势析出相。但275 ℃（即548 K）的吸热峰在亚稳相图上没有对应。

项目组在275 ℃的信号上卡了两周。后来回头查DSC曲线发现，这个峰的焓变只有0.18 J/g，比θ’的1.24 J/g小一个数量级。一个可能性是GP区（Guinier-Preston区）——它就是Cu原子在α-Al基体中形成的富Cu原子团簇，没有确定的晶体结构，但在DSC上确实会以微小吸热峰出现（对应GP区溶解）。

这一假说被HRTEM观测间接证实：在峰值温度下淬火的样品中观察到宽约2-3 nm的盘状GP区，与文献中Marion、Pollnow等人的工作吻合。

五、热力学一致性的回溯检验

论文投稿前，审稿人要求做热力学一致性检验。CALPHAD数据库中所有参数必须满足Redlich-Kister多项式在边界处的连续性，以及亚规则模型Gibbs-Duhem方程的相容性。验证结果是：Al-Cu二元系中L⁰_{Al,Cu}参数在800 K以上出现轻微的非物理跳跃（d²G/dx² < 0的”反常区”）。

这一发现被证明不可忽视——它意味着该参数在高温下预测的两相分离倾向会与实验数据冲突。把L⁰_{Al,Cu}从原有-32000 J/mol调整到-28700 J/mol，Chi²从2.41降到0.86，预测的两相区边界在500 ℃附近向高Cu方向移动约0.8 at%，与Matsuyama、Takeda等人2021年的最新实验评估一致。

六、危机：审稿人指出的”自由能符号争议”

第二轮审稿意见中最尖锐的一条是关于S相的ΔGf符号。审稿人引用2019年一项用PBE0杂化泛函的工作，报告S相的ΔGf = +0.8 kJ/mol（即热力学不稳定），与本工作的-23.1 kJ/mol差异巨大。

项目组回查发现，2019年那项工作的DFT计算里Cu的+U修正取5.5 eV，过高估计了Cu 3d的局域化，使得S相中Cu-Cu键被弱化。把U_d从4.5 eV调到3.0 eV后，DFT重新计算给出S相ΔGf = -21.6 kJ/mol，与CALPHAD的-22.4 kJ/mol差距缩小到0.8 kJ/mol。响应审稿人时附上的U_d参数敏感性分析表格，最终说服了审稿人接受”在不同U_d选择下S相均为稳定相”这一结论。

七、回顾性总结

这个项目最有价值的教训有三。第一，CALPHAD是热力学计算的基础工具，但DFT+热力学联合方法在亚稳相、未知相预测中具有不可替代的价值——前者提供外推能力，后者提供端基值的物理锚定。第二，DFT+U方法中U参数的选取对含d电子金属的化合物形成焓有显著影响，应做参数敏感性分析而非单点取值。第三，DSC、HRTEM等实验数据是热力学计算结论的”现实约束”，任何与实验对不上的计算结论都应优先怀疑参数设置而非数据本身。

整个Al-Cu-Mg项目的热力学计算链条，从CALPHAD外推到DFT回灌，从DSC对标到U_d敏感性检验，最终给出的是一份”自洽、可复现、与实验一致”的三元相图与亚稳相析出序列。这套方法后来被项目组用于Al-Zn-Mg-Cu四元系的析出行为研究，把整个研发周期缩短了约40%。