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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】Hse06计算时KPOINTS出错
311/4返回列表1234下一页
查看: 7694  |  回复: 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下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 让我中一个面上吧！ +12 大萍1987 2026-08-20 15/750 2026-08-23 20:44 by 大萍1987
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 什么时候开奖？ +9 CrisMessi 2026-08-18 10/500 2026-08-23 12:03 by 丶昵称占用
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 炎黄贵胄 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] filecode，4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗？20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
关闭
请填处理意见
关闭