| 查看: 4907 | 回复: 6 | |||
[求助]
用vasp计算，neb方法找过渡态，能垒高的构像老是跑偏
|
|
用vasp自带的neb方法找过渡态，乙腈水解（为的是节省计算时间），目的是熟悉一下vasp找过渡态的方法，以便于用vasp研究目前课题的反应机理。 第一个INCAR： SYSTEM=O the first job of majun ISMAER=0 SIGMA=0.01 ENCUT=400 NSW=500 IBRION=3 POTIM=0.05 SPRING=-5 IMAGES=8 EDIFF=1E-05 EDIFFG=-0.05 ISPIN=1 计算开始后第四个image很快突出，并且越来越突出，可是后来第三个点的能量起来了，超过了第四个点，到了后来第三个点的能量相对于第二个点也不突出了。感觉不对，杀掉文件。 在第二个INCAR在incar-1的基础上修改参数IMAGES=4（提供了新的构象文件）; EDIFFG=-0.1; 增加参数LCLIMB=.TRUE. 计算开始后第三个点很快突出，并且越来越突出，可是后来第三个点的能量下降，第二个点能量增加，并超过了第三个点，然后一直处于领先优势。为啥？！用vesta看发现时第三个点自己构想明显朝着末态方向变化了，堕落了。能量也是越来越接近最后一个image点的能量。这个时候发现第二个点的能量也在下降，构像也有变化。最后第二个image点的能量越来越低，结构越来越接近初始态结构，失败了。 第三个INCAR在第二个incar-2的基础上，增加参数：IOPT=3 吃饭回来后，发现nebdf.pl后，输出的能量都接近零了，失败了。 一直在用Gaussian，这是第一次用vasp找过渡态，请高手帮忙看看是怎么回事啊？谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP and MS |
» 猜你喜欢
广东省 085601 329分求调剂
已经有6人回复
-
求调剂
已经有18人回复
-
调剂 化学 307
已经有17人回复
-
297求调剂
已经有8人回复
-
材料类284调剂
已经有45人回复
-
368求调剂
已经有3人回复
-
295分求调剂
已经有5人回复
-
263能源动力专硕求调剂
已经有3人回复
-
272分材料子求调剂
已经有28人回复
-
一志愿西北工业大学289 085602
已经有30人回复
hakuna
木虫 (知名作家)
- 1ST强帖: 3
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
majun042楼2017-11-07 17:22:22
3楼2017-11-07 17:32:33
|
以前是用gaussian计算均想反应机理的。因为gaussian计算块体不是强项，周期性的块体表面也解决不了。现在在学vasp，期望解决块体表面催化反应机理。这两天很困惑，我选择得体系都是比较普通的，难道是有什么问题？还有因为刚刚开始用vasp，服务器上没有进行任何编译，现在想着先走通一个过渡态，感觉一下vasp计算反应机理是怎么回事。 现请教大家： 第一， 用vasp找过渡态，如果在比较明确过渡态的结构了，是不是用dimer方法比neb方法要节省时间？ 第二， 那么用dimer方法时候，需要MODECAR，这个modecar又需要两个POSCAR, 这两个POSCAR分别对应与一个是反应初始态构象，一个是你设计的过渡态构象？ 第三， 对于vasp而言如果自带的neb不好用，是不是用cl-neb是必须要编译，这个cl-neb不是vasp开发的，是vtst开发的，这个需要版权吗？是不是网上随便下载就可以拿来用。 第四，如果用dimer方法找过渡态，是不是也得先对vsap进行编译？ 第五，前天把以前用gaussian算过得过渡态，有机小分子得，用vasp自带得neb方法找，images=8，都找不到。考虑到vasp不善于计算有机小分子，昨天用vasp自带的neb计算金属Pd分解氢气得过渡态NSW=250，都没有给出结果。 |
4楼2017-11-09 09:54:06
5楼2017-11-09 15:40:41
hakuna
木虫 (知名作家)
- 1ST强帖: 3
- 应助: 82 (初中生)
- 金币: 8735.9
- 散金: 158
- 红花: 126
- 帖子: 5599
- 在线: 1367.2小时
- 虫号: 167622
- 注册: 2006-01-15
- 专业: 考古理论
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
majun04(franch代发): 金币+15, 谢谢回帖交流。 2018-01-04 11:08:01
majun04(franch代发): 金币+15, 谢谢回帖交流。 2018-01-04 11:08:01
|
1.不建议用VASP自带的NEB，找TS命中率不理想（VASP自带的TS搜索方法只有NEB和DIMER，其它都是第三方的） 2.建议CINEB，或VASP新版自带的DIMER（和VTST的dimer不同，这个手册里有详细使用说明，不需要MODCAR！） 3.是否加入VTST，你可以查一下编译时源代码里有没有VTST里的文件，或问一下编译的人，或按CINEB跑一下，查OUTCAR里关于NEB的信息，是否有CINEB 4.NEB类的方法需要并行，且计算核心数一定是插入image的倍数；dimer、lanczos(VTST)不需要并行，单个核心就可以算 5.用VTST的TS搜索方法，要好好消化一下Henkelmann网站上的说明，采用力收敛标准时，不同TS优化方法有不同参数组合 6.Henkelmann也分享了许多python脚本，很实用，他网站上可直接下载 7.网上有许多CINEB使用的例子，可以下来参考 8.Henkelmann主页上有一个论坛链接，注册一下，可以看到许多VTST的讨论 |
6楼2017-11-19 20:43:27
7楼2017-12-26 15:20:23
