VASP计算在锂离子电池正极材料研究中的关键参数与经验

VASP计算在锂离子电池正极材料领域几乎成了标配工具。前年帮一个课题组做NCM三元正极材料的电压曲线预测，当时他们用GGA算出来的电压比实验值系统偏低0.3-0.5V。这不是第一次遇到GGA对过渡金属氧化物电压预测偏低的案例，但这次客户要求把偏差控制在0.1V以内——因为他们在做材料组分筛选，0.1V的偏差直接影响候选材料的排序。

ENCUT与k点：说烂了但每次还会踩的坑

VASP计算里ENCUT和k点是最基础的参数，但正因为基础，容易被当成”按惯例设置”就跳过去。NCM材料跑结构优化时，按POSCAR手册的建议设了ENCUT=520 eV和500 eV的对比测试——两者算出来的晶格参数差了约0.01 Å，但总能量差了近30 meV/f.u.。

对于电压预测这个量级不够友好。锂离子电池的工作电压由Li脱嵌前后的总能差决定：V = -(E[脱锂后] – E[脱锂前]) / nF。30 meV的总能差如果传播到电压上，会带来约0.03-0.05V的偏差。单看不大，但加上后面泛函带来的偏差，很容易就超出0.1V的误差上限。

k点密度也有类似的情况。Gamma-centered的4×4×2网格和6×6×3网格，算出来的Li扩散势垒能差0.05 eV。对扩散系数评估来说，这个量级会通过Arrhenius关系放大一个数量级。如果时间允许，k点密度收敛性测试是最值得花时间的预处理步骤。

磁矩：NCM体系绕不开的变量

过渡金属氧化物的磁基态是VASP计算里最容易翻车的地方。NCM里Ni以Ni²⁺/Ni³⁺/Ni⁴⁺混合价态存在，Co³⁺是低自旋，Mn⁴⁺没有磁矩。初始磁矩设得不对，电子步跑几百步也未必收敛到正确的基态。

在这个项目里，Ni的初始磁矩设了5 μB（高自旋的2.0 μB太少），Co设0 μB（低自旋基态）或者0.5 μB试探。测试了几种初始磁构型组合之后发现，对富Ni的NCM组分，初始磁矩的选择对总能影响达到了0.1 eV/f.u.——和泛函选择的影响同级。

还有一个容易被”默认值”覆盖的参数是LMAXMIX。对含d电子的过渡金属氧化物，INCAR里设LMAXMIX=4能让电荷密度混合收敛速度明显加快，不设的话电子步迭代经常在这个体系里卡到60步才勉强收敛。

电压计算：从总能到电压，中间有一整套校正链

有了结构优化和电子结构的基础数据之后，电压计算本身是后处理。但这里最容易丢失精度的地方不是VASP本身，而是参考态的选择。

金属锂的体相能量作为Li⁺/Li参考电极的DFT能量，跑的时候需要注意：锂是bcc结构，k点要足够密（至少12×12×12），否则算出来的参考能量有0.1-0.2 eV的偏差。这个偏差直接加到每一组电压计算里，导致所有候选材料的电压被系统地抬高或压低。

另一个重要的校正项是范德华相互作用。NCM这类层状材料，层间有非键相互作用，不加范德华修正（IVDW=11或DFT-D3）的话，层间距会高估，直接影响到Li在层间的迁移势垒，进而影响倍率性能的评估。

团队在这个项目上最后给出的电压预测，偏差控制在了0.08V以内。方案上线后过了一段时间，在同领域其他课题组的文章里看到了对照结果——实验测试数据和计算预测之间的差距大约在0.05-0.1V，基本在所承诺的精度范围之内。在科研学术网首页上能找到更多电池材料DFT计算的实操案例和参数对照。

回头看

VASP本身是一个工具，工具怎么用带来的结果差异可能比工具本身的版本差异还大。这个项目带来的最大收获不是算出了哪几组具体的数据，而是弄清楚了电池正极材料DFT计算中每一个参数对最终电压预测的定量贡献——从ENCUT到磁矩，从参考态到范德华修正——形成了一个内部校验清单。这个清单后来成了每个电池材料项目启动时的第一轮自检步骤。