| 查看: 6404 | 回复: 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。 请问高手这样算的原因是什么呢？ |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
~~~点云数据求取k近邻并求取整个点云的特征点~~~
已经有17人回复- HSE06计算突然停止？
已经有6人回复
- vasp自洽计算K点不匹配问题，急！！
已经有15人回复
- 2011-06-04 PETS-5考试归来，听力题目还原，希望大家参与
已经有14人回复
- 【求助】GW0-Band structure erro
已经有10人回复
- 【求助】能带图横坐标中为什么总有特殊的K点？
已经有8人回复
- 【求助】vasp 5.2中用HSE06计算不能进入主循环
已经有7人回复
- 【求助】用castep计算BaTiO3的自发极化
已经有8人回复
- 【重点讨论】计算能带时如何确定高对称点和k-points path?
已经有40人回复
- 【求助】PBE0 计算Si、Ge能带结构
已经有9人回复
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
8楼2010-10-20 12:17:36
|
这是我在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) 想请高手讲讲这样算的原因~  |
9楼2010-10-20 19:01:58
29楼2021-01-22 09:37:42
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
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2010-10-20 10:51:36
|
我用了静态计算。 这是我计算能带时的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
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
7楼2010-10-20 12:05:04
 |
10楼2010-10-21 09:02:23