应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】怎么在VASP中设置NEB方法的参数
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
161/2返回列表12下一页
查看: 3155  |  回复: 15
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

阿勇

铜虫 (小有名气)

[交流] 【求助】怎么在VASP中设置NEB方法的参数

想做一个表面扩散的势垒，想用NEB方法计算，但不知怎么设置INCAR和POSCAR 参数。。。
回复此楼

» 收录本帖的淘帖专辑推荐

vasp计算总结 计算资料收集byFrank

» 猜你喜欢
1楼2009-06-05 16:01:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
spur(金币+2,VIP+0): 3Q！！谢谢帮助！ 6-5 18:58
6.55 Elastic band method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83

VASP说明书上 6.55 83页
赞一下(11人) 回复此楼
2楼2009-06-05 16:37:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cbFeng200881

金虫 (著名写手)
★
csfn(金币+1,VIP+0):3x 6-6 18:53
vasptst 官方网站
赞一下(20人) 回复此楼
3楼2009-06-06 00:20:01
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿勇

铜虫 (小有名气)

谢谢了

谢谢了。。。
回复此楼
4楼2009-06-06 08:22:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fegg7502

荣誉版主 (知名作家)

无

优秀版主
★ ★ ★ ★ ★
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
赞一下(20人) 回复此楼
心诚意正，方可始终；身无苦灭，永登光明。
5楼2009-06-06 16:41:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)
★ ★ ★ ★ ★ ★
小木虫(金币+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这个小节。
赞一下(21人) 回复此楼
6楼2009-06-06 18:32:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zeoliters(金币+2,VIP+0):thank you very much! 6-7 12:53
fegg兄弟的INCAR中IBRION不该取2，应该为1或者3，vasp manual 里面说了。PREC 取Medium就够了。
赞一下(7人) 回复此楼
7楼2009-06-06 19:06:28
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zeoliters

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
wuli8(金币+2,VIP+0):感谢交流，(*^__^*) …… 6-8 19:13
引用回帖:
Originally posted by spur at 2009-6-6 12:06:
fegg兄弟的INCAR中IBRION不该取2，应该为1或者3，vasp manual 里面说了。PREC 取Medium就够了。

对neb,取IBRION2，计算完neb 之后，IBRION取1结构优化，然后计算频率验证
赞一下(21人) 回复此楼
8楼2009-06-07 20:02:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿勇

铜虫 (小有名气)
好久没上小木虫了，谢谢大家的回帖
赞一下 回复此楼
9楼2009-06-12 21:37:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿勇

铜虫 (小有名气)
★ ★
fegg7502(金币+2,VIP+0):感谢交流 6-14 17:13
引用回帖:
Originally posted by zeoliters at 2009-6-7 20:02:

对neb,取IBRION2，计算完neb 之后，IBRION取1结构优化，然后计算频率验证

还想问大家一句，NEB计算时，应该也是在算通过线性插值得到的中间态IMAGE的结构优化吧？ 那算完之后还要做IBRION=1的结构优化吗？
赞一下(10人) 回复此楼
10楼2009-06-12 21:40:29
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿勇 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭