| 查看: 6407 | 回复: 30 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】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。 请问高手这样算的原因是什么呢？ |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
may-sucess
金虫 (正式写手)
- 应助: 21 (小学生)
- 金币: 1443.5
- 散金: 1043
- 红花: 10
- 帖子: 397
- 在线: 523.5小时
- 虫号: 621473
- 注册: 2008-10-09
- 专业: 凝聚态物性 II ：电子结构
★
jttjtl(金币+1):谢谢参与
jttjtl(金币+1):谢谢参与
|
我也是出现这样的问题，INCAR设置和你一样，KPOINTS是按照VASP说明书上的，发现没有相同的k点。 给出如下错误： LDA part: xc-table for Pade appr. of Perdew 9 error in IBZKPT_HF: two k-points are equivalent 11 73 K点设置如下： Automatically generated mesh 2 62 3 Reciprocal lattice 4 0.00000000000000 0.00000000000000 0.00000000000000 1 5 0.00000000000000 0.00000000000000 0.16666666666667 1 6 0.00000000000000 0.00000000000000 0.33333333333333 1 7 0.00000000000000 0.00000000000000 0.50000000000000 1 8 0.0000000E+00 0.0000000E+00 0.5555556E-01 1 9 0.0000000E+00 0.0000000E+00 0.1111111E+00 1 10 0.0000000E+00 0.0000000E+00 0.2222222E+00 1 11 0.0000000E+00 0.0000000E+00 0.2777778E+00 1 12 0.0000000E+00 0.0000000E+00 0.3888889E+00 1 13 0.0000000E+00 0.0000000E+00 0.4444444E+00 1 14 -0.3703333E-01 0.7407444E-01 0.5000000E+00 1 15 -0.7406667E-01 0.1481489E+00 0.5000000E+00 1 16 -0.1111000E+00 0.2222233E+00 0.5000000E+00 1 17 -0.1481333E+00 0.2962978E+00 0.5000000E+00 1 18 -0.1851667E+00 0.3703722E+00 0.5000000E+00 1 19 -0.2222000E+00 0.4444467E+00 0.5000000E+0 1 20 -0.2592333E+00 0.5185211E+00 0.5000000E+00 1 21 -0.2962667E+00 0.5925956E+00 0.5000000E+00 1 22 -0.3333000E+00 0.6666700E+00 0.5000000E+00 1 23 -0.3333000E+00 0.6666700E+00 0.4444444E+00 1 24 -0.3333000E+00 0.6666700E+00 0.3888889E+00 1 25 -0.3333000E+00 0.6666700E+00 0.3333333E+00 1 26 -0.3333000E+00 0.6666700E+00 0.2777778E+00 1 27 -0.3333000E+00 0.6666700E+00 0.2222222E+00 1 28 -0.3333000E+00 0.6666700E+00 0.1666667E+00 1 29 -0.3333000E+00 0.6666700E+00 0.1111111E+00 1 30 -0.3333000E+00 0.6666700E+00 0.5555556E-01 1 31 -0.3333000E+00 0.6666700E+00 0.0000000E+00 1 -0.2962667E+00 0.5925956E+00 0.0000000E+00 1 33 -0.2592333E+00 0.5185211E+00 0.0000000E+00 1 34 -0.2222000E+00 0.4444467E+00 0.0000000E+00 1 35 -0.1851667E+00 0.3703722E+00 0.0000000E+00 1 36 -0.1481333E+00 0.2962978E+00 0.0000000E+00 1 37 -0.1111000E+00 0.2222233E+00 0.0000000E+00 1 38 -0.7406667E-01 0.1481489E+00 0.0000000E+00 1 39 -0.3703333E-01 0.7407444E-01 0.0000000E+00 1 40 0.0000000E+00 0.5555556E-01 0.0000000E+00 1 41 0.0000000E+00 0.1111111E+00 0.0000000E+00 1 42 0.0000000E+00 0.1666667E+00 0.0000000E+00 1 43 0.0000000E+00 0.2222222E+00 0.0000000E+00 1 44 0.0000000E+00 0.2777778E+00 0.0000000E+00 1 45 0.0000000E+00 0.3333333E+00 0.0000000E+00 1 46 0.0000000E+00 0.3888889E+00 0.0000000E+00 1 47 0.0000000E+00 0.4444444E+00 0.0000000E+00 1 48 0.0000000E+00 0.5000000E+00 0.0000000E+00 1 49 0.0000000E+00 0.5000000E+00 0.5555556E-01 1 50 0.0000000E+00 0.5000000E+00 0.1111111E+00 1 51 0.0000000E+00 0.5000000E+00 0.1666667E+00 1 52 0.0000000E+00 0.5000000E+00 0.2222222E+00 1 53 0.0000000E+00 0.5000000E+00 0.2777778E+00 1 54 0.0000000E+00 0.5000000E+00 0.3333333E+00 1 55 0.0000000E+00 0.5000000E+00 0.3888889E+00 1 56 0.0000000E+00 0.5000000E+00 0.4444444E+00 1 57 0.0000000E+00 0.5000000E+00 0.5000000E+00 1 58 -0.3703333E-01 0.5185189E+00 0.5000000E+00 1 59 -0.7406667E-01 0.5370378E+00 0.5000000E+00 1 60 -0.1111000E+00 0.5555567E+00 0.5000000E+00 1 61 -0.1481333E+00 0.5740756E+00 0.5000000E+00 1 62 -0.1851667E+00 0.5925944E+00 0.5000000E+00 1 63 -0.2222000E+00 0.6111133E+00 0.5000000E+00 1 64 -0.2592333E+00 0.6296322E+00 0.5000000E+00 1 65 -0.2962667E+00 0.6481511E+00 0.5000000E+00 1 |
16楼2011-11-17 22:29:02
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★
ym23(金币+2):欢迎常来交流 2010-10-19 20:21:53
jttjtl(金币+5): 2010-10-20 10:08:36
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 ] |
2楼2010-10-19 19:36:04
3楼2010-10-20 10:11:24
wusongz
金虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 1077.3
- 红花: 1
- 帖子: 152
- 在线: 34.8小时
- 虫号: 839044
- 注册: 2009-09-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2010-10-20 10:13:02