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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求助： VASP Relaxation & Band structure calculation of alpha-Al2O3
51/1返回列表
查看: 1162  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

klwang

木虫 (著名写手)

[求助] 求助： VASP Relaxation & Band structure calculation of alpha-Al2O3 已有2人参与

正在学习建模，用VASP进行结构relaxation和band structure的计算。我将稀土元素（Nd，Er)分布掺入alpha-Al2O3的空隙中，结果收敛得很不好（下面的模型是我做的最好的结果了），请大侠给予建议，如何改进。
另外，试图计算电子结构（包括DOS和能带结构），不知道如何做，请给予详细的解释。不胜感激。
我的INCAR如下
System =  Al12O18Nd

ISMEAR = 0        ! Use 0 for  KPOINTS less than 4 otherwise -5.

PREC  = Accurate       ! low, medium, normal are other options. Use suitable one.
ENCUT = 600 eV    ! Decide considering the crystal size and accuracy you want.
EDIFF = 1.0E-5    ! Enegy difference covergence limit for electronic optimization.
EDIFFG = -1.0E-3   ! Enegy difference covergence limit for ionic optimization.

IBRION = 2        ! 0 for MD, 1 best, 2 for diff relaxation problems. Othere numbers for diffrent cases.
NSW    = 1000      ! Total number of ionic steps.
ISIF   = 3        ! 2 and 4 ionic, 7 volume and 3 both.

LREAL  = Auto     ! Projection on real space. use FALSE (default) for reciprocal space.
NPAR   = 4       ! Best sqrt of NCPUs used. should be >= NCPUs/32.
ALGO   = Fast     ! default is Normal.

++++++++++++++++++++++++++++++++++++++++++++++++
用grep E0 OSZICAR得到如下结果
  1 F= -.22657211E+03 E0= -.22637830E+03  d E =-.226572E+03
   2 F= -.22643608E+03 E0= -.22624023E+03  d E =0.136036E+00
   3 F= -.22643396E+03 E0= -.22623968E+03  d E =0.138149E+00
   4 F= -.22649364E+03 E0= -.22629965E+03  d E =0.784727E-01
   5 F= -.22651988E+03 E0= -.22632516E+03  d E =0.522362E-01
   6 F= -.22651796E+03 E0= -.22632327E+03  d E =0.541501E-01
   7 F= -.22656292E+03 E0= -.22637028E+03  d E =-.449578E-01
   8 F= -.22657929E+03 E0= -.22638289E+03  d E =-.613257E-01
   9 F= -.22658334E+03 E0= -.22638786E+03  d E =-.405303E-02
  10 F= -.22660753E+03 E0= -.22641141E+03  d E =-.282378E-01
  11 F= -.22651116E+03 E0= -.22631442E+03  d E =0.963654E-01
  12 F= -.22655311E+03 E0= -.22635636E+03  d E =0.544181E-01
  13 F= -.22652968E+03 E0= -.22633326E+03  d E =0.778444E-01
  14 F= -.22654235E+03 E0= -.22634586E+03  d E =-.126660E-01
  15 F= -.22655874E+03 E0= -.22636264E+03  d E =-.290588E-01
  16 F= -.22646476E+03 E0= -.22625910E+03  d E =0.939837E-01
  17 F= -.22650256E+03 E0= -.22629580E+03  d E =0.561816E-01
  18 F= -.22648921E+03 E0= -.22628608E+03  d E =0.695343E-01
  19 F= -.22645195E+03 E0= -.22625323E+03  d E =0.106795E+00
  20 F= -.22644833E+03 E0= -.22625366E+03  d E =0.110412E+00
  21 F= -.22655066E+03 E0= -.22635187E+03  d E =0.807837E-02
  22 F= -.22655919E+03 E0= -.22636233E+03  d E =-.448879E-03
  23 F= -.22655441E+03 E0= -.22635730E+03  d E =0.433572E-02
  24 F= -.22655623E+03 E0= -.22635913E+03  d E =0.250788E-02
  25 F= -.22655690E+03 E0= -.22636003E+03  d E =0.184416E-02
  26 F= -.22655513E+03 E0= -.22635853E+03  d E =0.361137E-02
  27 F= -.22655583E+03 E0= -.22635889E+03  d E =0.291569E-02
  28 F= -.22655755E+03 E0= -.22636045E+03  d E =0.119494E-02
  29 F= -.22655605E+03 E0= -.22635897E+03  d E =0.269185E-02
  30 F= -.22655354E+03 E0= -.22635651E+03  d E =0.519719E-02
  31 F= -.22655527E+03 E0= -.22635896E+03  d E =0.346698E-02
  32 F= -.22655615E+03 E0= -.22635914E+03  d E =0.259426E-02
  33 F= -.22655440E+03 E0= -.22635778E+03  d E =0.433895E-02
  34 F= -.22656165E+03 E0= -.22636538E+03  d E =-.291174E-02


如果做relaxation，INCAR 应该如何改进？其他建议不胜感激。
如果做电子结构计算，INCAR应该如何改动？其他建议不胜感激。
如果有相关的电子书或链接，请不吝分享。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-02-02 08:35:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skydark

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2015-02-15 07:49:28
klwang: 金币+20, ★★★很有帮助 2015-02-16 00:48:27
你一開始精確度別調太高  分次relax看看
也就是說第一跑的時候 把所有跟精確度有關的設定都調低一點
讓第一次relax能收斂後  再把第一次跑出?淼腃ONTCAR檔轉成POSCAR檔
再把精確度調高 跑第2次relax
另外在relax那一步
我ISMEAR都調成 -5  且 Kpoint 別設太少
你Kpoint檔也貼上?斫o大家看看吧

PS: 你這OSZICAR看起?碛虚_始收斂 怎不讓他跑完?
赞一下 回复此楼
4楼2015-02-14 22:38:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
klwang: 金币+5, ★★★很有帮助, 能否具体说明改NELM还是NSW或其他参数？谢谢！ 2015-02-02 14:48:34
难收敛很可能是结构不合理引起的，建议把步长适当调小一点，试试
赞一下 回复此楼
2楼2015-02-02 10:23:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
klwang: 金币+15, ★★★很有帮助, 谢谢！请问DOS 计算有任何建议吗？新手，请多指教。如果有操作步骤详细介绍将不胜感激。我搜到一个版本介绍如何算DOS，但是很多步骤不熟悉，不是很清楚如何做。比如先做自洽，然后算BAND STRUCTUE. 2015-02-02 22:56:48
liliangfang: 金币+1, 谢谢交流 2015-02-15 07:49:23
POTIM。主要还是调整模型，使结构更合理一些，并不是任何一种结构它都能帮你优化出最优结构的
赞一下 回复此楼
3楼2015-02-02 15:27:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

klwang

木虫 (著名写手)
引用回帖:
4楼: Originally posted by skydark at 2015-02-14 02:38:52
你一開始精確度別調太高  分次relax看看
也就是說第一跑的時候 把所有跟精確度有關的設定都調低一點
讓第一次relax能收斂後  再把第一次跑出?淼腃ONTCAR檔轉成POSCAR檔
再把精確度調高 跑第2次relax
另外在rela ...

您好！谢谢指点。其实我已经把CONTCAR拷贝到POSCAR了，精度都降低了，已经算了有近2000步了，还是不收敛，在10e-2和10e-3之间徘徊。
AMIX=0.5 BMIX=1. 我调低了AMIX 和BMIX，收敛变坏，所以在尝试调整这些参数，如果还是不行，可能要重新考虑其他参数了，比如BMIX and AMIN.
这些参数的调试我不是很有经验，以前的体系很简单。如果您在这些参数调试上有什么好的建议，请多多指教。谢谢！
赞一下(1人) 回复此楼
5楼2015-02-16 01:24:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭