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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 紧急求助vasp5.2 计算氧化铈表面吸附遇到的若干问题
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 2282  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 312求调剂 +10 gtw1 2026-04-06 10/500 2026-04-06 14:15 by 蒋皓禹
[考研] 调剂 +9 不逢春 2026-04-05 10/500 2026-04-06 11:58 by 考研学校招点人
[考研] 第一志愿东南大学物理313，有科研竞赛获奖经历，希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 专硕0854初试考材科基，求调剂 +6 3220548044 2026-04-06 9/450 2026-04-06 10:26 by barlinike
[考研] 化学调剂 +17 艾志恒 2026-04-03 18/900 2026-04-06 07:10 by jj987
[考研] 308求调剂 +4 maverick^_^ 2026-04-03 4/200 2026-04-05 19:08 by 蓝云思雨
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 270求调剂 +9 小杰pp 2026-03-31 11/550 2026-04-05 11:02 by 风雨无晴
[考研] 271分求调剂学校 +12 zph158488！ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 303求调剂 +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by ls刘帅
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 英一数一408，总分284，二战真诚求调剂 +13 12.27 2026-03-30 15/750 2026-04-03 14:41 by 氮气气气
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 08工科，295，接受跨专业调剂 +8 lmnlzy 2026-03-30 8/400 2026-04-03 13:08 by nalakaiqi
[考研] 318求调剂，计算材料方向 +10 吸喵有害笙命 2026-04-01 11/550 2026-04-02 16:29 by oooqiao
[考研] 生物学327，求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
信息提示
关闭
请填处理意见
关闭