| 4 | 1/1 | 返回列表 |
| 查看: 842 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題,作者socchu將贈(zèng)送您 5 個(gè)金幣 | ||
socchu金蟲 (正式寫手)
|
[求助]
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,開啟了自旋,設(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 ... 不知道是什么 原因。 |
木蟲 (正式寫手)
|
這個(gè)錯(cuò)誤通常出現(xiàn)在Quantum ESPRESSO的DFT+U計(jì)算中,尤其是后處理步驟如bands.x時(shí),pw.x的SCF和nscf/bands計(jì)算正常,但后處理檢測(cè)到Hubbard manifold(這里是Mn-3d軌道)有零占據(jù),這不符合代碼要求;陬愃瓢咐治,以下是可能原因和解決方案建議。你可以逐一嘗試,優(yōu)先從簡單調(diào)整開始。如果問題持續(xù),提供完整輸入文件(包括ATOMIC_SPECIES、HUBBARD卡和偽勢(shì)細(xì)節(jié))有助于進(jìn)一步診斷。 可能原因 軟件bug或兼容性問題:在QE 7.1及早期版本中,后處理工具(如bands.x、dos.x、projwfc.x)對(duì)Hubbard處理的bug會(huì)導(dǎo)致零占據(jù)錯(cuò)誤,尤其當(dāng)原子類型順序或多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è)類型2是Mn,第一個(gè)manifold是3d,占據(jù)5電子),但如果Mn偽勢(shì)有多個(gè)Hubbard通道(e.g., 3d外還有3p或其它背景manifold),未設(shè)置的可能默認(rèn)為零;蛟趎scf計(jì)算中,占據(jù)未正確傳播。 偽勢(shì)(PP)問題:如果Mn偽勢(shì)(UPF格式)缺少Hubbard manifold標(biāo)簽(e.g., 3d的principal/orbital quantum numbers),代碼無法正確構(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的自旋),但如果沒有per-atom starting_magnetization,初始猜想可能導(dǎo)致manifold不穩(wěn)定。 |
金蟲 (正式寫手)
木蟲 (正式寫手)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266求調(diào)劑 +5 | 哇塞王帥 2026-03-03 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +4 | 青科11 2026-03-02 | 4/200 |
|
|
[考研] 291求調(diào)劑 +3 | MuoLuo1312 2026-03-02 | 6/300 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 材料類考研調(diào)劑 +6 | gemmgemm 2026-03-01 | 7/350 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考研]
|
好好好1233 2026-02-28 | 16/800 |
|
|
[考研] 290分材料工程085601求調(diào)劑 數(shù)二英一 +8 | llx0610 2026-03-02 | 9/450 |
|
|
[考研] 272求調(diào)劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
[考研] 材料284求調(diào)劑,一志愿鄭州大學(xué)英一數(shù)二專碩 +15 | 想上岸的土撥鼠 2026-02-28 | 15/750 |
|
|
[考研] 材料復(fù)試調(diào)劑 +5 | 學(xué)材料的點(diǎn) 2026-03-01 | 6/300 |
|
|
[考研] 化學(xué),材料,環(huán)境類求調(diào)劑 +7 | 考研版棒棒 2026-03-02 | 7/350 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 322求調(diào)劑 +3 | 熊境喆 2026-03-01 | 3/150 |
|
|
[考研] 272求調(diào)劑 +6 | 田智友 2026-02-28 | 6/300 |
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|