| 查看: 5178 | 回复: 18 | ||||||||||
[交流]
【求助】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人回复
-
上海电力大学材料防护与新材料重点实验室招收调剂研究生（材料、化学、电化学，环境）
已经有4人回复
-
材料学求调剂
已经有6人回复
-
303求调剂
已经有5人回复
-
一志愿武理085500机械专业总分300求调剂
已经有7人回复
-
考研调剂
已经有4人回复
-
281求调剂
已经有4人回复
-
0805 316求调剂
已经有6人回复
-
085601求调剂总分293英一数二
已经有3人回复
-
08工学调剂
已经有17人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
» 抢金币啦！回帖就可以得到:
-
中国散裂中子源诚聘“AI材料计算”方向博士后/职工
+1/427
-
香港中文大学（深圳）靳羽华教授交叉实验室招募2026年材料学博士生（光致变色）
+2/134
-
材料、化学类硕士研究生调剂机会，学硕专硕均可，有意请联系，非诚勿扰！
+1/84
-
江南大学招收具有3年海外博后经历的青年教师
+1/82
-
中山大学环境材料团队招聘
+1/81
-
【全额奖学金+英国博士学位】国科大杭高院Nathan课题组招收2026级联培博士生1名！
+1/47
-
武汉高校接受研究生调剂，材料、化学、纺织方向
+2/44
-
哈尔滨工业大学（深圳）-何自开教授团队诚招化学工程与技术专业博士研究生
+1/34
-
Alicat 层流压差质量流量控制器在新型硅碳负极流化床的应用- 艾里卡特（Alicat）
+3/19
-
同济大学环境学院 肖倩特聘研究员课题组 招聘硕士/博士（长期有效）
+1/17
-
浙江中医药大学副校长团队招聘药学/生物信息学等方向高校教师7名，报名通道已开启
+1/12
-
南京林业大学化学工程学院柏惺峰课题组诚聘科研助理、水杉英才E类/F类青年教师或博后
+1/9
-
香港大学 化学系 博后招聘
+1/9
-
中国中医科学院中医药健康产业研究所（江西中医药健康产业研究院）2026年高层次人才等
+1/6
-
哈尔滨工业大学(深圳)国家级青年人才机械工程、热能与能源等专业26年入学博士招生
+1/6
-
中国人民大学物理学系蔡鹏课题组招收优秀硕士博士研究生
+1/6
-
宁波东方理工大学余鹏课题组2026年上海交大联培博士生、博士后、科研助理招聘
+1/5
-
中国地质大学（武汉）清洁能源与资源循环团队诚招科研助理/2026级硕士、博士生/博士后
+1/5
-
浙江大学高纯分离团队诚聘博士后
+1/1
-
11自命题 294分求调剂
+1/1
★ ★
zzy870720z(金币+2):谢谢指教 2010-12-04 21:36:29
wuli8(金币+100):感谢。我再试一试 2010-12-04 22:51:40
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点的数目是两部分的总和，画能带的时候只取下面高对称部分的本征值即可 |
4楼2010-12-04 21:15:41
2楼2010-12-02 00:00:21
3楼2010-12-02 08:37:51
5楼2011-03-11 00:08:02
6楼2011-03-11 00:08:47
7楼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
9楼2011-03-14 20:22:26
10楼2011-03-14 21:58:20
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
楼主找到例子了吗?我也需要关于这种算法计算的例子,楼主得到了例子计算文件的话,希望楼主给发一份,不胜感激,万分感谢!elijia2010@hotmail.com 此致, 敬礼! |
11楼2011-03-24 15:55:29
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
gzqdyouxia(金币+1): 鼓励交流 2011-04-01 10:26:04
小木虫(金币+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. ），那我是不是直接运行一遍自洽运算就可以得出结果，不需要像一楼所说的那样进行两遍运算，因为毕竟这个运算是相当耗时间的。 有哪位高人可以指点一下？ |
12楼2011-03-31 21:31:26
13楼2011-04-01 15:46:24
14楼2011-04-02 09:35:22
继续求助 15楼2011-09-26 15:58:46
16楼2012-07-25 20:59:00
17楼2013-03-18 09:37:05
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
老是有这样的提示： | | | 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
19楼2017-04-03 22:15:01
