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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
51/1返回列表
查看: 2581  |  回复: 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的回帖

贺仪

铁杆木虫 (著名写手)
★ ★
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的回帖
查看全部 5 个回答

贺仪

铁杆木虫 (著名写手)

【答案】应助回帖

★
感谢参与，应助指数 +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的回帖

zhuzisummer

银虫 (小有名气)
我不懂VASP，也很想知道。。。:
赞一下 回复此楼
5楼2011-12-30 19:09:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 16:33 by Z9YWQ5EAO3qp
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 cqwQDCxMcL3I 2026-09-29 6/300 2026-09-30 16:14 by Z9YWQ5EAO3qp
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
信息提示
关闭
请填处理意见
关闭