24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1046  |  回复: 3
【悬赏金币】回答本帖问题,作者socchu将赠送您 5 个金币

socchu

金虫 (正式写手)

[求助] Hubbard manifold with zero occupations is not allowed 已有1人参与

各位老师好:
在用qe计算能带的时候报错了,请老师们帮忙看看。
工作时在进行一个有机金属mn的杂化材料的能带计算,由于有mn,开启了自旋,设置了u值,自旋和u值参数如下:
/
&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
出现的状况是,进行scf和bands的时候(pw.x)的时候是合适的,但是在进行能带解析的时候(bands.x)的时候报错了,错误如下:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     error in routine offset_atom_wfc (1):
     hubbard manifold with zero occupations is not allowed
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
不知道是什么 原因。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
socchu: 金币+5, ★★★★★最佳答案, 感谢老师解疑,我先按照您提到的建议,再摸索一下,解决不了我再请教您,感谢感谢。 2025-10-14 08:44:52
这个错误通常出现在Quantum ESPRESSO的DFT+U计算中,尤其是后处理步骤如bands.x时,pw.x的SCF和nscf/bands计算正常,但后处理检测到Hubbard manifold(这里是Mn-3d轨道)有零占据,这不符合代码要求。基于类似案例分析,以下是可能原因和解决方案建议。你可以逐一尝试,优先从简单调整开始。如果问题持续,提供完整输入文件(包括ATOMIC_SPECIES、HUBBARD卡和伪势细节)有助于进一步诊断。
可能原因
软件bug或兼容性问题:在QE 7.1及早期版本中,后处理工具(如bands.x、dos.x、projwfc.x)对Hubbard处理的bug会导致零占据错误,尤其当原子类型顺序或多manifold时。你的设置(nspin=2, Hubbard ortho-atomic, U Mn-3d=2.5)是标准的collinear DFT+U,但后处理可能未正确继承占据。
Hubbard_occ设置不完整:你设置了hubbard_occ(2,1)=5(假设类型2是Mn,第一个manifold是3d,占据5电子),但如果Mn伪势有多个Hubbard通道(e.g., 3d外还有3p或其它背景manifold),未设置的可能默认为零。或在nscf计算中,占据未正确传播。
伪势(PP)问题:如果Mn伪势(UPF格式)缺少Hubbard manifold标签(e.g., 3d的principal/orbital quantum numbers),代码无法正确构建manifold,导致零占据。
占据计算模式:occupations='fixed'适合有带隙系统,但对于金属杂化材料或能带计算,如果费米能附近状态不准,可能导致有效占据为零。特に在bands.x中。
其它:tot_magnetization=20表示总磁矩20 μB(可能对应Mn的自旋),但如果没有per-atom starting_magnetization,初始猜想可能导致manifold不稳定。
2楼2025-10-13 23:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

socchu

金虫 (正式写手)

引用回帖:
2楼: Originally posted by rlafite at 2025-10-13 23:05:08
这个错误通常出现在Quantum ESPRESSO的DFT+U计算中,尤其是后处理步骤如bands.x时,pw.x的SCF和nscf/bands计算正常,但后处理检测到Hubbard manifold(这里是Mn-3d轨道)有零占据,这不符合代码要求。基于类似案例分 ...

感谢前面老师对Hubbard manifold with zero occupations is not allowed的解答,目前已经解决了,但是现在出现了新的错误:Error in routine write_p_avg (1):     write_p_avg not working with LDA+U,麻烦老师帮忙看看
3楼2025-10-20 19:33:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rlafite

木虫 (正式写手)

引用回帖:
3楼: Originally posted by socchu at 2025-10-19 23:33:16
感谢前面老师对Hubbard manifold with zero occupations is not allowed的解答,目前已经解决了,但是现在出现了新的错误:Error in routine write_p_avg (1):     write_p_avg not working with LDA+U,麻烦老师帮 ...

错误根本原因:write_p_avg 目前不支持与 LDA+U 同时使用。

解决方案:去掉 write_p_avg = .true.,或暂时禁用 LDA+U。

推荐做法:先完成 LDA+U 自洽计算,再在相同结构上执行非-U 计算输出平均势。
4楼2025-10-25 12:05:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 socchu 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 机械专硕274求调剂,不挑专业学校 +6 泛泛2333 2026-04-05 7/350 2026-04-05 23:20 by chyhaha
[考研] 327求调剂 +4 拾光任染 2026-04-05 4/200 2026-04-05 20:16 by 南航~万老师
[考研] 求调剂到0856材料工程 +3 程9915 2026-04-05 3/150 2026-04-05 18:15 by 蓝云思雨
[考研] 调剂 +3 好好读书。 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考博] 申博 +7 IQwQl 2026-04-04 7/350 2026-04-04 23:32 by mumin1990
[考研] 材料调剂 +11 一样YWY 2026-04-02 13/650 2026-04-04 23:10 by 无际的草原
[考研] 求调剂:085600材料与化工,考材科基,总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 一志愿沪985,326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
[考研] 283分材料与化工求调剂 +29 罗KAKA 2026-04-02 29/1450 2026-04-03 23:56 by userper
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 化学070300-总分378-求调剂 +5 挪椅子的泡泡糖 2026-04-02 5/250 2026-04-02 22:20 by ZXlzxl0425
[考研] 296求调剂 +4 sdhu 2026-04-02 4/200 2026-04-02 21:29 by baoball
[考研] 一志愿山东大学,085600,344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 一志愿北京科技材料科学与工程288分,求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 土木304求调剂 +5 顶级擦擦 2026-03-31 5/250 2026-04-01 08:15 by fdcxdystjk¥
信息提示
请填处理意见