 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (3734)
- >考研 (864)
- >導(dǎo)師招生 (727)
- >蟲(chóng)友互識(shí) (232)
- >文獻(xiàn)求助 (106)
- >休閑灌水 (87)
- >碩博家園 (58)
- >招聘信息布告欄 (41)
- >考博 (40)
- >綠色求助(高懸賞) (24)
- >博后之家 (22)
- >論文投稿 (18)
- >基金申請(qǐng) (15)
- >教師之家 (9)
- >攻關(guān)文獻(xiàn)(高獎(jiǎng)勵(lì)) (8)
- >外文書(shū)籍求助 (8)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4005 | 回復(fù): 6 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
zyqh鐵桿木蟲(chóng) (著名寫(xiě)手)
|
[求助]
HSE06做自洽計(jì)算時(shí)提示增大NBANDS，怎么做？ 已有1人參與
|
|
|
用HSE06計(jì)算時(shí)，輸出結(jié)果里出現(xiàn)了如下的錯(cuò)誤，提示增加空帶NBANDS。 我的INCAR里沒(méi)有設(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 貼出來(lái)，請(qǐng)大家看下有啥子問(wèn)題？ ----------------------------------------------------------------------------- | | | 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ù)此樓» 猜你喜歡
311求調(diào)劑
已經(jīng)有6人回復(fù)
-
289求調(diào)劑
已經(jīng)有25人回復(fù)
-
329求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料考研調(diào)劑
已經(jīng)有9人回復(fù)
-
327求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有18人回復(fù)
-
337求調(diào)劑
已經(jīng)有11人回復(fù)
-
301求調(diào)劑
已經(jīng)有5人回復(fù)
-
275求調(diào)劑
已經(jīng)有11人回復(fù)
-
292求調(diào)劑
已經(jīng)有14人回復(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ù)
2120160150
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6.3
- 帖子: 9
- 在線: 8.8小時(shí)
- 蟲(chóng)號(hào): 5161786
- 注冊(cè): 2016-10-25
- 性別: MM
7樓2018-05-09 11:42:37
valenhou001
至尊木蟲(chóng) (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲(chóng)號(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來(lái)計(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
zyqh
鐵桿木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 33 (小學(xué)生)
- 金幣: 6708
- 散金: 1454
- 紅花: 22
- 帖子: 2377
- 在線: 593.4小時(shí)
- 蟲(chóng)號(hào): 2762045
- 注冊(cè): 2013-10-29
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
萬(wàn)分榮幸得到候老師的回答，感謝候老師的回答。候老師，我有兩個(gè)問(wèn)題請(qǐng)教您： 1,我后來(lái)重新算了，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)問(wèn)候老師；這正常嗎？ 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
至尊木蟲(chóng) (職業(yè)作家)
- 1ST強(qiáng)帖: 13
- 應(yīng)助: 241 (大學(xué)生)
- 金幣: 25701.7
- 散金: 602
- 紅花: 166
- 帖子: 3782
- 在線: 873.8小時(shí)
- 蟲(chóng)號(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)問(wèn)候老師；這正常嗎？” 這種問(wèn)題，沒(méi)有辦法回答。盡量不要問(wèn)這種問(wèn)題。 “自洽計(jì)算不是重已有的結(jié)構(gòu)出發(fā)直接算電荷和波函數(shù)嗎？為何需要設(shè)置精度�！� 這說(shuō)明你根本沒(méi)有理解“結(jié)構(gòu)優(yōu)化”和“自洽計(jì)算”的差別。 自洽計(jì)算---》 對(duì)某個(gè)結(jié)構(gòu)進(jìn)行電子波函數(shù)的自洽計(jì)算，即根據(jù)初始電荷密度（這個(gè)可以是已自洽過(guò)的，或未自洽的，或原子電荷密度），構(gòu)造Kohn-Sham哈密頓量，求解本征波函數(shù)和本征值，根據(jù)本征波函數(shù)得到新的電荷密度，再重新構(gòu)造Kohn-Sham哈密頓量，重復(fù)前兩步直到所指定的物理量（總能或電密度等或波函數(shù)）的變化小于所設(shè)置的精度�？纯从嘘P(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
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 324分 085600材料與化工 +20 | 呆鵝oor 2026-03-27 | 20/1000 |
|
|---|---|---|---|---|
|
|
[考研] 考研調(diào)劑 +12 | Amber00 2026-03-31 | 12/600 |
|
|
[考研] 生物學(xué)327，求調(diào)劑 +5 | 書(shū)上的梅子 2026-04-01 | 6/300 |
|
|
|
[考博] 26年申博 +3 | staryer 2026-03-30 | 4/200 |
|
|
|
[考研] 11408 321分求調(diào)劑 +3 | huchun12138 2026-03-30 | 4/200 |
|
|
|
[考研] 材料科學(xué)與工程339求調(diào)劑 +11 | hyz0119 2026-03-31 | 12/600 |
|
|
|
[考研] 362求調(diào)劑 +13 | 西南交材料專碩3 2026-03-31 | 13/650 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +13 | 大饞小子 2026-03-27 | 14/700 |
|
|
|
[考研] 291求調(diào)劑 +3 | 迷蒙木木 2026-04-01 | 4/200 |
|
|
|
[考研] 調(diào)劑申請(qǐng) +8 | 張張張張zy 2026-03-31 | 9/450 |
|
|
|
[考研] 08開(kāi)頭看過(guò)來(lái)！��！ +3 | wwwwffffff 2026-03-31 | 5/250 |
|
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分，前三科223分 +11 | dongfang59 2026-03-27 | 11/550 |
|
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +6 | hanamiko 2026-03-29 | 6/300 |
|
|
|
[基金申請(qǐng)] 面上5B能上會(huì)嗎？ +8 | redcom 2026-03-29 | 8/400 |
|
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +11 | 多久上課 2026-03-27 | 12/600 |
|
|
|
[考研] 279求調(diào)劑 +12 | j的立方 2026-03-29 | 12/600 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂(lè)樂(lè)樂(lè) 2026-03-26 | 4/200 |
|