| 4 | 1/1 | 返回列表 |
| 查看: 868 | 回復(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)為零;蛟趎scf計(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ā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工,291,求調(diào)劑 +12 | 咕嚕咕嚕123123 2026-03-05 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑,307分 +10 | 張泳銘1 2026-03-09 | 10/500 |
|
|
[考研] 一志愿清華深研院材料專(zhuān)碩294分,專(zhuān)業(yè)課111分,本科中南大學(xué)材料,有六級(jí),有工作經(jīng)驗(yàn) +4 | H14528 2026-03-04 | 4/200 |
|
|
[考研] 一志愿鄭州大學(xué),學(xué)碩,物理化學(xué), 333求調(diào)劑 +8 | 李魔女斗篷 2026-03-04 | 9/450 |
|
|
[考研] 313分生物學(xué)求調(diào)劑 +4 | Yyt楊1 2026-03-09 | 6/300 |
|
|
[考研] 【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工,服從調(diào)劑,抗壓能力強(qiáng)! +13 | xiiiia 2026-03-04 | 14/700 |
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
[考研] 085701環(huán)境工程專(zhuān)碩求調(diào)劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
[考研] 一志愿中科大能動(dòng)297求調(diào)劑,本科川大 +4 | 邵11 2026-03-03 | 4/200 |
|
|
[考研] 301求調(diào)劑 +5 | 一二LV 2026-03-07 | 5/250 |
|
|
[考研] 08工科 +4 | li李樂(lè)成 2026-03-06 | 4/200 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 276求調(diào)劑 本科天津工業(yè)大學(xué) 一志愿北京化工大學(xué)0855 +4 | 路lyh123 2026-03-06 | 6/300 |
|
|
[考研] 085600材料與化工(高分子)290分求調(diào)劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
[考研] 0856材料與化工求調(diào)劑! +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|