进行 VASP 计算,尤其是进行过渡态搜索(如 NEB、Dimer 方法)时,VTST Scripts 绝对是不可或缺的效率神器。本文系统整理了 VTST 脚本库中所有命令的功能与用法。
一、通用脚本
主要用于日常的 VASP 文件管理与基础结构分析。
Vasp.pm:
核心 Perl 模块,提供读写 POSCAR、向量点积及求模等底层函数支持。
vef.pl:
提取并打印 VASP 运行过程中每一步的受力与能量。
vfin.pl
读取 OUTCAR 标签,自动打包拷贝结果文件到指定输出目录,并将运行目录中的 CONTCAR 覆盖 POSCAR,为下一步计算做好准备。
boxset.pl
根据指定的晶格常数重新缩放并生成新的 POSCAR。
posinterp.pl
在两个 POSCAR 结构之间按指定比例进行线性插值,生成中间结构。
pos2rdf.pl
计算特定原子的径向分布函数 (RDF)。
neighbors.pl
输出特定原子的近邻原子距离。
diffcon.pl
计算两个结构之间原子的位移差异(标量及矢量和)。
dist.pl
计算两个结构配置之间的均方根位移距离。
modemake.pl
生成两个 POSCAR 之间的单位向量,并输出至 MODECAR。
二、格式转换
用于在不同晶体建模软件和可视化工具间无缝切换。
pos2con.pl / con2xyz.pl
实现 POSCAR、CON 以及 XYZ 坐标文件的互相转换。
xdat2pos.pl
将 XDATCAR 轨迹中的指定步数提取并保存为 POSCAR 文件。
xdat2xyz.pl
将整个 XDATCAR 转换为可用 Jmol 播放观看的 XYZ 动画轨迹文件。
xdat2vdat.pl
(限 VASP 5.2) 使用前向差分法计算原子速度,并输出至 VDATCAR 文件。
三、 NEB 过渡态计算
专门针对微弹性带方法的强大辅助工具集。
nebmake.pl
核心制备工具!输入初态、末态和插点数,自动线性插值并生成从 00 到 NI+1 的所有镜像计算目录。
neb2dim.pl / neb2lan.pl
从 NEB 结果中自动提取最高点结构与曲率,无缝建立 Dimer 或 Lanczos 鞍点搜索计算。
nebef.pl
获取当前 NEB 链上各个镜像点的能量与受力。
nebbarrier.pl
提取反应路径数据,生成包含距离、能量和受力的 neb.dat 文件。
nebspline.pl
基于受力和能量数据绘制三次样条插值曲线,输出反应路径图 mep.eps 以及极值点信息。
nebmovie.pl
将反应路径生成动态 XYZ 动画。
nebconverge.pl
在作业运行时,实时监控每个镜像点的能量受力收敛情况。
nebresults.pl
一键分析工具!任务结束后自动连串运行提取、绘图、动画及收敛分析。
nebfreeze.pl / nebavoid.pl
几何微调工具,可用于冻结特定原子,或自动推开距离过近的原子以防结构崩溃。
四、态密度分析
快速处理能带与 DOS 数据,特别是 d 带中心理论分析。
split_dos
关键脚本!将总 DOSCAR 切割拆分为各个原子的分波态密度文件 (DOS1, DOS2 ...)。
dosanalyze.pl
计算特定原子、指定能带的态密度中心(常用于 d 带中心计算)。支持通过半峰宽截断计算加权平均。
doslplot.pl
快速绘制所选原子的局域态密度与体系总态密度对比图 (输出 EPS 格式)。
五、动力学矩阵与频率
用于计算振动频率、零点能及反应速率前因子。
dymmatrix.pl
处理 OUTCAR,构建质量缩放的动力学矩阵,并输出频率和简正振动模式。
dymseldsp.pl / dymselsph.pl
根据原子位移或几何中心,智能筛选需要考虑振动的原子,生成对应的 DISPLACECAR 位移控制文件。
dymeffbar.pl / dymzpbar.pl
计算零点能 (ZPE) 修正,以及包含量子隧穿效应的有效势垒。
dymcmpdisp.pl / dymfit.pl / dymextract.pl / dymreorder.pl
提供矩阵的对比、拼接、提取、重排序功能,支持分批次计算动力学矩阵并整合。
dymprefactor.pl
计算基态到过渡态的指前因子。
dymmodes2xyz.pl
将算得的振动模式转换为 XYZ 动画,直观检视虚频对应的振动方向。
六、 Dimer 与电荷密度
diminit.pl
初始化 Dimer 鞍点搜索目录。
dimplot.pl
绘制 Dimer 优化过程中的受力、能量及曲率曲线。
dimmins.pl
从收敛的 Dimer 结果出发,沿极小化模式生成用于回到初末态的起始构型。
chgavg.pl / chgparavg.pl
计算并输出平均电荷密度 (CHGCAR) 或部分电荷密度 (PARCHG)。
chgsum.pl / chgdiff.pl
计算多个电荷密度的加权和或电荷差分密度。
chg2cube.pl
将 VASP 电荷密度转换为常见的 CUBE 格式,方便不同软件可视化。