 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 3968 | 回復(fù): 6 | ||
[求助]
HSE06做自洽計算時提示增大NBANDS，怎么做？ 已有1人參與
|
|
用HSE06計算時，輸出結(jié)果里出現(xiàn)了如下的錯誤，提示增加空帶NBANDS。 我的INCAR里沒有設(shè)置NBANDS。 用的是默認(rèn)值， 我看了下OUTCAR里 NBANDS=144，根據(jù)手冊NBANDS = NION/2 + ELECTRON/2 ,我計算的提醒，NIONS=48 ELECTRON =200 ，帶入公式計算 NBANDS =124 ，為什么OUTCAR 里會顯示144呢？ 再者，如果我要在NBANDS=144的基礎(chǔ)上增大，我需要增大到多少較為合適呢？ 請各位大俠知道，金幣奉上，感激不盡，在錯誤的后面我會把我的INCAR,KPOINTS 貼出來，請大家看下有啥子問題？ ----------------------------------------------------------------------------- | | | 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ù)此樓» 猜你喜歡
0703化學(xué)調(diào)劑
已經(jīng)有12人回復(fù)
-
0703化學(xué) 305求調(diào)劑
已經(jīng)有4人回復(fù)
-
0703化學(xué)調(diào)劑，求各位老師收留
已經(jīng)有10人回復(fù)
-
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 關(guān)于（VASP）Fe100表面的SCF收斂
已經(jīng)有19人回復(fù)
- GW準(zhǔn)粒子計算會高估能隙？
已經(jīng)有21人回復(fù)
- 相同的vasp版本，不同機(jī)器上計算同一個體系，能量相差近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小時
- 蟲號: 1007127
- 注冊: 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并不是簡單地按這個公式NION/2 + ELECTRON/2來計算的，它在OUTCAR里輸出的值還取決其他計算參數(shù)的設(shè)置。 比如是否非共線磁性參數(shù)(LNONCOLLINEAR)、NPAR、ISPIN的值 你可以簡單的讀讀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），增加的太大的話，會增加計算量而已。 |
2樓2014-11-30 08:37:10
|
萬分榮幸得到候老師的回答，感謝候老師的回答。候老師，我有兩個問題請教您： 1,我后來重新算了，NBANDS 設(shè)置為默認(rèn)值的2.5倍，默認(rèn)144，設(shè)置成了300（根據(jù)網(wǎng)上經(jīng)驗(yàn)帖子），發(fā)現(xiàn)計算機(jī)運(yùn)行很慢了。我是在長沙超算上算的。10個節(jié)點(diǎn)，共120個核做的HSE06自洽計算，24個小時只跑了兩步，我想請問候老師；這正常嗎？ 2，我看別人HSE06自洽的INCAR中，設(shè)置了力的收斂精度和能量的收斂精度，自洽計算不是重已有的結(jié)構(gòu)出發(fā)直接算電荷和波函數(shù)嗎？為何需要設(shè)置精度。 3，候老師，我想確認(rèn)下HSE06計算能帶步驟 DFT自洽 + HSE自洽+ 能帶計算 對嗎？ 謝謝您！ |
3樓2014-11-30 13:12:36
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
300設(shè)置的似乎太多了，會顯著得增加計算量。除非你想得到導(dǎo)帶很大能量范圍內(nèi)的信息，否則的話，建議減少它。 “10個節(jié)點(diǎn)，共120個核做的HSE06自洽計算，24個小時只跑了兩步，我想請問候老師；這正常嗎？” 這種問題，沒有辦法回答。盡量不要問這種問題。 “自洽計算不是重已有的結(jié)構(gòu)出發(fā)直接算電荷和波函數(shù)嗎？為何需要設(shè)置精度�！� 這說明你根本沒有理解“結(jié)構(gòu)優(yōu)化”和“自洽計算”的差別。 自洽計算---》 對某個結(jié)構(gòu)進(jìn)行電子波函數(shù)的自洽計算，即根據(jù)初始電荷密度（這個可以是已自洽過的，或未自洽的，或原子電荷密度），構(gòu)造Kohn-Sham哈密頓量，求解本征波函數(shù)和本征值，根據(jù)本征波函數(shù)得到新的電荷密度，再重新構(gòu)造Kohn-Sham哈密頓量，重復(fù)前兩步直到所指定的物理量（總能或電密度等或波函數(shù)）的變化小于所設(shè)置的精度�？纯从嘘P(guān)Kohn-Sham方程自洽求解的簡單介紹就了解。這個可以在很多資料中找到。 做自洽計算因此當(dāng)然要設(shè)置收斂精度了。 結(jié)構(gòu)優(yōu)化是將對某個結(jié)構(gòu)的原子位置或晶格基矢進(jìn)行調(diào)整使得原子結(jié)構(gòu)中原子所受的力小于某個值，晶格的應(yīng)力小于某個值，總能變化小于某個值。在結(jié)構(gòu)優(yōu)化中不斷調(diào)整結(jié)構(gòu)，對其中的每個結(jié)構(gòu)都會有電子波函數(shù)的自洽計算后，再計算原子所受的力的步驟。 “我想確認(rèn)下HSE06計算能帶步驟 DFT自洽 + HSE自洽+ 能帶計算 對嗎？” 這種確認(rèn)從手冊里就可以找到，手冊上也明確建議了計算步驟，你順著上面的就對了。 |
4樓2014-11-30 13:42:03
valenhou001
至尊木蟲 (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時
- 蟲號: 1007127
- 注冊: 2010-04-27
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
5樓2014-11-30 13:53:19
6樓2014-11-30 15:52:44
7樓2018-05-09 11:42:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zyqh 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 281求調(diào)劑（0805） +6 | 煙汐憶海 2026-03-16 | 16/800 |
|
|
|
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4 | 守候夕陽CF 2026-03-18 | 4/200 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 9/450 |
|
|
|
[考研] 0854可跨調(diào)劑，一作一項(xiàng)核心論文五項(xiàng)專利，省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 總分331 +3 | ZY-05 2026-03-13 | 3/150 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +6 | 孔志浩 2026-03-12 | 11/550 |
|
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對?/a> 2026-03-13 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy！ 2026-03-13 | 6/300 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] 學(xué)碩285求調(diào)劑 +13 | Wisjxn 2026-03-12 | 46/2300 |
|
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +4 | 芬達(dá)46 2026-03-12 | 4/200 |
|