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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » GGA+U算出DOS图分析
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 2704  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dwd0826

铁杆木虫 (著名写手)

[交流] GGA+U算出DOS图分析

最近算了下两种不同的MnO2的TDOS，如下:


比较不解的是，两种结构算出来都为导体？而MnO2按文献为半导体，如下图：

不知是我的分析有误还是计算问题？求解答。以下为DOS的INCAR：
--------------------------------------
Start parameter for this run:
-----------------------------------------
   ISTART = 1
   ICHARG = 11

------------------------------------------
Electronic Relaxation:
------------------------------------------
   LREAL  = F
   PREC   = A            precision:Medium/High/Low
   ENCUT = 520
   EDIFF  = 1e-4              converge criterion: default=1e-4
   EDIFFG = -0.02
   MAXMIX = 50
#  ALGO  = N
   NELM = 300
   NELMIN = 5                 the minimum number of electronic SC steps
   IBRION = -1
   NSW    = 0
   ISIF   = 0
   ISYM   = 2                 symmetry (2-PAW on, 1-US-PP's on, 0-off)
   GGA    = PE                Perdew-Wang 91
   LORBIT = 11
   ISPIN = 2
   MAGMOM = n*4 m*0
--------------------------------------------
Output Control:
--------------------------------------------
   LWAVE  = .FALSE.
   LCHARG = .FALSE.
!   LVTOT  = .TRUE.
   NWRITE = 2
--------------------------------------------
DOS Related:
--------------------------------------------
   ISMEAR = -5
   LDAU = .TRUE.   ！Switches on the L(S)DA+U.
   LDAUTYPE = 2    ！1|2|4 Type of L(S)DA+U (Default: LDAUTYPE = 2)
   LDAUL = 2 -1
   LDAUU =3.9 0
   LDAUJ =0 0
   LDAUPRINT = 2
!   EMIN = -10
!   EMAX = 6
!   NEDOS = 3001

[ Last edited by dwd0826 on 2012-7-25 at 09:42 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2012-07-25 09:38:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dwd0826

铁杆木虫 (著名写手)
引用回帖:
12楼: Originally posted by Cony女孩 at 2013-12-12 22:23:20
您好！这是我做出来的金红石MnO2的态密度图好奇怪。 用vasp算的，LDA+U, U=6.7 J=1.2,算出来再用vasp_band减去费米能，再画图，这样的出来的图 0处是费米能吗？可是在0的位置没有间隙。 打扰您了！

看了下你的INCAR，主要问题出在以下四个参数上
LDAUTYDE=4
LDAUL=-1 -1 -1 -1 2 2
LDAUU=0 0 0 0 6.7 6.7
LDAUJ=0 0 0 0 1.2 1.2
首先，J的值的意义与LDAUTYPE有关的。当LDAUTYPE为2时，U-J的差值才有有意义，即有效的U参数。因此建议LDAUTYDE=2
第二就是LDAUL，这里是针对POTCAR当中的元素种类进行设置，如你INCAR当中的设置，POTCAR顺序为O、Mn，
LDAUL=-1 2
LDAUU=0 3.9
LDAUJ=0 0
其中U与J值是根据诸多已有文献中取值，
Phys. Chem. Chem. Phys., 2013, 15, 9075.
Phys. Rev. Lett., 96, 116405.
Phys. Rev. B, 75, 195128.
Chem. Mater., 24, 3943.
以供参考，祝好！
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

Cony女孩
13楼2013-12-13 20:49:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

dwd0826

铁杆木虫 (著名写手)
自己回复下，等解释
赞一下 回复此楼
2楼2012-07-25 10:11:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

itpfeng

木虫 (著名写手)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流。 2012-07-25 18:03:32
dwd0826: 金币+1, thx for your guidance！ 2012-07-26 08:35:19
U 值你取的多大呢？很明显你的gap没有文献中的大，有可能是这个原因吧。再有你算过LDA+U么，以及不加U么？GGA虽然号称精确，但有多人是不用的，因为它不是那么道地的first principal。像这种计算的问题要多算，然后明白不同的近似、不同的势、不同的U值以及不同的输入参数对结果的影响。经验在于总结，只是这么简单的问是没有意义的。
赞一下 回复此楼
3楼2012-07-25 12:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

quby7503

银虫 (著名写手)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流。 2012-07-25 18:03:43
dwd0826: 金币+1, thx for your guidance！ 2012-07-26 08:36:41
对于强关联体系，普通的DFT计算是得不到合理的结果的，GGA+U中用经验的U值去调体系的能隙，U值是可以计算的，但目前版本的VASP还做不到，所以可以将其当做一个经验参数，根据计算结果去调这个U值，直到得到合理的结果。
赞一下 回复此楼
4楼2012-07-25 13:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿河北工业大学0817化工278分求调剂 +10 jhybd 2026-03-23 15/750 2026-03-29 21:49 by nanaliuyun
[考博] 26申博自荐 +6 whh869393 2026-03-24 6/300 2026-03-29 21:05 by nxgogo
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +4 小小机灵虫 2026-03-29 10/500 2026-03-29 18:20 by 无际的草原
[考研] 085600 材料与化工 329分求调剂 +14 Mr. Z 2026-03-25 14/700 2026-03-29 17:27 by 学员8dgXkO
[考研] 一志愿郑州大学，080500学硕，总分317分求调剂 +8 举个栗子oi 2026-03-24 9/450 2026-03-29 13:08 by peike
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +7 WW.' 2026-03-23 9/450 2026-03-29 12:42 by fmesaito
[考研] 348求调剂 +5 小懒虫不懒了 2026-03-28 5/250 2026-03-29 10:34 by 唐沐儿
[考研] 352分-085602-一志愿985 +5 海纳百川Ly 2026-03-29 5/250 2026-03-29 09:57 by Sjndkwm
[考研] 材料学硕333求调剂 +11 北道巷 2026-03-24 11/550 2026-03-29 08:31 by fmesaito
[考研] 压国家一区线，求导师收留，有恩必谢！ +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 086502化学工程342求调剂 +6 阿姨复古不过 2026-03-27 6/300 2026-03-28 07:06 by wangy0907
[考研] 333求调剂 +3 question挽风 2026-03-23 3/150 2026-03-27 11:29 by 不吃魚的貓
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-25 3/150 2026-03-25 19:03 by Ainin_
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 一志愿吉林大学材料与化工303分求调剂 +4 为学666 2026-03-24 4/200 2026-03-25 11:27 by BruceLiu320
[考研] 材料考研调剂生 +3 黄粱一梦千年 2026-03-24 3/150 2026-03-24 17:00 by barlinike
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
信息提示
关闭
请填处理意见
关闭