24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1718  |  回复: 2
【悬赏金币】回答本帖问题,作者fo3mt2将赠送您 10 个金币

fo3mt2

金虫 (小有名气)

[求助] HSE06杂化泛函计算能带问题已有1人参与

大家好,

向大家请教一个问题,我在使用VASP代码里的HSE06杂化泛函计算体系的能带结构,我把自洽计算步骤产生的IBZKPT文件拷贝到了非自洽计算步骤里,并改名为KPOINTS文件,随后我又把能带的k路径坐标粘贴在KPOINTS文件的后面。(这个能带的k路径坐标是我之前利用VASP的PBE泛函计算相同体系的能带时,在OUTCAR文件里产生的k点坐标)。新的KPOINTS文件内容如下。

Automatically generated mesh
61
Reciprocal lattice
0.00000000000000 0.00000000000000 0.00000000000000 1
0.09090909090909 0.00000000000000 0.00000000000000 2
...
-0.09090909090909 0.45454545454546 0.00000000000000 2
0.00000000 0.00000000 0.00000000 0.000
...
0.50000000 0.50000000 0.00000000 0.000

非自洽计算步骤里的INCAR文件内容如下。

Syetem = C_based
ISTART = 1
LWAVE = .FALSE.
LCHARGE = .TRUE.
ICHARG = 11
IVDW = 12
ISPIN = 2
MAGMOM = 0.156 -0.155 -0.096 0.096 -0.156 0.155 -0.156 0.155 0.096\
-0.096 0.156 -0.155 -0.156 0.155 0.096 -0.096 0.156 -0.155\
0.156 -0.155 -0.096 0.096 -0.156 0.155
LORBIT = 11
ENCUT = 500
ISMEAR = 0
SIGMA = 0.03
PREC = Accurate
LHFCALC = .TRUE.
AEXX = 0.25
HFSCREEN = 0.2
ALGO = Damped
TIME = 0.4
PRECFOCK = Fast
ISYM = 3
LREAL = Auto
NELMIN = 5
NELM = 700
NELMDL = -10
EDIFF = 1e-5
NPAR = 32
KPAR= 16

计算结束后,我利用p4v软件读取vasprun.xml文件,但画出的能带完全错误。而且输出文件也提示警告信息,说k点与晶格的对称性不匹配。警告信息如下。

-----------------------------------------------------------------------------
| |
| W W AA RRRRR N N II N N GGGG !!! |
| W W A A R R NN N II NN N G G !!! |
| W W A A R R N N N II N N N G !!! |
| W WW W AAAAAA RRRRR N N N II N N N G GGG ! |
| WW WW A A R R N NN II N NN G G |
| W W A A R R N N II N N GGGG !!! |
| |
| Your generating k-point grid is not commensurate to the symmetry |
| of the lattice. This can cause slow convergence with respect |
| to k-points for HF type calculations |
| suggested SOLUTIONS: |
| ) if not already the case, use automatic k-point generation |
| ) shift your grid to Gamma (G) (e.g. required for hex or fcc lattice) |
| |
-----------------------------------------------------------------------------

我不知道是哪里出了问题。能麻烦大家给些建议,解决这个问题吗?

谢谢大家啦。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenzikun

铜虫 (著名写手)

K点文件的第二行,数字对应的是全部的K点数目吗?

发自小木虫Android客户端
低头要勇气,抬头要智慧
2楼2021-10-08 10:45:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

轻轻飘过之

金虫 (正式写手)

【答案】应助回帖

首先,如楼上所说,计算的时候K点文件的第二行,数字对应的是全部的K点数目;
其次,在画能带图时,需要把EIGENVAL中对应IBZKPT的K点下面的本征值删除,同时修改EIGENVAL中相应K点数目的总和
3楼2021-10-08 11:10:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fo3mt2 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见