| 查看: 2270 | 回复: 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 ... 不知道是什么 原因。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
rlafite
至尊木虫 (正式写手)
- 应助: 199 (高中生)
- 金币: 11352.5
- 红花: 29
- 沙发: 1
- 帖子: 956
- 在线: 361.9小时
- 虫号: 4563202
- 注册: 2016-04-02
- 性别: GG
- 专业: 同步辐射技术及其应用
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
socchu: 金币+5, ★★★★★最佳答案, 感谢老师解疑，我先按照您提到的建议，再摸索一下，解决不了我再请教您，感谢感谢。 2025-10-14 08:44:52
感谢参与，应助指数 +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
socchu
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1005.7
- 散金: 108
- 红花: 5
- 帖子: 611
- 在线: 237.9小时
- 虫号: 762844
- 注册: 2009-05-04
- 专业: 无机非金属类光电信息与功
3楼2025-10-20 19:33:16
rlafite
至尊木虫 (正式写手)
- 应助: 199 (高中生)
- 金币: 11352.5
- 红花: 29
- 沙发: 1
- 帖子: 956
- 在线: 361.9小时
- 虫号: 4563202
- 注册: 2016-04-02
- 性别: GG
- 专业: 同步辐射技术及其应用
4楼2025-10-25 12:05:00
