版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wuli8

荣誉版主 (知名作家)

[交流] 【求助】band structures using hybrid functionals 的Kpoints问题

http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?4.6633.10
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)

Note: A Hartee-Fock calculation can NOT be continued from an existing CHGCAR file, since the non-local exchange is not determined by the charge density but by the density matrix and/or the KS-orbitals.

以上是官方给出的解释。

我的问题是：
1、-) 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.
按照这样的方法产生的Kpoints运行后提示Kpoints错误，不能读取。
2、这一句不知所云：(to the number of the k-points used in the standard mesh PLUS the number of the k-points along the lines), 没有理解怎么回事？

请问谁用hse的方法计算过能带结构，请show一下下你的Kpoints。
谢谢。

[ Last edited by wuli8 on 2010-12-1 at 21:04 ]

回复此楼

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

VASP	第一性原理相关文档	表界面计算	vasp光学
杂化泛函方面的研究	VASP相关	VASP报错建模问题集

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有26人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有4人回复
磺酰氟产物，毕不了业了！已经有6人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

HF type calculations 这里的HF是Hatree-Fock 还是hybrid fuanctional的意思啊！已经有8人回复
使用多个势函数的问题，关于pair_style hybrid tersoff lj/cut 5.0 已经有5人回复
GW计算的错误提示【问题已解决】已经有4人回复
HSE能带计算中ISYM参数的问题已经有7人回复
求助：p4vasp为什么无法打开vasprun.xml？已经有4人回复
3-phase PMSM Motor using TMS320F240 已经有12人回复
【求助】VASP5.2 optical properties calculations 已经有20人回复
【求助】VASP计算的KPOINTS 设置，求助已经有10人回复
【求助】关于IBZKPT 已经有5人回复
【求助】Hse06计算时KPOINTS出错已经有30人回复
【求助】NBANDS和d band center 已经有6人回复
【求助成功】k-points 已经有4人回复
【求助】PBE0 计算Si、Ge能带结构已经有9人回复
【求助】diamond问题求助已经有11人回复
【转帖】KPOINTS 已经有11人回复
【求助】GW计算错误已经有6人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2010-12-01 21:03:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

282926214

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.1
帖子: 76
在线: 104.5小时
虫号: 886740

引用回帖:

Originally posted by vimtex at 2010-12-02 00:00:21:
show your KPOINTS instead.

IBZKPT
是单机产生的？

赞一下

回复此楼

3楼2010-12-02 08:37:51

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 19 个回答

vimtex

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 443.2
帖子: 34
在线: 17.5小时
虫号: 481120

show your KPOINTS instead.

赞一下

回复此楼

2楼2010-12-02 00:00:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jghe

金虫 (著名写手)

★ ★
zzy870720z(金币+2):谢谢指教 2010-12-04 21:36:29
wuli8(金币+100):感谢。我再试一试 2010-12-04 22:51:40

Automatically generated mesh
   57
Reciprocal lattice
  0.000000  0.000000  0.000000    1.000000
  0.333333  0.000000  0.000000    1.000000
-0.333333  0.000000  0.000000    1.000000
  0.000000  0.333333  0.000000    1.000000
  0.333333  0.333333  0.000000    1.000000
-0.333333  0.333333  0.000000    1.000000
  0.000000 -0.333333  0.000000    1.000000
  0.333333 -0.333333  0.000000    1.000000
-0.333333 -0.333333  0.000000    1.000000
  0.000000  0.000000  0.333333    1.000000
  0.333333  0.000000  0.333333    1.000000
-0.333333  0.000000  0.333333    1.000000
  0.000000  0.333333  0.333333    1.000000
  0.333333  0.333333  0.333333    1.000000
-0.333333  0.333333  0.333333    1.000000
  0.000000 -0.333333  0.333333    1.000000
  0.333333 -0.333333  0.333333    1.000000
-0.333333 -0.333333  0.333333    1.000000
  0.000000  0.000000 -0.333333    1.000000
  0.333333  0.000000 -0.333333    1.000000
-0.333333  0.000000 -0.333333    1.000000
  0.000000  0.333333 -0.333333    1.000000
  0.333333  0.333333 -0.333333    1.000000
-0.333333  0.333333 -0.333333    1.000000
  0.000000 -0.333333 -0.333333    1.000000
  0.333333 -0.333333 -0.333333    1.000000
-0.333333 -0.333333 -0.333333    1.000000
  0.00000000  0.00000000  0.50000000    0
  0.00000000  0.00000000  0.44444444    0
  0.00000000  0.00000000  0.38888889    0
  0.00000000  0.00000000  0.33333333    0
  0.00000000  0.00000000  0.27777778    0
  0.00000000  0.00000000  0.22222222    0
  0.00000000  0.00000000  0.16666667    0
  0.00000000  0.00000000  0.11111111    0
  0.00000000  0.00000000  0.05555556    0
  0.00000000  0.00000000  0.00000000    0
  0.00000000  0.00000000  0.00000000    0
  0.05555556  0.05555556  0.00000000    0
  0.11111111  0.11111111  0.00000000    0
  0.16666667  0.16666667  0.00000000    0
  0.22222222  0.22222222  0.00000000    0
  0.27777778  0.27777778  0.00000000    0
  0.33333333  0.33333333  0.00000000    0
  0.38888889  0.38888889  0.00000000    0
  0.44444444  0.44444444  0.00000000    0
  0.50000000  0.50000000  0.00000000    0
  0.50000000  0.50000000  0.00000000    0
  0.50000000  0.50000000  0.05555556    0
  0.50000000  0.50000000  0.11111111    0
  0.50000000  0.50000000  0.16666667    0
  0.50000000  0.50000000  0.22222222    0
  0.50000000  0.50000000  0.27777778    0
  0.50000000  0.50000000  0.33333333    0
  0.50000000  0.50000000  0.38888889    0
  0.50000000  0.50000000  0.44444444    0
  0.50000000  0.50000000  0.50000000    0

例子，上半部分是Monkhorst Pack产生的K点，下面是高对称的K点，第二行总K点的数目是两部分的总和，画能带的时候只取下面高对称部分的本征值即可

赞一下(2人)

回复此楼

4楼2010-12-04 21:15:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖