VASPKIT：波函数导出与宇称分析
WAVECAR必须由自己的VASP计算生成并保留。实部、虚部与模平方含义不同，波函数整体相位可以改变，不能直接用实部正负判断电荷转移。
阅读本篇目录
WAVECAR必须由自己的VASP计算生成并保留。实部、虚部与模平方含义不同，波函数整体相位可以改变，不能直接用实部正负判断电荷转移。
先准备示例#
下载输入、后处理数据与绘图脚本，解压后进入目录：
unzip vaspkit-examples-*.zip
cd vaspkit-examples
python3 prepare.py
module load vaspkit/1.5.1
prepare.py从集群安装的官方examples复制需要的原始示例；附件本身不转发这些分发文件。每个功能的目录独立，下面的cd cases/...均从附件根目录开始。完成一个功能后可用cd ../..回到根目录。
主要例子的文件#
本站VASP教程已在集群计算的 advanced-20261024/nacl-density 数据。
POSCAR#
Cl Na
1.0000000000000000
0.0000000000000000 2.8500000000000001 2.8500000000000001
2.8500000000000001 0.0000000000000000 2.8500000000000001
2.8500000000000001 2.8500000000000001 0.0000000000000000
Cl Na
1 1
Direct
-0.5000000000000000 0.5000000000000000 0.5000000000000000
0.0000000000000000 0.0000000000000000 -0.0000000000000000
INCAR#
SYSTEM = nacl-density
ENCUT = 550
PREC = Accurate
EDIFF = 1E-6
EDIFFG = -0.02
IBRION = -1
NSW = 0
ISIF = 2
ISMEAR = 0
SIGMA = 0.05
LREAL = .FALSE.
LASPH = .TRUE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 150
ALGO = Normal
ISYM = 2
LAECHG = .TRUE.
LELF = .TRUE.
NPAR = 1
NGXF = 120
NGYF = 120
NGZF = 120
KPOINTS#
Gamma mesh
0
Gamma
7 7 7
0 0 0
完整的配套文件通过附件和prepare.py准备。若本功能需要WAVECAR/WAVEDER，先按新计算与波函数准备运行自己的VASP，不使用其他体系的波函数代替。
各功能操作#
511 · 波函数实虚部#
输入k点1、能带1，分别导出实部和虚部；不同版本WAVECAR格式需匹配。
实例来源：本站VASP教程已在集群计算的 advanced-20261024/nacl-density 数据。
在附件根目录执行：
cd cases/511
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
511
1
1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、OUTCAR、POSCAR。
本次得到：WFN_IMAG_B0001_K0001.vasp、WFN_REAL_B0001_K0001.vasp。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
Which kpoint-index do you want to plot? (1<= ikpt <=20)
Which band-index do you want to plot? (1<= iband <=8)
SPIN = 1
NKPTS = 20
NBANDS = 8
ENCUT = 550.00
Maximum number of G values: GX = 8, GY = 8, GZ = 8
Estimated maximum number of plane-waves (NPWS): 1642
Meshgrid Size: 81 81 81
512 · 波函数实虚部Cube#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 advanced-20261024/nacl-density 数据。
在附件根目录执行：
cd cases/512
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
512
1
1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、OUTCAR、POSCAR。
本次得到：WFN_IMAG_B0001_K0001.vasp.cube、WFN_REAL_B0001_K0001_UP.vasp.cube。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
Which kpoint-index do you want to plot? (1<= ikpt <=20)
Which band-index do you want to plot? (1<= iband <=8)
SPIN = 1
NKPTS = 20
NBANDS = 8
ENCUT = 550.00
Maximum number of G values: GX = 8, GY = 8, GZ = 8
Estimated maximum number of plane-waves (NPWS): 1642
Meshgrid Size: 81 81 81
513 · 反演宇称#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
这是实验性入口；本例波函数在当前程序中出现平面波数不匹配，不能据此判断宇称或拓扑。请参考本站VASP能带不可约表示教程，用经过波函数格式核对的IrRep流程。
本例未提供可直接采用的原生结果。
515 · 波函数模平方#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 advanced-20261024/nacl-density 数据。
在附件根目录执行：
cd cases/515
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
515
1
1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、OUTCAR、POSCAR。
本次得到：WFN_SQUARED_B0001_K0001.vasp。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
Which kpoint-index do you want to plot? (1<= ikpt <=20)
Which band-index do you want to plot? (1<= iband <=8)
SPIN = 1
NKPTS = 20
NBANDS = 8
ENCUT = 550.00
Maximum number of G values: GX = 8, GY = 8, GZ = 8
Estimated maximum number of plane-waves (NPWS): 1642
Meshgrid Size: 81 81 81
516 · 波函数模平方Cube#
按下面顺序输入；原子编号从1开始，元素名与当前POSCAR一致。更换体系时先核对所需的文件和编号。
实例来源：本站VASP教程已在集群计算的 advanced-20261024/nacl-density 数据。
在附件根目录执行：
cd cases/516
vaspkit < answers.in
answers.in实际内容如下，每行对应一次输入；也可运行vaspkit后逐行手动键入。
516
1
1
目录中的主要配套文件：CONTCAR、DOSCAR、EIGENVAL、INCAR、KPOINTS、OUTCAR、POSCAR。
本次得到：WFN_SQUARED_B0001_K0001.vasp.cube。原始输入与结果保持在同一例子目录；表格文件先查看首行列名，再选择横纵轴。
本例输出摘要（数值仅对应上述输入）：
Which kpoint-index do you want to plot? (1<= ikpt <=20)
Which band-index do you want to plot? (1<= iband <=8)
SPIN = 1
NKPTS = 20
NBANDS = 8
ENCUT = 550.00
Maximum number of G values: GX = 8, GY = 8, GZ = 8
Estimated maximum number of plane-waves (NPWS): 1642
Meshgrid Size: 81 81 81
换成自己的体系#
先复制一份计算目录，确认POSCAR、INCAR、KPOINTS和各结果来自同一计算，再运行所需菜单。涉及静态、优化或声子等前置VASP任务，先完成计算再后处理；VASPKIT不会替代这些计算。保留结构、输入参数、工具版本和提取后的数据，后续改图时不必反复重算。
方法参考：VASPKIT官方教程；菜单与交互输入以本专题实际使用的安装版本为准。使用结果开展研究时引用VASPKIT论文。