 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2308 | 回復(fù): 9 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
nkleof銀蟲 (正式寫手)
|
[交流]
【原創(chuàng)】ICHARG=11和+U的計(jì)算 已有4人參與
|
||||
|
vasp計(jì)算DOS有一種方法，就是以較低k網(wǎng)格達(dá)到收斂，再以很高的k網(wǎng)格通過(guò)ICHARG=11讀取自洽CHGCAR進(jìn)行DOS計(jì)算。 我昨天進(jìn)行了這樣的計(jì)算，但是在計(jì)算中進(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ì)算進(jìn)行對(duì)比。 如下是三次的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é)果和其他兩個(gè)有點(diǎn)不同，費(fèi)米面以上的部分有所變化。 所以，我猜測(cè)可能是因?yàn)槭荌NCAR里面的+U設(shè)置在后來(lái)的非自洽計(jì)算中重復(fù)起了作用，使得實(shí)際作用的U值成了4eV。即原來(lái)以U=2進(jìn)行了自洽計(jì)算，得到的CHGCAR已經(jīng)包含2eV 的Coulomb repulsion，后來(lái)進(jìn)行非自洽計(jì)算時(shí)INCAR中的+U設(shè)置又起了作用，再加上了2eV的Coulomb repulsion，所以出來(lái)的DOS結(jié)果和原來(lái)U=2時(shí)的不同。 希望有高手來(lái)證實(shí)一下，或者解釋一下k21的非自洽結(jié)果與原自洽結(jié)果的DOS不同的原因。 [ Last edited by nkleof on 2010-4-29 at 10:51 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 材料計(jì)算模擬實(shí)用技巧 | VASP |
» 猜你喜歡
0703化學(xué)調(diào)劑
已經(jīng)有7人回復(fù)
-
0854，計(jì)算機(jī)類招收調(diào)劑
已經(jīng)有6人回復(fù)
-
【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級(jí)475
已經(jīng)有6人回復(fù)
-
0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國(guó)級(jí)證書40+數(shù)一英一287
已經(jīng)有4人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有8人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有3人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有11人回復(fù)
-
299求調(diào)劑
已經(jīng)有5人回復(fù)
-
收復(fù)試調(diào)劑生
已經(jīng)有3人回復(fù)
-
280求調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 靜態(tài)計(jì)算能收斂，同樣的參數(shù)弛豫不能收斂，這是為什么？
已經(jīng)有4人回復(fù)
- vasp5.2計(jì)算半導(dǎo)體Si光學(xué)性質(zhì)的總結(jié)
已經(jīng)有106人回復(fù)
- 態(tài)密度圖跟文獻(xiàn)里面的怎么都不一樣，求救
已經(jīng)有36人回復(fù)
- ISPIN=1時(shí)態(tài)密度有為負(fù)的區(qū)域，為嘛呢？
已經(jīng)有6人回復(fù)
- 【求助成功】DOS計(jì)算ICHARG=11出現(xiàn)警告，請(qǐng)高手幫忙
已經(jīng)有9人回復(fù)
- 【求助】VASP計(jì)算體系性能沒(méi)有f能帶電子
已經(jīng)有24人回復(fù)
- 【求助】態(tài)密度困惑
已經(jīng)有8人回復(fù)
- 【求助】關(guān)于擬合得到體彈性模量&自旋極化&LDA+U等問(wèn)題，望大家指點(diǎn)指點(diǎn)
已經(jīng)有4人回復(fù)
- 【求助】vasp靜態(tài)計(jì)算能量值錯(cuò)誤
已經(jīng)有21人回復(fù)
- 【求助】LDA+U優(yōu)化問(wèn)題
已經(jīng)有10人回復(fù)
- 【求助】VASP進(jìn)行DFT+HF計(jì)算的問(wèn)題？
已經(jīng)有7人回復(fù)
- 【求助】計(jì)算單個(gè)Fe原子，無(wú)法收斂，請(qǐng)教
已經(jīng)有19人回復(fù)
nkleof
銀蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時(shí)
- 蟲號(hào): 152581
- 注冊(cè): 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
7樓2010-04-30 08:54:07
nkleof
銀蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時(shí)
- 蟲號(hào): 152581
- 注冊(cè): 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
U=4的計(jì)算結(jié)果出來(lái)了，證明我前面提到的猜測(cè)不正確，U=4的DOS和U=2，ICHARG=11的DOS并不相同。其實(shí)從上個(gè)帖子的DOS圖可以看出來(lái)，當(dāng)U值增大到4時(shí)，費(fèi)米面以上的那部分峰肯定會(huì)向著高能量方向移動(dòng)一段距離，而不是幾乎在原地不動(dòng)。 因此還有待高手來(lái)說(shuō)明一下當(dāng)在+U=2的自洽CHGCAR基礎(chǔ)上以ICHARG=11并附加+U計(jì)算標(biāo)簽的情況下，DOS結(jié)果為什么會(huì)發(fā)生改變？ |
2樓2010-04-29 20:14:05
銳利的碎片
木蟲 (正式寫手)
star watcher
- 1ST強(qiáng)帖: 1
- 應(yīng)助: 136 (高中生)
- 金幣: 3637.1
- 散金: 252
- 紅花: 22
- 帖子: 988
- 在線: 1224.9小時(shí)
- 蟲號(hào): 961933
- 注冊(cè): 2010-03-05
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-04-29 20:32:49
nkleof
銀蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 2093.2
- 散金: 649
- 紅花: 5
- 帖子: 813
- 在線: 165.9小時(shí)
- 蟲號(hào): 152581
- 注冊(cè): 2006-01-02
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
4樓2010-04-29 23:00:53
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +4 | 煙汐憶海 2026-03-16 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過(guò)，有科研經(jīng)歷 +9 | 曦熙兮 2026-03-15 | 9/450 |
|
|
|
[考研] 297求調(diào)劑 +4 | 戲精丹丹丹 2026-03-17 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[基金申請(qǐng)] 被我言中：新模板不強(qiáng)調(diào)格式了，假專家開(kāi)始管格式了 +4 | beefly 2026-03-14 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運(yùn)連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 080500，材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12 | 11啦啦啦 2026-03-11 | 35/1750 |
|
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
|
[考研] 工科278分求調(diào)劑 +5 | 周慢熱啊 2026-03-12 | 7/350 |
|
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|