 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 7 | 1/1 | 返回列表 |
| 查看: 3970 | 回復(fù): 6 | |||
[求助]
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ù)此樓» 猜你喜歡
288求調(diào)劑，一志愿華南理工大學(xué)071005
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有5人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有12人回復(fù)
-
0856調(diào)劑，是學(xué)校就去
已經(jīng)有5人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有13人回復(fù)
-
梁成偉老師課題組歡迎你的加入
已經(jīng)有11人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑 ，六級(jí)已過，有科研經(jīng)歷
已經(jīng)有11人回復(fù)
-
321求調(diào)劑
已經(jīng)有10人回復(fù)
-
東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng)
已經(jīng)有8人回復(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)
【答案】應(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)
|
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)
5樓2014-11-30 13:53:19
6樓2014-11-30 15:52:44
2120160150
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6.3
- 帖子: 9
- 在線: 8.8小時(shí)
- 蟲號(hào): 5161786
- 注冊(cè): 2016-10-25
- 性別: MM
7樓2018-05-09 11:42:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zyqh 的主題更新
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +12 | 陌の森林 2026-03-18 | 12/600 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿 西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +3 | 晨昏線與星海 2026-03-19 | 3/150 |
|
|
|
[考研] 346求調(diào)劑[0856] +3 | WayneLim327 2026-03-16 | 6/300 |
|
|
|
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +10 | SUICHILD 2026-03-12 | 10/500 |
|
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
|
[考研] 求材料調(diào)劑 +10 | 隔壁陳先生 2026-03-12 | 10/500 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7 | ZY-05 2026-03-13 | 8/400 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 268求調(diào)劑 +8 | 一定有學(xué)上- 2026-03-14 | 9/450 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿211 0703方向310分求調(diào)劑 +3 | 努力奮斗112 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|