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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時(shí)
蟲號(hào): 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【原創(chuàng)】ICHARG=11和+U的計(jì)算已有4人參與

vasp計(jì)算DOS有一種方法，就是以較低k網(wǎng)格達(dá)到收斂，再以很高的k網(wǎng)格通過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)行對比。
如下是三次的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)米面以上的部分有所變化。
所以，我猜測可能是因?yàn)槭荌NCAR里面的+U設(shè)置在后來的非自洽計(jì)算中重復(fù)起了作用，使得實(shí)際作用的U值成了4eV。即原來以U=2進(jìn)行了自洽計(jì)算，得到的CHGCAR已經(jīng)包含2eV 的Coulomb repulsion，后來進(jìn)行非自洽計(jì)算時(shí)INCAR中的+U設(shè)置又起了作用，再加上了2eV的Coulomb repulsion，所以出來的DOS結(jié)果和原來U=2時(shí)的不同。
希望有高手來證實(shí)一下，或者解釋一下k21的非自洽結(jié)果與原自洽結(jié)果的DOS不同的原因。

[ Last edited by nkleof on 2010-4-29 at 10:51 ]

回復(fù)此樓

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

材料計(jì)算模擬實(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)價(jià)值貼推薦，對您同樣有幫助:

靜態(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)警告，請高手幫忙已經(jīng)有9人回復(fù)
【求助】VASP計(jì)算體系性能沒有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ì)算能量值錯(cuò)誤已經(jīng)有21人回復(fù)
【求助】LDA+U優(yōu)化問題已經(jīng)有10人回復(fù)
【求助】VASP進(jìn)行DFT+HF計(jì)算的問題？已經(jīng)有7人回復(fù)
【求助】計(jì)算單個(gè)Fe原子，無法收斂，請教已經(jīng)有19人回復(fù)

1樓 2010-04-29 10:49:09

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

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時(shí)
蟲號(hào): 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by nkleof at 2010-04-29 10:49:09:
vasp計(jì)算DOS有一種方法，就是以較低k網(wǎng)格達(dá)到收斂，再以很高的k網(wǎng)格通過ICHARG=11讀取自洽CHGCAR進(jìn)行DOS計(jì)算。
我昨天進(jìn)行了這樣的計(jì)算，但是在計(jì)算中進(jìn)行了+U，如下是INCAR
PREC = Accurate
LREAL = A
IALGO ...

U=4的計(jì)算結(jié)果出來了，證明我前面提到的猜測不正確，U=4的DOS和U=2，ICHARG=11的DOS并不相同。其實(shí)從上個(gè)帖子的DOS圖可以看出來，當(dāng)U值增大到4時(shí)，費(fèi)米面以上的那部分峰肯定會(huì)向著高能量方向移動(dòng)一段距離，而不是幾乎在原地不動(dòng)。
因此還有待高手來說明一下當(dāng)在+U=2的自洽CHGCAR基礎(chǔ)上以ICHARG=11并附加+U計(jì)算標(biāo)簽的情況下，DOS結(jié)果為什么會(huì)發(fā)生改變？

贊一下

回復(fù)此樓

2樓2010-04-29 20:14:05

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

銳利的碎片

木蟲 (正式寫手)

star watcher

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時(shí)
蟲號(hào): 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

理論上來說U應(yīng)該保持一致。

贊一下(1人)

回復(fù)此樓

3樓2010-04-29 20:32:49

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

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時(shí)
蟲號(hào): 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by 銳利的碎片 at 2010-04-29 20:32:49:
理論上來說U應(yīng)該保持一致。

意思是說在自洽CHGCAR的基礎(chǔ)上進(jìn)行ICHARG=11的計(jì)算也需要原來自洽計(jì)算時(shí)的+U設(shè)置？

贊一下

回復(fù)此樓

4樓2010-04-29 23:00:53

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

銳利的碎片

木蟲 (正式寫手)

star watcher

1ST強(qiáng)帖: 1
應(yīng)助: 136 (高中生)
金幣: 3637.1
散金: 252
紅花: 22
帖子: 988
在線: 1224.9小時(shí)
蟲號(hào): 961933
注冊: 2010-03-05
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by nkleof at 2010-04-29 23:00:53:

意思是說在自洽CHGCAR的基礎(chǔ)上進(jìn)行ICHARG=11的計(jì)算也需要原來自洽計(jì)算時(shí)的+U設(shè)置？

對啊，U值最好先算下帶隙看下，不同關(guān)聯(lián)勢都不一樣。一般要求整個(gè)過程都應(yīng)該+U。

贊一下(1人)

回復(fù)此樓

5樓2010-04-29 23:03:42

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

wuchenwf

榮譽(yù)版主 (職業(yè)作家)

應(yīng)助: 0 (幼兒園)
貴賓: 3.433
金幣: 19419.2
散金: 10
紅花: 4
帖子: 3560
在線: 1035.7小時(shí)
蟲號(hào): 398569
注冊: 2007-06-10
性別: GG
專業(yè): 凝聚態(tài)物性I:結(jié)構(gòu)、力學(xué)和
管轄: 第一性原理

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

U目前來講多半還沒有什么理論確定，也就說得我們自己試

贊一下(1人)

回復(fù)此樓

6樓2010-04-29 23:59:59

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

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時(shí)
蟲號(hào): 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by 銳利的碎片 at 2010-04-29 23:03:42:

對啊，U值最好先算下帶隙看下，不同關(guān)聯(lián)勢都不一樣。一般要求整個(gè)過程都應(yīng)該+U。

嗯，不過我想知道的是在+U得到的CHGCAR基礎(chǔ)上做非自洽計(jì)算時(shí)是否還需要+U設(shè)置。

贊一下

回復(fù)此樓

7樓2010-04-30 08:54:07

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

wujiandongde

木蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 4265.4
散金: 30
紅花: 5
帖子: 182
在線: 704.8小時(shí)
蟲號(hào): 941604
注冊: 2010-01-13
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

你好！你在設(shè)置了ISTART=1 ICHARG=11以后，會(huì)不會(huì)有這個(gè)Warning出現(xiàn)呢？WARNING: stress and forces are not correct，我最近計(jì)算遇到了這個(gè)問題，我想確定一下是不是我設(shè)置參數(shù)的問題，謝謝！

贊一下(1人)

回復(fù)此樓

8樓2011-04-11 12:12:43

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

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時(shí)
蟲號(hào): 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by wujiandongde at 2011-04-11 12:12:43:
你好！你在設(shè)置了ISTART=1 ICHARG=11以后，會(huì)不會(huì)有這個(gè)Warning出現(xiàn)呢？WARNING: stress and forces are not correct，我最近計(jì)算遇到了這個(gè)問題，我想確定一下是不是我設(shè)置參數(shù)的問題，謝謝！

我查了一下發(fā)現(xiàn)也有的，這應(yīng)該沒什么問題，因?yàn)槟阍O(shè)置了ICHARG=11后，做的是非自洽計(jì)算，k網(wǎng)格也改掉了，計(jì)算出來的力和能量什么自然也沒什么意義。

贊一下

回復(fù)此樓

9樓2011-04-15 09:49:56

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

wujiandongde

木蟲 (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 4265.4
散金: 30
紅花: 5
帖子: 182
在線: 704.8小時(shí)
蟲號(hào): 941604
注冊: 2010-01-13
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

引用回帖:

Originally posted by nkleof at 2011-04-15 09:49:56:
我查了一下發(fā)現(xiàn)也有的，這應(yīng)該沒什么問題，因?yàn)槟阍O(shè)置了ICHARG=11后，做的是非自洽計(jì)算，k網(wǎng)格也改掉了，計(jì)算出來的力和能量什么自然也沒什么意義。

多謝！我也咨詢了高手，說是可以忽略的。。。

贊一下(1人)

回復(fù)此樓

10樓2011-04-15 14:28:28

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

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

返回列表

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

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

[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 281求調(diào)劑（0805） +3	煙汐憶海 2026-03-16	9/450	2026-03-18 07:23 by 煙汐憶海
[考研] 材料與化工求調(diào)劑 +6	為學(xué)666 2026-03-16	6/300	2026-03-17 20:15 by peike
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 275求調(diào)劑 +4	太陽花天天開心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 考研調(diào)劑 +3	淇ya_~ 2026-03-17	5/250	2026-03-17 09:25 by Winj1e
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8	秋有木北 2026-03-14	8/400	2026-03-16 15:21 by 哦哦123
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 材料工程調(diào)劑 +9	咪咪空空 2026-03-12	9/450	2026-03-13 22:05 by 星空星月
[考研] 求材料調(diào)劑 +5	隔壁陳先生 2026-03-12	5/250	2026-03-13 22:03 by 星空星月
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 0703化學(xué)求調(diào)劑 +7	綠豆芹菜湯 2026-03-12	7/350	2026-03-13 17:25 by njzyff
[考研] 0856化學(xué)工程280分求調(diào)劑 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor
[考研] 270求調(diào)劑 085600材料與化工專碩 +3	YXCT 2026-03-11	3/150	2026-03-13 10:13 by houyaoxu
[考博] 2026年博士申請 +3	QwQwQW10 2026-03-11	3/150	2026-03-12 17:58 by gxch43

24小時(shí)熱門版塊排行榜

[交流] 【原創(chuàng)】ICHARG=11和+U的計(jì)算 已有4人參與

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

» 猜你喜歡

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

[交流] 【原創(chuàng)】ICHARG=11和+U的計(jì)算已有4人參與

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