电极材料的结构优化，在 DFT 计算流程的第一页，却是最容易埋下隐患的地方。一个没有收敛到位的结构，后面算的电子结构、Li 扩散能垒、插层电位——都会带着系统误差往下游传导。偏偏结构优化又是初学者最不重视的一步——”跑一遍 ISIF = 3 就行了”——事情没这么简单。
ISIF = 3 允许晶格参数和原子位置同时弛豫。对于电极材料来说，这听上去很合理——毕竟实验上的晶格参数有测量误差，理论上应该让 DFT 自己找到能量最低的晶格。但全自由度弛豫在过渡金属氧化物正极（如 NCM、LCO、LFP）中有个实际问题：GGA 泛函对过渡金属氧化物的晶格参数普遍高估 1-2%。
PBE 算 NCM811 的 c 轴晶格参数通常比实验值大 0.2-0.3 Å。这个偏差本身不是致命问题，但如果你用 ISIF = 3 得到了一个比实验大了 2% 的晶格体积，然后在这个晶格上计算 Li 扩散能垒——能垒对晶格体积非常敏感，1% 的体积膨胀可能让扩散能垒降 0.1-0.2 eV。这意味着你得到的是一个”在 DFT 偏大的晶格里 Li 更容易扩散”的假象。
一个更稳的做法是：先固定晶格参数为实验值，只做离子位置弛豫（ISIF = 2）；然后在弛豫后的结构中，对晶格参数做 ±1%、±2% 的扫描，拟合状态方程得到 DFT 最优晶格参数。 多了一步拟合的工作量，但能让结构优化过程和扩散能垒计算解耦——晶格参数用什么值是一个明确的决策，而不是 ISIF = 3 自动给出的”黑箱结果”。
NCM 材料里，Ni、Co、Mn 三种过渡金属在 3b 位点上的分布是”准无序”的——没有长程有序，但短程有偏好。DFT 需要对这种部分占位做近似。常用的办法是用特殊准随机结构（SQS）生成一个有代表性的超胞构型——但这意味着需要 3×2×1 以上的超胞（~100 个原子），计算量翻倍。
在计算资源有限的情况下，很多课题组直接用有序排布做近似：Ni 和 Co 交替排列、Mn 均匀分散。这个近似对于晶格参数和电子结构的计算结果偏差不大（<1% 的晶格偏差），但对 Li 扩散路径的影响就不能忽略了——Li 在 Ni 富集区域的扩散行为和 Mn 富集区域不同，有序排布强行把环境”平均化”了，得到的扩散能垒可能比真实材料更”均匀”。
一个折中方案是生成 2-3 种不同的有序排布构型，分别算扩散能垒，取平均和方差。三个构型之间成能垒差 0.1 eV 以内，说明离子占位对扩散的影响不大；方差超过 0.2 eV，就得考虑是否要用更大的超胞或 SQS 方法来降低近似的误差。
石墨负极和过渡金属氧化物正极的层间相互作用性质完全不同。石墨的层间以范德华力为主，PBE 在石墨层间距上偏差可以大到 0.5-1 Å，不用 vdW 修正直接算石墨嵌锂的插层电压，结果基本不能用。DFT-D3 和 DFT-D2 在石墨体系里都可以把层间距修正到实验值附近，但 D3-BJ 对长程 vdW 的描述更好，层间距偏差更小（~0.02 Å）。
过渡金属氧化物正极就不同了。NCM 的层间除了 vdW，还有 Li-O 的静电作用——vdW 修正对整个晶格的贡献比例比石墨小得多。用 D3 修正之后晶格参数可能反而偏小，因为 D3 在过渡金属氧化物中对晶格有轻微过收缩倾向。这种情况下，vaspkit 辅助的状态方程拟合可以帮助判断加了 vdW 之后的晶格参数是否还在合理范围。
一个经验法则是：石墨类负极 → D3-BJ；过渡金属氧化物 → 可以不加 vdW，或仅用 optB86b-vdW 这类对氧化物体系参数更好的 vdW 泛函。
EDIFFG = -0.01（力收敛标准 0.01 eV/Å）是大多数计算的工作标准。但对于电极材料，特别是含 Li 的正极，单个 Li 离子的受力收敛到这个量级之后，总能量可能还在以 1-2 meV/步的速度往下滑——因为过渡金属氧化物里 Li 离子的势能面非常平缓。力收敛了，能量没稳，这种结构在后续 NEB 计算里容易出现虚假的低能垒。
更严格的做法是追加一步高精度收敛：EDIFF 从默认的 1E-5 降到 1E-6，PREC 从 Normal 升到 Accurate，再跑 20-30 个离子步。往往这一步能让总能量再降 5-10 meV，对于插层电位的计算而言，这 5-10 meV 能直接影响电压平台的排序。
但这些精细收敛操作能不能真正改善结果，最终取决于一个问题：你跑的结构优化到底是为了获得一个”合理”的构型，还是为了做后续的定量分析？ 如果只是定性讨论电子结构趋势，0.01 eV/Å 足够了。如果接下来要跑 NEB 算能垒、算插层电压、和实验做定量对比——把收敛门槛抬高一级，前期的投入会在后期以”数据一致性”的形式收回来。
