| 查看: 3888 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
gaofenglili银虫 (小有名气)
|
[求助]
新手请教用NEB方法计算energy barrier.
|
|
|
想用NEB方法算体系的energy barrier, 请教各位大侠： 1. 我用nebmake.pl脚本产生中间状态的Image, 但是结构优化后都变成了几乎同一种最终结构，好像中间插入的各个结构之间没有相互作用和约束一样，得不知道是怎么回事，我的INCAR如下，请不吝指教： PREC = Mediate ENCUT = 500 EDIFF = 1E-5 EDIFFG = -0.05 ISMEAR = 0 SIGMA = 0.1 ISPIN = 2 LPLANE = .True. LWAVE = .False. LCHARG = .False. NPAR = 4 LSCALU = .False. ICHAIN = 0 IMAGES = 6 LCLIMB= .TRUE SPRING = -5.0 IBRION = 2 NSW = 500 POTIM = 0.1 NSIM = 4 ALGO = Fast LREAL = Auto 另外我的POSCAR1和POSCAR2中都写了selective dynamics, 然后让原子可移动，这样对吗？ 2. 对于初始结构应该单独进行计算吗？ 3. 我计算的体系是一个H原子从无限远接近吸附在表面的一个O，然后形成OH后又远离表面，这好像是两个过程啊，这种情况直接给出初态和末态位置用NEB方法计算好像不合理啊，是不是应该分两个过程算呢？ 谢谢了！  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有9人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有10人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- castep计算完能带之后如何挑选有用的数据画出清晰的能带图
已经有16人回复
- 请教FT-IR测试硅薄膜数据中基线处理的问题
已经有6人回复
- 反转录之后PCR扩增片段问题的请教，急急急
已经有5人回复
- 请教NEB收敛判据中能量的设置
已经有6人回复
- 【求助】请教CNEB
已经有15人回复
- 【求助】有人用dl_poly计算过solvation free energy吗？请教几个问题。。
已经有10人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【其他】vasp neb 计算讨论
已经有20人回复
- 【求助成功】NEB计算的一个stdout输出怎么解读?
已经有6人回复
- 【求助成功】怎么样提高NEB计算速度?
已经有19人回复
- 【求助/交流】请教NEB的pmeI的酶切体系
已经有7人回复
gaofenglili
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 126
- 散金: 8
- 红花: 1
- 帖子: 295
- 在线: 18小时
- 虫号: 307162
- 注册: 2006-12-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2011-12-10 00:00:31
likai0106
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 38 (小学生)
- 金币: 1352.4
- 帖子: 233
- 在线: 68.5小时
- 虫号: 634278
- 注册: 2008-10-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-12-09 08:12:02
likai0106
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 38 (小学生)
- 金币: 1352.4
- 帖子: 233
- 在线: 68.5小时
- 虫号: 634278
- 注册: 2008-10-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★
贺仪(金币+1): 多谢你的热心帮助 2011-12-14 09:37:08
贺仪(金币+1): 多谢你的热心帮助 2011-12-14 09:37:08
|
不好意思，我可能理解错误了。我认为你的过程应该是两个，H吸附形成OH是一个，OH脱离又是一个，分开做就好了。如果你有催化剂的话必须考虑到H的物理吸附和O成键的过程，所以在催化剂下可能就是三个了，第一个是H吸附，之后H和O成键形成OH。最后OH脱离。 |
» 本帖已获得的红花（最新10朵）
gaofenglili5楼2011-12-10 09:16:19
gaofenglili
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 126
- 散金: 8
- 红花: 1
- 帖子: 295
- 在线: 18小时
- 虫号: 307162
- 注册: 2006-12-16
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2011-12-11 11:15:29
