| 查看: 2133 | 回复: 2 | ||
zw_tju木虫 (小有名气)
|
[求助]
【求助】GW或HSE杂化泛函带结构计算的问题
|
|
请教，我在使用HSE06或GW计算能带时出现问题：error in IBZKPT_HF: two k-points are equivalent 1 60 this will cause problems in the HF routine。请问这是什么问题，如何解决啊。。。。 使用的计算过程是网上给出的计算Si的模版： GW Band Structure Plots To obtain band structure plots from the GW method, first you have to perform a DFT calculation. The INCAR file I used for DFT on Si looks like this: ISTART = 0 ISYM = 0 ICHARG = 2 ISMEAR = 0 SIGMA = 0.1 LORBIT = 11 NSW = 200 IBRION = 2 ISIF = 3 EDIFF = 0.0001 EDIFFG = -0.01 KPOINTS file for the run is: Monkhorst Pack 0 Gamma 4 4 4 0 0 0 POSCAR file is: Si 5.38936000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 2 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 POTCAR file is supplied with VASP 5.2, remember to use paw pseudopotentials so that in the band structure calculation the Wigner-Seitz radius for Si need not be directly specified in INCAR. After this run completes you are ready to perform the HF-type (HSE06) calculation. Change the INCAR file to: ISTART = 1 ICHARG = 2 EDIFF = 0.0001 EDIFFG = -0.01 ENCUT = 400 ENAUG = 800 LREAL = .FALSE. LWAVE = .TRUE. LCHARG = .TRUE. NELM = 200 NSW = 0 IBRION = -1 LMAXMIX = 4 ISMEAR = 0 SIGMA = 0.1 NSIM = 4 IALGO = 48 ISYM = 0 LHFCALC = .TRUE. HFSCREEN = 0.2 ALGO = Damped TIME = 0.4 ENCUTFOCK = 0 NKRED = 2 AEXX = 0.25 Keep the other files constant for this run (e.g. KPOINTS, POSCAR, POTCAR need not be changed). Once this run has completed, add to the bottom of the INCAR file from the HSE06 calculation: LOPTICS = .TRUE. This step calculates the dielectric matrix Once this run has completed, change ALGO = Damped to ALGO = GWo and ICHARG = 2 to ICHARG = 11 and add: NOMEGA = 50 ENCUTGW = 100 LORBIT = 11 and change the KPOINTS file to: k-points along high symmetry lines 20 ! 20 intersections Line-mode rec 0 0 0 ! gamma 0.5 0.0 0 ! X 0.5 0.0 0 ! X 0.5 0.5 0.0 ! M 0.5 0.5 0.0 ! M 0.0 0.0 0.0 ! gamma 0.0 0.0 0.0 ! gamma 0.5 0.5 0.0 ! M 0.5 0.5 0.0 ! M 0.5 0.5 0.5 ! R 0.5 0.5 0.5 ! R 0.0 0.0 0.0 ! gamma These values are the points of high symmetry along the lines of high symmetry within the first Brillouin Zone of the Si unit cell structure in reciprocal space. And that’s it! Once complete you can plot the PROCAR, EIGENVAL, and DOSCAR files from the output of this run as described in the previous posts. The band structure plot produced from the plotting utility is shown below. Quantum Quibbles » GW Band Structure Plots [ Last edited by zw_tju on 2011-5-1 at 09:02 ] |
回复此楼» 猜你喜欢
fileCode有新解读？
已经有13人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 【文献学习】Li,Na掺杂ZnO空穴的杂化密度泛函计算(PRB，2009)
已经有21人回复
- 【专家在线答疑贴之三】杂化泛函计算FAQ
已经有28人回复
- vasp 杂化泛函HSE 它怎么就不run涅？
已经有12人回复
- HSE（GW类似）计算能带总结
已经有51人回复
- 用杂化泛函计算VO2，总是算不对，可有高手指导下啊，
已经有7人回复
- [悬赏交流]大家讨论一下vasp中HF以及杂化泛函的应用吧
已经有34人回复
- 【求助】求助关于各杂化密度泛函在计算上适用性 计算结果优劣分析
已经有3人回复
- 【求助】求助关于各杂化密度泛函在计算上适用性 计算结果优劣分析
已经有8人回复
- 【其他】VASP5.2的杂化泛函计算速度如何？
已经有3人回复
- 【求助】PWSCF是否支持杂化泛函？
已经有4人回复
- 【求助】GW /GW0/G0W0 计算是不是一定要基于HSE 计算的结果？
已经有4人回复
- 【求助】关于杂化泛函HSE计算
已经有20人回复
2楼2012-04-28 11:39:31
xueht987
木虫 (正式写手)
- 应助: 112 (高中生)
- 贵宾: 0.081
- 金币: 3514.2
- 散金: 106
- 红花: 6
- 帖子: 593
- 在线: 273.4小时
- 虫号: 1059317
- 注册: 2010-07-17
- 性别: GG
- 专业: 金属材料表面科学与工程
3楼2012-04-29 13:51:21
100