 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3995 | 回復(fù): 6 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
HSE06做自洽計(jì)算時(shí)提示增大NBANDS，怎么做？ 已有1人參與
|
|||
|
用HSE06計(jì)算時(shí)，輸出結(jié)果里出現(xiàn)了如下的錯(cuò)誤，提示增加空帶NBANDS。 我的INCAR里沒有設(shè)置NBANDS。 用的是默認(rèn)值， 我看了下OUTCAR里 NBANDS=144，根據(jù)手冊(cè)NBANDS = NION/2 + ELECTRON/2 ,我計(jì)算的提醒，NIONS=48 ELECTRON =200 ，帶入公式計(jì)算 NBANDS =124 ，為什么OUTCAR 里會(huì)顯示144呢？ 再者，如果我要在NBANDS=144的基礎(chǔ)上增大，我需要增大到多少較為合適呢？ 請(qǐng)各位大俠知道，金幣奉上，感激不盡，在錯(cuò)誤的后面我會(huì)把我的INCAR,KPOINTS 貼出來，請(qǐng)大家看下有啥子問題？ ----------------------------------------------------------------------------- | | | W W AA RRRRR N N II N N GGGG !!! | | W W A A R R NN N II NN N G G !!! | | W W A A R R N N N II N N N G !!! | | W WW W AAAAAA RRRRR N N N II N N N G GGG ! | | WW WW A A R R N NN II N NN G G | | W W A A R R N N II N N GGGG !!! | | | | Your highest band is occupied at some k-points! Unless you are | | performing a calculation for an insulator or semiconductor, without | | unoccupied bands, you have included TOO FEW BANDS!! Please increase | | the parameter NBANDS in file 'INCAR' to ensure that the highest band | | is unoccupied at all k-points. It is always recommended to | | include a few unoccupied bands to accelerate the convergence of | | molecular dynamics runs (even for insulators or semiconductors). | | Because the presence of unoccupied bands improves wavefunction | | prediction, and helps to suppress 'band-crossings.' | | | ----------------------------------------------------------------------------- INCAR： ISTART = 1 ICHARG = 2 EDIFF = 0.00001 EDIFFG = -0.001 ENCUT = 400 LREAL = .FALSE. LWAVE = .TRUE. LCHARG = .TRUE. NELM = 200 NSW = 0 IBRION = -1 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.1 NSIM = 4 ISYM = 0 LHFCALC = .TRUE. HFSCREEN = 0.2 ALGO = D TIME = 0.4 ENCUTFOCK = 0 AEXX = 0.25 KPOINTS: Automatic generation 0 Monhkorst-Pack 4 4 4 0.0 0.0 0.0 |
回復(fù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有13人回復(fù)
-
071010 323 分求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料調(diào)劑
已經(jīng)有9人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有3人回復(fù)
-
一志愿南昌大學(xué)324求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 關(guān)于（VASP）Fe100表面的SCF收斂
已經(jīng)有19人回復(fù)
- GW準(zhǔn)粒子計(jì)算會(huì)高估能隙？
已經(jīng)有21人回復(fù)
- 相同的vasp版本，不同機(jī)器上計(jì)算同一個(gè)體系，能量相差近10ev
已經(jīng)有17人回復(fù)
- 【已經(jīng)完成】Optical properties and density functional perturbation theory 的翻譯
已經(jīng)有50人回復(fù)
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
300設(shè)置的似乎太多了，會(huì)顯著得增加計(jì)算量。除非你想得到導(dǎo)帶很大能量范圍內(nèi)的信息，否則的話，建議減少它。 “10個(gè)節(jié)點(diǎn)，共120個(gè)核做的HSE06自洽計(jì)算，24個(gè)小時(shí)只跑了兩步，我想請(qǐng)問候老師；這正常嗎？” 這種問題，沒有辦法回答。盡量不要問這種問題。 “自洽計(jì)算不是重已有的結(jié)構(gòu)出發(fā)直接算電荷和波函數(shù)嗎？為何需要設(shè)置精度。” 這說明你根本沒有理解“結(jié)構(gòu)優(yōu)化”和“自洽計(jì)算”的差別。 自洽計(jì)算---》 對(duì)某個(gè)結(jié)構(gòu)進(jìn)行電子波函數(shù)的自洽計(jì)算，即根據(jù)初始電荷密度（這個(gè)可以是已自洽過的，或未自洽的，或原子電荷密度），構(gòu)造Kohn-Sham哈密頓量，求解本征波函數(shù)和本征值，根據(jù)本征波函數(shù)得到新的電荷密度，再重新構(gòu)造Kohn-Sham哈密頓量，重復(fù)前兩步直到所指定的物理量（總能或電密度等或波函數(shù)）的變化小于所設(shè)置的精度。看看有關(guān)Kohn-Sham方程自洽求解的簡(jiǎn)單介紹就了解。這個(gè)可以在很多資料中找到。 做自洽計(jì)算因此當(dāng)然要設(shè)置收斂精度了。 結(jié)構(gòu)優(yōu)化是將對(duì)某個(gè)結(jié)構(gòu)的原子位置或晶格基矢進(jìn)行調(diào)整使得原子結(jié)構(gòu)中原子所受的力小于某個(gè)值，晶格的應(yīng)力小于某個(gè)值，總能變化小于某個(gè)值。在結(jié)構(gòu)優(yōu)化中不斷調(diào)整結(jié)構(gòu)，對(duì)其中的每個(gè)結(jié)構(gòu)都會(huì)有電子波函數(shù)的自洽計(jì)算后，再計(jì)算原子所受的力的步驟。 “我想確認(rèn)下HSE06計(jì)算能帶步驟 DFT自洽 + HSE自洽+ 能帶計(jì)算 對(duì)嗎？” 這種確認(rèn)從手冊(cè)里就可以找到，手冊(cè)上也明確建議了計(jì)算步驟，你順著上面的就對(duì)了。 |
4樓2014-11-30 13:42:03
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
zyqh: 金幣+5 2014-11-30 12:52:58
感謝參與，應(yīng)助指數(shù) +1
zyqh: 金幣+5 2014-11-30 12:52:58
|
默認(rèn)的NBANDS并不是簡(jiǎn)單地按這個(gè)公式NION/2 + ELECTRON/2來計(jì)算的，它在OUTCAR里輸出的值還取決其他計(jì)算參數(shù)的設(shè)置。 比如是否非共線磁性參數(shù)(LNONCOLLINEAR)、NPAR、ISPIN的值 你可以簡(jiǎn)單的讀讀main.F里有關(guān)NBANDS變量的賦值. .... NBANDS=MAX(NINT(INFO%NELECT+2)/2+MAX(T_INFO%NIONS/2,3),INT(0.6*INFO%NELECT))+NMAG IF (WDES%LNONCOLLINEAR) NBANDS = NBANDS*2 NBANDS=((NBANDS+NPAR-1)/NPAR)*NPAR ... 你可在NBANDS=144的基礎(chǔ)在適當(dāng)?shù)脑黾右恍?0~20），增加的太大的話，會(huì)增加計(jì)算量而已。 |
2樓2014-11-30 08:37:10
|
萬分榮幸得到候老師的回答，感謝候老師的回答。候老師，我有兩個(gè)問題請(qǐng)教您： 1,我后來重新算了，NBANDS 設(shè)置為默認(rèn)值的2.5倍，默認(rèn)144，設(shè)置成了300（根據(jù)網(wǎng)上經(jīng)驗(yàn)帖子），發(fā)現(xiàn)計(jì)算機(jī)運(yùn)行很慢了。我是在長(zhǎng)沙超算上算的。10個(gè)節(jié)點(diǎn)，共120個(gè)核做的HSE06自洽計(jì)算，24個(gè)小時(shí)只跑了兩步，我想請(qǐng)問候老師；這正常嗎？ 2，我看別人HSE06自洽的INCAR中，設(shè)置了力的收斂精度和能量的收斂精度，自洽計(jì)算不是重已有的結(jié)構(gòu)出發(fā)直接算電荷和波函數(shù)嗎？為何需要設(shè)置精度。 3，候老師，我想確認(rèn)下HSE06計(jì)算能帶步驟 DFT自洽 + HSE自洽+ 能帶計(jì)算 對(duì)嗎？ 謝謝您！ |
3樓2014-11-30 13:12:36
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲號(hào): 1007127
- 注冊(cè): 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2014-11-30 13:53:19
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程調(diào)劑 +5 | 深V宿舍吧 2026-03-30 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 332求調(diào)劑 +12 | 032500 2026-03-25 | 12/600 |
|
|
|
[考研] 材料工程專碩求調(diào)劑 +6 | hyl3153942 2026-03-29 | 6/300 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +12 | 李李不服輸 2026-03-26 | 12/600 |
|
|
|
[考研] 085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富 +6 | 來財(cái)aa 2026-03-25 | 6/300 |
|
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 320分，材料與化工專業(yè)，求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
|
[考研]
 304材料求調(diào)劑
+4
|
鐘llll 2026-03-26 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪！？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|