| 4 | 1/1 | 返回列表 |
| 查看: 844 | 回復(fù): 3 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者socchu將贈(zèng)送您 5 個(gè)金幣 | |||
socchu金蟲(chóng) (正式寫(xiě)手)
|
[求助]
Hubbard manifold with zero occupations is not allowed 已有1人參與
|
||
|
各位老師好: 在用qe計(jì)算能帶的時(shí)候報(bào)錯(cuò)了,請(qǐng)老師們幫忙看看。 工作時(shí)在進(jìn)行一個(gè)有機(jī)金屬mn的雜化材料的能帶計(jì)算,由于有mn,開(kāi)啟了自旋,設(shè)置了u值,自旋和u值參數(shù)如下: / &system ibrav=0, nat= 132, ntyp= 6, ecutwfc = 50, ecutrho = 200, nbnd= 428 nspin=2 tot_magnetization=20 occupations='fixed', hubbard_occ(2,1)=5 。。。。。。 hubbard (ortho-atomic) !加u的一些操作 u mn-3d 2.5 出現(xiàn)的狀況是,進(jìn)行scf和bands的時(shí)候(pw.x)的時(shí)候是合適的,但是在進(jìn)行能帶解析的時(shí)候(bands.x)的時(shí)候報(bào)錯(cuò)了,錯(cuò)誤如下: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% error in routine offset_atom_wfc (1): hubbard manifold with zero occupations is not allowed %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% stopping ... 不知道是什么 原因。 |
木蟲(chóng) (正式寫(xiě)手)
|
這個(gè)錯(cuò)誤通常出現(xiàn)在Quantum ESPRESSO的DFT+U計(jì)算中,尤其是后處理步驟如bands.x時(shí),pw.x的SCF和nscf/bands計(jì)算正常,但后處理檢測(cè)到Hubbard manifold(這里是Mn-3d軌道)有零占據(jù),這不符合代碼要求。基于類(lèi)似案例分析,以下是可能原因和解決方案建議。你可以逐一嘗試,優(yōu)先從簡(jiǎn)單調(diào)整開(kāi)始。如果問(wèn)題持續(xù),提供完整輸入文件(包括ATOMIC_SPECIES、HUBBARD卡和偽勢(shì)細(xì)節(jié))有助于進(jìn)一步診斷。 可能原因 軟件bug或兼容性問(wèn)題:在QE 7.1及早期版本中,后處理工具(如bands.x、dos.x、projwfc.x)對(duì)Hubbard處理的bug會(huì)導(dǎo)致零占據(jù)錯(cuò)誤,尤其當(dāng)原子類(lèi)型順序或多manifold時(shí)。你的設(shè)置(nspin=2, Hubbard ortho-atomic, U Mn-3d=2.5)是標(biāo)準(zhǔn)的collinear DFT+U,但后處理可能未正確繼承占據(jù)。 Hubbard_occ設(shè)置不完整:你設(shè)置了hubbard_occ(2,1)=5(假設(shè)類(lèi)型2是Mn,第一個(gè)manifold是3d,占據(jù)5電子),但如果Mn偽勢(shì)有多個(gè)Hubbard通道(e.g., 3d外還有3p或其它背景manifold),未設(shè)置的可能默認(rèn)為零。或在nscf計(jì)算中,占據(jù)未正確傳播。 偽勢(shì)(PP)問(wèn)題:如果Mn偽勢(shì)(UPF格式)缺少Hubbard manifold標(biāo)簽(e.g., 3d的principal/orbital quantum numbers),代碼無(wú)法正確構(gòu)建manifold,導(dǎo)致零占據(jù)。 占據(jù)計(jì)算模式:occupations='fixed'適合有帶隙系統(tǒng),但對(duì)于金屬雜化材料或能帶計(jì)算,如果費(fèi)米能附近狀態(tài)不準(zhǔn),可能導(dǎo)致有效占據(jù)為零。特に在bands.x中。 其它:tot_magnetization=20表示總磁矩20 μB(可能對(duì)應(yīng)Mn的自旋),但如果沒(méi)有per-atom starting_magnetization,初始猜想可能導(dǎo)致manifold不穩(wěn)定。 |
金蟲(chóng) (正式寫(xiě)手)
木蟲(chóng) (正式寫(xiě)手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 材料328求調(diào)劑 +4 | 一個(gè)蘿卜02 2026-03-03 | 4/200 |
|
|
[考研] 085602化學(xué)工程350,調(diào)劑,有沒(méi)有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
[考研] 0805總分292,求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
[考研] 一志愿華中科技大學(xué),化學(xué)專(zhuān)業(yè)344分,求調(diào)劑 +6 | 邢x(chóng)ing1 2026-03-02 | 6/300 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 278求調(diào)劑 +3 | 滿天星11_22 2026-03-02 | 3/150 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 281電子信息求調(diào)劑 +3 | jhtfeybgj 2026-03-02 | 6/300 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 11408,學(xué)碩276求調(diào)劑 +3 | 崔wj 2026-03-02 | 5/250 |
|
|
[考研] 285求調(diào)劑 +9 | 滿頭大汗的學(xué)生 2026-02-28 | 9/450 |
|
|
[考研] 085600材料工程一志愿中科大總分312求調(diào)劑 +9 | 吃宵夜1 2026-02-28 | 11/550 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-02 | 3/150 |
|
|
[考研] 材料學(xué)調(diào)劑 +10 | 提神豆沙包 2026-02-28 | 12/600 |
|
|
[考研] 調(diào)劑 +3 | 簡(jiǎn)木ChuFront 2026-02-28 | 3/150 |
|
|
[考研] 311求調(diào)劑 +9 | 南迦720 2026-02-28 | 10/500 |
|
|
[論文投稿]
求助coordination chemistry reviews 的寫(xiě)作模板
10+3
|
ljplijiapeng 2026-02-27 | 4/200 |
|
|
[論文投稿]
Optics letters投稿被拒求助
30+3
|
luckyry 2026-02-26 | 4/200 |
|
|
[高分子] 求環(huán)氧樹(shù)脂研發(fā)1名 +3 | 孫xc 2026-02-25 | 11/550 |
|