NEB+dimer计算过渡态经验-无VTST篇

夜夜衣

有时间更新vtst篇的 自带的容易被质疑 就不误人子弟了

下午茶

您好！我刚接触vasp,看了您的帖子，有几个问题想向您请教1. 这是我优化的初态和末态，算NEB时需要初态和末态的坐标文件和OUTCAR，原子要一一对应，是需要算初态时就把H原子加上吗？我在这个地方有点混乱，希望您能给解答一下，谢谢！

ghifi37

下午茶 发表于 2019-12-13 15:19
您好！我刚接触vasp,看了您的帖子，有几个问题想向您请教1.这是我优化的初态和末态，算NEB时需要初态和 ...

原子数目肯定需要一致的，不然能量都没法比较，至少也得分开算

阿宋啦

你好，我想请问一下，如果我算的是表面的反应，那基地在算dimer的时候要用F关闭，只打开参与反应的原子吗？

夜夜衣

阿宋啦 发表于 2020-1-6 11:37
你好，我想请问一下，如果我算的是表面的反应，那基地在算dimer的时候要用F关闭，只打开参与反应的原子吗？

这得看你计算频率固定的哪些原子，那这些固定就没有震动，xyz三个方向都是0，dimer不固定应该也没事

小枫叶

楼主，您好，我看了您的帖子帮助蛮大，就是我不知道为何用您的方法，参数设置，却总是在第一步刚要跑完时却出现这种报错ERROR: missing or invalid vector defining dimer,一直没解决这个问题，比较郁闷，这个为什么啊？谢谢您

夜夜衣

小枫叶 发表于 2020-2-24 21:41
楼主，您好，我看了您的帖子帮助蛮大，就是我不知道为何用您的方法，参数设置，却总是在第一步刚要跑完时却 ...

应该是你那三列频率没弄好，安装了vtst和没有vtst是有差别的，没有vtst是放在POSCAR坐标下面，安装了vtst是尤文专门的文件放频率

Siling_S

有个问题想请教一下，能不能在粗算完NEB之后，用modemake.pl算出来的初始MODECAR作为初始的搜索方向放到POSCAR里面算DIMER呢？？这样可以省了算频率的时间。（体系有点大，算频率的时间可能比较长……）另外楼主算TS的INCAR里的ICHAIN = 0语句应该是编译了VTST之后才能调用的吧？放在这里是否直接会在计算过程中被忽略掉？

夜夜衣

下午茶 发表于 2019-12-13 15:19
您好！我刚接触vasp,看了您的帖子，有几个问题想向您请教1.这是我优化的初态和末态，算NEB时需要初态和 ...

不好意思，现在回复你不知道还能不能帮助到你，IS和FS的原子数类型必须是一模一样的，针对你的体系，IS应该在加一个H原子，尽量靠近和slab接触的O，这样能让H由最简单的路径到FS，我猜你是做电催化，barrier的计算我们是这样处理的，是你在neb的过程中去除TS的能量减去 *CO2减去1/2H2就是barrier，所以说IS我们只需要用它的构型，不要他的能量，

下午茶

夜夜衣 发表于 2020-4-12 15:51
不好意思，现在回复你不知道还能不能帮助到你，IS和FS的原子数类型必须是一模一样的，针对你的体系，IS应 ...

感谢您的回复，我不是太明白您说的“是你在neb的过程中去除TS的能量减去 *CO2减去1/2H2就是barrier，”这句是什么意思？你能再详细的讲一下吗？万分感谢

夜夜衣

下午茶 发表于 2020-4-14 22:27
感谢您的回复，我不是太明白您说的“是你在neb的过程中去除TS的能量减去 *CO2减去1/2H2就是barrier，”这 ...

我打字打错了，正常算barrier是Eb=E（TS）-E(IS)是吧？ 但是电催化里面有两种机制 一种是H吸附在表面之后再加氢，这样计算能垒就是上面这个公式。但是电催化里面，H离子可以认为是游离的，那么能垒就成了Eb=E（TS）-E（catalyst）-1/2H2    因为标准氢电极电势下H+ + e- = 1/2H2的能量，所以我们只用了IS的构型来找过渡态 但是不用他的能量

下午茶

夜夜衣 发表于 2020-4-17 22:08
我打字打错了，正常算barrier是Eb=E（TS）-E(IS)是吧？ 但是电催化里面有两种机制 一种是H吸附在表面之后 ...

好的，谢谢您

小枫叶

夜夜衣 发表于 2020-2-26 10:36
应该是你那三列频率没弄好，安装了vtst和没有vtst是有差别的，没有vtst是放在POSCAR坐标下面，安装了vtst ...

感谢您的回复，我还有两个问题想请教一下您。
1。就是我想问一下我把MODECAR单独拿出来之后，第一步结束倒是没有报错，但是却不生成dimcar,并且一直不收敛，但是cl-neb却可以用，这种情况也是VTST没有编译好吗？
2.我在用dimer法时，并没有等到整体的插点收敛，就直接选取切向力接近0的结构计算频率，发现只有一个虚频，接着就用此结构跑dimer，这种方法可以吗？谢谢您了。

夜夜衣

小枫叶 发表于 2020-5-22 23:17
感谢您的回复，我还有两个问题想请教一下您。
1。就是我想问一下我把MODECAR单独拿出来之后，第一步结束 ...

回复你这个问题我其实有点纠结，我算dimer其实也很不好收敛，我都是在计算多少步之后把结构挑出来算频率，看震动，如果符合一个大虚频 振动方向也对，不管dimer计算收敛不收敛我就停止计算了，但是现在我一般都是用CINEB计算过渡态，我觉得还比较好用，现在比较少用dime了。你第一个问题，你既然能用neb.ef想来你应该是安装了VTST的，你用nebef脚输出的各列都应有数，特别是切线力，而不是都为零。至于你说dimercar没出来，肯定是你的输入文件的问题。第二个问题，我以前写的这个还是比教片面，近期我会根据我的经验重新写一部分vtst计算过渡态的部分，可不用dimer，也能较为精准的定位过渡态，最重要的是nebspline.pl这个绘制出来的势能面图得好看 也就是一个漂亮的山峰型，并且力收敛，才能真正确定是过渡态，若是那种波浪型，那肯定是错的，。基元反应的势能面应该是很好看的，要么INACAR不行，要么 你的IS和FS没设计好。

小枫叶

夜夜衣 发表于 2020-5-29 21:20
回复你这个问题我其实有点纠结，我算dimer其实也很不好收敛，我都是在计算多少步之后把结构挑出来算频率 ...

感谢楼主回复，经过测试，发现的确是因为INCAR参数丢失一个，目前DIMER可以跑起来了。搬好板凳蹲一波楼主的新作，嘿嘿。