24小时热门版块排行榜    

查看: 443  |  回复: 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个字符以上)
信息提示
请填处理意见