3. 电池材料热力学与相变

在相场模拟中，热力学不是背景知识，而是决定“相如何演化”的驱动力来源。锂离子电池正极材料（如 LiFePO₄、LiCoO₂、NCM）在充放电过程中，锂的嵌入与脱出会引发浓度变化，进而驱动固溶体分解、两相共存或有序—无序转变。相场模型中的自由能密度函数 f(c, η, T) 正是把这些热力学信息编码进控制方程的核心。本节围绕锂化/脱锂相变、固溶体与两相共存、自由能曲线构建、相图与平衡态计算四个要点展开，目标是让读者能够从热力学数据出发，构建可用于相场模拟的自由能表达式，并判断平衡态下的相组成。

一、锂化/脱锂相变的热力学本质

以 Li_xFePO₄ 为例，锂含量 x 从 0 变化到 1 的过程中，材料并非均匀连续地改变锂浓度，而是在很大范围内表现为 LiFePO₄（富锂相）与 FePO₄（贫锂相）两相共存。这种行为的根源在于自由能曲线 f(x) 在中间浓度区间出现“非凸”形状。当自由能曲线存在凹向下的区域时，均匀固溶体不再稳定，系统会自发分解为两个浓度不同但自由能总和更低的相。

判断相变类型的关键是自由能对浓度的二阶导数 ∂²f/∂x²。若在整个浓度范围内 ∂²f/∂x² > 0，则自由能曲线凸向下，任意浓度下均匀固溶体都是稳定的，锂化/脱锂表现为连续固溶体行为。若存在 ∂²f/∂x² < 0 的区间，则发生两相分离。相场模拟中，这一信息直接决定是否需要在自由能中引入双阱势或正则溶体模型。

二、固溶体与两相共存的自由能描述

描述固溶体与两相共存最常用的自由能模型是正则溶体模型：

f(c) = Ω·c·(1-c) + k_B·T·[c·ln c + (1-c)·ln(1-c)]

其中 c 为归一化锂浓度，Ω 为相互作用参数（单位 J/mol 或 eV/atom），k_B 为玻尔兹曼常数，T 为温度。第一项为焓贡献，Ω > 0 表示同类原子相互排斥，倾向于相分离；第二项为熵贡献，倾向于均匀混合。两者竞争决定了相行为：当 Ω 足够大或温度足够低时，熵项不足以抵消焓项，自由能曲线出现双阱，两相共存。

在相场模拟中，通常将自由能密度写为 f(c) = W·g(c) + κ|∇c|²，其中 g(c) 为双阱势，W 为势垒高度，κ 为梯度能系数。双阱势的极小值位置对应平衡相浓度，势垒高度控制相变驱动力，梯度能系数控制界面宽度和界面能。这些参数需要与实验或第一性原理计算得到的自由能曲线匹配。

三、自由能曲线构建方法

构建自由能曲线通常有三条路径：第一性原理计算、CALPHAD 热力学数据库、以及实验数据拟合。对于电池材料，第一性原理计算（DFT）可以给出不同锂浓度下的形成能，再通过拟合得到 Ω 和 k_BT 项。CALPHAD 方法则利用已有的热力学数据库，通过优化参数得到多组元体系的自由能表达式。实验数据拟合则基于开路电压（OCV）曲线，因为 OCV 与锂化学势直接相关：V(x) = −(1/F)·∂f/∂x。

实际操作中，一个典型流程是：先获取不同锂含量下的形成能数据点，再用正则溶体模型拟合。例如，假设某材料在 x = 0.1、0.3、0.5、0.7、0.9 处的形成能分别为 −0.02、−0.05、−0.06、−0.05、−0.02 eV/atom，通过最小二乘拟合可得到 Ω ≈ 0.08 eV/atom，T = 300 K 时 k_BT ≈ 0.026 eV。此时 Ω/(k_BT) ≈ 3.1，大于临界值 2，说明存在两相共存区。这些数值属于教学示例，实际材料需用真实计算或实验数据替换。

四、相图与平衡态计算

相图是自由能曲线的几何表达。对于二元体系，相图横轴为浓度，纵轴为温度。两相共存区由公切线构造确定：在自由能曲线上作一条同时切于两个极小值点的直线，切点对应的浓度即为平衡相浓度 c_α 和 c_β。公切线斜率对应化学势，截距对应各相的化学势。相场模拟中，平衡态计算就是求解 ∂f/∂c = μ 的根，得到平衡浓度分布。

对于 LiFePO₄ 体系，实验 OCV 曲线在约 3.4 V 附近出现平坦平台，对应两相共存。平台宽度对应两相区浓度范围，平台高度对应化学势。相场模拟中，若自由能模型正确，模拟得到的平衡浓度应落在实验相图的两相区内。若模拟结果与相图不符，通常需要检查 Ω 值是否合理、梯度能系数是否过大导致界面过宽、或温度设置是否偏离实际工况。

五、从热力学到相场模拟的衔接

在相场模拟中，自由能密度函数需要与 Cahn-Hilliard 方程耦合：

∂c/∂t = ∇·(M ∇μ)
μ = ∂f/∂c − κ∇²c

其中 M 为迁移率，μ 为化学势。若自由能曲线存在双阱，则 μ 在双阱之间会出现非单调变化，驱动浓度场自发分离为富锂和贫锂区域。模拟时，初始均匀浓度场会因微小扰动而放大，最终形成两相共存结构。判断模拟是否达到平衡态的标准是：总自由能不再下降，且各相浓度稳定在公切线切点附近。

实际计算中，常用无量纲化处理：令 c* = c/c_max，f* = f/(k_BT)，x* = x/L，t* = t·M·k_BT/L²。这样可以将方程简化为只含 Ω/(k_BT) 和 κ/(k_BT·L²) 两个无量纲参数。例如，取 Ω/(k_BT) = 3.0，κ/(k_BT·L²) = 0.5，网格数 128×128，时间步长 Δt* = 0.01，可观察到典型的旋节分解形貌。

六、平衡态计算的数值实现

平衡态计算的核心是求解公切线条件：f′(c_α) = f′(c_β) = [f(c_β) − f(c_α)]/(c_β − c_α)。对于正则溶体模型，可解析求解：令 f′(c) = Ω(1−2c) + k_BT·ln[c/(1−c)]，通过数值迭代（如牛顿-拉夫森法）可得到 c_α 和 c_β。例如，Ω/(k_BT) = 3.0 时，计算得到 c_α ≈ 0.07，c_β ≈ 0.93。这些值可作为相场模拟的初始条件或验证基准。

对于多组元体系（如 NCM 三元材料），自由能曲线需要扩展为多元函数，公切线变为公切面。此时通常借助 CALPHAD 软件（如 Thermo-Calc）计算平衡相组成，再将结果导入相场模型。相场模拟中，多组元扩散耦合会引入交叉扩散项，需要额外考虑 Onsager 系数矩阵。这部分内容将在后续章节展开。

总结来说，本节的核心逻辑是：自由能曲线决定相变类型，公切线构造确定平衡相浓度，相图给出温度-浓度全局视图，而相场模拟则是在这些热力学约束下求解浓度场的时空演化。掌握自由能曲线的构建与解读，是连接热力学数据与相场模拟的关键桥梁。