 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2301 | 回復(fù): 9 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
nkleof銀蟲 (正式寫手)
|
[交流]
【原創(chuàng)】ICHARG=11和+U的計算 已有4人參與
|
||||
|
vasp計算DOS有一種方法，就是以較低k網(wǎng)格達(dá)到收斂，再以很高的k網(wǎng)格通過ICHARG=11讀取自洽CHGCAR進(jìn)行DOS計算。 我昨天進(jìn)行了這樣的計算，但是在計算中進(jìn)行了+U，如下是INCAR PREC = Accurate LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 ISTART = 1 ICHARG = 11 NELMIN = 4 LORBIT = 12 EDIFF = 1E-5 ISMEAR = -5 EMIN = -20 EMAX = 20 NEDOS = 11000 ISPIN = 2 LDAU = .TRUE. LDAUTYPE = 2 LDAUL = -1 2 -1 LDAUU = 0 2 0 LDAUJ = 0 0 0 LDAUPRINT = 2 其中自洽收斂是以11 11 11網(wǎng)格，非自洽DOS是以21 21 21網(wǎng)格，另外也進(jìn)行了21 21 21網(wǎng)格的自洽計算進(jìn)行對比。 如下是三次的DOS 圖，分別為k11網(wǎng)格自洽的DOS，在k11網(wǎng)格自洽基礎(chǔ)上k21的DOS結(jié)果，k21網(wǎng)格自洽的DOS    從圖上可以看到，在k11自洽基礎(chǔ)上進(jìn)行k21的ICHARG=11的DOS結(jié)果和其他兩個有點(diǎn)不同，費(fèi)米面以上的部分有所變化。 所以，我猜測可能是因為是INCAR里面的+U設(shè)置在后來的非自洽計算中重復(fù)起了作用，使得實(shí)際作用的U值成了4eV。即原來以U=2進(jìn)行了自洽計算，得到的CHGCAR已經(jīng)包含2eV 的Coulomb repulsion，后來進(jìn)行非自洽計算時INCAR中的+U設(shè)置又起了作用，再加上了2eV的Coulomb repulsion，所以出來的DOS結(jié)果和原來U=2時的不同。 希望有高手來證實(shí)一下，或者解釋一下k21的非自洽結(jié)果與原自洽結(jié)果的DOS不同的原因。 [ Last edited by nkleof on 2010-4-29 at 10:51 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 材料計算模擬實(shí)用技巧 | VASP |
» 猜你喜歡
焦慮
已經(jīng)有10人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
297求調(diào)劑
已經(jīng)有3人回復(fù)
-
302求調(diào)劑
已經(jīng)有8人回復(fù)
-
環(huán)境工程調(diào)劑
已經(jīng)有7人回復(fù)
-
268求調(diào)劑
已經(jīng)有6人回復(fù)
-
311求調(diào)劑
已經(jīng)有11人回復(fù)
-
工科材料085601 279求調(diào)劑
已經(jīng)有4人回復(fù)
-
301求調(diào)劑
已經(jīng)有9人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 靜態(tài)計算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經(jīng)有4人回復(fù)
- vasp5.2計算半導(dǎo)體Si光學(xué)性質(zhì)的總結(jié)
已經(jīng)有106人回復(fù)
- 態(tài)密度圖跟文獻(xiàn)里面的怎么都不一樣，求救
已經(jīng)有36人回復(fù)
- ISPIN=1時態(tài)密度有為負(fù)的區(qū)域，為嘛呢？
已經(jīng)有6人回復(fù)
- 【求助成功】DOS計算ICHARG=11出現(xiàn)警告，請高手幫忙
已經(jīng)有9人回復(fù)
- 【求助】VASP計算體系性能沒有f能帶電子
已經(jīng)有24人回復(fù)
- 【求助】態(tài)密度困惑
已經(jīng)有8人回復(fù)
- 【求助】關(guān)于擬合得到體彈性模量&自旋極化&LDA+U等問題，望大家指點(diǎn)指點(diǎn)
已經(jīng)有4人回復(fù)
- 【求助】vasp靜態(tài)計算能量值錯誤
已經(jīng)有21人回復(fù)
- 【求助】LDA+U優(yōu)化問題
已經(jīng)有10人回復(fù)
- 【求助】VASP進(jìn)行DFT+HF計算的問題？
已經(jīng)有7人回復(fù)
- 【求助】計算單個Fe原子，無法收斂，請教
已經(jīng)有19人回復(fù)
nkleof
銀蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時
- 蟲號: 152581
- 注冊: 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2010-04-29 23:00:53
nkleof
銀蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時
- 蟲號: 152581
- 注冊: 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2010-04-29 20:14:05
銳利的碎片
木蟲 (正式寫手)
star watcher
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 136 (高中生)
- 金幣: 3637.1
- 散金: 252
- 紅花: 22
- 帖子: 988
- 在線: 1224.9小時
- 蟲號: 961933
- 注冊: 2010-03-05
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-29 20:32:49
銳利的碎片
木蟲 (正式寫手)
star watcher
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 136 (高中生)
- 金幣: 3637.1
- 散金: 252
- 紅花: 22
- 帖子: 988
- 在線: 1224.9小時
- 蟲號: 961933
- 注冊: 2010-03-05
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2010-04-29 23:03:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +6 | 簡單點(diǎn)0 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 工科材料085601 279求調(diào)劑 +4 | 困于星晨 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +9 | 小熊joy 2026-03-14 | 9/450 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 321求調(diào)劑 +5 | 大米飯！ 2026-03-15 | 5/250 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +3 | 一顆meteor 2026-03-13 | 3/150 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +7 | 綠豆芹菜湯 2026-03-12 | 7/350 |
|
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|