博文
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://blog.sciencenet.cn/blog-671981-634791.html
上一篇：一些有用的博客
下一篇：处理vasp的WAVECAR工具
