| 查看: 7695 | 回复: 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 |
» 猜你喜欢
建议基金发布提前给出明确的时间点
已经有14人回复
-
范进中举一文的中心思想
已经有5人回复
-
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- ~~~点云数据求取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人回复
|
这是我在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
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
