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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】band structures using hybrid functionals 的Kpoints问题
191/1返回列表
查看: 5382  |  回复: 18
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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

» 猜你喜欢

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

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

查看全部散金贴
1楼 2010-12-01 21:03:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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

vimtex

银虫 (初入文坛)
show your KPOINTS instead.
赞一下 回复此楼
2楼2010-12-02 00:00:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

282926214

金虫 (小有名气)
引用回帖:
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的回帖

wuli8

荣誉版主 (知名作家)
我算的时候还是有提示K点设置有问题
赞一下 回复此楼
5楼2011-03-11 00:08:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
引用回帖:
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
...

能不能给我个具体的例子任何体系都可以。谢谢
赞一下 回复此楼
6楼2011-03-11 00:08:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jghe

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by wuli8 at 2011-03-11 00:08:47:
能不能给我个具体的例子任何体系都可以。谢谢

什么提示啊，我给你的就是个例子
赞一下(1人) 回复此楼
7楼2011-03-14 02:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
引用回帖:
Originally posted by jghe at 2011-03-14 02:37:00:
什么提示啊，我给你的就是个例子

老是有这样的提示：

| |
| W W AA RRRRR N N II N N GGGG !!! |
| W W A A R R NN N II NN N G G !!! |
| W W A A R R N N N II N N N G !!! |
| W WW W AAAAAA RRRRR N N N II N N N G GGG ! |
| WW WW A A R R N NN II N NN G G |
| W W A A R R N N II N N GGGG !!! |
| |
| Your generating k-point grid is not commensurate to the symmetry |
| of the lattice. This can cause slow convergence with respect |
| to k-points for HF type calculations |
| suggested SOLUTIONS: |
| ) if not already the case, use automatic k-point generation |
| ) shift your grid to Gamma (G) (e.g. required for hex or fcc lattice) |
赞一下 回复此楼
8楼2011-03-14 20:18:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
引用回帖:
Originally posted by jghe at 2011-03-14 02:37:00:
什么提示啊，我给你的就是个例子

我想要个计算好的例子。包含所有的输入和输出文件。
如果没有算个小的也可以。例如：纯diamond的能带结构，用HSE方法。
谢谢
赞一下 回复此楼
9楼2011-03-14 20:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TIGERYZZ

金虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
ellsaking(金币+1): 感谢回帖交流 2011-07-12 15:02:29
关注下，我做GW计算，能带也不知道怎么做，线性撒点用不了。
期待大牛解答。
lz这样加K点的话ISMEAR用-4，-5k点简单加上应该会出错吧，只能用0了。还有就是为什么静态计算的时候不先把需要的K点加上呢，这样的话HF计算时的k点与之前一致，是否有可能消除一些影响。

[ Last edited by TIGERYZZ on 2011-3-14 at 22:02 ]
赞一下(2人) 回复此楼
10楼2011-03-14 21:58:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

215215

银虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by wuli8 at 2011-03-11 00:08:47:
能不能给我个具体的例子任何体系都可以。谢谢

楼主找到例子了吗?我也需要关于这种算法计算的例子,楼主得到了例子计算文件的话,希望楼主给发一份,不胜感激,万分感谢!elijia2010@hotmail.com
此致,
                                                 敬礼!
赞一下(1人) 回复此楼
11楼2011-03-24 15:55:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江湖之远

银虫 (小有名气)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+1): 鼓励交流 2011-04-01 10:26:04
一楼的那个官方回复我有两点不太明白。
1）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,
这个是说把我们想要计算的沿高对称方向的K点自己加到KPOINTS的末尾部分，问题是自己计算这些K点的坐标是一件很麻烦的事情，像是有点能带计算可能需要计算八九十个K点，这要是自己一个一个算坐标，那可费劲了。我们之所以在VASP5.2之前的版本中使用的KPOINTS文件中使用Line-mode模式，就是因为vasp可以根据我们给出的高对称点，自行算出这些K点的坐标。
所以不知道哪位朋友知道有没有比较简单的方法可以得出沿高对称点的K点的坐标？
2) Perform a second VASP run: 在这个方法中自洽预算需要运行两遍，那我想知道第一遍自洽运算的作用是什么？我的想法是：如果我的KPOINTS文件事先就按照他上面所说的程序都设置好了（即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. ），那我是不是直接运行一遍自洽运算就可以得出结果，不需要像一楼所说的那样进行两遍运算，因为毕竟这个运算是相当耗时间的。
有哪位高人可以指点一下？
赞一下(2人) 回复此楼
12楼2011-03-31 21:31:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
★ ★
zzy870720z(金币+2): 谢谢指教 2011-04-01 16:31:44
引用回帖:
Originally posted by 江湖之远 at 2011-03-31 21:31:26:
一楼的那个官方回复我有两点不太明白。
1）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 KPOINT ...

回复你的问题2.
第一次运算是普通的DFT得到一个初始的波函数，第二个SCF是加上HSE的计算。这个HSE运算需要初始的波函数
赞一下(1人) 回复此楼
13楼2011-04-01 15:46:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江湖之远

银虫 (小有名气)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 谢谢指教 2011-04-02 14:14:00
引用回帖:
Originally posted by wuli8 at 2011-04-01 15:46:24:
回复你的问题2.
第一次运算是普通的DFT得到一个初始的波函数，第二个SCF是加上HSE的计算。这个HSE运算需要初始的波函数

在一楼你贴的那个官方回答里边说的很明白
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
第一遍运行就是HSE自洽运算，而非普通的DFT运算。

[ Last edited by 江湖之远 on 2011-4-2 at 09:45 ]
赞一下(2人) 回复此楼
14楼2011-04-02 09:35:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)
继续求助
回复此楼
15楼2011-09-26 15:58:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

robinson8607

铜虫 (初入文坛)
引用回帖:
1882049楼: Originally posted by 江湖之远 at 2011-04-02 09:35:22
在一楼你贴的那个官方回答里边说的很明白
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
第一遍运行就是HSE自洽运算，而非普通的DFT运算。
...

同问
回复此楼
16楼2012-07-25 20:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
1882049楼: Originally posted by 江湖之远 at 2011-04-02 09:35:22
在一楼你贴的那个官方回答里边说的很明白
-) First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file.
第一遍运行就是HSE自洽运算，而非普通的DFT运算。
...

楼主的第一步进行DFT运算的意思是结构优化中的普通DFT运算吧
官方那句话中的First perform a selfconsistent Hartree-Fock/HSE calculation using a conventional KPOINTS file应该也是在结构优化之后进行的静态自洽，官方的回答忽略了之前一步的结构优化步骤。
不知道这样理解对不对
赞一下 回复此楼
17楼2013-03-18 09:37:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gaozhibin

铜虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
8楼: Originally posted by wuli8 at 2011-03-14 20:18:40
老是有这样的提示：

| |
| W W AA RRRRR N N II N N GGGG !!! |
| W W A A R R NN N II NN N G G !!! |
| W W A A R R N N N II N N N G !!! |
| W WW W AAAAAA RRRRR N N N II N N N G GGG ! |
| WW ...

老是有这样的提示：

| |
| W W AA RRRRR N N II N N GGGG !!! |
| W W A A R R NN N II NN N G G !!! |
| W W A A R R N N N II N N N G !!! |
| W WW W AAAAAA RRRRR N N N II N N N G GGG ! |
| WW WW A A R R N NN II N NN G G |
| W W A A R R N N II N N GGGG !!! |
| |
| Your generating k-point grid is not commensurate to the symmetry |
| of the lattice. This can cause slow convergence with respect |
| to k-points for HF type calculations |
| suggested SOLUTIONS: |
| ) if not already the case, use automatic k-point generation |
| ) shift your grid to Gamma (G) (e.g. required for hex or fcc lattice) |



这个提示应该不是问题，只是一个 warning ， 继续算应该没问题吧？  您最后解决了吗？
赞一下 回复此楼
18楼2015-11-22 15:54:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xujing15

新虫 (初入文坛)
★
小木虫: 金币+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中的信息   但是好像错了   求指点
赞一下 回复此楼
19楼2017-04-03 22:15:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli8 的主题更新
191/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭