形成能DFT计算案例：LiFePO₄中Li空位形成能与电压平台预测

形成能算不准，整个后续分析都建立在错误基础上。这个项目在评估LiFePO₄正极材料的Li脱嵌电压时，沿着缺陷形成能→化学势参考态→电压预测这条链路逐级排查，最终看到的问题是：DFT+U参数怎么选，比选什么超胞尺寸更重要。

缺陷形成能的三块拼图

LiFePO₄中Li空位的形成能定义为：

E_form = E_defect − E_perfect + μ_Li

三块拼图：含Li空位的超胞能量、完美晶胞能量、Li的化学势——前两块是VASP直接算出来的，问题不大。真正的变量在μ_Li。

对于LiFePO₄正极，Li的化学势不能直接用体心立方Li金属的能量。这是因为电池工作条件下Li处于脱嵌态，化学势对应的是Li⁺在电解质中的电化学势，而不是金属Li的标准态。项目采用了两种化学势参考态分别计算：极端贫锂条件下μ_Li取Li金属能量减去实验脱锂电压对应的能量偏移；富锂条件下取与FePO₄平衡时的化学势。两种条件下Li空位形成能分别为2.35 eV和1.87 eV——差了0.48 eV，这个差异最终会传导到电压预测上，不可忽视。

DFT+U：U值选多大，电压就偏多远

Fe是3d过渡金属，GGA/PBE对Fe 3d电子的自相互作用误差会导致d带位置不对，进而影响氧化还原电位。项目对Fe的3d轨道加Hubbard U修正，测试了U=3.0、4.0、4.5、5.0 eV四个值。

U=4.0 eV时计算电压3.42 V，最接近实验值3.45 V。U=3.0给出3.18 V（偏低），U=5.0给出3.61 V（偏高）。这个趋势不是线性的——U从3.0到4.0改变0.23 V，从4.0到5.0改变0.19 V，说明Fe的3d电子在U=4.0附近趋于合理描述。

有意思的是U值对晶格参数的影响比对电压的影响小：U从3.0到5.0，a轴变化不到0.3%，b轴不到0.5%。这意味着如果只用晶格参数来”标定”U值，很可能会选到一个对电子结构描述不够准确的U。

超胞尺寸：2×1×1够不够

项目用2×1×1超胞（共112个原子，含8个Li位）计算单个Li空位。更大的2×2×1超胞（224个原子，16个Li位）形成能变化不到0.03 eV，说明2×1×1的Li空位间相互作用已在可接受范围。但2×1×1的问题是倒空间采样不对称——a轴方向K点折合后比b/c轴稀疏，对能带色散沿a轴方向的体系可能引入额外误差。这个项目最终保留2×1×1但用Gamma-centered K点弥补了部分不对称误差。

整套流程下来，LiFePO₄的形成能计算真正被参数选择左右的部分占比约70%——不是VASP算不准，是参数怎么给决定了结果往哪个方向走。化学参考态和U值是两个最容易被低估的变量。

参考文献：Zhou et al., Phys. Rev. B, 2004, 69, 201101；Ong et al., Energy Environ. Sci., 2011, 4, 3680-3688.