| 查看: 742 | 回复: 4 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】大家帮忙看一看这个NEB的INCAR有啥问题么？ 已有4人参与
|
|||
|
ISTART = 0; ICHARG = 2 ENCUT = 400 ISMEAR = 2 * LREAL=FLASE NSW=60 IBRION=2 EDIFFG=-0.01 *ISIF =3 PREC=High NGX = 58 NGY = 58 NGZ = 56 NBANDS = 50 IMAGES=8 SPRING=-5 想算单空位的迁移能 算出来是40多电子伏 我手动让空位走直线算出来的才0.5电子伏。。。 0 0.00000000 -127.53892500 0.00000000 1 0.00000000 -117.35116000 10.18776500 2 0.00000000 -105.60393500 21.93499000 3 0.00000000 -90.23796600 37.30095900 4 0.00000000 -85.14998300 42.38894200 5 0.00000000 -96.81678600 30.72213900 6 0.00000000 -113.40630700 14.13261800 7 0.00000000 -127.51704500 0.02188000 |
回复此楼» 猜你喜欢
博士延得我，科研能力直往上蹿
已经有7人回复
-
退学或坚持读
已经有27人回复
-
面上基金申报没有其他的参与者成吗
已经有5人回复
-
有70后还继续奋斗在职场上的吗？
已经有5人回复
-
遇见不省心的家人很难过
已经有22人回复
-
多组分精馏求助
已经有6人回复
4楼2010-04-18 21:51:33
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
2楼2010-04-17 14:59:53
zengxiangm
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1981.2
- 帖子: 221
- 在线: 96.9小时
- 虫号: 676900
- 注册: 2008-12-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
dazao(金币+1):谢谢参与
qli0218(金币+1):谢谢 2010-04-17 15:21
dazao(金币+3):内牛满面 2010-04-18 22:01
dazao(金币+1):谢谢参与
qli0218(金币+1):谢谢 2010-04-17 15:21
dazao(金币+3):内牛满面 2010-04-18 22:01
|
IBRION=2 这个设置好像有问题 The NEB and min-mode following (dimer/Lanczos) saddle point finding methods use a force projection in order to direct the optimizers towards minimum energy paths and saddle points. These modifications to the force mean that the energy is no longer consistent with the force being optimized. Because of this, we can only use optimizers that are solely based upon the force (and not the energy). The quasi-Newton and quick-min (IBRION=1 and 3 respectively) optimizers that are built into VASP are both force-based, but the conjugate-gradient method (IBRION=2) is not. IMAGES最好设置成奇数，我最近也在学NEB，有时间多交流啊 |
3楼2010-04-17 15:19:23
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
dazao(金币+3):谢谢指点 内牛满面啊 2010-04-19 16:20
qasd(金币+3):谢谢 2010-04-19 23:47
小木虫(金币+0.5):给个红包，谢谢回帖交流
dazao(金币+3):谢谢指点 内牛满面啊 2010-04-19 16:20
qasd(金币+3):谢谢 2010-04-19 23:47
|
1.楼主计算的INCAR中IMAGE=8,但是楼主的给出的结果中间的IMAGE只有6个，这好想有点问题吧！！呵呵！ 2.楼主计算NEB时最好不要用CG弛豫（即IBRION=2），请参考三楼的说法。呵呵！！ 3.还有一点就是最好不要使用全弛豫（楼主似乎没有使用啊），如果要使用全弛豫的话，楼主必须保证中间的IMAGE的构型是连续变化，不能产生突变，否则要考虑自由度造成的影响。楼主计算的如果是金属体系的话，使用全弛豫应该是没有问题的。呵呵！！ 4.三楼说的把IMAGES设置成奇数的问题，这个好像也不一定，如果楼主单纯使用NEB的话，设成奇数还是要好些，如果使用CI-NEB，那就不用考虑这方面的问题了。呵呵！ 以上拙见供楼主参考！ NEB作者Graeme Henkelman教授组织创办的论坛：http://theory.cm.utexas.edu/forum/index.php |
5楼2010-04-19 13:04:16