| 查看: 5237 | 回复: 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报错建模问题集 |
» 猜你喜欢
» 本主题相关价值贴推荐,对您同样有帮助:
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人回复
» 抢金币啦!回帖就可以得到:
科研需求-硫化镍钴精矿样品
+1/399
国家“千人计划”青年项目获得者团队最后2人07、08等学科研究生调剂
+1/284
伦敦布鲁内尔大学诚招国家公派读博,联合培养博士生,物流供应链方向
+1/91
广西高校3个材料与化工专业的硕士名额,计划做光谱、光电器件和机器学习相结合的工作
+1/88
物理学 调剂
+1/87
大湾区大学李红庚与中国科学技术大学王柳教授联合招收博士后
+1/85
平顶山学院 化工材料专业研究生报考调剂名额充足
+1/40
招生物专业考研生,有4个名额
+1/36
招收调剂生(一志愿07开头专业),海洋生物专业4个名额
+1/36
B区调剂 兰州交通大学材料学院调剂
+1/20
浙江理工大学蒋仲庆教授课题组招收物理或者材料背景研究生(调剂1名)
+1/16
粤港澳大湾区量子中心与南方科技大学联合招收2026级博士生(材料、凝聚态物理方向)
+1/16
求助请教各位:如果您的算法被顶尖院校学术不端怎么办?
+1/13
复试调剂-招收调剂生(一志愿07开头专业),海洋生物专业4个名额
+1/10
C13、N15标记植物,助力土壤碳氮循环与地学研究实操
+1/10
河南师大全国重点实验室—0780药学、08600生物与医药仍有调剂指标!
+1/7
最后调剂名额,丽水学院材料与能源工程学院
+1/5
江西水利电力大学市政工程接受调剂
+1/4
江西水利电力大学土木水利第二轮调剂系统已开
+1/4
活性环氧稀释剂
+1/3
6楼2011-03-11 00:08:47
2楼2010-12-02 00:00:21
3楼2010-12-02 08:37:51
★ ★
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














回复此楼