| 查看: 5202 | 回复: 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 |
» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
5楼2017-11-09 15:40:41
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
