VASPKIT：静电势、平面平均与真空参照
沿真空方向平均LOCPOT，可用于检查真空平台和电势变化。应先说明使用LVHAR的离子加Hartree势还是LVTOT总局域势，不混用两种定义。
阅读本篇目录
沿真空方向平均LOCPOT，可用于检查真空平台和电势变化。应先说明使用LVHAR的离子加Hartree势还是LVTOT总局域势，不混用两种定义。
先准备示例#
下载输入、后处理数据与绘图脚本，解压后进入目录：
unzip vaspkit-examples-*.zip
cd vaspkit-examples
python3 prepare.py
module load vaspkit/1.5.1
prepare.py从集群安装的官方examples复制需要的原始示例；附件本身不转发这些分发文件。每个功能的目录独立，下面的cd cases/...均从附件根目录开始。完成一个功能后可用cd ../..回到根目录。
本例结果图#
图由实际导出的数据绘制。附件根目录执行module load phonopy/4.8.3，再执行python plot.py可重画；此模块提供绘图所需的Python环境。
主要例子的文件#
本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
POSCAR#
B N
1.0000000000000000
2.5000000000000000 0.0000000000000000 0.0000000000000000
-1.2500000000000000 2.1650635094610964 0.0000000000000000
0.0000000000000000 0.0000000000000000 20.0000000000000000
B N
1 1
Direct
0.0000000000000000 0.0000000000000000 0.5000000000000000
0.3333333333333333 0.6666666666666666 0.5000000000000000
INCAR#
SYSTEM = hbn-field-0.00
ENCUT = 550
PREC = Accurate
EDIFF = 1E-6
IBRION = -1
NSW = 0
ISMEAR = 0
SIGMA = 0.05
LREAL = .FALSE.
LASPH = .TRUE.
LWAVE = .FALSE.
LCHARG = .FALSE.
LVTOT = .FALSE.
LVHAR = .TRUE.
NELM = 200
ALGO = Normal
ISYM = 0
NBANDS = 32
LDIPOL = .TRUE.
IDIPOL = 3
DIPOL = 0.5 0.5 0.5
EFIELD = 0
LORBIT = 11
KPOINTS#
Gamma
0
Gamma
9 9 1
0 0 0
完整的配套文件通过附件和prepare.py准备。若本功能需要WAVECAR/WAVEDER，先按新计算与波函数准备运行自己的VASP，不使用其他体系的波函数代替。
各功能操作#
420 · 电势切片#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/420
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
420
1
0 0 1
0.5
1 1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：SLICE_001.grd、X.grd、Y.grd。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
1) LOCPOT
2) Specified Potential-Related File
Select One Atom to Determine the Postion of Plane [1 <= & <= 2]
(Tip: MUST be two integers, e.g., 2 3)
1 0 0
0 1 0
0 0 1
Lattice Constants in New Cell: 2.500 2.500 20.000
Lattice Angles in New Cell: 90.00 90.00 120.00
Total Atoms in New Cell: 2
Interpolation Grid: 62 62 500
422 · 电势线平均#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/422
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
422
1
1 1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：POTLAVG.grd、X.grd、Y.grd。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
1) ab Plane
2) ac Plane
3) bc Plane
(Tip: MUST be two integers, e.g., 2 3)
Maximum value : 1.072
Minimum value : -1.750
425 · 沿路径取电势#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/425
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
425
1
0 0 0
0 0 1
101
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：POTENTIAL_1D.dat。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
1) LOCPOT
2) Specified Potential-Related File
426 · 电势平面平均#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/426
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
426
3
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：PLANAR_AVERAGE.dat。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
1) Lattice a Direction
2) Lattice b Direction
3) Lattice c Direction
Vacuum-Level (eV): 1.347
Work Function (eV): 5.524
427 · 电势宏观平均#
宏观窗口2.5 Å。真空平台判断以未任意平移的势能及同一费米参照为准。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/427
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
427
3
2.5
1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：MACROSCOPIC_AVERAGE.dat。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
1) Lattice a Direction
2) Lattice b Direction
3) Lattice c Direction
(Typical value: interlayer spacing, ~1-3 A)
(Must be integer. You can gradually increase from 1 to 3.)
Maximum of macroscopic average: 1.347
Minimum of macroscopic average: -8.976
428 · 电势导出XSF#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/428
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
428
LOCPOT
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：LOCPOT、LOCPOT.xsf。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
可用head 文件名查看表头，结构结果可在VESTA中打开，曲面数据按对应网格绘制。
429 · 电势导出Cube#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/429
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
429
LOCPOT
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：LOCPOT、LOCPOT.cube。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
可用head 文件名查看表头，结构结果可在VESTA中打开，曲面数据按对应网格绘制。
430 · 电势超胞#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 band-alignment-20261025/hbn-monolayer-hartree 数据。
在附件根目录执行：
cd cases/430
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
430
1
2 0 0
0 1 0
0 0 1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、LOCPOT、OUTCAR、POSCAR、PROCAR。
本次得到：LOCPOT、SLOCPOT.vasp。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
(MUST be three integers, e.g., 1 2 3)
Enter the 2nd row of Transformation Matrix to redefine cell:
(MUST be three integers, e.g., 1 2 3)
Enter the 3th row of Transformation Matrix to redefine cell:
(MUST be three integers, e.g., 1 2 3)
2 0 0
0 1 0
0 0 1
Lattice Constants in New Cell: 5.000 2.500 20.000
Lattice Angles in New Cell: 90.00 90.00 120.00
Total Atoms in New Cell: 4
Interpolation Grid: 99 50 399
换成自己的体系#
先复制一份计算目录，确认POSCAR、INCAR、KPOINTS和各结果来自同一计算，再运行所需菜单。涉及静态、优化或声子等前置VASP任务，先完成计算再后处理；VASPKIT不会替代这些计算。保留结构、输入参数、工具版本和提取后的数据，后续改图时不必反复重算。
方法参考：VASPKIT官方教程；菜单与交互输入以本专题实际使用的安装版本为准。使用结果开展研究时引用VASPKIT论文。