计算化学公社

标题: 关于VASP找过渡态一直以来很苦恼的问题，求助！ [打印本页]

作者
Author: nosugar 时间: 2023-12-24 15:56
标题: 关于VASP找过渡态一直以来很苦恼的问题，求助！
各位老师，关于找过渡态有一些问题想请教，困扰本人很久了。在用VASP找过渡态的过程中，用CI-NEB的方法，先进行粗收敛，计算结束后进行精收敛，然后此过程发现力一直不收敛。于是采用Dimer的方法，把ci-neb粗收敛结构中像过渡态的结构拿出来做频率分析，发现有一个虚频，且震动方向也是朝向断键的方向，然后再拿着这个结构去跑dimer，就出现力难收敛的状况；尝试用IDM方法跑，结构则是出现了散架的情况。请问有没有大神能够解答！

举个例子下面为最近计算过渡态，结构为ZnO,五层，上面两层放开，下面三层固定，反应为：H20→H+OH。初末态dist.pl的返还值为4.86，插点数为四个，根据返还值/0.8少插了两个(这个会对结构造成影响吗？)，用的脚本idpp.py非线性插点。
nebef.pl结果：
(, 下载次数 Times of downloads: 14)
INCAR
(, 下载次数 Times of downloads: 19)

(, 下载次数 Times of downloads: 17)
IS
(, 下载次数 Times of downloads: 17)
FS
(, 下载次数 Times of downloads: 15)

作者
Author: vasp新新手1 时间: 2024-1-5 16:10
你好，请问IDM方法计算，结构出现了散架的情况怎么解决？

作者
Author: glb88 时间: 2024-1-6 11:07
可以固定一些远离反应区域的原子，另外多插一些点

作者
Author: nosugar 时间: 2024-1-31 14:37

glb88 发表于 2024-1-6 11:07
可以固定一些远离反应区域的原子，另外多插一些点

感谢您的回复我尝试一下

作者
Author: wypkdhd 时间: 2024-1-31 16:25
NCORE=你使用的核数/2
KPAR=2
PREC=Normal
ADDGRID注释掉
使用vtst脚本生成mep.eps看看cineb的曲线如何
最高点拿出来做dimer
做dimer使用IOPT=3
TIMESTEP=0.02
具体参数大小写去vtst官网看，我记不清了
http://theory.cm.utexas.edu/vtsttools/scripts.html

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