| 查看: 3584 | 回复: 15 | |||||
| 当前主题已经存档。 | |||||
阿勇铜虫 (小有名气)
|
[交流]
【求助】怎么在VASP中设置NEB方法的参数
|
||||
| 想做一个表面扩散的势垒，想用NEB方法计算，但不知怎么设置INCAR和POSCAR 参数。。。 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp计算总结 | 计算资料收集byFrank |
» 猜你喜欢
时间戳变了，能看出什么问题？
已经有10人回复
-
重要消息，中午系统在维护
已经有8人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有105人回复
-
明天放榜？
已经有4人回复
-
投稿咨询
已经有6人回复
-
两姐妹，一个辛辛苦苦30岁工作当大学教师，另一个23岁就工作当开车教练
已经有5人回复
-
时间戳又变了8-15
已经有26人回复
-
今天维护系统维护 祝所有人 高中
已经有9人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有16人回复
-
2027广东省杰青
已经有8人回复
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
2楼2009-06-05 16:37:26
cbFeng200881
金虫 (著名写手)
- 应助: 0 (幼儿园)
- 金币: 262.8
- 散金: 3414
- 红花: 2
- 沙发: 12
- 帖子: 2772
- 在线: 699.4小时
- 虫号: 524662
- 注册: 2008-03-14
- 专业: 凝聚态物性 II ：电子结构
3楼2009-06-06 00:20:01
阿勇
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184.9
- 散金: 61
- 帖子: 86
- 在线: 16.5小时
- 虫号: 33779
- 注册: 2004-01-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2009-06-06 08:22:49
fegg7502
荣誉版主 (知名作家)
无
- 应助: 20 (小学生)
- 贵宾: 13.133
- 金币: 18822.8
- 散金: 2753
- 红花: 31
- 沙发: 13
- 帖子: 7489
- 在线: 591.3小时
- 虫号: 352334
- 注册: 2007-04-22
- 性别: GG
- 专业: 催化化学
- 管辖: 第一性原理
★ ★ ★ ★ ★
csfn(金币+5,VIP+0):3x 6-6 18:54
csfn(金币+5,VIP+0):3x 6-6 18:54
|
SYSTEM = Al2O3 # Startparameter for this Run: # NWRITE = 2 ! how much will be written to the file OUTCAR ('verbosity flag')verbosity write-flag ISTART = 0 ! 0-new 1-cont 2-samecut job,/if WAVECAR exists 1 restart with constant energy cut-off'/2 restart with constant basis set': Continuation job -- read wave functions from the file WAVECAR # ISPIN = 2 ! 1-no, 2-yes INIWAV = 1 ! 0-jellium WF; 1-random numbers (recomm.) ICHARG = 2 ! ISTART = 0 2, else 0 # WEIMIN = 0.001 ! for dynamic calculation IBRION >=0/for static calculation IBRION =-1 dynamic: LPLANE = .TRUE. NPAR = 8 # Electronic Relaxation ENCUT = 350 ! Ecut (eV) PREC = High ! precission: #Low for MD/Medium for optimization/Accurate for TS and frequence EDIFF = 1.E-05 ! stopping-criterion for electronic upd. LREAL = Auto ! .FALSE. projection done in reciprocal space/.TRUE. real space/On or O real space, projection operators re-optimized/Auto real space,fuly auto opti of projection operators no user interference required NELMIN = 5 ! the minimum number of electronic SC steps # NELM = 150 ! max. # of electronic SC steps # NELMDL = 10 ! # of non-consistent steps at the beggining:>0 or <0 # LDIAG = .TURE. ! LDIAG = .TRUE. or .FALSE. (perform sub space rotation) ALGO = Fast # IALGO = 48 ! algorithm: use only 8 (CG) or 48 (RMM-DIIS) GGA = 91 ! xc-type: PB, PW (86), LM, 91 (Perdew-Wang 91) VOSKOWN = 1 ! it is desirable to use this interpolation for single atom or whenever the PW91 functional is applied. ISMEAR = 0 ! part. occup.:-5 Blochl,-4 tet,-3,-2,-1-fermi,0-Gauss,>0MP broadening in eV -4 tet -1-fermi 0-gauss SIGMA = 0.2 ! # IWAVPR = 12 ! if IBRION=0 (MD) 2,12/if IBRION=1,2 (relaxation) 1,11/else (static calculation) 0 # Ionic relaxation EDIFFG = -0.05 ! break condition for the ionic relaxation loop/Default :EDIFF*10/negative ,will stop if all forces are smaller than [EDIFFG] IBRION = 2 ! ionic relax: -1:no move; 0-MD;1-Newton;2-CG;3-damped;4-? NSW = 400 ! steps for ionic update (def:0) ISIF = 2 ! calc.stress:0-no;1-tr;2-7 yes; 3-change vol.4-change shape POTIM = 0.3 # TS flags IMAGES = 8 # NFREE = 2 SPRING = -5 ENERGY1 = -589.41939 ENERGY2 = -590.393218 # Write flags LWAVE = .False. ! write WAVECAR LCHARG = .False. ! write CHGCAR and CHG LVTOT = .False. ! write the local potential LOCPOT LELF = .False. !create ELFCAR file LORBIT = .F. !create PROOUT # Others # APACO = 10.0 !distance for P.C. AMIX = 0.2 BMIX = 0.0001 AMIX_MAG = 0.8 BMIX_MAG = 0.0001 |
5楼2009-06-06 16:41:00
stractor
金虫 (著名写手)
- 1ST强帖: 6
- 应助: 91 (初中生)
- 金币: 1487.6
- 散金: 666
- 红花: 30
- 沙发: 1
- 帖子: 2224
- 在线: 426.8小时
- 虫号: 394087
- 注册: 2007-06-06
- 专业: 半导体材料
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
csfn(金币+5,VIP+0):3x 6-6 18:54
小木虫(金币+0.5):给个红包，谢谢回帖交流
csfn(金币+5,VIP+0):3x 6-6 18:54
|
用VASP做NBE计算，主要有这几个参数： image=N, 初态与末态之间的image个数，SPRING=弹性恢复系数，用默认值-5就可以了，单位是ev/A^2，IBRION=1或3。在当前工作目录下要建N+2个子目录00，01，...，N+1.每个目录包括一个POSCAR，初态,末态的POSCAR分别放在00，N+1目录下面。中间的N个POSCAR需要在初态和末态之间通过线性插值得到。 VASP手册有Elastic band method这个小节。 |
6楼2009-06-06 18:32:18
spur
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 贵宾: 1.21
- 金币: 4994.7
- 散金: 559
- 红花: 5
- 帖子: 866
- 在线: 374.9小时
- 虫号: 215683
- 注册: 2006-03-11
- 性别: GG
- 专业: 第一性原理计算
7楼2009-06-06 19:06:28
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.9小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
8楼2009-06-07 20:02:56
阿勇
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184.9
- 散金: 61
- 帖子: 86
- 在线: 16.5小时
- 虫号: 33779
- 注册: 2004-01-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2009-06-12 21:37:25
阿勇
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 184.9
- 散金: 61
- 帖子: 86
- 在线: 16.5小时
- 虫号: 33779
- 注册: 2004-01-03
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
10楼2009-06-12 21:40:29
