简介：这是面向Vasp使用者的vtstscripts 0.91脚本工具包，适用于Linux环境下开展结构优化、能带计算、分子动力学及第一性原理计算等任务，可帮助材料科学、凝聚态物理方向的研究者减少重复操作并提升计算效率。压缩包内含117个文件，主要包含93个Perl脚本、8个Shell脚本、7个Perl模块及4个GNU绘图脚本，涵盖自动生成输入文件、解析输出数据、设置并行任务、检查收敛误差及后处理绘图等功能，整体包体仅257KB，轻量易部署。资源已吸引578人学习浏览，适合具备一定Vasp基础、希望借助脚本自动化完成多步骤模拟的计算用户。通过掌握该脚本集，读者可快速复用成熟计算流程，完成KPOINTS生成、能带解析、态密度与电荷密度后处理等操作，同时也可依据自身研究需求对脚本进行定制扩展，较为适合日常科研中的重复性计算场景。
1. vtstscripts_0.91.tar.gz 是什么：把 NEB 从“玄学”变成可重复的过渡态搜索流水线
过渡态搜索是材料模拟里最让人头疼的一类任务。一个反应路径上鞍点找不准，后续的速率常数、能垒全是空中楼阁。vtstscripts_0.91.tar.gz 就是这样一套专为 VASP 准备的过渡态脚本工具包：它把初末态插值、NEB 力投影、收敛判断、Dimer 搜索这些环节封装成可复现的命令行工具，让我在多个模拟项目里不再靠手改 POSCAR 碰运气。0.91 这个版本分支结构简单、依赖少，适合刚接触过渡态计算的人上手，也够老手做批量筛选。这篇笔记会把我的用法、参数和踩过的坑写清楚。
2. 装进环境：解包、依赖与三分钟跑通最小命令
拿到 vtstscripts_0.91.tar.gz 的第一件事不是急着跑 NEB，而是把环境收拾利索。这包脚本的定位是“VASP 过渡态计算的瑞士军刀”，脚本本身不负责算力，所有数值还是交给 VASP 完成。这意味着两件事：一是脚本要在你能登录计算节点的 shell 里能直接调用；二是你编译的 VASP 可执行文件必须带 VTST 补丁，否则脚本生成的 INCAR 项会被 VASP 静默忽略。很多翻车现场就出在这两步没做干净。
2.1 解包与目录结构：0.91 包里到底有什么
常见的做法是先给这包脚本找个固定住处，不要解压在临时目录里。我会把它放到~/tools/下面，然后用软链接的方式把脚本目录挂到 PATH，这样升级或切换版本时只要动一个链接。
mkdir -p ~/tools tar -xzvf vtstscripts_0.91.tar.gz -C ~/tools cd ~/tools/vtstscripts ls -1 | head -30解压后你会看到一批以pl结尾的 Perl 脚本，还有少量 Python 脚本和示例目录。0.91 版本分支的特点是脚本数量不多，核心功能基本靠 Perl 完成，Python 脚本集中在势能面插值、spline 拟合这些后处理环节。目录里通常还会带几个示例，我的习惯是把示例目录完整留着，等配置完环境后逐个跑一遍，确认脚本行为正常再用于正式任务。
这里有个容易忽略的点：这类脚本包没有 make install 步骤，解压即用。但正因为如此，脚本对“当前工作目录”极其敏感。大部分脚本会假设你是在某个 NEB 任务目录下运行它，而不是在自己的家目录里。所以解包后第一件事不是直接跑，而是先把 PATH 配好。
2.2 依赖检查与 PATH 配置：别让脚本卡在 Perl 模块上
0.91 的脚本主体是 Perl，所以系统里必须有一个可用的 Perl 解释器，同时还需要几个常见的 Perl 模块。如果你用的计算节点是精简安装，很可能缺这些东西。我的习惯是先跑一条命令验证脚本能否被正确加载，然后再配 PATH。
export PATH=~/tools/vtstscripts:$PATH which nebmake.pl nebef.pl perl -e "use strict; use Getopt::Long; print qq(perl-ok\n)"which命令用来确认脚本目录已经被 shell 找到，避免后面在并行任务里因为找不到脚本而报错。第二行那条perl -e是快速检查 Perl 核心模块是否可用，Getopt::Long是解析命令行参数的基础模块，几乎每个脚本都会用到。如果这一步报错，先装系统自带的 Perl 模块包。
PATH 配置通常写进~/.bashrc，但注意计算节点和登录节点可能是两套环境。我见过的问题是：登录节点配了 PATH，提交到计算节点后脚本找不到。这时候需要在作业脚本或者~/.bashrc里重新 export 一次。0.91 版本分支对 Python 的依赖通常只在 spline 拟合和频率后处理阶段，如果只是跑最小 NEB 验证，可以先不管 Python 版本。
2.3 最小 NEB 命令：从两个 POSCAR 到可检查的中间构型
跑通最小流程只需要三样东西：初态结构、末态结构、一份能跑的 VASP 输入。初态和末态各自放在00和02目录里，目录名是脚本约定的，别随意改。中间插值点由脚本生成，不需要手动写。
mkdir -p ~/test_neb/{00,02} cp ~/my_structures/reactant/POSCAR ~/test_neb/00/POSCAR cp ~/my_structures/product/POSCAR ~/test_neb/02/POSCAR cd ~/test_neb nebmake.pl 00 02 2nebmake.pl的第三个参数是插入的中间点数，这里写 2 表示在初末态之间插入 2 个构型。运行后当前目录下会多出01目录。加上00和02，整个 NEB 带上共有 4 个 image。你可以用nepos.pl之类的脚本抽查生成的插值构型是不是存在原子重叠。如果初始结构差异很大，插值点数太少会让中间构型离谱，甚至出现原子挨得过近的情况。
生成目录后，还要自己补齐 INCAR、POTCAR、KPOINTS。注意 INCAR 里IMAGES这个参数必须和插值点数一致，这里就是 2。如果我这里写 4 而插值只有 2，VASP 会多算两个不可控的 image，结果完全没法解释。这也是 0.91 时代最常见的低级错误之一。
cp ~/my_templates/INCAR . cp ~/my_templates/POTCAR . cp ~/my_templates/KPOINTS . echo "IMAGES = 2" >> INCARINCAR 里其余关键项建议保持 VTST 补丁的约定，后面第 3 章会展开。配置完成后直接提交 VASP 任务。等到算完后，第一时间用nebef.pl检查每个 image 的受力，而不是关心能量曲线。
nebef.plnebef.pl会打印每个 image 的最大原子受力，单位是 eV/Å，还会给出力的最大值。我一般判断是否收敛，就看这个最大受力有没有压到 0.05 eV/Å 以下。如果没收敛，就回去调整插值点数、弹簧常数，或者干脆换 Dimer 方案。
3. 过渡态计算流程：从初末态插值到 NEB 收敛参数怎么设
最小命令能跑通之后，真正的问题才浮出来：怎么让过渡态搜索可靠地收敛。这一章讲我固定用的那一套流程，每一步都有对应的参数和判断标准。初末态准备、插值设置、爬坡选择、收敛判据，这四个环节里任何一个偷懒，最后都会变成能垒数值上的谜之偏差。
3.1 初末态准备：结构优化做到什么程度才不算白跑
NEB 的初末态是整条路径的两个锚点。锚点不稳，整条带子都是漂的。很多计算失败的根本原因不在 NEB 参数，而在初末态压根没优化好。我会把初末态分别做结构优化，收敛标准压到每个原子受力最大不超过 0.02 eV/Å，而且能量波动的变化要小于 1e-6 eV 量级。
这里说的“优化好”不仅指力收敛，还要确认结构没有停留在明显不合理的高能构型。比如表面吸附体系，吸附物可能有两个能量很接近的局域极小，初态取错一个，整个过渡态路径都会绕远。我的习惯是先做一遍从头到尾的几何优化，再用优化后的结构重新检查配位数、键长是否合理，最后才作为 NEB 端点。
端点还有一个隐含要求：初末态必须用同一套晶胞、同一个截断能、同一种泛函。如果初态是 PBE 算的，末态换成了别的泛函，VASP 读入结构后能量基线不一致，能垒数值会非常难解释。0.91 时代没有自动检查这些的工具，全靠人工盯住 INCAR 里的ENMAX和EDIFF。我一般是直接复制同一份 INCAR 模板给初末态优化，只改结构文件。
3.2 插值与“爬坡”设置：nebmake 参数和 INCAR 关键项
插值点数量直接影响路径分辨率。点数太少，MEP 会漏掉关键拐点；点数太多，计算成本线性上升，而且容易出现相邻 image 间构型变化过小、力投影数值抖动的现象。我的经验是：简单的气相反应或表面单原子扩散，插 4 到 6 个点够用；涉及分子重排、键断裂的复杂路径，先插 6 到 8 个点，跑通后再加密验证。
nebmake.pl 的做法是线性插值笛卡尔坐标。这个近似在端点结构相似时非常好用，但端点结构差别很大时，线性插值可能让中间构型出现原子碰撞。我一般会在生成后检查最短原子间距，如果明显小于该元素应有的共价半径之和，就改用 Image Dependent Pair Potential 插值，或者直接减少单次插值步长。0.91 脚本包里的生成工具支持不同插值方式，但最简单的检查方法是用可视化软件打开 01 目录的 POSCAR。
INCAR 的设置才是 NEB 收敛的核心。下面是我跑 CI-NEB 的基准模板，注意LCLIMB = .TRUE.表示开启爬坡弹性带，SPRING是弹簧常数。
IBRION = 3 POTIM = 0 NSW = 120 LCLIMB = .TRUE. IMAGES = 4 SPRING = -5 EDIFF = 1E-6 EDIFFG = -0.05IBRION = 3表示使用 VTST 补丁内置的优化器，而不是 VASP 自带的共轭梯度。POTIM = 0是必须的，因为 NEB 的力是投影后的力，步长由优化器自行控制。EDIFFG这里是负值，表示按力的阈值判断收敛，我一般取 -0.05，也就是最大受力小于 0.05 eV/Å 就停止。SPRING = -5表示弹簧常数 5 eV/Ų，负号只是约定。
如果初末态能量差很大，或者路径中存在明显的亚稳态，我不建议一上来就开LCLIMB。常见的做法是先跑几轮普通 NEB，把路径大致拉出来，再把LCLIMB打开做精细收敛。0.91 版本的脚本包里，你可以直接改 INCAR 重启，不需要重新插值。这个策略能减少爬坡带在早期震荡的几率。
3.3 收敛判断：nebef 力曲线怎么看
很多人把注意力放在能量上，看初末态能量差就算完事，这是本末倒置。NEB 的收敛判据是每个 image 上投影后的原子力，而不是能量。我会用nebef.pl拉出所有 image 的受力，然后找那个受力最大的点，判断它离收敛线有多远。
nebef.pl nebconverge.plnebconverge.pl会把所有 image 的最大受力和均方根受力汇总成一个简要列表。看的时候重点看两个值：最大受力MAX和均方根受力RMS。我一般要求最大受力小于 0.05 eV/Å，均方根小于 0.03 eV/Å。只盯着最大力还不够，要看受力最大的 image 是不是出现在路径中部。如果最大受力始终集中在端点附近，说明初末态没有固定好或没优化好，而不是插值点的问题。
力曲线如果呈现锯齿状，比如相邻 image 受力方向交替变化，通常是弹簧常数太大或插值点太密。这时把SPRING往小调，比如从 -5 改成 -3，或者减少插值点数重新跑。相反，受力曲线平缓但整体偏高，收敛很慢，说明弹簧常数偏小或步长控制有问题，可以适当调大SPRING。
另外，nebjob.pl这类任务管理脚本能帮你在作业结束后自动调用nebef.pl，但我习惯手动跑一次，因为有时候作业脚本里的环境变量没带全，自动调用会静默失败。
4. 避坑指南：vtstscripts_0.91 使用中的五个常见翻车现场
任何脚本包用久了都会攒下一堆血泪经验。这一章我挑了五个频率最高的翻车现场，每条按现象、原因、解决三步来写。前三条和参数配置有关，后两条和环境、方法选择有关。新手照着排查能少走弯路，熟手也可以对照看看自己的流程里有没有隐藏的坑。
4.1 现象：nebmake.pl 报错，提示目录已存在
第一次跑 NEB 的人最容易遇到这个。插值点目录被上一次计算占用，脚本为了不覆盖已有数据，直接拒绝执行。原因通常是任务中断后你想重新插值，但忘了清理中间目录。解决方法是把旧的中间目录清干净，或者把整个任务目录复制一份再插值。
cd ~/test_neb rm -rf 01 nebmake.pl 00 02 2这里有个更隐蔽的情况：如果你之前用 4 个插值点跑过，现在想改成 2 个，只清理 01 目录是不够的，因为旧的 01、02、03、04 都在，新脚本生成的文件会和旧目录混在一起。我会把任务目录里除了00、02之外的所有 image 目录和OUTCAR备份后删除，再重新生成。
4.2 现象：NEB 跑完力还是不收敛，力曲线像锯齿
力曲线锯齿状通常是插值点数或者弹簧常数不合适。原因我来回踩过几次：插值点太少导致相邻构型差异大，弹簧力分量会振荡；插值点太多又会让优化器在细碎路径上反复调整。解决方法是先看最大受力集中在哪个 image，如果是中间几个点轮流高，就把插值点减少到 4 附近，同时把SPRING调成 -3 再跑。
不要一上来就无限增大NSW。NEB 优化步数给到 300 已经是上限，超过 300 仍不收敛，问题几乎都在路径初值或弹簧设置上，而不是时间不够。我现在的习惯是固定NSW = 120，如果 120 步没收敛，就检查路径而不是加步数。
4.3 现象：INCAR 里明明写了 LCLIMB = .TRUE. 却不起作用
这是最隐蔽的坑，提示乱用的那种乏力感。原因不是脚本的问题，而是 VASP 可执行文件没有打 VTST 补丁。VASP 原生版本不认识LCLIMB、SPRING、ICHAIN这些关键词，遇到后要么忽略要么报错。解决方法是重新编译 VASP，把 VTST 补丁打进去，生成带过渡态功能的二进制。
怎么确认自己的 VASP 打了补丁？跑一个只有两三个离子的极小体系，在 INCAR 里加上ICHAIN = 0，看能否正常启动且输出文件里出现 NEB 相关信息。或者直接在编译目录里查看二进制和源码补丁的时间戳。0.91 脚本包本身不会检查这一点，所以别指望脚本给你报错。另外还要注意，同一集群上不同节点可能装了不同版本的 VASP，提交作业前确认作业脚本加载的哪个二进制。
4.4 现象：脚本报错，提示找不到某个 Perl 模块
这类错误通常出现在你从源码编译 Perl 模块的环境里，或者计算节点的系统 Perl 路径和登录节点不一致。原因就是@INC路径里没有包含脚本需要的模块目录。解决方法是先定位模块位置，再设置PERL5LIB环境变量。
perl -e "print qq(@INC), qq(\n);" export PERL5LIB=/path/to/perl/modules:$PERL5LIB把PERL5LIB写进作业脚本里，而不是只放在登录节点的~/.bashrc。因为很多作业管理系统在计算节点上不会加载登录节点的环境变量。这个坑折腾过我一整天，后来我把环境检查写成了作业脚本的前几行，每次提交自动打印关键环境变量，出了问题一眼就能定位。
4.5 现象：初末态能量比中间构型还高，MEP 曲线倒挂
这是方法选择问题。初末态没有被定位到足够的局域极小，插值出来的路径会越过一个不存在的能垒，导致能量曲线看起来像倒过来的山峰。另一个原因是初末态里的原子编号不一致，同一原子在初态和末态映射错位，线性插值会制造出明显不合理的高能结构。解决方法是回到初末态优化，确认它们各自确实是邻近的局域极小，还要检查 POSCAR 里原子排序是否一致。
如果初末态本身没问题，但 MEP 依然倒挂，我会放弃 NEB，改用 Dimer 方法从附近的鞍点初猜搜索。0.91 脚本包里带了 dimer 相关工具，先用小体系粗搜出鞍点，再用这个鞍点附近的结构作为 NEB 端点之一。这个策略比硬调 NEB 参数更省时间，后面第 5 章我会再展开。
5. 进阶用法：Dimer 搜索、频率验证与批量筛选技巧
过渡态脚本包真正的价值不在单点跑通，而在批量筛选和验证。这一章讲我常用的三个进阶套路：当 NEB 初猜不可靠时用 Dimer 兜底；搜到鞍点后用频率计算验证是不是真鞍点；最后用脚本循环批量跑多条反应路径。
5.1 Dimer 方法：初猜不好时的兜底方案
当两条反应路径纠缠不清时，我会先用 Dimer 方法粗搜。做法是把脚本包里的 dimer 工具链挂到任务目录，INCAR 里设置ICHAIN = 2之类的 Dimer 关键词，从一个中等大小的结构扰动出发，让 VASP 自动旋转 Dimer 寻找曲率反转方向。这个方法的成本比 NEB 高，但它不依赖初末态插值，适合处理困难路径。
我的习惯是先用小基组或低精度快速粗搜，收敛到受力 0.1 eV/Å 左右，拿到鞍点初猜再过渡到 CI-NEB 精细收敛。这样能把最耗时的精细搜索压到很小范围。
5.2 频率验证：别把鞍点找成极小点
搜到鞍点后，必须做振动频率验证。做法是把鞍点构型拿去算一次频率，确认虚频数量为 1，且虚频对应的振动方向连接反应物和产物。这个步骤能筛掉大量假鞍点。0.91 脚本包里提供的后处理脚本可以帮助生成频率计算的输入结构，但频率验证本身由 VASP 的离子响应计算完成。
我一般会用脚本把鞍点构型复制出来，INCAR 换成频率计算模板，只做单点二阶导计算。虚频方向判断最稳妥的方法是可视化振动模式：把OUTCAR里的虚频位移向量叠加到结构上，看两端各指向哪个极小点。如果虚频不止一个，说明结构处在更高阶鞍点，需要沿最低虚频方向扰动后重新优化。
5.3 一个习惯：用循环脚本批量筛选
做表面反应网络时，十几个初始构型要逐个跑，我固定用一个循环脚本来管理目录生成、提交和初步判断。下面是一个最小可用的批量循环示例，它遍历初末态目录名，自动调用 nEB 脚本生成插值目录。
for reac in r1 r2 r3; do prod=${reac}_prod mkdir -p batch/$reac cd batch/$reac nebmake.pl ../$reac ../../products/$prod 4 cp ../../template/INCAR . echo "IMAGES = 4" >> INCAR cd ../.. done这个循环的核心价值是保证所有任务用同一份模板、同一个插值点数，避免手动复制时漏改参数。批量跑完后，我用grep从每个目录的 OUTCAR 里提取最大受力，做一次全局排序，优先处理那些离收敛最远的任务。这个排查逻辑让我在上周的模拟项目 X 里一天筛完了 12 条反应路径，比手动逐个跑节省了大半时间。
最后说个教训：有一次在模拟项目 X 里，我偷懒把初态优化收敛标准放到了 0.05，直接扔给 NEB，结果能垒偏差了 0.3 eV，白跑了一周。现在我的习惯是所有端点统一压到 0.02，无论看上去多接近。宁可多花半天优化端点，也不要在过渡态搜索里猜结构。这套工作流我已经跑了几十个体系，希望帮到你。
本文还有配套的精品资源，点击获取
