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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
51/1返回列表
查看: 2252  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 304求调剂 +5 小熊joy 2026-03-14 5/250 2026-03-14 21:07 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 0856材料与化工309分求调剂 +6 ZyZy…… 2026-03-10 6/300 2026-03-14 00:38 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 5/250 2026-03-13 22:06 by 星空星月
[考研] 285学硕求调剂 +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 282分材料专业求调剂院校 +18 枫桥ZL 2026-03-09 25/1250 2026-03-13 10:47 by 白夜悠长
[考研] 296求调剂 +3 大口吃饭 身体健 2026-03-13 3/150 2026-03-13 10:31 by 学员8dgXkO
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3 我是一个兵， 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考博] 2026年博士申请 +3 QwQwQW10 2026-03-11 3/150 2026-03-12 17:58 by gxch43
[考研] 研究生招生 +3 徐海涛11 2026-03-10 7/350 2026-03-12 14:26 by 徐海涛11
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
关闭
请填处理意见
关闭