应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】Hse06计算时KPOINTS出错
北京石油化工学院2026年研究生招生接收调剂公告
311/4返回列表1234下一页
查看: 6858  |  回复: 30
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jttjtl

银虫 (小有名气)

[交流] 【求助】Hse06计算时KPOINTS出错 已有4人参与

本人在使用HSE06计算能带时出现问题：
error in IBZKPT_HF: two k-points are equivalent           
  this will cause problems in the HF routine
我的KPOINTS文件为
K-points
5
Line-mode
rec
0.0 0.0 0.0
0.0 0.0 0.5

0.0 0.0 0.5
0.5 0.5 0.0

0.5 0.5 0.0
0.5 0.5 0.5

0.5 0.5 0.5
0.0 0.0 0.0

错误原因好像是说不能有k点相同，我于是将相同的点做了一下偏移
如第一个点
0 0 0
变为
0 0 0.01
其他重复的点也是偏移了原来位置，但是还是和原来位置在同一直线上，最后算出来了能带结构图。请问这种方法可用吗？
此外我还在论坛上看到了一下做法：
First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.-) Copy the IBZKPT file to KPOINTS, and explicitely add all desired k-points along high-symmetry lines of the BZ that are needed for the bandstructure plot. Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0。
请问高手这样算的原因是什么呢？
回复此楼

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

VASP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-19 15:00:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

y1ding

铁杆木虫 (著名写手)
★ ★
jttjtl(金币+1):谢谢参与
bingmou(金币+1):感谢分享 2010-10-20 13:09:15
jttjtl(金币+5):请问：我计算能带时用了之前自洽得到的WAVECAR+CHGCAR ，这样做对吗？ 2010-10-20 18:59:11
HSE计算能带是不能只靠CHGCAR做非自洽计算的
赞一下(3人) 回复此楼
8楼2010-10-20 12:17:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)
这是我在vasp官方论坛看到的HSE计算的帖子：
band structures using hybrid functionals have to be calculated the following way:
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
-) Copy the IBZKPT file to KPOINTS, and explicitely add all desired k-points along high-symmetry lines of the BZ that are needed for the bandstructure plot. Add the points at the end of the KPOINTS file, but set the weights of these added k-points to 0.
-) Do not forget to set the number of k-points in KPOINTS correctly (to the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines), such that all k-points are used for the calculations.
-) NKRED can not be used.
-) Perform a second VASP run:
It is recommended to use the Davidson algorithm, since it converges that eigen energies at the new k-points fastest. Since VASP terminates when the total energy is converged to a certain threshold, it is important to force VASP to do a minimum number of steps, so that the orbitals at the new k-points are fully converged (note: since their weight is zero, they do not contribute to the total energy).
This can be done using e.g.

ALGO = N ; NELMIN = 5 ! Davidson, minimum 5 scf-steps
IBRION = 1 ! Use simple charge mixer, since
Pulay might blow up

-) The KS-eigenvalues of the states along the high-symmetry lines are written in OUTCAR, EIGENVAL, vasprun.xml; please cut the k-points required for the bandstructure from one of these files and proceed as usual (using p4vasp or any other graphics tool you usually use to produce bandstructure-plots)

想请高手讲讲这样算的原因~
赞一下(1人) 回复此楼
9楼2010-10-20 19:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

大懒猫王浩

金虫 (著名写手)

学渣
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
挖坑来了，帮同学算了一个HSE的能带，大概就是HSE算能带这一步的K点路径取起来比较复杂，可以直接使用vaspkit先生成常规DFT路径，再使用251功能生成HSE的KPOINTS，这样比手动加KPOINTS要快的多。。。
赞一下(1人) 回复此楼
一步一步往上走，，，，
29楼2021-01-22 09:37:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
ym23(金币+2):欢迎常来交流 2010-10-19 20:21:53
jttjtl(金币+5): 2010-10-20 10:08:36
第一个问题。你可以自己画个图看一下，按照你的K点路径有两条线是相交的，分别是（0 0 0.5）到（0.5 0.5 0）这条线与（0.5 0.5 0.5）到（0 0 0）这条线。你取了连线上的5个K点，正好有一个(0.25 0.25 0.25)是相同的，因此会有这个错误。抛开这个问题不谈，貌似你这个路径选取本身就不太正确，一般都取一个大概能走遍所有高对称点的路径，而你选的只是一个面
    第二个问题。由于自洽计算与能带计算取得不同的K点，一般都是做完自洽计算再修改KPOINTS取高对称K点算能带。这段话的意思是，你也可以在做自洽计算的KPOINTS里加入计算能带要用的K点，这样会计算出这些K点的能量本征值(用做画能带)，将这些K点的权重设为0是为了不影响自洽计算(这样新的电荷密度与这些K点对应的电子波函数无关)。这样就达到了在自洽计算的同时完成了能带计算的目的。

[ Last edited by sunyang1988 on 2010-10-19 at 20:29 ]
赞一下(1人) 回复此楼
心之所向无惧无悔
2楼2010-10-19 19:36:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)
不好意思，我写错了，应该是
K-points
5
Line-mode
rec
0.0 0.0 0.0
0.5 0.0 0.0

0.5 0.0 0.0
0.5 0.5 0.0

0.5 0.5 0.0
0.5 0.5 0.5

0.5 0.5 0.5
0.0 0.0 0.0
还是这个错误：
error in IBZKPT_HF: two k-points are equivalent           
  this will cause problems in the HF routine
赞一下 回复此楼
3楼2010-10-20 10:11:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wusongz

金虫 (小有名气)
引用回帖:
Originally posted by sunyang1988 at 2010-10-19 19:36:04:
第一个问题。你可以自己画个图看一下，按照你的K点路径有两条线是相交的，分别是（0 0 0.5）到（0.5 0.5 0）这条线与（0.5 0.5 0.5）到（0 0 0）这条线。你取了连线上的5个K点，正好有一个(0.25 0.25 0.25)是相同 ...

学习了
回复此楼
----Whathappenedtoyou?----Life.
4楼2010-10-20 10:13:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★
jttjtl(金币+2): 2010-10-20 11:22:36
zzy870720z(金币+1): 谢谢提示 2011-09-03 17:00:36
引用回帖:
Originally posted by jttjtl at 2010-10-20 10:11:24:
不好意思，我写错了，应该是
K-points
5
Line-mode
rec
0.0 0.0 0.0
0.5 0.0 0.0

0.5 0.0 0.0
0.5 0.5 0.0

0.5 0.5 0.0
0.5 0.5 0.5

0.5 0.5 0.5
0.0 0.0 0.0
还是这个错误：
error in IBZKP ...

你是不是做完自洽计算然后静态计算能带呢，INCAR修改了没？
赞一下(1人) 回复此楼
心之所向无惧无悔
5楼2010-10-20 10:51:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)
引用回帖:
Originally posted by sunyang1988 at 2010-10-20 10:51:36:

你是不是做完自洽计算然后静态计算能带呢，INCAR修改了没？

我用了静态计算。
这是我计算能带时的INCAR：
ISTART=1
GGA = PE
ENCUT = 400
ISMEAR = 0
SIGMA = 0.1
IBRION = -1
ISIF = 0
NSW = 0
EDIFF = 1E-05
LREAL = .FALSE.
LDIAG=.TRUE.
NELM = 100

LHFCALC = .TRUE.
HFSCREEN = 0.207
ALGO = Damped
TIME = 0.2
AEXX = 0.25
AGGAX = 0.75
AGGAC = 1.0
ALDAC = 1.0
ENCUTFOCK = 0

请问用静态计算有问题吗？
赞一下 回复此楼
6楼2010-10-20 11:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunyang1988

至尊木虫 (著名写手)
★ ★
bingmou(金币+2):感谢你的积极参与，有无兴趣集中给讲讲VASP中HSE03计算的过程呢？谢谢 2010-10-20 13:10:03
引用回帖:
Originally posted by jttjtl at 2010-10-20 11:21:20:


我用了静态计算。
这是我计算能带时的INCAR：
ISTART=1
GGA = PE
ENCUT = 400
ISMEAR = 0
SIGMA = 0.1
IBRION = -1
ISIF = 0
NSW = 0
EDIFF = 1E-05
LREAL = .FALSE.
LDIAG=.TRUE.
NELM = 100 ...

你加个ICHARGE = 11看看吧，另外要确定存在正确的WAVECAR与CHGCAR
赞一下(1人) 回复此楼
心之所向无惧无悔
7楼2010-10-20 12:05:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jttjtl

银虫 (小有名气)
10楼2010-10-21 09:02:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jttjtl 的主题更新
311/4返回列表1234下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703调剂，一志愿天津大学319分 +4 haaaabcd 2026-04-05 4/200 2026-04-05 21:09 by 学员8dgXkO
[考研] 307分材料专业求调剂 +7 Hll胡 2026-04-05 7/350 2026-04-05 18:47 by 无际的草原
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +13 lv945 2026-04-04 13/650 2026-04-05 18:14 by 猪会飞
[考研] 288求调剂 一志愿哈工大 材料与化工 +13 洛神哥哥 2026-04-03 13/650 2026-04-05 17:27 by zzx2138
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +8 哇呼哼呼哼 2026-04-02 9/450 2026-04-05 17:14 by lbsjt
[考研] 313求调剂 +5 海日海日 2026-04-04 7/350 2026-04-05 13:58 by imissbao
[考研] 材料工程专硕求调剂 +11 hyl3153942 2026-03-29 11/550 2026-04-05 11:12 by 风雨无晴
[考研] 材料化工306分找合适调剂 +14 沧海轻舟e 2026-04-04 14/700 2026-04-05 09:53 by 朱云虎202
[考研] 求生物学学硕调剂——364分 +7 云朵遛弯指南 2026-04-04 7/350 2026-04-04 22:49 by zhyzzh
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 297求调剂 +11 ljy20040718！ 2026-04-03 13/650 2026-04-04 09:23 by 来看流星雨10
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 0703一志愿南师大334求调剂 +4 seven7yu 2026-03-30 4/200 2026-04-01 16:10 by oooqiao
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
信息提示
关闭
请填处理意见
关闭