| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 850 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題,作者socchu將贈(zèng)送您 5 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
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ě)手)
木蟲(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ù),這不符合代碼要求。基于類似案例分析,以下是可能原因和解決方案建議。你可以逐一嘗試,優(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)原子類型順序或多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)定。 |
木蟲(chóng) (正式寫(xiě)手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 271求調(diào)劑 +7 | 月色c 2026-03-05 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|
[考研] 材料類考研調(diào)劑 +7 | gemmgemm 2026-03-01 | 8/400 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 267調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
[考研] 調(diào)劑材料學(xué)碩 +6 | 詞凝Y 2026-03-02 | 6/300 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研]
材料325求調(diào)劑
30+5
|
mariusuki 2026-03-02 | 10/500 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 0805總分292,求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +11 | February_Feb 2026-03-01 | 11/550 |
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,? 2026-03-02 | 10/500 |
|
|
[考研] 材料工程求調(diào)劑 +3 | 1431251 2026-03-03 | 3/150 |
|
|
[考研] 清華大學(xué) 材料與化工 353分求調(diào)劑 +5 | awaystay 2026-03-02 | 6/300 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
[考研] 化工專碩342,一志愿大連理工大學(xué),求調(diào)劑 +6 | kyf化工 2026-02-28 | 7/350 |
|