24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dwd0826

铁杆木虫 (著名写手)

应助: 7 (幼儿园)
金币: 6550.9
帖子: 1329
在线: 426.9小时
虫号: 1100436

[交流] 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 ]

回复此楼

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

查看全部散金贴

1楼 2012-07-25 09:38:15

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

itpfeng

木虫 (著名写手)

应助: 10 (幼儿园)
金币: 8159.5
帖子: 1023
在线: 757.5小时
虫号: 1869788

★ ★ ★ ★
小木虫: 金币+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值以及不同的输入参数对结果的影响。经验在于总结，只是这么简单的问是没有意义的。