24小时热门版块排行榜    

查看: 4264  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zyqh

铁杆木虫 (著名写手)

[求助] HSE06计算能带求组 已有1人参与

各位大虾: 我最近用VASP计算材料的能带,因为直接用LDA 和 GGA 方法忽略了强关联作用,造成了能带带隙偏小,所以我采用HSE泛函杂化的方法。

根据VASP手册和网上虫子的资料,我的计算步骤是一 ,正常的优化,二,DFT的自洽(SCF), 三,HSE06的自洽(用到上一步DFT自洽的WAVECAR),四,做HSE能带计算。
我做到第三步的时候,也就是HSE06自洽,提交任务后,显示了如下的提示:
running on   72 nodes
distr:  one band on    1 nodes,   72 groups
vasp.5.2.2 15Apr09 complex
POSCAR found type information on POSCAR  Br H  Ge C  N
POSCAR found :  5 types and   48 ions

-----------------------------------------------------------------------------
|                                                                             |
|  ADVICE TO THIS USER RUNNING 'VASP/VAMP'   (HEAR YOUR MASTER'S VOICE ...):  |
|                                                                             |
|      You have a (more or less) 'large supercell' and for larger cells       |
|      it might be more efficient to use real space projection operators      |
|      So try LREAL= Auto  in the INCAR   file.                               |
|      Mind: If you want to do a very accurate calculations keep the          |
|      reciprocal projection scheme          (i.e. LREAL=.FALSE.)             |
|                                                                             |
-----------------------------------------------------------------------------

LDA part: xc-table for Pade appr. of Perdew
POSCAR found type information on POSCAR  Br H  Ge C  N
POSCAR found :  5 types and   48 ions
found WAVECAR, reading the header
  number of bands has changed, file:     0 present:   144
  trying to continue reading WAVECAR, but it might fail
  number of k-points has changed, file:     0 present:    32
  trying to continue reading WAVECAR, but it might fail
WAVECAR: different cutoff or change in lattice found
POSCAR, INCAR and KPOINTS ok, starting setup
WARNING: small aliasing (wrap around) errors must be expected
FFT: planning ...(           1 )
reading WAVECAR
random initialization beyond band            1
the WAVECAR file was read sucessfully
initial charge from wavefunction

请问计算过程有错误吗?为什么说K 点和带数都改变了,当然能够一直算下去,计算不会中断,我就想知道K点改变的提示,带数改变的提示,会不会对最后自洽的结果有影响,进而影响能带的结果。 谢谢各位大侠了,我是第一次计算,很困惑,请懂的给小弟略为点拨一二。金币奉上
回复此楼

» 猜你喜欢

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

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

雨后方能见彩虹
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyqh

铁杆木虫 (著名写手)

引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-30 09:19:57
"因为直接用LDA 和 GGA 方法忽略了强关联作用,造成了能带带隙偏小,所以我采用HSE泛函杂化的方法。"

纠正你的理解:
HSE杂化泛函只是改进了交换作用项。对关联项并没有任何的改进,还是LDA/GGA的。
...

谢谢候老师的回答,懂了许些!

那么,候老师,我按照现有方法计算能带应该是可行的吧
雨后方能见彩虹
3楼2014-11-30 13:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
zyqh: 金币+5 2014-11-30 13:21:02
"因为直接用LDA 和 GGA 方法忽略了强关联作用,造成了能带带隙偏小,所以我采用HSE泛函杂化的方法。"

纠正你的理解:
HSE杂化泛函只是改进了交换作用项。对关联项并没有任何的改进,还是LDA/GGA的。

  至于另一个问题:计算过程中能带数目的变化,见我下面一个帖子中的应答。
http://muchong.com/bbs/viewthread.php?tid=8228013&fpage=1

另外你也可以读读fileio.F的代码理解为什么会提示能带数目和k点数目的改变。
WAVECAR里所包含的能带数目和k点数目信息与当前计算设置里的不一致。所以提示它们的值变化了。

”我就想知道K点改变的提示,带数改变的提示,会不会对最后自洽的结果有影响,进而影响能带的结果。“
不会影响计算结果。

» 本帖已获得的红花(最新10朵)

2楼2014-11-30 09:19:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

对马拉二胡

至尊木虫 (职业作家)

楼主,我想请教一下用HSE06计算能带是不是需要新的势函数才行(而不是LDA或者GGA-PBE)?
相互学习,共同进步。
4楼2015-07-01 14:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangxiaoli-

木虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by valenhou001 at 2014-11-30 09:19:57
"因为直接用LDA 和 GGA 方法忽略了强关联作用,造成了能带带隙偏小,所以我采用HSE泛函杂化的方法。"

纠正你的理解:
HSE杂化泛函只是改进了交换作用项。对关联项并没有任何的改进,还是LDA/GGA的。
...

“HSE杂化泛函只是改进了交换作用项。对关联项并没有任何的改进,还是LDA/GGA的。“那我可不可以交换项用HSE,关联相用GGA+U呢?如果可以该怎么设置参数呢?把U值得参数和杂化泛函相关参数都写进INCAR吗?
希望越来越好
5楼2017-01-19 02:18:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[基金申请] 申请科研项目纯靠实力能成功吗? 200+4 David1982 2026-10-13 7/350 2026-10-14 14:24 by 13108017953
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
信息提示
请填处理意见