小木虫

当前位置：首页 > 第一原理 >【讨论】vasp 单点能计算收敛

【讨论】vasp 单点能计算收敛

作者 zhangfan2192: 来源: 小木虫 700 14 举报帖子

用vasp计算单点能时，遇到不收敛的情况时，大家都是这么解决的？欢迎大家讨论。

我看到网上说可以调整INCAR 里面的tag 大家都是这么调的啊。

[ Last edited by zhangfan2192 on 2009-7-11 at 10:50 ] 返回小木虫查看更多

今日热帖

精华评论

bluce

要看你算得什么体系了。
六198

问题不清
zhangfan2192

引用回帖:
Originally posted by 六198 at 2009-7-11 08:33:

问题不清

体系也不是很复杂，就是Cr的表面加Hg原子。算单点能的时候就是不收敛，不知道调INCAR 里面的哪些参数可以增加收敛。谢谢
hakuna

首先，你得告诉大家你的INCAR，让别人猜不是解决问题的办法。
其次，虽然是单点不收敛，还是建议检查一下初始结构是否合理。
zhangfan2192

INCAR 如下
SYSTEM = Cu-222
  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  = 338.00 #Ecut (eV)
PREC = Medium  *  #precission: #Low/Medium/High
NELM = 100    max. # of electronic SC steps
NELMDL  = -12    # of non-consistent steps at the beggining:>0 or <0
EDIFF  = 1.0e-4 stopping-criterion for electronic upd.
VOSKOWN = 1
ISMEAR = -5
SIGMA = 0.1
AMIX = 0.2
BMIX = 0.0001
NBANDS = 249

  Ionic relaxation
EDIFFG =  1.0e-3 stopping-conditions for ionic relax.loop (EDIFF*10)
NSW = 1    *  # 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，
nkleof

引用回帖:
Originally posted by zhangfan2192 at 2009-7-14 20:43:
INCAR 如下
SYSTEM = Cu-222
  Startparameter for this Run:
NWRITE =    2    verbosity write-flag
LPETIM=    F    write-flag & timer #
ISTART =    0    0-new 1-cont 2-samec ...

你的INCAR设置很复杂啊，所有参数确认都经过调试是合适的？NSW为什么设置为1？这样的话既要做离子弛豫，又限制离子运动步数这么小，结构优化很奇怪啊。确认在OUTCAR里面Iteration还在继续进行吗？另外建议如果不是特殊必要的话，很多tag可以删掉，用默认的。
zhangfan2192

引用回帖:
Originally posted by nkleof at 2009-7-14 22:13:

你的INCAR设置很复杂啊，所有参数确认都经过调试是合适的？NSW为什么设置为1？这样的话既要做离子弛豫，又限制离子运动步数这么小，结构优化很奇怪啊。确认在OUTCAR里面Iteration还在继续进行吗？另外建议如果不 ...

我要算的是单点能
所以NSW 设为1
至于其他参数是因为默认的不管用我加上去的但是还是不收敛
是在不行的话我只有去改初始构型了。