计算化学公社

标题: NEB计算总是能量跳变 [打印本页]

作者
Author: ikea1984 时间: 2023-5-11 20:02
标题: NEB计算总是能量跳变
用VTST和VASP编译，然后做CI-NEB。发现某些image总是能量隔若干步就会突然升高2、3eV，检查结构发现此时有个别键长出现跳变，从而导致整个NEB过程不能收敛。这种问题应该怎么处理呢？如果用dimer，有些计算可以找到一个过渡态，但是怎么判断这个过渡态是不是想要的过程仍然有点困难。

作者
Author: Deepast 时间: 2023-5-12 00:26
跳变大的话估计是你插点原子顺序不对
找到过渡态的话，算个频率验证就行
根据我的经验，稍微复杂一点的化学反应neb收敛都很慢，甚至不收敛

作者
Author: ikea1984 时间: 2023-5-12 05:54

Deepast 发表于 2023-5-12 00:26
跳变大的话估计是你插点原子顺序不对
找到过渡态的话，算个频率验证就行
根据我的经验，稍微复杂一点的化 ...

嗯，您说的有道理。但是我这个反应是一个有机反应，不是那种简单的扩散或者单个键断裂，感觉原子在空间中存在一定的重排，这种很难认定原子顺序一定是怎么样。能否把同类原子按可能顺序都排列一下，然后选择ini和final的dist最小的映射关系？

作者
Author: Deepast 时间: 2023-5-12 08:41

ikea1984 发表于 2023-5-12 05:54
嗯，您说的有道理。但是我这个反应是一个有机反应，不是那种简单的扩散或者单个键断裂，感觉原子在空间中 ...

距离最短不一定是你想要的，而且好像也没现成调整原子顺序的脚本
原子顺序一样是建模的时候就需要考虑的问题，插完点后记得把结构都看一下
一般来讲的话，后面加的原子会出现在文件末尾，可以利用这一点

欢迎光临计算化学公社 (http://bbs.keinsci.com/)