24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 2112  |  回复: 4

流浪猫8386

银虫 (正式写手)

[求助] 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题

关于vasp5.2版本以及GGA+U的方法还在摸索中,所以有很多问题想向各位虫请教
      1. 对于CeO2(111)表面吸附的构型优化,在INCAR中设置了自旋ISPIN=2,那是否还需要设置MAGMOM,如果需要设置,应该怎么设置。与NUPDOWN一起使用又该考虑哪些问题。
    2. U值的选取及如何处理+U后不收敛的情况。
    3.对于优化结果,vasp4.5版本算出来正常,固定的原子受力为0,但是5.2版本计算的放开的原子达到了受力收敛标准,反而固定的原子受力很大,但是计算结果说已经达到所需精度,4.6版本的我也测试了,出现和5.2版本一样的问题。

附上我的INCAR
   LWAVE =  .FALSE.  do not write WAVECAR file
   LCHARG = .FALSE.  do not write CHGCAR file
   LVTOT =  .TRUE.  do not write LOCPOT file

# settings for parallel version (ignored in serial run)
  LPLANE=.TRUE.
  NPAR= 2
  NSIM= 1

Electronic Relaxation 1
#   ISTART = 3
   PREC = Medium       Medium=default, Low, High; affects ENMAX, mesh, pspot
   EDIFF  =  1E-05    energy stopping-criterion for electr. iterations
   EDIFFG =   -.01    force stopping-criterion for geometry steps
   IDIPOL = 3
#  MAXMIX = 80
#  AMIX = 0.2
#  BMIX = 0.0001
#  AMIX_MAG = 0.8
#  BMIX_MAG = 0.0001
#  NBANDS = 194
#   NELMDL = -120
Ionic Relaxation
    NSW    = 1000     max number of geometry steps
    IBRION =   2     ionic relax: 0-MD 1-quasi-New 2-CG
    ISYM   =   0     1=use symmetry, 0 = no symmetry
    ISPIN = 2
    ISIF = 0
    NELMIN  = 3
    POTIM = 0.5
#    MAGMOM=0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3
#    NUPDOWN = 5
DOS related values:
   ISMEAR =      0    -4-tet -1-fermi 1=Methfessel/Paxton 1.order
   SIGMA  =     .05    broadening in eV


Electronic Relaxation
   IALGO  =     48    algorithm (8=CG for small, 48=RMM for big systems)
   LREAL= Auto    real-space projection

#   ICORELEVEL = 1

     LDAU = .TRUE.
     LMAXMIX = 6
     LDAUTYPE = 2
     LDAUL = 3 -1 -1
     LDAUU = 4.5 0 0
     LDAUJ = 0 0 0
     LDAUPRINT = 2
回复此楼

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

吸附研究

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
fzx2008(金币+1): 谢谢指教 2011-12-28 23:09:59
流浪猫8386(金币+1): 2012-01-06 21:16:22
EDIFFG =   -.01    force stopping-criterion for geometry steps
这个值有问题 明确定义吧
EDIFFG =  0.0001或者0.001
嘟啦啦嘟啦啦嘟
2楼2011-12-28 22:27:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主


youzhizhe(金币+1): 谢谢交流。 2011-12-29 16:09:30
引用回帖:
2楼: Originally posted by 贺仪 at 2011-12-28 22:27:58:
EDIFFG =   -.01    force stopping-criterion for geometry steps
这个值有问题 明确定义吧
EDIFFG =  0.0001或者0.001

vasp是用fortran语言写的,所以写-.01应该就是-0.01,这个应该木有问题
3楼2011-12-28 23:09:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贺仪

铁杆木虫 (著名写手)

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-12-29 16:09:18
引用回帖:
3楼: Originally posted by fzx2008 at 2011-12-28 23:09:48:
vasp是用fortran语言写的,所以写-.01应该就是-0.01,这个应该木有问题

恩,把负号连在一起我都没明白过来
他说的不是原子受力的问题嘛, 不过这个值太大了 -0.001和-0.0001或许会好点,
嘟啦啦嘟啦啦嘟
4楼2011-12-29 08:37:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhuzisummer

银虫 (小有名气)

我不懂VASP,也很想知道。。。:
5楼2011-12-30 19:09:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 流浪猫8386 的主题更新
信息提示
请填处理意见