在表面吸附、非对称薄膜和孤立分子计算中，周期性边界条件可能让具有净偶极矩的体系与自身周期镜像发生非物理相互作用，从而影响能量、势和受力。处理“VASP怎么设置偶极修正，VASP偶极修正后体系总能量异常如何排查”时，不能只看修正前后能量是否发生变化，还要结合体系维度、真空层、偶极方向以及电子自洽情况判断结果是否合理。

　　一、VASP怎么设置偶极修正

　　VASP中的能量修正和势、受力修正由不同参数控制。对于常见的非对称slab，应先确定表面法向对应哪一个晶格矢量，再配置偶极方向；孤立分子的设置方式则有所不同。

　　1、为表面体系指定偶极方向

　　①确认POSCAR中真空层位于哪个晶格方向。

　　②如果真空层沿第三个晶格矢量，也就是常见的z方向，可在INCAR中设置：

　　LDIPOL=.TRUE.

　　IDIPOL=3

　　③如果表面法向对应第一或第二晶格矢量，则分别将【IDIPOL】设为【1】或【2】。

　　④重新计算后检查OUTCAR中的偶极矩，确认主要分量位于预期方向。

　　【IDIPOL=1、2、3】分别计算第一、第二和第三晶格矢量方向上的偶极修正，VASP官方建议slab体系让该方向与表面法向一致。【LDIPOL=.TRUE.】则进一步修正周期镜像导致的势和受力误差。

　　2、设置偶极中心

　　偶极矩的计算依赖参考中心。对于表面和孤立分子，官方建议让DIPOL大致对应体系原子的质心位置，并使用直接坐标表示。

　　①根据结构计算或估算原子质心在晶胞中的直接坐标。

　　②例如slab基本位于晶胞中心，可以先测试：

　　DIPOL=0.5 0.5 0.5

　　③如果slab整体偏离晶胞中心，则根据实际结构调整【DIPOL】。

　　④比较不同合理DIPOL位置下的偶极矩与总能量，确认结果没有对参考点表现出异常敏感。

　　如果电荷密度在真空方向没有充分衰减，偶极矩会对参考中心变得非常敏感，这通常提示真空层或晶胞设置还需要继续检查。

　　3、孤立分子采用全方向修正

　　①确认计算对象属于放置在真空超胞中的孤立分子。

　　②将【IDIPOL】设置为【4】，计算三个方向上的完整偶极修正。

　　③需要同时修正势时，再启用【LDIPOL=.TRUE.】。

　　④保证分子与周期镜像之间留有足够空间，再进行超胞尺寸收敛测试。

　　VASP官方将【IDIPOL=4】推荐用于孤立分子，而表面体系通常只处理表面法向上的偶极分量。

　　二、VASP偶极修正后体系总能量异常如何排查

　　启用偶极修正后总能量发生变化本身并不代表计算错误，因为修正项本来就会改变有限超胞中的静电能量。真正需要关注的是能量突然大幅跳动、自洽无法稳定、不同真空层下结果持续漂移，或者相近结构之间出现不合理的能量排序。

　　1、检查IDIPOL方向是否设置错误

　　①重新查看POSCAR三个晶格矢量，确认真正的表面法向。

　　②检查INCAR中的【IDIPOL】是否与该方向一致。

　　③如果真空沿z方向，却误设成【IDIPOL=1】或【IDIPOL=2】，改为【IDIPOL=3】重新计算。

　　④重新读取OUTCAR中的偶极矩，观察主要分量是否符合结构非对称方向。

　　对于slab，只有垂直表面的偶极分量具有明确的修正意义，方向选错会让能量修正对应错误的周期方向。

　　2、检查真空层是否足够

　　①保持slab结构和面内晶格参数不变，只增加真空层厚度。

　　②分别计算不同真空尺寸下的总能量、偶极矩和功函数等结果。

　　③如果扩大真空后能量仍持续明显变化，继续检查电荷密度是否已经在真空区域衰减。

　　④不要只比较单个真空尺寸下“开修正”和“不开修正”的能量差。

　　官方指出，晶胞过小时可能出现电荷向真空区域来回移动，既影响偶极矩精度，也会让电子自洽明显变慢；增加真空宽度通常有助于改善这一问题。

　　3、先排除LDIPOL引起的自洽问题

　　【LDIPOL】会修改局域势，因此开启后电子基态可能需要更多迭代才能稳定。如果计算没有真正收敛，最终总能量和受力都可能失去可靠性。

　　①先使用【LDIPOL=.FALSE.】完成一次稳定的电子自洽计算。

　　②保留生成的WAVECAR。

　　③再从该WAVECAR重新启动，并改为【LDIPOL=.TRUE.】。

　　④如果仍出现电荷振荡，可尝试把【AMIN】减小到【0.01】。

　　⑤必要时增加【NELM】，并将【EDIFF】收紧到【1E-6】后重新检查收敛。

　　VASP官方的收敛建议明确推荐先关闭LDIPOL预收敛，再利用已有WAVECAR启动带偶极修正的计算。

　　4、确认体系是否带净电荷

　　①检查INCAR中的【NELECT】，确认电子数是否人为改变。

　　②如果计算的是带电slab，不要直接把不同真空宽度下的总能量进行普通收敛比较。

　　③保持真空宽度一致时，再比较组成和电荷状态相同的相对能量。

　　④需要处理带电缺陷或其他带电周期体系时，应根据体系维度采用相应的静电有限尺寸修正方法。

　　VASP官方特别指出，带净电荷的表面体系总能量不会随着slab间真空距离增加而正常收敛，因此这种情况下“总能量随真空变化”并不能简单归咎于偶极修正参数错误。

　　三、怎样判断偶极修正后的结果是否可信

　　偶极修正是否有效，不能只根据一次计算的最终能量判断。更可靠的方法是保持研究对象不变，通过真空尺寸、参考中心和修正设置进行小范围对照，让物理量表现出稳定趋势。

　　1、用统一条件比较相对能量

　　①比较吸附构型、表面终止或结构优化结果时，统一使用相同的【IDIPOL】、【LDIPOL】和真空厚度。

　　②相近体系尽量保持一致的【DIPOL】定义原则。

　　③分别记录总能量、偶极矩和最终电子迭代情况。

　　④发现某个结构的修正能量明显偏离其他构型时，再检查电荷重新分布和偶极中心。

　　统一设置能够减少数值条件差异混入本来需要比较的结构能量中。

　　2、做真空层收敛验证

　　①选取一个代表性结构，准备多个不同真空层厚度。

　　②保持ENCUT、K点、结构和其他INCAR设置一致。

　　③分别完成带偶极修正的计算。

　　④比较能量差、偶极矩及相关表面物性是否逐渐稳定。

　　⑤达到研究所需精度后，再采用相同设置计算其他构型。

　　对于具有净偶极矩的表面，足够的真空空间和稳定的偶极定义是获得可靠有限尺寸修正的重要条件。

　　总结

　　“VASP怎么设置偶极修正，VASP偶极修正后体系总能量异常如何排查”的核心，在于判断周期镜像产生的静电误差是否被合理控制，并保证不同结构之间的能量比较建立在一致的计算条件上。对于非对称表面、孤立分子和带电体系，偶极修正的适用方式并不完全相同，尤其不能把带电slab的有限尺寸问题简单理解成参数设置错误。希望本文对大家判断VASP偶极修正结果有所帮助，如果在表面体系静电修正、总能量异常或真空层收敛判断方面还有疑问，欢迎联系咨询。