24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4648  |  回复: 14
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zhangfan2192

金虫 (著名写手)

[交流] 【讨论】vasp 单点能计算收敛

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


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

[ Last edited by zhangfan2192 on 2009-7-11 at 10:50 ]
回复此楼
相逢是缘!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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
相逢是缘!
6楼2009-07-14 20:43:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

bluce

木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuchenwf(金币+1,VIP+0):谢谢回帖交流 7-16 22:51
要看你算得什么体系了。
天性驽钝,唯勤学以补之。出生寒门,校非名校,师非名门,一切靠自己打拼。
2楼2009-07-11 05:17:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

六198

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流

问题不清
穷理求真真贯纵横皆入咏致知格物物含妙理总堪寻
3楼2009-07-11 08:33:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)

★ ★
fegg7502(金币+2,VIP+0):thank you very much! 7-13 03:47
引用回帖:
Originally posted by 六198 at 2009-7-11 08:33:

问题不清

体系也不是很复杂,就是Cr的表面加Hg原子。算单点能的时候就是不收敛,不知道调INCAR 里面的哪些参数可以增加收敛。谢谢
相逢是缘!
4楼2009-07-11 10:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见