| 查看: 4001 | 回复: 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方法计算好像不合理啊，是不是应该分两个过程算呢？ 谢谢了！  |
回复此楼» 猜你喜欢
申请科研项目纯靠实力能成功吗？
已经有9人回复
-
40+中年人失业后，还能找到好工作吗？
已经有4人回复
-
极限拉扯，如何抉择
已经有15人回复
-
中科院海洋所招收海洋材料相关领域博士研究生及研究助理
已经有4人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有4人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
wyf198166
至尊木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 19322.9
- 红花: 4
- 帖子: 1403
- 在线: 47.1小时
- 虫号: 1481340
- 注册: 2011-11-07
- 性别: GG
- 专业: 金属结构材料
★ ★
franch(金币+2): 谢谢回帖交流 2011-12-13 09:47:23
franch(金币+2): 谢谢回帖交流 2011-12-13 09:47:23
| 如果路径你知道的话，不需用分几步做，把中间的IMAGES多设几个，只要找到鞍点就可以了，但如果对称性不高，路径不清楚的话，建议使用CLNEB，你这个INCAR的参数也是写的CLNEB的，而你上面却写着NEB计算，都不知道你到底用的那个？而且我同意楼上的说法，你这个吸氧过程说不定是个稳态，NEB主要是研究非稳态的，你这个首先最好先验证下。 |
gaofenglili
8楼2011-12-12 10:49:17
likai0106
金虫 (小有名气)
- 1ST强帖: 2
- 应助: 38 (小学生)
- 金币: 1352.4
- 帖子: 233
- 在线: 68.5小时
- 虫号: 634278
- 注册: 2008-10-23
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-12-09 08:12:02
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 ：电子结构
★
贺仪(金币+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朵）
5楼2011-12-10 09:16:19
200
