计算化学公社

标题: 关于NEB+DIMER方法搜寻过渡态过程中一直没有收敛趋势的求助 [打印本页]

作者
Author: 奶盖盖 时间: 2023-9-19 13:34
标题: 关于NEB+DIMER方法搜寻过渡态过程中一直没有收敛趋势的求助
本帖最后由奶盖盖于 2023-9-20 11:07 编辑

各位老师好，最近我想要计算一下Zn基催化剂丙烷脱氢反应最后一步生成氢气的过渡态，在用CINEB粗收敛得到一个可能的结构后，转用VTST提供的DIMER方法进行精收敛搜索，但是在读取DIMCAR时发现Curvature在后期会突然变正。感觉没有任何能够正常收敛的趋势。NEB能量趋势已贴出，完全不知道从何改起。如图是我用来跑DIMER的POSCAR结构，INCAR以及DIMCAR的输出文件也已贴出，初态结构为POSCAR1，末态结构为POSCAR2，恳请老师们能够帮忙查看一下我的问题所在，万分感谢。
  0       0.071948    -2283.443800       0.000000
1       0.659017    -2283.265200       0.178600
2       0.218020    -2283.108400       0.335400
3       0.394007    -2282.986200       0.457600
4       0.204644    -2282.632100       0.811700
5       0.434662    -2281.890800       1.553000
6       0.256942    -2282.619800       0.824000
7       0.443042    -2282.845900       0.597900
8       0.159587    -2282.878800       0.565000

(, 下载次数 Times of downloads: 4) (, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 10) (, 下载次数 Times of downloads: 3)

作者
Author: ahxb 时间: 2023-9-19 17:34
你给的INCAR_DIMER里用不是DIMER方法，而是NEB。先给出正确的用DIMER的INCAR。

作者
Author: sobereva 时间: 2023-9-20 03:53
第一性原理程序的问题不要发到量子化学版块，这次给你移动了，以后注意

作者
Author: wypkdhd 时间: 2023-9-20 08:56
你的初态和末态难道不应该是 00和08吗？看你给的一个0.071948 的力，一个 0.159587 的力，这样的做丙烷脱氢是不行的。别说你incar错了，你结构优化力最高不能超过0.05，一般都是0.02。

作者
Author: 奶盖盖 时间: 2023-9-20 11:09

ahxb 发表于 2023-9-19 17:34
你给的INCAR_DIMER里用不是DIMER方法，而是NEB。先给出正确的用DIMER的INCAR。

好滴谢谢，我贴成NEB的了，已经更改啦，可以继续检查啦感谢

作者
Author: ahxb 时间: 2023-9-20 15:43

奶盖盖发表于 2023-9-20 11:09
好滴谢谢，我贴成NEB的了，已经更改啦，可以继续检查啦感谢

为何把DRotMax设置为1？Dimer方法需要尽量准确的Dimer方向用于估计曲率，设为1相当于不优化Dimer方向，所以完全估计不准。为了平衡精度和耗时，一般设为3~5。

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