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报错建模问题集

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有10人回复
2025年遐想已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
求个博导看看已经有18人回复

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

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的回帖

xujing15

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 9
帖子: 3
在线: 6.1小时
虫号: 4873211

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

4楼: Originally posted by jghe at 2010-12-04 21:15:41
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
  ...

高手我想问一下你的下半部分是从哪里得来的呢我原来以为是普通的band计算的IBZKPT中的信息但是好像错了求指点

赞一下