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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】Hse06计算时KPOINTS出错
311/4返回列表1234下一页
查看: 6776  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工 求调剂 +12 enenenhui 2026-03-13 13/650 2026-03-16 08:30 by Linda Hu
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 337一志愿华南理工材料求调剂（有希望2吗？） +3 mysdl 2026-03-09 3/150 2026-03-14 02:53 by JourneyLucky
[考研] 求调剂！ +4 朔朔话 2026-03-09 4/200 2026-03-14 01:38 by JourneyLucky
[考研] 0703化学调剂 +4 快乐的香蕉 2026-03-11 4/200 2026-03-13 22:41 by JourneyLucky
[考研] 材料与化工304求B区调剂 +5 邱gl 2026-03-11 6/300 2026-03-13 22:37 by JourneyLucky
[考研] 0703，333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 泣血叩求调剂恩，愿以丹心报师恩 +6 Iuruoh 2026-03-11 6/300 2026-03-13 22:06 by JourneyLucky
[考研] 求材料调剂 +5 隔壁陈先生 2026-03-12 5/250 2026-03-13 22:03 by 星空星月
[考研] 332求调剂 +3 Zz版 2026-03-13 3/150 2026-03-13 20:36 by 18595523086
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂！ +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 0857 资源与环境 285分 +6 未名考生 2026-03-09 6/300 2026-03-11 21:08 by 30660438
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[考研] 279求调剂 +3 莫xiao 2026-03-10 4/200 2026-03-11 08:06 by 斩魂滴兔子！
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 085602化工求调剂 +7 董boxing 2026-03-10 7/350 2026-03-10 17:07 by BruceLiu320
信息提示
关闭
请填处理意见
关闭