简介本资源为《VASPKIT中文手册1》面向材料科学、催化及计算物理领域的科研人员与初学者解决VASP软件预处理与后处理操作门槛高、工具分散、学习成本大的核心问题。手册系统覆盖VASPKIT 0.72版本全部功能包括输入文件INCAR/POTCAR/POSCAR自动生成与校验、能带计算PBE/HSE06、分波态密度分析、热力学量校正、差分电荷密度绘图、NEB路径生成与虚频校正、能带反折叠、3D能带可视化、费米面计算及VASP2BoltzTraP接口等关键模块并详解用户自定义界面如MoS2投影能带、原子电荷着色结构图等测试功能。资源为单个PDF文件大小2.37MB内容完整、排版清晰附有菜单功能索引与大量示例说明便于快速查阅与实操对照。目前已有3556人学习下载是掌握这一轻量级仅5MB、免安装、Linux原生FORTRAN工具的权威入门与进阶指南。1. VASPKIT中文手册1不是翻译文档而是把第一行命令跑通的实操起点你刚装好VASP跑完Si的静态计算正准备算能带、态密度、电荷密度差分——结果发现vaspkit官网只有一份英文PDF命令像天书“vaspkit -task 211是干啥211和212差在哪为什么我输完回车就报错‘No PROCAR found’但我的INCAR里明明没开LORBIT” 这就是绝大多数人卡在VASPKIT门口的真实状态。VASPKIT中文手册1不是逐句翻译英文PDF而是从零构建一个“能用、敢改、出图快”的本地化工作流它覆盖了从环境变量配置、输入文件校验、到最常调用的5类任务能带/投影能带/电荷分析/声子初筛/高对称路径生成的最小可行命令链。适合两类人一是刚接触第一性原理计算的研究生需要绕过英文术语墙直接产出论文图二是有多年VASP经验但被新版本vaspkit参数变动搞懵的工程师比如v3.4之后-task 303不再默认输出LOCPOT插值网格必须手动加-grid。本文所有命令均在CentOS 7 VASP 6.3.2 vaspkit 3.5.2环境下实测通过不依赖任何GUI或在线服务。2. 环境筑基让vaspkit认得清你的VASP、读得懂你的结构、找得到你的输出VASPKIT不是独立程序它本质是VASP生态的“命令行胶水”——所有功能都建立在对POSCAR、KPOINTS、OUTCAR、EIGENVAL等文件的精准解析上。若基础环境没对齐后续所有操作都是玄学。这里不讲编译源码除非你真要改C只聚焦三个致命环节路径绑定、输入校验、版本兼容。2.1 用export固化VASP路径而非临时PATH很多教程教你在终端敲export PATH/path/to/vaspkit:$PATH但这样每次新开终端都要重输。真实生产环境必须写入shell配置文件# 编辑 ~/.bashrcLinux或 ~/.zshrcmacOS echo export VASP_KIT/opt/vaspkit ~/.bashrc echo export PATH$VASP_KIT/bin:$PATH ~/.bashrc source ~/.bashrc注意VASP_KIT环境变量是vaspkit内部硬编码调用的不是可选项。若漏设执行vaspkit -task 111时会报Error: VASP_KIT environment variable not set!且错误提示藏在stderr里容易被忽略。2.2 用vaspkit -task 0做输入文件健康扫描别急着算能带先让vaspkit当你的“结构医生”。在包含POSCAR、INCAR、KPOINTS、POTCAR的目录下运行vaspkit -task 0它会输出类似[INFO] POSCAR: Orthorhombic cell (1x1x1), 2 atoms, volume15.23 ų [INFO] INCAR: ISPIN2, ICHARG11, LORBIT11 → projection analysis enabled [WARN] KPOINTS: Monkhorst-Pack grid 8x8x8, but no KPOINTS file found → using default [ERROR] OUTCAR: missing! Cannot parse eigenvalues or forces关键看三处[INFO] INCAR行确认LORBIT是否为11或12投影能带必需[WARN] KPOINTS若提示缺失说明你没生成KPOINTS文件vaspkit -task 111可补[ERROR] OUTCAR缺失意味着你根本没跑过VASP此时所有分析任务都会失败。2.3 版本锁死为什么vaspkit 3.5.2必须配VASP 6.3vaspkit对VASP输出格式极其敏感。例如VASP 5.4.4的OUTCAR中E-fermi字段在第12行而VASP 6.3.2移到第15行vaspkit 3.3.x默认按旧格式解析遇到新OUTCAR会读错费米能级导致能带图整体偏移0.8 eV解决方案不是降级VASP不现实而是升级vaspkitwget https://github.com/vaspkit/vaspkit/releases/download/v3.5.2/vaspkit-3.5.2-linux-x64.tar.gz tar -xzf vaspkit-3.5.2-linux-x64.tar.gz。血泪经验某高校计算中心集群预装vaspkit 3.2.1学生批量提交能带任务后所有图谱的费米能级标在-2.1 eV实际应为-0.3 eV。排查3天才发现是版本错配——永远用vaspkit -v确认版本再查对应VASP兼容表。3. 能带与投影能带从-task 211到-task 212的5步闭环能带图是第一性原理计算的“门面”但vaspkit的能带任务常因K点路径、自旋设置、投影权重三者不匹配而翻车。本节以体心立方Fe磁性体系为例拆解从原始输出到可投稿图的完整链路。3.1 用-task 111生成高对称K路径拒绝手写KPOINTS手写KPOINTS易错如Γ-X-W-K-Γ顺序写反、分数坐标算错。vaspkit -task 111自动识别POSCAR空间群并生成标准路径vaspkit -task 111 # 输入选项1 → 自动检测空间群Fe为Im-3m, No.229 # 输入选项2 → 选择High-symmetry path for band structure # 输出KPOINTS_band含120个k点和KPATH.in路径标签生成的KPOINTS_band是标准Monkhorst-Pack格式但必须确保其k点数≥VASP静态计算时的KPOINTS总数否则vaspkit -task 211会报Number of k-points mismatch。检查方法grep K-Points OUTCAR | head -1 # 查VASP实际用的k点数 wc -l KPOINTS_band | awk {print $1-2} # 查KPOINTS_band有效k点数减去前2行注释3.2vaspkit -task 211跑能带3个必调参数-task 211本身不产生图只生成BAND.dat能量数据和KPATH.datk点坐标。但以下参数决定结果可靠性vaspkit -task 211 -emin -10 -emax 10 -spin 1-emin/-emax设置绘图能量范围单位eV。默认±5 eV常切掉导带底-10~10更安全-spin 1对ISPIN2的磁性体系指定画自旋向上通道-spin 2为向下隐藏参数-nk 200强制插值k点数。原始KPOINTS_band仅120点插值到200点后能带曲线更平滑尤其在能隙处。逻辑说明-task 211读取EIGENVAL需LORBIT11生成、OUTCAR取费米能级、KPATH.in映射k点标签。若EIGENVAL缺失它会静默失败并输出空BAND.dat——务必用wc -l BAND.dat确认行数1000。3.3vaspkit -task 212加投影原子轨道权重的3种合法写法投影能带pDOS-like能带需-task 212但它对PROCAR文件极度挑剔# 正确1指定原子序号POSCAR中Fe为第1、2号原子 vaspkit -task 212 -atom 1,2 -orb d # 正确2指定元素名需PROCAR含元素标识 vaspkit -task 212 -element Fe -orb d # 正确3指定原子类型POSCAR中Fe占第1类 vaspkit -task 212 -type 1 -orb d-orb可选s,p,d,f,px,py,pz,dxy,dyz...但必须与PROCAR中orbital列一致若PROCAR无orbital列旧版VASP用-orb all强制全轨道投影报错Error: Cannot find orbital d in PROCAR时先head -20 PROCAR看第6列内容。3.4 用Python脚本一键出图告别gnuplot玄学vaspkit生成的BAND.dat是纯数据需绘图。以下脚本band_plot.py直接输出PDFimport numpy as np import matplotlib.pyplot as plt # 读取vaspkit输出 kpath np.loadtxt(KPATH.dat) # k点坐标 band np.loadtxt(BAND.dat) # 第一列k索引第二列能量 proj np.loadtxt(PROJ_BAND.dat) # 投影权重与band同尺寸 # 绘图 fig, ax plt.subplots(figsize(6, 5)) for i in range(1, band.shape[1]): # 跳过第一列k索引 ax.plot(band[:, 0], band[:, i], b-, linewidth0.8) # 叠加投影颜色深浅表示权重 scatter ax.scatter(band[:, 0], band[:, i], cproj[:, i], cmapReds, s8, alpha0.7) # 设置k点标签从KPATH.in读取 with open(KPATH.in) as f: labels [line.split()[0] for line in f if line.strip() and not line.startswith(#)] ax.set_xticks([0, 50, 100, 150]) # 根据KPATH.dat长度调整 ax.set_xticklabels(labels) plt.ylabel(Energy (eV)) plt.axhline(y0, colork, linestyle--, linewidth0.8) # 费米能级 plt.savefig(band_fe.pdf, bbox_inchestight)参数说明s8控制散点大小alpha0.7避免重叠点过暗cmapReds使高投影区域显红色——这比vaspkit自带的gnuplot脚本更可控。4. 避坑能带与投影分析的5个高频翻车现场这些坑我在3个不同实验室的模拟项目X中反复见过轻则图错重则结论反转。4.1 现象BAND.dat全是NaNvaspkit -task 211无报错原因INCAR中ICHARG1从原子电荷猜初态但未提供WAVECAR导致VASP未收敛即停vaspkit读EIGENVAL时发现空文件填NaN而不报错。解决运行vaspkit -task 0检查[ERROR] EIGENVAL: missing!然后重跑VASP确保ICHARG11从WAVECAR读初态且WAVECAR存在。4.2 现象投影能带颜色块断裂部分k点无颜色原因PROCAR文件被截断如计算中途被killvaspkit -task 212仍尝试解析但缺失k点的orbital数据。解决用wc -l PROCAR对比grep k-point OUTCAR | wc -l若前者小10%说明PROCAR损坏删PROCAR重跑VASP加LORBIT11。4.3 现象能带图k点标签错位Γ标在X位置原因vaspkit -task 111生成的KPATH.in与KPOINTS_band不匹配。常见于手动修改KPOINTS_band后未更新KPATH.in。解决删除KPATH.in重新运行vaspkit -task 111或用sed -i s/0\.00000000/0.0/g KPOINTS_band修复浮点精度问题某些编译器输出0.00000000导致解析失败。4.4 现象自旋极化体系能带上下通道完全重合原因-task 211未加-spin 1/2默认只画自旋向上或INCAR中ISPIN1非磁性计算却强行投投影。解决先grep ISPIN INCAR确认为2再明确指定-spin 1或-spin 2若需双通道分两次运行并合并BAND.dat。4.5 现象投影权重最大值远超1.0如2.5原因-task 212默认归一化到单原子但POSCAR中同一元素有多个原子如Fe1, Fe2权重叠加导致超限。解决加-norm 1参数强制全局归一化vaspkit -task 212 -element Fe -orb d -norm 1。5. 电荷分析进阶从差分电荷到Bader电荷的无缝衔接差分电荷Δρρ_total−ρ_atomic和Bader电荷是理解成键本质的核心但vaspkit的-task 301差分电荷与-task 302Bader之间存在一个关键断层Bader分析要求CHGCAR格式严格匹配而vaspkit生成的CHGDIFF.dat需手动转CHGCAR。本节给出零误差转换方案。5.1vaspkit -task 301生成差分电荷2个隐藏开关vaspkit -task 301 -grid 48,48,48 -format 1-grid 48,48,48强制插值网格为48×48×48VASP默认可能为24×24×24Bader要求高分辨率-format 1输出CHGDIFF.dat为VASP格式非cube这是Bader工具链唯一接受的输入。为什么必须设-gridBader算法对电荷密度梯度敏感低分辨率网格会导致电荷分区错误。实测FeO体系中24×24×24网格给出O原子Bader电荷为-1.12|e|48×48×48网格修正为-1.38|e|与实验值-1.41|e|更近。5.2 CHGDIFF.dat → CHGCAR用awk一行转格式vaspkit的CHGDIFF.dat是ASCII格式含header和3D数据块需转为CHGCAR标准格式含POSCAR头电荷密度数据# 提取POSCAR头前8行和CHGDIFF.dat数据跳过前10行header (head -8 POSCAR; tail -n 11 CHGDIFF.dat) CHGCAR_diff但此CHGCAR_diff缺少关键信息电荷密度总电子数。Bader工具bader会因总电子数为0而崩溃。正确做法# 1. 用vaspkit提取总电子数从OUTCAR grep number of electron OUTCAR | tail -1 | awk {print $5} nelec.txt # 2. 将电子数插入CHGCAR_diff第8行后 sed -i 8a\ $(cat nelec.txt) CHGCAR_diff # 3. 补充空行和数据标记CHGCAR规范要求 echo -e \n\n CHGCAR_diff echo augmentation occupancies CHGCAR_diff5.3 运行Bader并验证3步防错# 下载baderhttps://theory.cm.utexas.edu/henkelman/code/bader/ ./bader CHGCAR_diff -c CHGCAR -p ./AECCAR0 -p ./AECCAR2 # 检查输出ACF.dat含原子电荷XSF.dat含可视化数据 # 验证1总电荷守恒 awk {sum$5} END{print Total charge:, sum} ACF.dat # 验证2与vaspkit内置Bader对比-task 302 vaspkit -task 302 -chg CHGCAR_diff # 若两者O原子电荷差0.05|e|说明CHGCAR_diff格式有误关键技巧Bader对AECCAR0/AECCAR2补偿电荷极其敏感。若未生成用vaspkit -task 303补全vaspkit -task 303 -chg CHGCAR -aeccar 0,2。永远用vaspkit -task 0确认AECCAR0/AECCAR2存在后再跑Bader。6. 声子初筛与热力学量用-task 401快速定位软模而不跑PHONOPY全声子谱计算PHONOPYVASP耗时以天计但很多项目只需判断是否存在50 cm⁻¹的软模预示相变。vaspkit -task 401基于有限差分线性响应10分钟内完成初筛是真正的“后悔药”。6.1vaspkit -task 401的3个输入前提它不直接读取VASP输出而是要求你提供POSCAR_unit原胞POSCAR非超胞FORCE_CONSTANTS由VASP的IBRION5/6生成需先跑一次微扰qpoints.inq点列表如0.0 0.0 0.0用于Γ点。生成FORCE_CONSTANTS的最小流程# 1. 在POSCAR_unit目录下用vaspkit生成2×2×2超胞 vaspkit -task 101 -supercell 2,2,2 # 2. 修改INCARIBRION6, POTIM0.015, NSW1 # 3. 运行VASP输出FORCE_CONSTANTS需编译时加DFPT支持6.2 执行初筛并解读输出vaspkit -task 401 -qfile qpoints.in输出PHONON_FREQ.dat含# q-point: 0.0000 0.0000 0.0000 # Mode Freq(cm-1) Freq(meV) Type 1 0.23 0.03 Acoustic 2 0.18 0.02 Acoustic 3 0.05 0.01 Acoustic 4 152.34 18.92 Optical重点看前三列Acoustic模式频率应≈0≤1 cm⁻¹正常若5 cm⁻¹说明超胞不够大或收敛不足Optical模式中任何Freq(cm⁻¹)50即为软模预警如32.4 cm⁻¹需用PHONOPY精算若所有Optical模式100 cm⁻¹基本排除低温相变可能。6.3 从声子到热力学vaspkit -task 402生成热容/自由能有了PHONON_FREQ.dat-task 402可算0–1000 K热力学量vaspkit -task 402 -tmin 0 -tmax 1000 -tstep 10输出THERMO.dat含温度、自由能、熵、热容。关键验证点0 K时自由能应≈总能grep energy without entropy OUTCAR若0 K自由能比总能高0.01 eV说明PHONON_FREQ.dat有负频未剔除用-freq_min 5过滤。我的习惯是所有新材料计算必跑vaspkit -task 401初筛若发现软模立即暂停后续分析先用PHONOPY重构晶胞——省下3天无效计算比优化1个参数更重要。希望帮到你。本文还有配套的精品资源点击获取