24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2490  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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

» 猜你喜欢

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

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

查看全部散金贴

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

Cony女孩

银虫 (初入文坛)



小木虫: 金币+0.5, 给个红包,谢谢回帖
您好!这是我做出来的金红石MnO2的态密度图好奇怪。 用vasp算的,LDA+U, U=6.7 J=1.2,算出来再用vasp_band减去费米能,再画图,这样的出来的图 0处是费米能吗?可是在0的位置没有间隙。 打扰您了!

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : DOSCAR
  • 2013-12-12 22:19:42, 420.4 K
  • 附件 2 : band.dat
  • 2013-12-12 22:20:47, 1.22 M
  • 附件 3 : INCAR
  • 2013-12-12 22:21:43, 370 bytes
12楼2013-12-12 22:23:20
已阅   回复此楼   关注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的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见