LOCPOT 怎么用于表面电势分析？

引言

LOCPOT 不是只用来画一条漂亮的电势曲线。它能帮助判断真空能级、功函数、界面电势降和偶极修正是否可靠。真正要盯的是三件事：平面平均后有没有稳定平台，平台读取是否避开原子层振荡，吸附或界面引起的电势变化是否能被结构和电荷分布一起解释。只要这三件事没说清，一个功函数或电势降数字就很容易被写过头。

LOCPOT 里到底存了什么？

为什么不能直接看三维电势？

LOCPOT 里的局域电势本来是三维网格数据，原子核附近会有很尖的振荡。直接在三维图里挑一个点读数，读到的往往是局部起伏，而不是 slab 外侧真正可作为能级参照的真空电势。做功函数、真空能级或界面电势降时，第一步应把电势沿表面平面平均，先得到随 z 方向变化的一维曲线。

在 VASP 结果里，LOCPOT 提供局域电势，OUTCAR 提供 E_F。只有把平面平均后的真空平台和 E_F 放到同一套计算里比较，功函数Φ 才有清楚来源。这里要注意，平面平均不是为了把图画得更平滑，而是为了把原子层附近的快速振荡和真空区的宏观平台分开。

式(1) V(z) =

真空平台为什么是关键？

图1应该先看平台，而不是先看曲线高低。比较不同设置时，重点看真空区是否出现足够长的平坦区，左右两侧平台是否一致，以及吸附或偶极修正是否让平台整体倾斜。若平台本身不稳，后面再报功函数小数点后三位也没有意义。

如果真空区没有平台，不能强行取最高点或最低点当 E_vac。这种读法会把原子层振荡、真空层不足或偶极未修正带来的漂移误当成物理差异。更稳的做法是先加厚真空层、检查偶极修正，再重新判断平台区间是否收敛。

实际保存结果时，建议把 LOCPOT、OUTCAR、平面平均脚本、平台读取区间和最终数值放在同一目录里。后面若功函数或电势降出现异常，能直接回头检查到底是平台没取稳，还是模型本身确实引起了电势变化。

图1 不同计算设置下沿 z 方向的静电势曲线对比。DOI: 10.1021/acsomega.0c04420。

平面平均曲线怎样读？

平台位置应该怎样选？

平面平均后先看真空层厚度，再看左右平台是否对称，最后看吸附或界面导致的电势偏移。这一步最怕把默认设置当作真理。不同材料、不同吸附态、不同电荷边界下，后处理图的可读性和能量趋势都会变。

只改变一个变量，其余参数保持一致。平台读取要避开原子层振荡和边界噪声。

式(2) Φ = E_vac – E_F

电势波动说明什么问题？

带边相对真空电势的图能说明参照的重要性。图2能帮助说明这个判断，因为它把理论量和可讨论的材料差异放在同一幅图里。

这里不要急着追求一句万能结论。电势曲线漂移可能来自真实偶极，也可能来自真空层太薄或偶极修正缺失。电势波动可能意味着真空层或偶极处理不足。

如果两个模型给出相反趋势，先检查结构是否真的可比，再检查参照态、磁性、电荷、自旋和收敛阈值。很多“机制反转”其实来自输入条件不一致，而不是材料发生了神秘变化。

图2 相对于真空电势与 NHE 的能带边位置计算示例。DOI: 10.1039/d4ra08172b。

偶极修正什么时候必须检查？

LDIPOL 和 IDIPOL 怎么配合？

不对称 slab 通常需要检查 LDIPOL = .TRUE.、IDIPOL = 3，并确认 DIPOL 位置合适。建议把这个流程写成固定清单，而不是每次凭记忆操作。POSCAR、INCAR、KPOINTS、POTCAR、OUTCAR、CHGCAR、LOCPOT 或 vasprun.xml 各自承担不同角色，混在一个目录里后期很难复查。

对关键物理量，最好保留原始文件、提取脚本和最终数值三层证据。不对称 slab 应把偶极修正作为复查项。

式(3) ΔV = V_right – V_left

吸附前后怎样比较功函数？

界面模型的电势剖面能展示不同边界条件。图3适合提醒我们，图像解释必须回到计算对象本身：材料、位点、吸附物、覆盖度和边界条件都要说清楚。

若开启偶极修正后Φ、E_vac 或ΔV 大幅变化，说明原模型对长程电势敏感。如果只截一张好看的图，而没有能量、结构或电子态旁证，这张图很难独自撑住结论。功函数比较必须使用同一套 E_vac 读取规则。

交付前建议逐项复核上下标。像 E_F、E_ads、ΔG_*H、O₂、H⁺、e^– 这样的符号必须用真实上下标，不能混用普通字符或残留 HTML 标记。

图3 slab 与纳米液滴模型中的静电势剖面对比。DOI: 10.1038/s41467-025-60298-z。

表面电势怎样写进论文？

怎样避免把数值漂移写成机制？

论文中应说明平台读取方式、偶极修正和真空层收敛，而不是只报一个功函数数字。比较稳的写法是先描述观察，再给出支撑，再限定适用条件。不要把一个后处理图直接翻译成“活性提高”，也不要用单个指标替代整个反应路径。

在某模型、某参照和某收敛标准下，某指标支持某种趋势。论文中不要把数值漂移直接解释为电子转移。

式(4) E_vac = plateau[V(z)]

界面电势降需要哪些旁证？

界面电势降不是从曲线里随手截两个点就结束。先看结构旁证：水层厚度、分子取向、氢键网络、吸附物位置和界面间距是否可比。若两个模型的水层排列已经完全不同，ΔV 的变化可能来自结构重排，而不是单纯的电荷转移。

再看电荷旁证。差分电荷密度可以告诉电子在哪一侧累积、在哪一侧耗尽；平面平均电荷差可以把这种重排投影到 z 方向；Bader 电荷只能作为净得失的旁证，不能替代电势曲线本身。

还要看边界条件。真空平台是否稳定、LDIPOL/IDIPOL 是否一致、DIPOL 位置是否合理、真空层厚度改变后 ΔV 是否收敛，这些都会影响界面电势降。尤其是不对称 slab，偶极修正前后的曲线最好都保留。

图4的意义就在这里：界面水层诱导的电势降要和电荷密度一起读。更稳的写法不是“电势降变大，所以电荷转移增强”，而是先说明结构和水层取向，再给出 ΔV、差分电荷和平面平均电荷差，最后限定这个判断只对应当前界面模型。

图4 界面水层诱导的电势降与电荷密度关系。DOI: 10.1021/jacs.5c04728。

本章要点总结

要点1：LOCPOT 平面平均能帮助读取真空电势、功函数和界面电势变化。

要点2：真空区必须出现稳定电势平台，否则 E_vac 和功函数读数不可靠。

要点3：不对称 slab 应检查偶极修正前后的电势曲线。

要点4：界面电势降要结合电荷密度和模型边界条件解释。