24小时热门版块排行榜    

查看: 904  |  回复: 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个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +3 大火山小火山 2026-03-16 5/250 2026-03-16 16:54 by barlinike
[考研] 070303一志愿西北大学学硕310找调剂 +5 d如愿上岸 2026-03-12 8/400 2026-03-16 15:19 by peike
[考研] 材料与化工专硕调剂 +3 heming3743 2026-03-16 3/150 2026-03-16 15:05 by peike
[考研] 一志愿华中师范071000,325求调剂 +6 RuitingC 2026-03-12 6/300 2026-03-16 14:50 by 可淡不可忘
[基金申请] NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗 20+5 NSFC2026我来了 2026-03-10 14/700 2026-03-15 23:53 by 不负韶华的虎
[考研] 294求调剂 +3 Zys010410@ 2026-03-13 4/200 2026-03-15 10:59 by zhq0425
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[基金申请] 现在如何回避去年的某一个专家,不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 考研调剂 +4 芬达46 2026-03-12 4/200 2026-03-13 16:04 by ruiyingmiao
[考研] 268求调剂 +4 好运连绵不绝 2026-03-12 4/200 2026-03-13 10:45 by hyswxzs
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子!
[考研] 一志愿:武汉理工,材料工程,英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
信息提示
请填处理意见