版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

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

最近算了下兩種不同的MnO2的TDOS，如下:

比較不解的是，兩種結(jié)構(gòu)算出來都為導(dǎo)體？而MnO2按文獻為半導(dǎo)體，如下圖：

不知是我的分析有誤還是計算問題？求解答。以下為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 ]

回復(fù)此樓

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

引用回帖:

12樓: Originally posted by Cony女孩 at 2013-12-12 22:23:20
您好！這是我做出來的金紅石MnO2的態(tài)密度圖好奇怪。用vasp算的，LDA+U, U=6.7 J=1.2,算出來再用vasp_band減去費米能，再畫圖，這樣的出來的圖 0處是費米能嗎？可是在0的位置沒有間隙。打擾您了！

看了下你的INCAR，主要問題出在以下四個參數(shù)上
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有關(guān)的。當LDAUTYPE為2時，U-J的差值才有有意義，即有效的U參數(shù)。因此建議LDAUTYDE=2
第二就是LDAUL，這里是針對POTCAR當中的元素種類進行設(shè)置，如你INCAR當中的設(shè)置，POTCAR順序為O、Mn，
LDAUL=-1 2
LDAUU=0 3.9
LDAUJ=0 0
其中U與J值是根據(jù)諸多已有文獻中取值，
Phys. Chem. Chem. Phys., 2013, 15, 9075.
Phys. Rev. Lett., 96, 116405.
Phys. Rev. B, 75, 195128.
Chem. Mater., 24, 3943.
以供參考，祝好！

贊一下(2人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

Cony女孩

13樓2013-12-13 20:49:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

自己回復(fù)下，等解釋

贊一下

回復(fù)此樓

2樓2012-07-25 10:11:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

itpfeng

木蟲 (著名寫手)

應(yīng)助: 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值以及不同的輸入?yún)?shù)對結(jié)果的影響。經(jīng)驗在于總結(jié)，只是這么簡單的問是沒有意義的。

贊一下

回復(fù)此樓

3樓2012-07-25 12:58:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
帖子: 1488
在線: 263.7小時
蟲號: 638688

★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
franch: 金幣+2, 謝謝回帖交流。 2012-07-25 18:03:43
dwd0826: 金幣+1, thx for your guidance！ 2012-07-26 08:36:41

對于強關(guān)聯(lián)體系，普通的DFT計算是得不到合理的結(jié)果的，GGA+U中用經(jīng)驗的U值去調(diào)體系的能隙，U值是可以計算的，但目前版本的VASP還做不到，所以可以將其當做一個經(jīng)驗參數(shù)，根據(jù)計算結(jié)果去調(diào)這個U值，直到得到合理的結(jié)果。

贊一下

回復(fù)此樓

4樓2012-07-25 13:44:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

鳳凰-涅磐

金蟲 (著名寫手)

應(yīng)助: 50 (小學生)
金幣: 779.3
帖子: 2963
在線: 462小時
蟲號: 1903667

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

這個圖有毛病吧，峰應(yīng)該向下呀？郁悶

贊一下

回復(fù)此樓

5樓2012-07-25 14:16:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

引用回帖:

3樓: Originally posted by itpfeng at 2012-07-25 12:58:55
U 值你取的多大呢？很明顯你的gap沒有文獻中的大，有可能是這個原因吧。再有你算過LDA+U么，以及不加U么？GGA雖然號稱精確，但有多人是不用的，因為它不是那么道地的first principal。像這種計算的問題要多算，然后 ...

這個U值與pseudopotential是根據(jù)文獻U值確定，只是對比起來計算結(jié)果不是很能理解

贊一下

回復(fù)此樓

6樓2012-07-26 08:34:27

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

引用回帖:

4樓: Originally posted by quby7503 at 2012-07-25 13:44:04
對于強關(guān)聯(lián)體系，普通的DFT計算是得不到合理的結(jié)果的，GGA+U中用經(jīng)驗的U值去調(diào)體系的能隙，U值是可以計算的，但目前版本的VASP還做不到，所以可以將其當做一個經(jīng)驗參數(shù)，根據(jù)計算結(jié)果去調(diào)這個U值，直到得到合理的結(jié) ...

這些都是根據(jù)文獻取值再做對比的，所以不存在U值取值范圍問題，謝謝你的交流

贊一下

回復(fù)此樓

7樓2012-07-26 08:36:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dwd0826

鐵桿木蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 6550.9
帖子: 1329
在線: 426.9小時
蟲號: 1100436

引用回帖:

5樓: Originally posted by 鳳凰-涅磐 at 2012-07-25 14:16:07
這個圖有毛病吧，峰應(yīng)該向下呀？郁悶

能說得更明白些嗎？

贊一下

回復(fù)此樓

8樓2012-07-26 08:37:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

dujunhui

金蟲 (正式寫手)

應(yīng)助: 16 (小學生)
金幣: 1769.3
帖子: 403
在線: 136.3小時
蟲號: 739582

★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
franch: 金幣+2, 謝謝回帖交流。 2012-07-26 19:40:07
dwd0826: 金幣+1, thx for your guidance！ 2012-07-30 10:08:45
dwd0826: 金幣+2 2012-12-29 10:23:19

引用回帖:

7樓: Originally posted by dwd0826 at 2012-07-26 08:36:24
這些都是根據(jù)文獻取值再做對比的，所以不存在U值取值范圍問題，謝謝你的交流

...

文獻上的數(shù)據(jù)也不一定可靠啊，而且計算結(jié)果應(yīng)該是和很多參數(shù)設(shè)置都有關(guān)系的，就像itpfeng說的。

贊一下

回復(fù)此樓

9樓2012-07-26 08:52:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

quby7503

銀蟲 (著名寫手)

應(yīng)助: 7 (幼兒園)
金幣: 485.7
帖子: 1488
在線: 263.7小時
蟲號: 638688

★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖
franch: 金幣+2, 謝謝回帖交流。 2012-07-26 19:40:14
dwd0826: 金幣+1, thx for your guidance！ 2012-07-30 10:09:07
dwd0826: 金幣+3 2012-12-29 10:23:04

引用回帖:

7樓: Originally posted by dwd0826 at 2012-07-26 08:36:24
這些都是根據(jù)文獻取值再做對比的，所以不存在U值取值范圍問題，謝謝你的交流

...

文獻也是根據(jù)計算結(jié)果調(diào)的，文獻可以調(diào)，樓主同樣可以調(diào)，還有不同的軟件，U值也不一定一樣，所以沒有必要非和文獻一致，和實驗結(jié)果一樣才是王道

贊一下

回復(fù)此樓

10樓2012-07-26 10:29:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

umklapp

銀蟲 (小有名氣)

應(yīng)助: 9 (幼兒園)
金幣: 463.6
帖子: 123
在線: 59.1小時
蟲號: 1181441

★
dwd0826: 金幣+1, thx for your guidance！ 2012-07-30 10:09:24

能否把你參考文獻的圖貼上來；
能否把你的 density matrix 貼出來；

贊一下

回復(fù)此樓

11樓2012-07-28 13:45:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Cony女孩

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 269.8
帖子: 46
在線: 105.5小時
蟲號: 2366315

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

您好！這是我做出來的金紅石MnO2的態(tài)密度圖好奇怪。用vasp算的，LDA+U, U=6.7 J=1.2,算出來再用vasp_band減去費米能，再畫圖，這樣的出來的圖 0處是費米能嗎？可是在0的位置沒有間隙。打擾您了！

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負責。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

Cony女孩

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 269.8
帖子: 46
在線: 105.5小時
蟲號: 2366315

送紅花一朵

引用回帖:

13樓: Originally posted by dwd0826 at 2013-12-13 20:49:51
看了下你的INCAR，主要問題出在以下四個參數(shù)上
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有關(guān)的。當LDAUTYPE為2時，U-J的差值才有有意義 ...

非常感謝樓主�。。�

回復(fù)此樓

14樓2013-12-16 15:59:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 dwd0826 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 276求調(diào)劑 +3	YNRYG 2026-03-21	4/200	2026-03-23 08:31 by 醉在風里
[考研] 287求調(diào)劑 +8	晨昏線與星海 2026-03-19	9/450	2026-03-22 17:01 by i_cooler
[考研] 生物學071000 329分求調(diào)劑 +5	我愛生物生物愛?/a> 2026-03-17	5/250	2026-03-22 16:42 by tcx007
[考研] 269專碩求調(diào)劑 +6	金恩貝 2026-03-21	6/300	2026-03-22 14:31 by ColorlessPI
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學員8dgXkO
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學員8dgXkO
[考研] 316求調(diào)劑 +6	梁茜雯 2026-03-19	6/300	2026-03-21 06:32 by Ecowxq666！
[考研] 華東師范大學-071000生物學-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 271材料工程求調(diào)劑 +8	.6lL 2026-03-18	8/400	2026-03-21 00:58 by JourneyLucky
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-18	6/300	2026-03-21 00:32 by JourneyLucky
[考研] 22408 344分求調(diào)劑一志愿華電計算機技術(shù) +4	solanXXX 2026-03-20	4/200	2026-03-20 23:49 by alg094825
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 304求調(diào)劑 +7	司空. 2026-03-18	7/350	2026-03-20 23:08 by JourneyLucky
[考研] 材料與化工求調(diào)劑 +7	為學666 2026-03-16	7/350	2026-03-19 14:48 by 盡舜堯1
[考研] 085600材料與化工調(diào)劑 324分 +10	llllkkkhh 2026-03-18	12/600	2026-03-19 14:33 by llllkkkhh
[考研] 0703化學調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3	2026eternal 2026-03-18	3/150	2026-03-18 19:09 by 搏擊518
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao

24小時熱門版塊排行榜

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

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

» 本帖已獲得的紅花（最新10朵）

» 本帖附件資源列表

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

» 搶金幣啦！回帖就可以得到:

查看全部散金貼