| 5 | 1/1 | 返回列表 |
| 查看: 843 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者socchu將贈(zèng)送您 5 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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)先從簡(jiǎn)單調(diào)整開始。如果問(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)原子類型順序或多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)問(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)定。 |
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +4 | 是樂(lè)渝哇 2026-03-03 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
[考研] 環(huán)境調(diào)劑 +8 | chenhanheng 2026-03-02 | 8/400 |
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè),申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
[考研] 299求調(diào)劑 +5 | kkcoco25 2026-03-02 | 9/450 |
|
|
[考研] 26考研報(bào)考西工大材料308分求調(diào)劑 +4 | weizhong123 2026-03-01 | 5/250 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,。 2026-03-02 | 10/500 |
|
|
[考研] 調(diào)劑 +5 | 13853210211 2026-03-02 | 7/350 |
|
|
[考研] 0857調(diào)劑 +6 | 一ll半 2026-02-28 | 8/400 |
|
|
[考研] 材料化工調(diào)劑 +13 | 今夏不夏 2026-03-01 | 16/800 |
|
|
[考研] 085600求調(diào)劑 +4 | LRZZZZZZ 2026-03-02 | 5/250 |
|
|
[考研] 0856材料求調(diào)劑 +12 | hyf hyf hyf 2026-02-28 | 13/650 |
|
|
[考研] 302材料工程求調(diào)劑 +5 | Doleres 2026-03-01 | 6/300 |
|
|
[考研] 306分材料調(diào)劑 +5 | chuanzhu川燭 2026-03-01 | 6/300 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 材料調(diào)劑 +3 | 恒順自然 2026-03-02 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | 熬夜的貓頭鷹 2026-03-02 | 3/150 |
|
|
[考研] 279求調(diào)劑 +3 | dua1 2026-03-01 | 4/200 |
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|