简介一份面向VASP使用者的中文功能手册系统讲解VASPKIT在材料科学、催化计算中的预处理与后处理应用适合从入门到进阶的科研人员。手册内容按菜单功能组织涵盖VASPKIT简介与配置、VASP输入文件生成和检查、PBE与HSE06杂化泛函能带计算、分波态密度与局域态密度、热力学量校正、差分电荷密度绘制以及NEB路径生成与可视化、虚频校正、能带反折叠、3D能带、费米面、VASP2BoltzTraP接口等催化与电子结构分析模块。此外还介绍了用户自定义界面可用于绘制原子电荷着色结构图、二维材料投影能带等扩展场景并附有引用方式与菜单功能附录。资源共1个文件为PDF格式包体仅2.37MB轻量易保存章节目录清晰便于按需定位。已有3564人学习下载。借助书本示例与逐步说明读者可快速掌握从输入文件准备、电子结构分析到催化反应路径处理的关键操作提升日常VASP计算效率。1. 从 GLIBC 报错说起先把 VASPKIT 0.72 装进 VASP 工作流第一次跑 VASPKIT 0.72 时我差点被一行/lib64/libc.so.6: version GLIBC_2.14 not found劝退。当时只是想用它的能带反折叠功能处理一个超胞结构的投影能带结果./vaspkit敲下去屏幕上先蹦出来的是库版本报错。查了一圈才明白工具没问题是运行环境的 glibc 太老。VASPKIT 是 FORTRAN 写成、在 LINUX 下运行的纯二进制程序体积只有 5M几乎不依赖第三方库但对系统 glibc 版本有硬要求。这个名为 VASPKIT 的 VASP 前处理/后处理工具箱在材料计算和催化计算里早就成了标配——从自动生成 INCAR、POTCAR、KPOINTS 并做格式检查到能带、态密度、差分电荷密度、NEB 路径和虚频校正全在一个交互菜单里完成。适合两类人刚入坑 VASP、不想满世界找脚本的新手以及在各种第三方脚本版本打架里消耗过大量时间的熟手。2. 环境准备下载 0.72 免安装包先看 bin 目录和 EXAMPLES2.1 版本怎么选为什么是 0.72装任何工具之前都要先定版本VASPKIT 也不例外。0.72 是目前功能覆盖最全的 release 版本这个版本的菜单做了一次大调整把功能相似的选项归类成“结构 Options”“电子 Options”“电荷与势与波函数 Options”“Misc Utilities”四个大组同时补上了很多催化相关的新功能。我特别关注的是几个新增任务task 28 加入了 DFT、hybrid-DFT 和 Meta-GGA 级别的能带反折叠Band-Structure Unfoldingtask 25 支持在杂化泛函能带计算里沿不同 k 路径做均匀 k 点采样task 505 可以线性插点生成 NEB 路径task 507 可以做虚频校正还有 task 112 到 115 这组投影态密度功能能一次处理单原子、单元素或指定原子轨道的投影态密度和投影能带。这些功能对做催化、做半导体材料的人来说非常实用。比如 NEB 路径可视化以前要把多个中间态结构导成 PDB 文件得自己写脚本处理 OUTCAR 里的坐标现在 task 405 直接转再比如 HSE06 杂化泛函能带老版本处理投影能带数据时会有 bug0.72 里专门修了提取 hybrid-DFT 级别投影能带数据的问题。所以如果你还在用 0.7x 早期版本建议直接换 0.72菜单结构变了功能入口和参数模板都按新版本对齐。另外要注意VASPKIT 是纯二进制分发没有图形界面安装向导。它几乎不依赖第三方库但依赖系统运行时库尤其是 glibc。装之前先确认你的发行版不太老ldd --version | head -n1常见的 CentOS 7、Ubuntu 16.04 以上都能满足 0.72 的运行要求。如果输出显示 glibc 版本低于 2.14后面运行 vaspkit 就大概率会遇到第 4 章里写的 GLIBC_2.14 not found 报错这点提前有个心理准备。2.2 解压、放置、赋权三步把二进制安顿好拿到 vaspkit.0.72 的压缩包后我一般先解压到临时目录再整体移动到 /home 或 /opt 下统一管理避免直接散落在下载目录里tar -zxvf vaspkit.0.72.tar.gz mv vaspkit.0.72 /home/ cd /home/vaspkit.0.72 ls -l解压后你会看到几个关键内容bin 目录里是 vaspkit 可执行文件how_to_set_environment_variable 是环境变量配置模板EXAMPLES 目录下是各主要功能的测试例子。建议先看一眼文件的权限ls -l /home/vaspkit.0.72/bin/vaspkit如果权限位里没有 x后面直接运行会报 Permission denied。常规操作是先赋执行权限再测试是否可用chmod ux /home/vaspkit.0.72/bin/vaspkit /home/vaspkit.0.72/bin/vaspkit -h-vaspkit -h 能打印出命令行帮助说明这个二进制在你的系统上能正常加载动态库。若这里就报错基本可以断定是 glibc 版本问题直接跳到第 4 章的排查方案。2.3 EXAMPLES 目录比任何文档都快的上手途径VASPKIT 的 EXAMPLES 目录里放了很多功能的测试例子这是我建议新手最先翻的地方。比如 examples/hybrid_DFT_band 目录下有一个用 HSE06 杂化泛函计算磷化镓能带结构的完整案例从 POSCAR 到最终能带数据都有examples/band_unfolding 目录下则是能带反折叠的多个例子对应 task 28 的功能。每个例子目录里基本都备好了输入文件和说明你可以直接在里面运行 vaspkit 看某个功能怎么工作不需要自己从头搭结构。我自己的习惯是装完一个新工具先跑一遍自带 examples不是为了算科学结果而是为了确认“工具本身能跑通”。如果 examples 里的案例跑不出合理数据那多半是配置有问题而不是你手头的研究体系有问题。这样排查范围就一下缩小了。3. 配置 ~/.vaspkit赝势路径、INCAR 写入模式与 PATH 三件事3.1 环境变量文件从哪里来VASPKIT 运行时需要读一个用户级配置文件默认位置是~/.vaspkit。首次使用前要把软件包自带的模板复制过去cp -f how_to_set_environment_variable ~/.vaspkit vi ~/.vaspkit这个文件不复制的话VASPKIT 也能启动但很多功能会因为找不到赝势路径而报错尤其是生成 POTCAR 的功能。复制时要特别注意-f参数它会强制覆盖已存在的 ~/.vaspkit。如果你之前配过旧版本先备份再复制cp ~/.vaspkit ~/.vaspkit.bak3.2 配置项逐条说明打开 ~/.vaspkit 后核心配置项大致如下不同版本字段略有差异0.72 是以这组为准VASP5 .TRUE. # .TRUE. or .FALSE.; Set .FALSE. if using vasp.4.x GGA_PATH ~/POTCAR/GGA # Path of GGA potential PBE_PATH ~/POTCAR/PBE # Path of PBE potential LDA_PATH ~/POTCAR/LDA # Path of LDA potential POTCAR_TYPE PBE # PBE, GGA or LDA GW_POTCAR .FALSE. # .TRUE. or .FALSE. RECOMMENDED_POTCAR .TRUE. # .TRUE. or .FALSE. MINI_INCAR .FALSE. # .TRUE. or .FALSE. SET_FERMI_ENERGY_ZERO .TRUE. # .TRUE. or .FALSE. SET_INCAR_WRITE_MODE OVERRIDE # OVERRIDE, APPEND, BACK-UP-OLD, BACK-UP-NEW这里我按实际使用经验拆开说参数作用我的建议VASP5是否使用 VASP 5.x 语法的 POSCAR/POTCAR现在绝大多数人用 VASP 5.4保持 .TRUE.只有老项目还在用 4.x 才改 .FALSE.GGA_PATH / PBE_PATH / LDA_PATH三种赝势库的根目录一定改成你机器上真实存在的路径建议用绝对路径POTCAR_TYPE生成 POTCAR 时默认从哪个赝势库取默认 PBE 即可别随手改成 GGAGW_POTCAR是否使用 GW 专用赝势不跑 GW 计算就 .FALSE.RECOMMENDED_POTCAR是否按 VASP 官网推荐自动选赝势新手保持 .TRUE.需要特殊后缀时再手动干预MINI_INCAR是否生成精简版 INCAR建议 .FALSE.让模板把该有的参数都写全SET_FERMI_ENERGY_ZERO处理能带/态密度时是否把费米能级平移到零保持 .TRUE.否则画图时费米能级位置对不齐SET_INCAR_WRITE_MODE已有 INCAR 时怎么处理默认 OVERRIDE 直接覆盖想留旧文件就设 BACK-UP-OLD这里最容易被忽略的是 PATH 类配置项里的~。在单引号里~不会展开如果 VASPKIT 读取时不做展开处理就会拿着字符串~/POTCAR/PBE去拼路径结果就是找不到赝势。我一般直接写成/home/username/POTCAR/PBE这种绝对路径省得在不同 shell、不同登录方式下出现“刚才还能用换个终端就找不到赝势”的玄学问题。3.3 把 wanjaspkit 加进 PATH配置文件弄好后为了方便在任意目录启动把 bin 目录加进 PATHecho export PATH/home/vaspkit.0.72/bin/:$PATH ~/.bashrc source ~/.bashrc如果你用的是 zsh就写到 ~/.zshrc别只记 bashrc。加完之后在终端直接输入vaspkit就能启动也可以继续用绝对路径/home/vaspkit.0.72/bin/vaspkit运行。这样做的意义是后面在 examples 里切换目录测试时不用每次找路径。3.4 启动界面长什么样配置正确后启动 vaspkit 会看到一个用 ASCII 字符画的菜单界面。很多新手看到这个“很萌的界面”会愣一下不确定有没有启动成功。其实看到下面这个主菜单就说明启动正常 Structural Options 1) VASP Input Files Generator 2) Elastic-Properties Kit 3) K-Path Generator 4) Structure Editor 5) Catalysis-ElectroChemi Kit 6) Symmetry Search Electronic Options 11) Density-of-States 21) DFT Band-Structure 23) 3D Band-Structure 25) Hybrid-DFT Band-Structure 26) Fermi-Surface 28) Band-Structure Unfolding Charge Potential Wavefunction Options 31) Charge Spin Density 42) Potential-Related 51) Wave-Function Analysis Misc Utilities 71) Linear Optics 72) Molecular-Dynamics Kit 73) VASP2BoltzTraP Interface 91) Semiconductor Calculator 92) 2D-Materials Kit 0) Quit菜单按数字分组输入对应编号回车进入子菜单。这里我提醒一句VASPKIT 是命令行交互式操作不支持鼠标点选。如果你是在远程服务器上通过 SSH 用完全没问题终端字符界面跑得稳。4. 避坑手册装完 VASPKIT 最常见的五个问题4.1 vaspkit: /lib64/libc.so.6: version GLIBC_2.14 not found这是我遇到最多的启动报错现象就是第 1 章那段描述执行 vaspkit 后系统直接抛出 “version GLIBC_2.14 not found (required by vaspkit)”。原因是二进制文件在编译时链接了较新版本的 glibc而运行机器上的 glibc 低于 2.14。常见于老旧的 CentOS 6.x 或同类发行版。解决有两个方向。第一个是升级系统 glibc但不建议因为 glibc 是系统最底层的基础库贸然升级可能连带破坏系统里其他依赖旧版本的应用属于“为了一个工具掀翻整个系统”的做法。第二个是找开发者要一份在低版本 LINUX 环境下编译的版本这个更稳妥相当于让二进制来适配你的系统而不是让系统去冒险适配它。如果不想走这一步也可以考虑换一台 glibc 版本更高的机器比如 CentOS 7 以上或 Ubuntu 16.04 以上直接绕开这个问题。4.2 Permission denied可执行权限没给现象是在 bash 里执行/home/vaspkit.0.72/bin/vaspkit或直接敲vaspkit提示-bash: line 7: ./vaspkit: Permission denied。原因是压缩包解压后文件默认没有执行权限。很多新手会先去查 PATH 配置、查环境变量其实问题根本没到那一步就是最简单的权限位没有 x。解决是赋权chmod ux /home/vaspkit.0.72/bin/vaspkit我建议解压后立刻执行这一步不要等报错再处理。赋权之后运行vaspkit -h如果能看到帮助信息说明二进制本身没问题后续再报错就集中在配置层面。4.3 POSCAR 里藏着制表符VASP 只读出了第一个元素这是一个提交 VASP 任务后才会暴露的坑。现象是任务秒失败日志里写着ERROR: the type information is not consistent with the number of typesVASP 只读出了 POSCAR 里的第一种元素。原因是从某些晶体数据库下载的 POSCAR 文件里分隔符用了制表符\t而不是 VASP 期望的空格。VASP 在读字符串时不能正确处理\t于是元素类型列表没有完整解析。这个问题在 Windows 系统下编辑生成的 POSCAR 或 INCAR 里也会出现表现形式可能是各种奇怪的读入错误。解决是用 VASPKIT 的功能 109 做格式纠正。运行 vaspkit进入功能 1VASP Input Files Generator再选 109程序会自动纠正 INCAR 和 POSCAR 的格式问题并检查赝势元素是否一致。执行完之后重新提交任务就能正常跑起来。如果你不想用 VASPKIT也可以用 sed 把制表符批量替换成空格但最省事的还是 109一步到位。4.4 POSCAR 元素类型和 POTCAR 不一致不报错但结果全错这个坑最隐蔽VASP 不会自动检查 POSCAR 里的元素类型是否与 POTCAR 一致。意味着你拿一份石墨烯的 POSCAR配一个 H 的 POTCARVASP 照样能跑但它算出来的结果是错的而且不会给你任何报错提示。原因是 VASP 的设计里POSCAR 只提供晶格和坐标POTCAR 提供势函数二者是否匹配全靠用户自觉。很多新手在做“替换元素”操作时会顺手把 POSCAR 改了而忘记重新生成 POTCAR或者反过来导致整个计算的物理意义就是错的。解决是靠 vaspkit 的 109 功能做检查。它会读取 POSCAR 里的元素顺序和 POTCAR 里的元素逐一比对发现不一致就明确提示。提交任务前跑一遍 109比自己盯着一行行核对靠谱得多。这一条我已经写进了自己的固定流程拿到任何新结构都先 109 再交 VASP。4.5 赝势路径写~换个终端就找不到 POTCAR现象是这个终端里 vaspkit 能正常生成 POTCAR换个终端或者换一种登录方式再运行就提示找不到指定赝势目录。原因是配置文件里写的是~/POTCAR/PBE这类带波浪号的路径。~在 shell 里会展开成家目录但 VASPKIT 是 FORTRAN 程序它读取配置文件时未必会做同样的展开处理。不同环境下表现不一致于是就成了“时好时坏”的玄学问题。解决是打开 ~/.vaspkit把 GGA_PATH、PBE_PATH、LDA_PATH 全部改成绝对路径例如/home/username/POTCAR/PBE。改完保存重启 vaspkit问题消失。这条经验我后来写进了自己的部署笔记所有配置文件里涉及路径的参数一律用绝对路径不给~留发挥空间。5. 生成输入文件与能带计算把 ZnO 从 POSCAR 带到能带反折叠5.1 第一步准备一份能用的 POSCAR配置好环境、看完避坑之后就该用真实体系验证 vaspkit 的完整流程了。我一般拿纤锌矿 ZnO 做测试因为这个结构简单、空间群明确而且能带和态密度的参考数据容易找。先从公开晶体数据库里检索 ZnO空间群选 P63mc 的六方相下载 POSCAR 文件典型文件名类似ZnO-1811117.vasp。把它复制成 VASP 标准文件名cp -f ZnO-1811117.vasp POSCAR如果你拿到的是 CIF 文件需要先用 VESTA 或 VTST 脚本转成 POSCAR注意原子位置的分数坐标是否处理正确。这一步不能省后续 vaspkit 生成的 KPOINTS、POTCAR 都以 POSCAR 里的元素和晶格为准。5.2 功能 101按任务类型生成 INCAR在存放 POSCAR 的目录下运行 vaspkit输入1进入 VASP Input Files Generator会看到功能子菜单101) Customize INCAR File 102) Generate KPOINTS File for SCF Calculation 103) Generate POTCAR File with Default Setting 104) Generate POTCAR File with User Specified Potential 107) Reformat POSCAR File by Specified Order of Elements 108) Successive Procedure to Generate VASP Files and Check 109) Check All VASP Files输入101会进入 INCAR 模板选择界面里面列出了针对不同计算任务的模板ST) Static-Calculation SR) Standard Relaxation MG) Magnetic Properties SO) Spin-Orbit Coupling D3) DFT-D3 no-damping Correction H6) HSE06 Calculation PU) DFTU Calculation MD) Molecular Dynamics GW) GW0 Calculation BS) BSE Calculation DC) Elastic Constant EL) ELF Calculation BD) Bader Charge Analysis OP) Optical Properties EC) Static Dielectric Constant PC) Decomposed Charge Density FD) Phonon-Finite-Displacement DT) Phonon-DFPT NE) Nudged Elastic Band (NEB) DM) The Dimer Method FQ) Frequence Calculations LR) Lattice Relaxation输入LR就会生成一个用于晶格弛豫的 INCAR。新手看到这么多模板可能会慌其实这里的选择逻辑很简单你要算什么任务就选对应的首字母组合。比如跑 Bader 电荷分析就选 BD跑 NEB 就选 NE跑 HSE06 静态计算就选 H6。每个模板生成后部分参数还需要手动调整比如 DFTU 的 U 值、NEB 的 IMAGES 数目vaspkit 会在生成时用一句 “You MUST Know What You Are Doing” 提示你别无脑跑。这里关于 INCAR 里的 ENCUT我想多说一句。默认生成的 INCAR 里 ENCUT400 往往是注释掉的程序保留的是PREC NormalVASP 会自动把 ENCUT 设为所有元素 ENMAX 的最大值。如果你用 LR 模板优化晶格常数vaspkit 会生成PREC HighVASP 此时会把 ENCUT 自动设为 max(ENMAX)*1.3这样对晶格常数优化的收敛更有利。不过较新的官方建议已经不推荐用 PRECHigh而是设 PRECAccurate 并手动定 ENCUT你可以在模板基础上改成这种写法确保不同体系间 ENCUT 一致。5.3 功能 102KPOINTS 与 POTCAR 一次生成INCAR 就绪后回到子菜单输入102生成 SCF 计算的 KPOINTS。程序会先提示你选择撒点方式输入2选择 Gamma Scheme这是大多数体系下最稳妥的选择。接着提示输入 K 点密度精度值* Accuracy Levels: (1) Low: 0.06~0.04; (2) Medium: 0.04~0.03; (3) Fine: 0.02-0.01; (4) Gamma-Only: 0. Input KP-Resolved Value (unit: 2*PI/Angstrom):输入0.04一般是够用的。0.04 这个值的意思是倒空间每个 K 点间隔约 0.04 × 2π/Åvaspkit 会根据你 POSCAR 里的倒空间晶格矢量比例自动算出各个方向的 K 点数。生成结果看起来类似K-Mesh Generated with KP-Resolved Value : 0.0400 Gamma 9 9 5 0.0 0.0 0.0这里 9×9×5 是程序根据你晶格倒空间比例自动分配的不是固定值。对于正交体系有个简单的判断准则实空间晶格矢量乘以对应方向的 K 点数约大于 20Å 就比较稳妥ZnO 这个例子里 ka≈29.62Å满足经验规律。非正交体系因为倒格矢长度和晶格常数不满足简单的反比关系用 vaspkit 自动分配反而更合理。生成 KPOINTS 的同时vaspkit 会根据 POSCAR 里的元素类型从你配置的赝势库里提取并组合生成 POTCAR。前提是 ~/.vaspkit 里的 PBE_PATH 指向正确目录POTCAR_TYPE 保持 PBE。生成成功后用以下命令确认grep TIT POTCAR输出会列出 POTCAR 里实际包含的元素赝势。如果显示的是 O 和 Zn说明 POTCAR 和 POSCAR 的元素顺序是对上的。5.4 功能 104手动指定赝势类型比如 Zn_pv默认的 103 功能是按 VASP 官网推荐自动选赝势。但有些场景需要手动指定比如做 Bader 电荷分析时建议使用带半芯态的赝势。PBE 赝势库里常见后缀含义如下后缀含义无后缀只把价电子当作价态处理_pv次外层的 p 半芯态也当价电子_sv次外层的 s 半芯态也当价电子_d次外层的 d 半芯态也当价电子_h针对短键二聚体的硬赝势如 O2、CO、N2、F2 等常见的做法是先测试两种赝势对自己关心的性质影响大不大。如果影响不明显优先用无后缀版本因为包含更多价电子意味着截断能上升、计算量明显增大。如果体系里有 O2 这类短键分子官方建议用 _h 势。手动干预时用功能 104按提示依次输入元素和赝势类型比如给 Zn 输入Zn_pv。如果指定目录下没有对应赝势文件程序会提示重新输入。生成后同样用grep TIT POTCAR验证元素顺序和类型。5.5 功能 109/108提交 VASP 前最后一道检查如果这时直接把任务交出去你可能会遇到我在第 4 章描述过的场景从数据库下载的 POSCAR 里藏着制表符VASP 只读出了 O 没有读出 Zn任务秒失败。vaspkit 从 0.71 开始加入了功能 109专门做 INCAR 和 POSCAR 的格式纠正、检查赝势元素一致性。输入109程序会告诉你所有必要文件是否存在并比对 POSCAR 和 POTCAR 的元素类型是否一致。这个检查非常关键因为 VASP 本身不检查你拿石墨烯配 H 的赝势它也照跑不误只是结果毫无意义。所以更省心的做法是直接用功能 108依次完成 INCAR 设置、KPOINTS 生成、POTCAR 生成最后自动调用 109 检查所有文件。整套流程一句话概括就是“VASP 之前先 KIT 一下”。我后来所有新体系的第一次计算都走这个流程没有再被制表符和元素不匹配坑过。5.6 能带路径与能带计算功能 3、21、25、28能带计算的 KPOINTS 和自洽计算不一样需要沿着布里渊区高对称点路径取点。传统做法是去 SeeK-Path 网站或者用 Materials Studio 找高对称点再写脚本插值。vaspkit 把这一步集成了输入3进入 K-Path Generator选择 3D bulk structure程序会分析晶体对称性并给出建议路径。对纤锌矿 ZnOvaspkit 会生成类似GAMMA-X-U|K-GAMMA-L-W-X的两条路径同时生成两个文件KPATH.in和PRIMCELL.vasp。KPATH.in 的第二行是每段路径插值的 K 点数目如果计算量紧张可以调小PRIMCELL.vasp 是程序生成的单胞能带路径只针对 primitive cell所以要把这个文件作为后续计算的结构cp -f PRIMCELL.vasp POSCAR接下来按计算层级选择后处理功能。普通 PBE 泛函能带用功能 21HSE06 杂化泛函能带用功能 25想直观展示三维能带走势可以用功能 23超胞或者表面体系想看原始晶格能带投影则用功能 28 做能带反折叠。0.72 的能带反折叠支持 DFT、hybrid-DFT 和 Meta-GGA 级别examples/band_unfolding 目录下给了多个测试例子跑之前先去那里看一眼用法比自己盲试快得多。功能 73 还能把 VASP 计算结果转成 BoltzTraP 需要的格式做输运性质计算时用得上。6. 装完怎么验证用 ZnO 能带数据检查整套环境6.1 三步验证运行链路装完不要直接上大体系先走一遍最小验证。第一步确认二进制可用运行vaspkit -h能打印帮助就说明动态库加载正常。第二步确认菜单能进直接启动 vaspkit在主菜单输入 0 退出这一步验证交互逻辑没坏。第三步确认配置有效进功能 1 然后选 109对一个已有 POSCAR 做格式检查看它能不能正确读出元素类型并和 POTCAR 比对。三步全过环境基本就稳了。6.2 结构层验证第三步完成后看具体输出是否合理。对 ZnO 这个例子执行grep TIT POTCAR确认元素是 O 和 Zn_pv执行cat KPOINTS确认 Gamma 网格和晶格匹配。重点看 POSCAR 和 POTCAR 的元素顺序是否一致这是 108/109 检查的核心价值。6.3 能带层验证跑完 ZnO 的能带计算后用功能 21 或 25 提取能带数据把带隙数值和文献对一下。纤锌矿 ZnO 的实验带隙约在 0.7 eV 附近PBE 泛函通常会低估HSE06 会更接近实验值。如果你的结果和这个数量级差出一倍以上先回头查结构有没有收敛、K 点密度够不够、赝势选得对不对。想看三维能带走势功能 23 直接出 3D 能带图能快速判断整体色散关系是否合理。处理态密度时注意 SET_FERMI_ENERGY_ZERO 保持 .TRUE.这样费米能级会被统一平移到零点画图时不会出现能级对不齐的问题。6.4 反折叠与最终习惯如果你手上有超胞体系功能 28 是检验 vaspkit 是否真正跑通的最好试金石。选一个 examples/band_unfolding 里的例子走一遍输出的反折叠能带数据如果能在原始布里渊区路径上看到清晰的能带分支说明整个链路没有问题。从那以后我每次拿到新结构都会强制走一遍 108/109 再交 VASP这个习惯帮我省下了至少三次因元素不匹配而浪费的机时。希望这份从装到验的流程能让你少走几趟我走过的弯路也希望帮到你。本文还有配套的精品资源点击获取
简介一份面向VASP使用者的中文功能手册系统讲解VASPKIT在材料科学、催化计算中的预处理与后处理应用适合从入门到进阶的科研人员。手册内容按菜单功能组织涵盖VASPKIT简介与配置、VASP输入文件生成和检查、PBE与HSE06杂化泛函能带计算、分波态密度与局域态密度、热力学量校正、差分电荷密度绘制以及NEB路径生成与可视化、虚频校正、能带反折叠、3D能带、费米面、VASP2BoltzTraP接口等催化与电子结构分析模块。此外还介绍了用户自定义界面可用于绘制原子电荷着色结构图、二维材料投影能带等扩展场景并附有引用方式与菜单功能附录。资源共1个文件为PDF格式包体仅2.37MB轻量易保存章节目录清晰便于按需定位。已有3564人学习下载。借助书本示例与逐步说明读者可快速掌握从输入文件准备、电子结构分析到催化反应路径处理的关键操作提升日常VASP计算效率。1. 从 GLIBC 报错说起先把 VASPKIT 0.72 装进 VASP 工作流第一次跑 VASPKIT 0.72 时我差点被一行/lib64/libc.so.6: version GLIBC_2.14 not found劝退。当时只是想用它的能带反折叠功能处理一个超胞结构的投影能带结果./vaspkit敲下去屏幕上先蹦出来的是库版本报错。查了一圈才明白工具没问题是运行环境的 glibc 太老。VASPKIT 是 FORTRAN 写成、在 LINUX 下运行的纯二进制程序体积只有 5M几乎不依赖第三方库但对系统 glibc 版本有硬要求。这个名为 VASPKIT 的 VASP 前处理/后处理工具箱在材料计算和催化计算里早就成了标配——从自动生成 INCAR、POTCAR、KPOINTS 并做格式检查到能带、态密度、差分电荷密度、NEB 路径和虚频校正全在一个交互菜单里完成。适合两类人刚入坑 VASP、不想满世界找脚本的新手以及在各种第三方脚本版本打架里消耗过大量时间的熟手。2. 环境准备下载 0.72 免安装包先看 bin 目录和 EXAMPLES2.1 版本怎么选为什么是 0.72装任何工具之前都要先定版本VASPKIT 也不例外。0.72 是目前功能覆盖最全的 release 版本这个版本的菜单做了一次大调整把功能相似的选项归类成“结构 Options”“电子 Options”“电荷与势与波函数 Options”“Misc Utilities”四个大组同时补上了很多催化相关的新功能。我特别关注的是几个新增任务task 28 加入了 DFT、hybrid-DFT 和 Meta-GGA 级别的能带反折叠Band-Structure Unfoldingtask 25 支持在杂化泛函能带计算里沿不同 k 路径做均匀 k 点采样task 505 可以线性插点生成 NEB 路径task 507 可以做虚频校正还有 task 112 到 115 这组投影态密度功能能一次处理单原子、单元素或指定原子轨道的投影态密度和投影能带。这些功能对做催化、做半导体材料的人来说非常实用。比如 NEB 路径可视化以前要把多个中间态结构导成 PDB 文件得自己写脚本处理 OUTCAR 里的坐标现在 task 405 直接转再比如 HSE06 杂化泛函能带老版本处理投影能带数据时会有 bug0.72 里专门修了提取 hybrid-DFT 级别投影能带数据的问题。所以如果你还在用 0.7x 早期版本建议直接换 0.72菜单结构变了功能入口和参数模板都按新版本对齐。另外要注意VASPKIT 是纯二进制分发没有图形界面安装向导。它几乎不依赖第三方库但依赖系统运行时库尤其是 glibc。装之前先确认你的发行版不太老ldd --version | head -n1常见的 CentOS 7、Ubuntu 16.04 以上都能满足 0.72 的运行要求。如果输出显示 glibc 版本低于 2.14后面运行 vaspkit 就大概率会遇到第 4 章里写的 GLIBC_2.14 not found 报错这点提前有个心理准备。2.2 解压、放置、赋权三步把二进制安顿好拿到 vaspkit.0.72 的压缩包后我一般先解压到临时目录再整体移动到 /home 或 /opt 下统一管理避免直接散落在下载目录里tar -zxvf vaspkit.0.72.tar.gz mv vaspkit.0.72 /home/ cd /home/vaspkit.0.72 ls -l解压后你会看到几个关键内容bin 目录里是 vaspkit 可执行文件how_to_set_environment_variable 是环境变量配置模板EXAMPLES 目录下是各主要功能的测试例子。建议先看一眼文件的权限ls -l /home/vaspkit.0.72/bin/vaspkit如果权限位里没有 x后面直接运行会报 Permission denied。常规操作是先赋执行权限再测试是否可用chmod ux /home/vaspkit.0.72/bin/vaspkit /home/vaspkit.0.72/bin/vaspkit -h-vaspkit -h 能打印出命令行帮助说明这个二进制在你的系统上能正常加载动态库。若这里就报错基本可以断定是 glibc 版本问题直接跳到第 4 章的排查方案。2.3 EXAMPLES 目录比任何文档都快的上手途径VASPKIT 的 EXAMPLES 目录里放了很多功能的测试例子这是我建议新手最先翻的地方。比如 examples/hybrid_DFT_band 目录下有一个用 HSE06 杂化泛函计算磷化镓能带结构的完整案例从 POSCAR 到最终能带数据都有examples/band_unfolding 目录下则是能带反折叠的多个例子对应 task 28 的功能。每个例子目录里基本都备好了输入文件和说明你可以直接在里面运行 vaspkit 看某个功能怎么工作不需要自己从头搭结构。我自己的习惯是装完一个新工具先跑一遍自带 examples不是为了算科学结果而是为了确认“工具本身能跑通”。如果 examples 里的案例跑不出合理数据那多半是配置有问题而不是你手头的研究体系有问题。这样排查范围就一下缩小了。3. 配置 ~/.vaspkit赝势路径、INCAR 写入模式与 PATH 三件事3.1 环境变量文件从哪里来VASPKIT 运行时需要读一个用户级配置文件默认位置是~/.vaspkit。首次使用前要把软件包自带的模板复制过去cp -f how_to_set_environment_variable ~/.vaspkit vi ~/.vaspkit这个文件不复制的话VASPKIT 也能启动但很多功能会因为找不到赝势路径而报错尤其是生成 POTCAR 的功能。复制时要特别注意-f参数它会强制覆盖已存在的 ~/.vaspkit。如果你之前配过旧版本先备份再复制cp ~/.vaspkit ~/.vaspkit.bak3.2 配置项逐条说明打开 ~/.vaspkit 后核心配置项大致如下不同版本字段略有差异0.72 是以这组为准VASP5 .TRUE. # .TRUE. or .FALSE.; Set .FALSE. if using vasp.4.x GGA_PATH ~/POTCAR/GGA # Path of GGA potential PBE_PATH ~/POTCAR/PBE # Path of PBE potential LDA_PATH ~/POTCAR/LDA # Path of LDA potential POTCAR_TYPE PBE # PBE, GGA or LDA GW_POTCAR .FALSE. # .TRUE. or .FALSE. RECOMMENDED_POTCAR .TRUE. # .TRUE. or .FALSE. MINI_INCAR .FALSE. # .TRUE. or .FALSE. SET_FERMI_ENERGY_ZERO .TRUE. # .TRUE. or .FALSE. SET_INCAR_WRITE_MODE OVERRIDE # OVERRIDE, APPEND, BACK-UP-OLD, BACK-UP-NEW这里我按实际使用经验拆开说参数作用我的建议VASP5是否使用 VASP 5.x 语法的 POSCAR/POTCAR现在绝大多数人用 VASP 5.4保持 .TRUE.只有老项目还在用 4.x 才改 .FALSE.GGA_PATH / PBE_PATH / LDA_PATH三种赝势库的根目录一定改成你机器上真实存在的路径建议用绝对路径POTCAR_TYPE生成 POTCAR 时默认从哪个赝势库取默认 PBE 即可别随手改成 GGAGW_POTCAR是否使用 GW 专用赝势不跑 GW 计算就 .FALSE.RECOMMENDED_POTCAR是否按 VASP 官网推荐自动选赝势新手保持 .TRUE.需要特殊后缀时再手动干预MINI_INCAR是否生成精简版 INCAR建议 .FALSE.让模板把该有的参数都写全SET_FERMI_ENERGY_ZERO处理能带/态密度时是否把费米能级平移到零保持 .TRUE.否则画图时费米能级位置对不齐SET_INCAR_WRITE_MODE已有 INCAR 时怎么处理默认 OVERRIDE 直接覆盖想留旧文件就设 BACK-UP-OLD这里最容易被忽略的是 PATH 类配置项里的~。在单引号里~不会展开如果 VASPKIT 读取时不做展开处理就会拿着字符串~/POTCAR/PBE去拼路径结果就是找不到赝势。我一般直接写成/home/username/POTCAR/PBE这种绝对路径省得在不同 shell、不同登录方式下出现“刚才还能用换个终端就找不到赝势”的玄学问题。3.3 把 wanjaspkit 加进 PATH配置文件弄好后为了方便在任意目录启动把 bin 目录加进 PATHecho export PATH/home/vaspkit.0.72/bin/:$PATH ~/.bashrc source ~/.bashrc如果你用的是 zsh就写到 ~/.zshrc别只记 bashrc。加完之后在终端直接输入vaspkit就能启动也可以继续用绝对路径/home/vaspkit.0.72/bin/vaspkit运行。这样做的意义是后面在 examples 里切换目录测试时不用每次找路径。3.4 启动界面长什么样配置正确后启动 vaspkit 会看到一个用 ASCII 字符画的菜单界面。很多新手看到这个“很萌的界面”会愣一下不确定有没有启动成功。其实看到下面这个主菜单就说明启动正常 Structural Options 1) VASP Input Files Generator 2) Elastic-Properties Kit 3) K-Path Generator 4) Structure Editor 5) Catalysis-ElectroChemi Kit 6) Symmetry Search Electronic Options 11) Density-of-States 21) DFT Band-Structure 23) 3D Band-Structure 25) Hybrid-DFT Band-Structure 26) Fermi-Surface 28) Band-Structure Unfolding Charge Potential Wavefunction Options 31) Charge Spin Density 42) Potential-Related 51) Wave-Function Analysis Misc Utilities 71) Linear Optics 72) Molecular-Dynamics Kit 73) VASP2BoltzTraP Interface 91) Semiconductor Calculator 92) 2D-Materials Kit 0) Quit菜单按数字分组输入对应编号回车进入子菜单。这里我提醒一句VASPKIT 是命令行交互式操作不支持鼠标点选。如果你是在远程服务器上通过 SSH 用完全没问题终端字符界面跑得稳。4. 避坑手册装完 VASPKIT 最常见的五个问题4.1 vaspkit: /lib64/libc.so.6: version GLIBC_2.14 not found这是我遇到最多的启动报错现象就是第 1 章那段描述执行 vaspkit 后系统直接抛出 “version GLIBC_2.14 not found (required by vaspkit)”。原因是二进制文件在编译时链接了较新版本的 glibc而运行机器上的 glibc 低于 2.14。常见于老旧的 CentOS 6.x 或同类发行版。解决有两个方向。第一个是升级系统 glibc但不建议因为 glibc 是系统最底层的基础库贸然升级可能连带破坏系统里其他依赖旧版本的应用属于“为了一个工具掀翻整个系统”的做法。第二个是找开发者要一份在低版本 LINUX 环境下编译的版本这个更稳妥相当于让二进制来适配你的系统而不是让系统去冒险适配它。如果不想走这一步也可以考虑换一台 glibc 版本更高的机器比如 CentOS 7 以上或 Ubuntu 16.04 以上直接绕开这个问题。4.2 Permission denied可执行权限没给现象是在 bash 里执行/home/vaspkit.0.72/bin/vaspkit或直接敲vaspkit提示-bash: line 7: ./vaspkit: Permission denied。原因是压缩包解压后文件默认没有执行权限。很多新手会先去查 PATH 配置、查环境变量其实问题根本没到那一步就是最简单的权限位没有 x。解决是赋权chmod ux /home/vaspkit.0.72/bin/vaspkit我建议解压后立刻执行这一步不要等报错再处理。赋权之后运行vaspkit -h如果能看到帮助信息说明二进制本身没问题后续再报错就集中在配置层面。4.3 POSCAR 里藏着制表符VASP 只读出了第一个元素这是一个提交 VASP 任务后才会暴露的坑。现象是任务秒失败日志里写着ERROR: the type information is not consistent with the number of typesVASP 只读出了 POSCAR 里的第一种元素。原因是从某些晶体数据库下载的 POSCAR 文件里分隔符用了制表符\t而不是 VASP 期望的空格。VASP 在读字符串时不能正确处理\t于是元素类型列表没有完整解析。这个问题在 Windows 系统下编辑生成的 POSCAR 或 INCAR 里也会出现表现形式可能是各种奇怪的读入错误。解决是用 VASPKIT 的功能 109 做格式纠正。运行 vaspkit进入功能 1VASP Input Files Generator再选 109程序会自动纠正 INCAR 和 POSCAR 的格式问题并检查赝势元素是否一致。执行完之后重新提交任务就能正常跑起来。如果你不想用 VASPKIT也可以用 sed 把制表符批量替换成空格但最省事的还是 109一步到位。4.4 POSCAR 元素类型和 POTCAR 不一致不报错但结果全错这个坑最隐蔽VASP 不会自动检查 POSCAR 里的元素类型是否与 POTCAR 一致。意味着你拿一份石墨烯的 POSCAR配一个 H 的 POTCARVASP 照样能跑但它算出来的结果是错的而且不会给你任何报错提示。原因是 VASP 的设计里POSCAR 只提供晶格和坐标POTCAR 提供势函数二者是否匹配全靠用户自觉。很多新手在做“替换元素”操作时会顺手把 POSCAR 改了而忘记重新生成 POTCAR或者反过来导致整个计算的物理意义就是错的。解决是靠 vaspkit 的 109 功能做检查。它会读取 POSCAR 里的元素顺序和 POTCAR 里的元素逐一比对发现不一致就明确提示。提交任务前跑一遍 109比自己盯着一行行核对靠谱得多。这一条我已经写进了自己的固定流程拿到任何新结构都先 109 再交 VASP。4.5 赝势路径写~换个终端就找不到 POTCAR现象是这个终端里 vaspkit 能正常生成 POTCAR换个终端或者换一种登录方式再运行就提示找不到指定赝势目录。原因是配置文件里写的是~/POTCAR/PBE这类带波浪号的路径。~在 shell 里会展开成家目录但 VASPKIT 是 FORTRAN 程序它读取配置文件时未必会做同样的展开处理。不同环境下表现不一致于是就成了“时好时坏”的玄学问题。解决是打开 ~/.vaspkit把 GGA_PATH、PBE_PATH、LDA_PATH 全部改成绝对路径例如/home/username/POTCAR/PBE。改完保存重启 vaspkit问题消失。这条经验我后来写进了自己的部署笔记所有配置文件里涉及路径的参数一律用绝对路径不给~留发挥空间。5. 生成输入文件与能带计算把 ZnO 从 POSCAR 带到能带反折叠5.1 第一步准备一份能用的 POSCAR配置好环境、看完避坑之后就该用真实体系验证 vaspkit 的完整流程了。我一般拿纤锌矿 ZnO 做测试因为这个结构简单、空间群明确而且能带和态密度的参考数据容易找。先从公开晶体数据库里检索 ZnO空间群选 P63mc 的六方相下载 POSCAR 文件典型文件名类似ZnO-1811117.vasp。把它复制成 VASP 标准文件名cp -f ZnO-1811117.vasp POSCAR如果你拿到的是 CIF 文件需要先用 VESTA 或 VTST 脚本转成 POSCAR注意原子位置的分数坐标是否处理正确。这一步不能省后续 vaspkit 生成的 KPOINTS、POTCAR 都以 POSCAR 里的元素和晶格为准。5.2 功能 101按任务类型生成 INCAR在存放 POSCAR 的目录下运行 vaspkit输入1进入 VASP Input Files Generator会看到功能子菜单101) Customize INCAR File 102) Generate KPOINTS File for SCF Calculation 103) Generate POTCAR File with Default Setting 104) Generate POTCAR File with User Specified Potential 107) Reformat POSCAR File by Specified Order of Elements 108) Successive Procedure to Generate VASP Files and Check 109) Check All VASP Files输入101会进入 INCAR 模板选择界面里面列出了针对不同计算任务的模板ST) Static-Calculation SR) Standard Relaxation MG) Magnetic Properties SO) Spin-Orbit Coupling D3) DFT-D3 no-damping Correction H6) HSE06 Calculation PU) DFTU Calculation MD) Molecular Dynamics GW) GW0 Calculation BS) BSE Calculation DC) Elastic Constant EL) ELF Calculation BD) Bader Charge Analysis OP) Optical Properties EC) Static Dielectric Constant PC) Decomposed Charge Density FD) Phonon-Finite-Displacement DT) Phonon-DFPT NE) Nudged Elastic Band (NEB) DM) The Dimer Method FQ) Frequence Calculations LR) Lattice Relaxation输入LR就会生成一个用于晶格弛豫的 INCAR。新手看到这么多模板可能会慌其实这里的选择逻辑很简单你要算什么任务就选对应的首字母组合。比如跑 Bader 电荷分析就选 BD跑 NEB 就选 NE跑 HSE06 静态计算就选 H6。每个模板生成后部分参数还需要手动调整比如 DFTU 的 U 值、NEB 的 IMAGES 数目vaspkit 会在生成时用一句 “You MUST Know What You Are Doing” 提示你别无脑跑。这里关于 INCAR 里的 ENCUT我想多说一句。默认生成的 INCAR 里 ENCUT400 往往是注释掉的程序保留的是PREC NormalVASP 会自动把 ENCUT 设为所有元素 ENMAX 的最大值。如果你用 LR 模板优化晶格常数vaspkit 会生成PREC HighVASP 此时会把 ENCUT 自动设为 max(ENMAX)*1.3这样对晶格常数优化的收敛更有利。不过较新的官方建议已经不推荐用 PRECHigh而是设 PRECAccurate 并手动定 ENCUT你可以在模板基础上改成这种写法确保不同体系间 ENCUT 一致。5.3 功能 102KPOINTS 与 POTCAR 一次生成INCAR 就绪后回到子菜单输入102生成 SCF 计算的 KPOINTS。程序会先提示你选择撒点方式输入2选择 Gamma Scheme这是大多数体系下最稳妥的选择。接着提示输入 K 点密度精度值* Accuracy Levels: (1) Low: 0.06~0.04; (2) Medium: 0.04~0.03; (3) Fine: 0.02-0.01; (4) Gamma-Only: 0. Input KP-Resolved Value (unit: 2*PI/Angstrom):输入0.04一般是够用的。0.04 这个值的意思是倒空间每个 K 点间隔约 0.04 × 2π/Åvaspkit 会根据你 POSCAR 里的倒空间晶格矢量比例自动算出各个方向的 K 点数。生成结果看起来类似K-Mesh Generated with KP-Resolved Value : 0.0400 Gamma 9 9 5 0.0 0.0 0.0这里 9×9×5 是程序根据你晶格倒空间比例自动分配的不是固定值。对于正交体系有个简单的判断准则实空间晶格矢量乘以对应方向的 K 点数约大于 20Å 就比较稳妥ZnO 这个例子里 ka≈29.62Å满足经验规律。非正交体系因为倒格矢长度和晶格常数不满足简单的反比关系用 vaspkit 自动分配反而更合理。生成 KPOINTS 的同时vaspkit 会根据 POSCAR 里的元素类型从你配置的赝势库里提取并组合生成 POTCAR。前提是 ~/.vaspkit 里的 PBE_PATH 指向正确目录POTCAR_TYPE 保持 PBE。生成成功后用以下命令确认grep TIT POTCAR输出会列出 POTCAR 里实际包含的元素赝势。如果显示的是 O 和 Zn说明 POTCAR 和 POSCAR 的元素顺序是对上的。5.4 功能 104手动指定赝势类型比如 Zn_pv默认的 103 功能是按 VASP 官网推荐自动选赝势。但有些场景需要手动指定比如做 Bader 电荷分析时建议使用带半芯态的赝势。PBE 赝势库里常见后缀含义如下后缀含义无后缀只把价电子当作价态处理_pv次外层的 p 半芯态也当价电子_sv次外层的 s 半芯态也当价电子_d次外层的 d 半芯态也当价电子_h针对短键二聚体的硬赝势如 O2、CO、N2、F2 等常见的做法是先测试两种赝势对自己关心的性质影响大不大。如果影响不明显优先用无后缀版本因为包含更多价电子意味着截断能上升、计算量明显增大。如果体系里有 O2 这类短键分子官方建议用 _h 势。手动干预时用功能 104按提示依次输入元素和赝势类型比如给 Zn 输入Zn_pv。如果指定目录下没有对应赝势文件程序会提示重新输入。生成后同样用grep TIT POTCAR验证元素顺序和类型。5.5 功能 109/108提交 VASP 前最后一道检查如果这时直接把任务交出去你可能会遇到我在第 4 章描述过的场景从数据库下载的 POSCAR 里藏着制表符VASP 只读出了 O 没有读出 Zn任务秒失败。vaspkit 从 0.71 开始加入了功能 109专门做 INCAR 和 POSCAR 的格式纠正、检查赝势元素一致性。输入109程序会告诉你所有必要文件是否存在并比对 POSCAR 和 POTCAR 的元素类型是否一致。这个检查非常关键因为 VASP 本身不检查你拿石墨烯配 H 的赝势它也照跑不误只是结果毫无意义。所以更省心的做法是直接用功能 108依次完成 INCAR 设置、KPOINTS 生成、POTCAR 生成最后自动调用 109 检查所有文件。整套流程一句话概括就是“VASP 之前先 KIT 一下”。我后来所有新体系的第一次计算都走这个流程没有再被制表符和元素不匹配坑过。5.6 能带路径与能带计算功能 3、21、25、28能带计算的 KPOINTS 和自洽计算不一样需要沿着布里渊区高对称点路径取点。传统做法是去 SeeK-Path 网站或者用 Materials Studio 找高对称点再写脚本插值。vaspkit 把这一步集成了输入3进入 K-Path Generator选择 3D bulk structure程序会分析晶体对称性并给出建议路径。对纤锌矿 ZnOvaspkit 会生成类似GAMMA-X-U|K-GAMMA-L-W-X的两条路径同时生成两个文件KPATH.in和PRIMCELL.vasp。KPATH.in 的第二行是每段路径插值的 K 点数目如果计算量紧张可以调小PRIMCELL.vasp 是程序生成的单胞能带路径只针对 primitive cell所以要把这个文件作为后续计算的结构cp -f PRIMCELL.vasp POSCAR接下来按计算层级选择后处理功能。普通 PBE 泛函能带用功能 21HSE06 杂化泛函能带用功能 25想直观展示三维能带走势可以用功能 23超胞或者表面体系想看原始晶格能带投影则用功能 28 做能带反折叠。0.72 的能带反折叠支持 DFT、hybrid-DFT 和 Meta-GGA 级别examples/band_unfolding 目录下给了多个测试例子跑之前先去那里看一眼用法比自己盲试快得多。功能 73 还能把 VASP 计算结果转成 BoltzTraP 需要的格式做输运性质计算时用得上。6. 装完怎么验证用 ZnO 能带数据检查整套环境6.1 三步验证运行链路装完不要直接上大体系先走一遍最小验证。第一步确认二进制可用运行vaspkit -h能打印帮助就说明动态库加载正常。第二步确认菜单能进直接启动 vaspkit在主菜单输入 0 退出这一步验证交互逻辑没坏。第三步确认配置有效进功能 1 然后选 109对一个已有 POSCAR 做格式检查看它能不能正确读出元素类型并和 POTCAR 比对。三步全过环境基本就稳了。6.2 结构层验证第三步完成后看具体输出是否合理。对 ZnO 这个例子执行grep TIT POTCAR确认元素是 O 和 Zn_pv执行cat KPOINTS确认 Gamma 网格和晶格匹配。重点看 POSCAR 和 POTCAR 的元素顺序是否一致这是 108/109 检查的核心价值。6.3 能带层验证跑完 ZnO 的能带计算后用功能 21 或 25 提取能带数据把带隙数值和文献对一下。纤锌矿 ZnO 的实验带隙约在 0.7 eV 附近PBE 泛函通常会低估HSE06 会更接近实验值。如果你的结果和这个数量级差出一倍以上先回头查结构有没有收敛、K 点密度够不够、赝势选得对不对。想看三维能带走势功能 23 直接出 3D 能带图能快速判断整体色散关系是否合理。处理态密度时注意 SET_FERMI_ENERGY_ZERO 保持 .TRUE.这样费米能级会被统一平移到零点画图时不会出现能级对不齐的问题。6.4 反折叠与最终习惯如果你手上有超胞体系功能 28 是检验 vaspkit 是否真正跑通的最好试金石。选一个 examples/band_unfolding 里的例子走一遍输出的反折叠能带数据如果能在原始布里渊区路径上看到清晰的能带分支说明整个链路没有问题。从那以后我每次拿到新结构都会强制走一遍 108/109 再交 VASP这个习惯帮我省下了至少三次因元素不匹配而浪费的机时。希望这份从装到验的流程能让你少走几趟我走过的弯路也希望帮到你。本文还有配套的精品资源点击获取
平方复杂度无法被常规优化消除｜TaoToken 统一 Key 通道实测)
