企业官网建设流程全解析

简介：本资源为《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配置文件：

# 编辑 ~/.bashrc（Linux）或 ~/.zshrc（macOS） 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, volume=15.23 ų [INFO] INCAR: ISPIN=2, ICHARG=11, LORBIT=11 → 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.3+

vaspkit对VASP输出格式极其敏感。例如：

• VASP 5.4.4的OUTCAR中E-fermi字段在第12行，而VASP 6.3.2移到第15行；
• vaspkit 3.3.x默认按旧格式解析，遇到新OUTCAR会读错费米能级，导致能带图整体偏移0.8 eV；
• 解决方案不是降级VASP（不现实），而是升级vaspkit：wget 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.dat（k点坐标）。但以下参数决定结果可靠性：

vaspkit -task 211 -emin -10 -emax 10 -spin 1

• -emin/-emax：设置绘图能量范围（单位eV）。默认±5 eV常切掉导带底，-10~10更安全；
• -spin 1：对ISPIN=2的磁性体系，指定画自旋向上通道（-spin 2为向下）；
• 隐藏参数-nk 200：强制插值k点数。原始KPOINTS_band仅120点，插值到200点后能带曲线更平滑（尤其在能隙处）。

逻辑说明：-task 211读取EIGENVAL（需LORBIT=11生成）、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）直接输出PDF：

import 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-', linewidth=0.8) # 叠加投影（颜色深浅表示权重） scatter = ax.scatter(band[:, 0], band[:, i], c=proj[:, i], cmap='Reds', s=8, alpha=0.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(y=0, color='k', linestyle='--', linewidth=0.8) # 费米能级 plt.savefig("band_fe.pdf", bbox_inches='tight')

参数说明：s=8控制散点大小，alpha=0.7避免重叠点过暗，cmap='Reds'使高投影区域显红色——这比vaspkit自带的gnuplot脚本更可控。

4. 避坑：能带与投影分析的5个高频翻车现场

这些坑我在3个不同实验室的模拟项目X中反复见过，轻则图错，重则结论反转。

4.1 现象：BAND.dat全是NaN，vaspkit -task 211无报错

原因：INCAR中ICHARG=1（从原子电荷猜初态），但未提供WAVECAR，导致VASP未收敛即停；vaspkit读EIGENVAL时发现空文件，填NaN而不报错。
解决：运行vaspkit -task 0检查[ERROR] EIGENVAL: missing!，然后重跑VASP，确保ICHARG=11（从WAVECAR读初态）且WAVECAR存在。

4.2 现象：投影能带颜色块断裂，部分k点无颜色

原因：PROCAR文件被截断（如计算中途被kill），vaspkit -task 212仍尝试解析，但缺失k点的orbital数据。
解决：用wc -l PROCAR对比grep "k-point" OUTCAR | wc -l，若前者小10%，说明PROCAR损坏，删PROCAR重跑VASP（加LORBIT=11）。

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中ISPIN=1（非磁性计算）却强行投投影。
解决：先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 302（Bader）之间存在一个关键断层：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×48（VASP默认可能为24×24×24，Bader要求高分辨率）；
• -format 1：输出CHGDIFF.dat为VASP格式（非cube），这是Bader工具链唯一接受的输入。

为什么必须设-grid？Bader算法对电荷密度梯度敏感，低分辨率网格会导致电荷分区错误。实测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_diff

5.3 运行Bader并验证：3步防错

# 下载bader（https://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

全声子谱计算（PHONOPY+VASP）耗时以天计，但很多项目只需判断是否存在<50 cm⁻¹的软模（预示相变）。vaspkit -task 401基于有限差分+线性响应，10分钟内完成初筛，是真正的“后悔药”。

6.1vaspkit -task 401的3个输入前提

它不直接读取VASP输出，而是要求你提供：

• POSCAR_unit：原胞POSCAR（非超胞）；
• FORCE_CONSTANTS：由VASP的IBRION=5/6生成（需先跑一次微扰）；
• qpoints.in：q点列表（如0.0 0.0 0.0用于Γ点）。

生成FORCE_CONSTANTS的最小流程：

# 1. 在POSCAR_unit目录下，用vaspkit生成2×2×2超胞 vaspkit -task 101 -supercell 2,2,2 # 2. 修改INCAR：IBRION=6, POTIM=0.015, NSW=1 # 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个参数更重要。希望帮到你。

本文还有配套的精品资源，点击获取