| 查看: 1391 | 回复: 6 | ||
[求助]
求助： 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应该如何改动？其他建议不胜感激。 如果有相关的电子书或链接，请不吝分享。 |
回复此楼» 猜你喜欢
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有23人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Too Small Band Gap of Silicon (~0.65 eV) Calculated by VASP
已经有3人回复
- Band structure of Si within GW and VASP
已经有63人回复
- vasp5.2如何计算DOS和band structure？
已经有10人回复
- VASP计算输出文件里面哪里可以看到band gap啊？
已经有4人回复
- [求助成功]VASP 计算 graphene 的错误
已经有15人回复
- 【求助】VASP计算的KPOINTS 设置，求助
已经有10人回复
- 【求助】band structures using hybrid functionals 的Kpoints问题
已经有18人回复
- 【求助】vasp计算的结果显示segmentation fault是什么错误
已经有5人回复
- 【求助】请教vasp如何使本征值按每条能带输出
已经有18人回复
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6521.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.8小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
2楼2015-02-02 10:23:34
未岸亡灵
铁杆木虫 (著名写手)
暂时离开小木虫
- 1ST强帖: 2
- 应助: 622 (博士)
- 金币: 6521.3
- 散金: 1221
- 红花: 111
- 帖子: 1345
- 在线: 517.8小时
- 虫号: 1730773
- 注册: 2012-04-01
- 性别: GG
- 专业: 理论和计算化学
3楼2015-02-02 15:27:14
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
liliangfang: 金币+1, 谢谢交流 2015-02-15 07:49:28
klwang: 金币+20, ★★★很有帮助 2015-02-16 00:48:27
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
5楼2015-02-16 01:24:53
skydark
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 44 (小学生)
- 金币: 3652.2
- 散金: 30
- 红花: 6
- 帖子: 302
- 在线: 249.8小时
- 虫号: 1927082
- 注册: 2012-08-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-02-16 12:24:49
7楼2015-02-17 14:12:14
28