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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » GGA+U算出DOS图分析
51/1返回列表
查看: 2597  |  回复: 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭