| 查看: 4010 | 回复: 20 | ||||||||
zhangfan2192金虫 (著名写手)
|
[交流]
【其他】vasp neb 计算讨论 已有11人参与
|
|||||||
|
希望各位用vasp的虫友能够讨论下neb计算中的注意事项： 不知道为什么，我在用neb计算过度态的时候总是感觉计算结果不是很理想，很难找到自己需要的结果，比如说计算的8个image的能量是单调上升或下降的，而不是抛物线形状，不知道大家在计算的时候遇到这样的情况没，如果遇到的话大家都是怎么解决的？相信现在用vasp软件的不在少数，欢迎大家积极讨论建议！ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP | first principle | 第一性原理和电化学 | vasp计算总结 |
| 计算资料收集byFrank |
» 猜你喜欢
材料求调剂
已经有5人回复
-
354求调剂
已经有9人回复
-
化学工程321分求调剂
已经有22人回复
-
考研调剂
已经有3人回复
-
326求调剂
已经有8人回复
-
333求调剂
已经有7人回复
-
一志愿东华大学控制学硕320求调剂
已经有3人回复
-
考研化学学硕调剂，一志愿985
已经有7人回复
-
0703化学调剂
已经有15人回复
-
0703化学调剂 ，六级已过，有科研经历
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp计算过渡态问题
已经有13人回复
- 新手请教用NEB方法计算energy barrier.
已经有15人回复
- vasp用neb计算迁移势垒
已经有11人回复
- vasp NEB 如何选取原胞大小 ？
已经有10人回复
- vasp 过渡态计算 线性插值方法。
已经有8人回复
- 【求助】vasp计算过渡态结果出错
已经有17人回复
- 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析
已经有30人回复
- 【求助】关于vasp NEB
已经有8人回复
- 【求助】NEB算过渡态出现致命错误
已经有5人回复
- 【求助成功】怎么样提高NEB计算速度?
已经有19人回复
- 【讨论】NEB寻找过渡态专题讨论
已经有49人回复
uuv2010
荣誉版主 (职业作家)
- 1ST强帖: 6
- 应助: 160 (高中生)
- 贵宾: 3.14
- 金币: 3681.7
- 散金: 5802
- 红花: 79
- 沙发: 35
- 帖子: 4560
- 在线: 916小时
- 虫号: 966216
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
- 管辖: 第一性原理
3楼2010-07-06 10:20:14
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3573.1
- 散金: 2141
- 红花: 37
- 帖子: 1432
- 在线: 681.4小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
8楼2010-07-07 12:21:50
18楼2010-11-20 21:29:59
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
2楼2010-07-05 22:16:51
zhangfan2192
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1007.3
- 散金: 62
- 红花: 3
- 帖子: 1170
- 在线: 171.3小时
- 虫号: 648132
- 注册: 2008-11-07
- 性别: GG
- 专业: 理论和计算化学
4楼2010-07-06 14:10:44
fanchen021
铁杆木虫 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 7994.9
- 散金: 1844
- 红花: 4
- 帖子: 776
- 在线: 569小时
- 虫号: 773188
- 注册: 2009-05-18
- 专业: 理论和计算化学
5楼2010-07-07 10:30:40
zhangfan2192
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1007.3
- 散金: 62
- 红花: 3
- 帖子: 1170
- 在线: 171.3小时
- 虫号: 648132
- 注册: 2008-11-07
- 性别: GG
- 专业: 理论和计算化学
|
INCAR文件如下： SYSTEM = Sul-Fe Startparameter for this Run: NWRITE = 2 verbosity write-flag LPETIM= F write-flag & timer # ISTART = 0 0-new 1-cont 2-samecut job ISPIN = 2 1-no, 2-yes INIWAV = 1 0-jellium WF; 1-random numbers (recomm.) Electronic Relaxation ! ENCUT = 400.00 #Ecut (eV) PREC = Medium * #precission: #Low/Medium/High NELM = 200 max. # of electronic SC steps NELMDL = 5 # of non-consistent steps at the beggining:>0 or <0 EDIFF = 1.0e-4 stopping-criterion for electronic upd. VOSKOWN = 1 ISMEAR = 1 SIGMA = 0.1 ! AMIX = 0.2 ! BMIX = 0.0001 Ionic relaxation EDIFFG = 1.0e-3 stopping-conditions for ionic relax.loop (EDIFF*10) NSW = 200 * # of steps for ionic update (def:0) IBRION = 2 * ionic relax: -1:no move; 0-MD;1-Newton;2-CG;3-damped;4-? ISIF = 2 * calc.stress:0-no;1-tr;2-7 yes; 3-change vol.4-change shape ISYM = 0 symmetry:0-OFF, 1-ON,2-ON+memory conserving for CD POTIM = 0.3 if IBRION=1,2,3:scaling constant forces;IBRION=0,time st. LPLANE = .TRUE. NPAR = 8 IALGO = 48 * algorithm: use only 8 (CG) or 48 (RMM-DIIS) LDIAG = .TRUE. sub-space diagonalisation GGA = 91 * xc-type: PB, PW (86), LM, 91 (Perdew-Wang 91) LWAVE = .F. * write WAVECAR LCHARG = .T. * write CHGCAR and CHG LVTOT = .F. * write the local potential LOCPOT LELF = .F. create ELFCAR file LORBIT = .F. create PROOUT LREAL = Auto ! NEB IMAGES = 8 ! Number of images in the band, excluding endpoints SPRING = -5.0 eV/A2 ! Spring force between images; negative value turns on nudging LCLIMB = .False. ! Flag to turn on the climbing image algorithm |
6楼2010-07-07 10:40:56
7楼2010-07-07 12:07:50
9楼2010-07-12 09:39:04
10楼2010-07-12 09:43:21
