锂镍钴锰热力学计算

结合案例图分析，第一张案例图横轴是Li含量x（在Li_xNi_{1/3}Co_{1/3}Mn_{1/3}O₂中），范围1.0到0.3；纵轴是电池电压（V vs Li/Li⁺）。曲线整体呈下降趋势，从x=1.0时的约4.45 V下降到x=0.3时的约3.9 V。图中标注”charge plateau ~3.9 V”，对应x≈0.55附近的电压平台。这个平台通常与层状结构到尖晶石或岩盐相变的两相共存区有关，是NCM材料在深度脱锂时的特征。继续脱锂到x<0.5，电压快速下降，说明材料热力学稳定性变差。箭头标注”delithiation Li_xNiCoMn → Li-poor”，表示脱锂方向。

第二张案例图横轴是Mn摩尔分数x（在Li(NiCoMn)O₂中），范围0到1；纵轴是形成能ΔH_f（eV/f.u.）。蓝色曲线是DFT计算的各成分形成能，红色虚线是连接稳定端点（x=0和x=1）的基态线（ground state line，也叫hull）。红色圆点是位于hull下方的稳定或亚稳组成。可以看到在x≈0.5附近，形成能达到最低值（约−0.09 eV/f.u.），说明这个Mn比例在热力学上最有利，对应常见的NCM111组成。其他成分如果落在hull上方，则是热力学不稳定的，会发生相分离。

从电压-锂含量曲线和形成能凸包图的组合来看，这是一个锂镍钴锰（NCM）三元正极材料的热力学计算项目——客户想知道不同SOC下的开路电压曲线以及Mn含量对相稳定性的影响。电压曲线回答了电化学性能问题，凸包图回答了组成稳定性问题。

项目计算方法上，我用VASP做了NCM不同Li含量和不同Mn比例的超胞计算。采用PBE+U方法，Ni、Co、Mn的U值分别取6.0、3.3、4.0 eV。Li脱嵌的电压用公式V=−[E(Li_{x2})−E(Li_{x1})−(x2−x1)E(Li)]/[(x2−x1)e]计算。形成能用凸包法，比较各组成相对于端点相的能量。

我重点关注了两个问题。第一个是Li/Ni混排：NCM合成中常出现Li占位被Ni占据的阳离子混排，这会改变局部配位环境和电压。我在计算中考虑了不同混排度，发现混排使平均电压降低约0.05 V。第二个是氧空位：深度脱锂时氧释放是安全隐患，我通过氧空位形成能评估了x<0.5时的热力学不稳定性。

踩坑经验有三条。第一个是U值选择：不同文献对Mn的U值差异很大，我先用氧化物标定，确保带隙和磁矩合理。第二个是构型熵：高温下Li和过渡金属的构型熵不可忽略，我用cluster expansion结合Monte Carlo方法做了有限温度修正。第三个是电压平台：单相DFT计算给出的是平滑电压曲线，实验中的平台来自两相共存，需要用常见切线法或相图计算才能得到。

这张图对客户的核心结论是：NCM111在x≈0.5时热力学最稳定，充电至x<0.5会出现约3.9 V的电压平台，提示相变风险。关于锂镍钴锰热力学计算的更多细节，我整理在锂镍钴锰热力学计算栏目里。