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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
51/1返回列表
查看: 2253  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

流浪猫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
回复此楼

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

吸附研究

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-12-28 20:06:47
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
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，这个应该木有问题
赞一下(1人) 回复此楼
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或许会好点，
赞一下(1人) 回复此楼
嘟啦啦嘟啦啦嘟
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 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 梁成伟老师课题组欢迎你的加入 +6 一鸭鸭哟 2026-03-14 7/350 2026-03-15 22:12 by Winj1e
[文学芳草园] 伙伴们，祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 283求调剂 +6 小楼。 2026-03-12 9/450 2026-03-14 22:15 by cfxczx2015
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 复试调剂 +3 呼呼？~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 求调剂，药学 +3 归零lbm 2026-03-09 5/250 2026-03-14 02:21 by JourneyLucky
[考研] 265求调剂 +9 小木虫085600 2026-03-09 12/600 2026-03-14 01:11 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
关闭
请填处理意见
关闭