24小时热门版块排行榜    

查看: 1554  |  回复: 1

lascqy

银虫 (小有名气)

[求助] 杂化泛函计算出能带需要的费米能级应该取哪一个?

在杂化泛函计算能带结构的时候,通常有以下顺序
针对杂化泛函特殊K点的PBE计算——〉
针对杂化泛函特殊K点的HSEscf计算,利用前一步的波函数(ISTART=1)——〉
针对杂化泛函特殊K点的HSEnoscf计算(ISTART=1,ICHARG=11)

经过测试发现,利用第二步或者第三步出的能带结构均可。。没啥太大改变。。
但是。出现与PBE计算不同的是,比如纯半导体的体系。在PBE计算出能带结构后,取PBE静态自洽的或者DOS的费米能级,那么价带顶正好是在能量零点处。。。

而用HSE方法计算能带发现无论取HSEscf还是HSEnoscf步的费米能级都达不到这样。。
后来却发现一个很奇怪的现象,就是利用
普通K点(自动生成)PBE的电荷和波函数,直接进行HSE的非自洽计算
(即ISTART=1,ICHARG=11),在前N步电子步都会出现不收敛的情况
entering main loop
       N       E                     dE             d eps       ncg     rms          ort
         N     E                     dE            rms(c)
   ROT:  1  -0.507082286949E+02   -0.50708E+02
   ROT:  2  -0.506938482460E+02    0.14380E-01
。。。。。。。
   ROT: **  -0.680892037469E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037469E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037469E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
   ROT: **  -0.769335096046E+02   -0.88443E+01
   ROT: **  -0.680892037469E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037469E+02    0.88443E+01
   ROT: **  -0.769335096047E+02   -0.88443E+01
   ROT: **  -0.680892037470E+02    0.88443E+01
gam= 0.000 g(H,U,f)=  0.949E+00-0.182E+01 0.157E-14 ort(H,U,f) = 0.000E+00 0.000E+00 0.000E+00
SDA:   1    -0.237059317507E+02   -0.23706E+02    0.34776E+00   864  -0.869E+00 0.000E+00
。。。。。。。
   ROT:  1  -0.107404641718E+03    0.47400E-01
   ROT:  2  -0.107404640662E+03    0.10556E-05
   ROT:  3  -0.107404641256E+03   -0.59404E-06
gam= 0.382 g(H,U,f)=  0.108E-05 0.878E-05 0.111E-16 ort(H,U,f) = 0.141E-05 0.598E-05 0.827E-16
DMP:  19    -0.243777224012E+02   -0.40814E-05   -0.50766E-05   864   0.987E-05 0.739E-05
   1 F= -.24377722E+02 E0= -.24377722E+02  d E =-.153305E-14  mag=     0.0000
writing wavefunctions
直到最后收敛,这样得到的费米能级却是合理的。。。但是觉得这样的计算不合理啊
HSE应该先用波函数进行自洽的啊。。不知道大家有什么看法?
回复此楼

» 收录本帖的淘帖专辑推荐

杂化泛函方面的研究

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

xueht987

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢交流,没错选2 2012-04-01 09:30:20
fzx2008: 金币+1, 谢谢交流! 2012-04-01 09:31:01
lascqy: 金币+1, 有帮助 2012-04-02 14:14:16
HSE06计算能带貌似不能把ICHARG设为11吧,会出问题的,楼主没遇到吗,在官网论坛也有说明。
2楼2012-04-01 09:27:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lascqy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +46 医学老男孩 2026-08-13 103/5150 2026-08-18 15:33 by dcqxinyang
[基金申请] 明天放榜? +3 Shxjjxjkx 2026-08-18 3/150 2026-08-18 15:32 by a441914426
[基金申请] 重要消息,中午系统在维护 +6 yuleib84 2026-08-18 7/350 2026-08-18 15:02 by sytp
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 8/400 2026-08-18 11:47 by gltch
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 5/250 2026-08-18 11:35 by Tide man
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见