| 查看: 3244 | 回复: 13 | ||
[交流]
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 ] |
回复此楼
|
看了下你的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. 以供参考，祝好！ |
» 本帖已获得的红花（最新10朵）
Cony女孩13楼2013-12-13 20:49:51
2楼2012-07-25 10:11:20
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
franch: 金币+2, 谢谢回帖交流。 2012-07-25 18:03:32
dwd0826: 金币+1, thx for your guidance！ 2012-07-26 08:35:19
小木虫: 金币+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
4楼2012-07-25 13:44:04
5楼2012-07-25 14:16:07
6楼2012-07-26 08:34:27
7楼2012-07-26 08:36:24
8楼2012-07-26 08:37:07
9楼2012-07-26 08:52:49
10楼2012-07-26 10:29:06
11楼2012-07-28 13:45:13
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
| 您好！这是我做出来的金红石MnO2的态密度图好奇怪。 用vasp算的，LDA+U, U=6.7 J=1.2,算出来再用vasp_band减去费米能，再画图，这样的出来的图 0处是费米能吗？可是在0的位置没有间隙。 打扰您了！ |
12楼2013-12-12 22:23:20
14楼2013-12-16 15:59:26