博文
HSE计算能带
||
http://emuch.net/bbs/viewthread.php?tid=4016318首先给出以前总结的HSE计算的链接，这段时间自己用HSE方法计算了能带，72个原子（掺杂体系），42个K点，现在给出自己的一点小经验，希望对大家有帮助。
1.hse band calculation之前不需要进行hse scf calculation,只需要进行pbe scf calculation,并且注意pbe scf calculation时就需要将KPOINTS设置为IBZKPT+inp.kpt,这样hse band calculation时读取的上一步的波函数才是正确的。
2.hse band calculation时，ISTART=1，ICHARG不能取11而是默认值。因为hse时哈密顿量变了，也就是CHARG变了，不能从上一步pbe自洽计算读。
3.能带计算完成后数据处理时，遇到的问题是能级出现交叉，如图中的红线标记，而实际上是不应该交叉的，这里图中交叉是由于显示的问题。本人解决方式是：在origin里每一列Sort Column--Ascending，这样重新排序之后再做图即可。注意：线性产生的点也就是权重不为零的点不能用来做图，只用那些权重为0的高对称点来画能带图。
http://emuch.net/html/201203/4232787.html
3.band计算
INCAR
SYSTEM = Si
ISTART = 1
ICHARG = 2
GGA = PS
EDIFF = 0.00001
EDIFFG = -0.001
ENCUT = 400
ENAUG = 800
LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
NELM = 200
NSW = 0
IBRION = 1
LMAXMIX = 4
ISMEAR = 0
SIGMA = 0.1
NSIM = 4
ISYM = 0
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = N
TIME = 0.4
ENCUTFOCK = 0
AEXX = 0.25
NELMIN = 5
KPOINTS
从HSE06文件夹的IBZKPT拷贝到KPOINTS,然后把高对称点拷贝在后面，权重为0
更新k点总数
Automatically generated mesh
127（36+91后面91为高对称）Reciprocal lattice
0.00000000000000 0.00000000000000 0.00000000000000 1
0.25000000000000 0.00000000000000 0.00000000000000 2
0.50000000000000 0.00000000000000 0.00000000000000 1
0.00000000000000 0.25000000000000 0.00000000000000 2
0.25000000000000 0.25000000000000 0.00000000000000 2
0.50000000000000 0.25000000000000 0.00000000000000 2
-0.25000000000000 0.25000000000000 0.00000000000000 2
0.00000000000000 0.50000000000000 0.00000000000000 1
0.25000000000000 0.50000000000000 0.00000000000000 2
0.50000000000000 0.50000000000000 0.00000000000000 1
0.00000000000000 0.00000000000000 0.25000000000000 2
0.25000000000000 0.00000000000000 0.25000000000000 2
0.50000000000000 0.00000000000000 0.25000000000000 2
-0.25000000000000 0.00000000000000 0.25000000000000 2
0.00000000000000 0.25000000000000 0.25000000000000 2
0.25000000000000 0.25000000000000 0.25000000000000 2
0.50000000000000 0.25000000000000 0.25000000000000 2
-0.25000000000000 0.25000000000000 0.25000000000000 2
0.00000000000000 0.50000000000000 0.25000000000000 2
0.25000000000000 0.50000000000000 0.25000000000000 2
0.50000000000000 0.50000000000000 0.25000000000000 2
-0.25000000000000 0.50000000000000 0.25000000000000 2
0.00000000000000 -0.25000000000000 0.25000000000000 2
0.25000000000000 -0.25000000000000 0.25000000000000 2
0.50000000000000 -0.25000000000000 0.25000000000000 2
-0.25000000000000 -0.25000000000000 0.25000000000000 2
0.00000000000000 0.00000000000000 0.50000000000000 1
0.25000000000000 0.00000000000000 0.50000000000000 2
0.50000000000000 0.00000000000000 0.50000000000000 1
0.00000000000000 0.25000000000000 0.50000000000000 2
0.25000000000000 0.25000000000000 0.50000000000000 2
0.50000000000000 0.25000000000000 0.50000000000000 2
-0.25000000000000 0.25000000000000 0.50000000000000 2
0.00000000000000 0.50000000000000 0.50000000000000 1
0.25000000000000 0.50000000000000 0.50000000000000 2
0.50000000000000 0.50000000000000 0.50000000000000 1
0.500000 0.000000 0.500000 0.00
0.475000 0.000000 0.475000 0.00
0.450000 0.000000 0.450000 0.00
。。。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。。。
0.093750 0.093750 0.187500 0.00
0.075000 0.075000 0.150000 0.00
0.056250 0.056250 0.112500 0.00
0.037500 0.037500 0.075000 0.00
0.018750 0.018750 0.037500 0.00
0.000000 0.000000 0.000000 0.00
其他文件从HSE06中拷贝
运行vasp
4.从EIGENVAL中取权重为0的高对称点做能带图
EIGENVAL：
2 2 1 1
0.2010146E+02 0.3845247E-09 0.3845247E-09 0.3845247E-09 0.5000000E-15
1.000000000000000E-004
CAR
Si-Diamond
8 127 8
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.1562500E-01
1 -7.830170
2 5.377682
3 5.377684
4 5.377686
5 8.798907
6 8.798909
7 8.798911
8 9.693314
0.2500000E+00 0.0000000E+00 0.0000000E+00 0.3125000E-01
1 -6.949546
2 1.136336
3 4.578605
4 4.578607
5 8.218893
6 9.797289
7 9.797292
8 13.389561
。。。。。。。。。。。。。。。
。。。。。。。。。。。。。。。
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 -7.830170
2 5.377682
3 5.377684
4 5.377686
5 8.798907
6 8.798909
7 8.798911
8 9.693314
去掉前面36个
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.1562500E-01
1 -7.830170
2 5.377682
3 5.377684
4 5.377686
5 8.798907
6 8.798909
7 8.798911
8 9.693314
留下后面91个高对称点，并改变高对称点数
2 2 1 1
0.2010146E+02 0.3845247E-09 0.3845247E-09 0.3845247E-09 0.5000000E-15
1.000000000000000E-004
CAR
Si-Diamond
8 91 8
。。。。。。。。。。。。
。。。。。。。。。。。。
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 -7.830170
2 5.377682
3 5.377684
4 5.377686
5 8.798907
6 8.798909
7 8.798911
8 9.693314
然后作图
https://wap.sciencenet.cn/blog-671981-634791.html
上一篇：一些有用的博客
下一篇：处理vasp的WAVECAR工具
全部精选博文导读
相关博文
- • Minerals线下恳谈会：履践致远、与时偕行——对话中国科学院广州地球化学研究所期刊合作学者
- • 聚英才 建高地 | 北京理工大学“特立青年学者”全球招聘开启
- • 700年后日本或濒临灭绝？日本学者推算预测：届时或仅剩1名15岁以下孩子
- • [转载]【同位素视角】非英语母语学者如何区分’e.g.’, ‘i.e.’, ‘namely’与‘such as’等混淆难题
- • 美国佐治亚大学等机构学者：刈割策略对Bulldog 805紫花苜蓿+Tifton 85狗牙根混播草地产量及品质的影响
- • 美国堪萨斯州立大学、密苏里大学等机构学者研究成果：土壤水分管理策略和品种多样性对紫花苜蓿产量、营养品质和农场盈利能力的影