vasp做过渡态经验<一>：NEB

sky

这算是一篇新手教学贴，简单谈一下做过渡态的经验，方便各种萌新入门。

NEB具体算法原理可见sob老师的帖子：过渡态、反应路径的计算方法及相关问题 - 思想家公社的门口：量子化学·分子模拟·二次元

最好使用vtst tool版本的vasp，里面有很多好用的算法，比原版vasp好用。编译方法见官网：Transition State Tools for VASP

以下NEB计算的方法只介绍vtst tool版本的vasp

NEB计算步骤如下：

1. 准备好初末构型的vasp文件: is.vasp，fs.vasp，以及脚本：nebmake.pl（来源：Transition State Tools for VASP），然后输入命令：

1. > nebmake.pl is.vasp fs.vasp 8

复制代码

最后一个8代表生成8个image，如果运行成功，会在当前文件夹生成编号为00-09的十个文件夹，每个文件夹包含一个POSCAR。编号00和09的文件夹对应is.vasp和fs.vasp。

之后把配置文件INCAR，KPOINTS，POTCAR放到当前文件夹。INCAR文件需要有以下字段：

1. IMAGES=8
   
2. LCLIMB = .TRUE.

复制代码

IMAGES代表有8个image，LCLIMB是是否开启climbing image方法，详见官网。

nebmake.pl会用线性插值的方法在is.vasp和fs.vasp之间插入若干个image，这会导致其中某些image的构型极度不合理，比如某两个原子重叠在了一起，需要调整。可以自己把01-08文件夹中的POSCAR文件下载下来，调整后再上传，但比较麻烦。推荐使用我写的脚本nebinfo（来源：Mabinogiysk/VASP-script），输入image的数量之后，会生成一个压缩包NEBINFO.tar.gz. 里面包含有每个image的vasp格式文件，gjf格式文件以及xyz格式的动画文件（可用VMD观看）。

调整好image之后在当前文件夹执行vasp命令就可以开始NEB计算。

刚开始算的时候，由于构型不合理，因此选用最速下降法会比较快，也不容易出错。此时选择IOPT=7，跑个二三十步之后再换用IOPT=1，当然也可以从头到尾只用IOPT=1.

对于复杂的过渡态，NEB计算需要非常多的image才能收敛。因此我一般只跑个五六十步，根据经验来看，后面基本不会有太大变化了。

之后可以采用vtst tool中的脚本来输出最后一步NEB的信息，如能量、切线力等，用法请到官网查找。也可以使用我脚本仓库中的nebvtst.pl脚本：

1. > nebvtst.pl 8

复制代码

该脚本将输出NEB计算每一步每个image的信息：

1. steps: 60
   
2. images: 1    -314.70152    -0.20180     0.07899
   
3. images: 2    -314.66766    -0.10107     0.08997
   
4. images: 3    -314.21920    -0.03095     0.06242
   
5. images: 4    -314.82366     0.36034     0.08089
   
6. images: 5    -315.02908     0.17184     0.10215
   
7. images: 6    -315.14834     0.10812     0.06171
   
8. images: 7    -315.21699     0.05799     0.05365
   
9. images: 8    -315.25374     0.02400     0.12696

复制代码

三列数字从左到右分别是能量，切线力和受力最大的原子所受的力（都小于EDIFFG时NEB计算就收敛了）。

其中最重要的是切线力，切线力为零的点就是过渡态。正常的NEB计算切线力应该是从负值到正值或者从正值到负值，并且能量应该先升高再下降。不满足这两个条件时意味着计算很可能出错了。

对于上面这个计算，可以看到image3的切线力很接近0了，也就意味着这个image跟过渡态很接近了，可以把它单独拿出来做个优化（IBRION=1），有可能就收敛到过渡态了。

如果所有的image的切线力都比0差很多，如下面这样：

1. steps: 42
   
2. images: 1    -314.71090    -0.19732     0.67034
   
3. images: 2    -314.68128    -0.10245     0.18674
   
4. images: 3    -314.22095    -0.23067     0.47518
   
5. images: 4    -314.80703     0.24911     0.40551
   
6. images: 5    -315.01767     0.16162     0.27819
   
7. images: 6    -315.13853     0.10266     0.38136
   
8. images: 7    -315.21270     0.04747     0.28035
   
9. images: 8    -315.25231     0.02491     0.24664

复制代码

就只能把零两边的image，即image3和image4拿出来重新用nebmake.pl插值再跑NEB，直到找到切线力接近零的点。这是原始版本的vasp做过渡态的方法，如果编译了vtst tool工具，采用climbing image方法，即设置LCLIMB = .TRUE. 基本上跑一遍NEB之后都会有一个image的切线力接近零，避免了重复插值。因此用vasp做过渡态一定要编译上vtst tool.

另外经常遇到的情况是切线力接近零的image拿出来做优化之后，发现力总是降不下去。此时NEB方法就无能为力了，需要另一个强大的过渡态搜索方法Dimer上场了。这就留做下一篇帖子来说了。

tyut-siggma

请问一下在用IOPT= 7跑了二三十步修改IOPT = 1是指，这个IOPT=7任务还没有计算完，不需要停止，直接修改INCAR文件中的IOPT=1吗？

ljb874722957

Leix 发表于 2024-7-6 10:31
**** 作者被禁止或删除 内容自动屏蔽 ****

您好，我最近在计算cneb路径也散了，这个问题您解决了吗？

Leix

追翔少年 发表于 2019-3-16 19:45
个人观点，我用VTST搜过渡态结果会很奇怪，后来经过人指点还是用VASP自带的CNEB工具才有所改善，我觉得VTST ...

我想问一下是什么工具呢，我在计算CNEB时，分子跑散了

2422759918

nebvtst.pl这个脚本不是VTST自带的啊，请问楼主怎么获取呢

Nora

您好，我想问一下为什么计算之后发现CI-NEB不按照原来制定的路径走？

lyzhang

您好，请问过渡态计算收敛，但是没有虚频是为什么呢？增加收敛精度到E-7和-0.001后，有的结构出现了较小虚频，但是有的结构依然没有虚频，而且对照文献了路径没问题，

QH1995

能不能简单介绍一下IDPP插点法呢？目前我就是在IDPP插点法插入8个点，然后精度比较低（EDIFF = 1.0E-5, EDIFFG = -0.8）的粗算，再能量最高的两个点之间再插入5个点进行精度稍微高一点的计算（EDIFF = 1.0E-5, EDIFFG = -0.5），得到能量最高的那个点基本上就是我想要的过渡态

yydy

多次拒绝斯佳丽 发表于 2020-12-30 22:06
麻烦请问一下，vasp自带的cneb就是不用vtst吧？那是不是就没必要设置LCLIMB=T呢？如果这样的话，还是CINE ...

设置CHAIN= 0时，CLINEB默认为T;要用NEB的话应该要改为F

fight

sky 发表于 2017-12-8 13:48
跟NEB一样就行了

请问楼主，这个image跟过渡态很接近了，可以把它拿出来跑一个IBRION=1，那么INCAR的其他参数和结构优化一样还是和跑NEB的参数一样呀，不胜感激！

夙雪夜归

Q-Chembio-llg 发表于 2021-4-28 17:59
请问如何调用VASP自带的NEB呢？我们的VASP在学校服务器上，没法改动谢谢

可以看一下官网给的这个例子，自带的NEB其实就多了images和spring两个参数https://www.vasp.at/wiki/index.p ... t_adatom_on_fcc-Pt_(001):_Nudged_Elastic_Band_Calculation
插值的话可以直接用vaspkit的505功能

Q-Chembio-llg

请问如何调用VASP自带的NEB呢？我们的VASP在学校服务器上，没法改动谢谢

CCCCCCLEAR

请问sky老师，切线力是使用nebbarrier.pl所得到力的数值，而原子所受的力是nebef.pl所得到的数值吗？因为我们实验室目前服务器暂时不能装python3，您的nebvtst.pl脚本无法使用，只能使用vtst自带脚本，想区分一下您在此脚本中得到的切线力与原子所收到的力这两个概念。若能得到您的解惑将给我很大帮助，感谢！

多次拒绝斯佳丽

追翔少年 发表于 2019-3-16 19:45
个人观点，我用VTST搜过渡态结果会很奇怪，后来经过人指点还是用VASP自带的CNEB工具才有所改善，我觉得VTST ...

麻烦请问一下，vasp自带的cneb就是不用vtst吧？那是不是就没必要设置LCLIMB=T呢？如果这样的话，还是CINEB吗？我感觉就成了NEB？还是您的意思只是优化方法用vasp优化器而不是vtst优化器呢？谢谢，因为我师弟之后会算到旋转的有机分子集团，特此问一下您

张丽

好帖，先码了